Papers

Every smooth first-order real planar elliptic system admits a universal complex form $w_{\bar z} - \mu w_z + \mathcal{A} w + \mathcal{B} \bar w = \mathcal{F}$, which we call the Beltrami-Vekua equation: the data $(\mu, \mathcal{A}, \mathcal{B}, \mathcal{F})$ are produced from the original system by algebraic operations and differentiations, with no auxiliary PDE. On this space we study the joint action of multiplicative gauges $w \mapsto \phi w$ and orientation-preserving diffeomorphisms. Our main result is that the 2-form $\Theta = |\mathcal{B}|^2 / (1 - |\mu|^2) \, dx \, dy$ is gauge-invariant and pulls back covariantly under diffeomorphisms; its form is forced, with $|\mathcal{B}|^2$ the unique $\mathcal{B}$-quadratic combination invariant under $\mathcal{B} \mapsto \mathcal{B}\phi/\bar\phi$ and $1 - |\mu|^2$ the conformal distortion factor from the diffeomorphism law for $\mu$. The total mass $\mathcal{M}(D) = \int_\Omega \Theta$, the pseudo-analytic mass, vanishes precisely on the analytic class $\mathcal{B} \equiv 0$ and separates a continuous family of pairwise inequivalent pseudo-analytic equations on the disk. As a by-product, Vekua's two-stage reduction - uniformization then gauge elimination - requires only one variable-coefficient PDE solve: the Beltrami diffeomorphism supplies the integrating factor for a flat $\bar\partial$-equation.

Motivated by the average partition function of c free bosons $($Afhkami-Jeddi et al. \cite{Afhk2021}$)$ and the average of the genus 1 partition function over the Narain moduli space $($Maloney-Witten \cite{Witten2020}$)$, we investigate ratios of theta functions. In this paper, we completely classify the minimizers (or maximizers) for ratios of theta and Epstein zeta functions. We find that the hexagonal lattice plays a pivotal role there. These results have direct applications in conformal and Liouville field theory via partition functions. Additionally, they yield the minima of differences of theta and Epstein zeta functions, which have implications for the mathematics of crystallization and interacting particle theory (\cite{Bet2016,Bet2019AMP}).

We consider the isoperimetric inequality involving the $s$-perimeter and the $t$-perimeter with $0<s<t<1$, and show that the ball is a local minimizer of the (scale-invariant) isoperimetric ratio $\mathcal{F}(E):=P_t(E)^{\frac{1}{n-t}}/ P_s(E)^{\frac{1}{n-s}}$ among sets $E$ that are nearly spherical. To this end, we rewrite $\mathcal{F}$ as a functional of $u$, where $u$ is a scalar function on the unit sphere in $\mathbb{R}^n$ that parametrizes the boundary of $E$, and prove a quantitative stability result for $\mathcal{F}$ around $u=0$ with respect to a suitable Sobolev norm. This parallels known results where the $s$-perimeter is replaced by the volume.

We establish the nonlinear stability threshold $O(\nu^{3/2})$ for the three-dimensional Couette flow governed by the compressible Navier–Stokes equations. While stability thresholds are well understood in two dimensions for both compressible and incompressible flows, and in three dimensions for incompressible flows, the three-dimensional compressible case remains open due to additional structural features, strong mode interactions, and wave coupling. The proof is based on a refined frequency-space approach. For zero modes, we improve upon two-dimensional methods by clearly separating and precisely estimating the main contributions from diffusion waves, acoustic waves, and the lift-up mechanism, leading to a systematic way to handle their nonlinear coupling. For the non-zero modes, we introduce new multiplier estimates and a decomposition based on the structure of the compressible system, which allows us to track the interaction between dissipation and acoustic effects.

It has been observed that the dynamics of the Toda lattice can be well described by solutions of the Korteweg-de Vries (KdV) equation in the continuum limit. We show that, under the KdV scaling and a suitable translation, the solution of the Toda lattice with H^1 initial data converges to that of the KdV equation globally in time. Our proof relies on tools from harmonic analysis and also on the construction and the conservation of mass and energy of the Toda lattice, the latter of which are derived from the completely integrable structure of the Toda lattice. As a consequence, we obtain long-wave KdV limits for the Toda lattice.

This paper studies singularity formation for conservation laws with a source term. Motivated by John (1974) and Barlin (2023), we prove finite-time blow-up under initial data conditions weaker than those in Barlin. Moreover, we show that a sufficiently small compact support length of the initial data promotes blow-up. Hence, global existence can only be achieved when the initial data have a large compact support length.

We present here the Robin boundary condition and its significance in mathematics and physics.

We study the 1D Brenner-Navier-Stokes-Fourier (BNSF) system, proposed as a refinement of the classical Navier–Stokes–Fourier model through the introduction of the volume velocity, distinct from the mass velocity describing convective transport. When formulated in the Lagrangian mass coordinates with the volume velocity, the discrepancy between the two velocities induces a dissipative structure in the mass conservation law. We prove the global existence of classical solutions for arbitrarily large initial data. More precisely, for initial data in $H^k(\mathbb{R})$ with $k\ge3$, with the specific volume and absolute temperature initially bounded away from zero, we construct global-in-time solutions that remain in the same regularity class. Our result accommodates arbitrarily large initial data. A major difficulty is to establish lower and upper bounds for the specific volume \(v\). The additional dissipation yields an $L_t^2 L_x^2$ bound for $v_x$, which is further improved to an $L_t^\infty L_x^\infty$ bound of $v$ and $1/v$ via the parabolic De Giorgi method. We also apply the maximum principle to obtain a positive lower bound for the absolute temperature.

The Vlasov-Poisson-BGK (VPBGK) model is a kinetic model for describing the dynamics of collisional plasmas. Although various numerical schemes have been developed for it, a corresponding convergence theory has been absent. This paper fills this gap by presenting the first convergence analysis for a non-splitting, finite-difference Eulerian scheme discretized on the full phase-space grid. A major theoretical obstacle is the mixing of velocity indices induced by the electric field, which hinders the derivation of a uniform lower bound for the discrete solution. To overcome this stability challenge, we propose a modified lower bound estimate suitable for ionized systems that incorporates the step-wise degradation. Under a truncated velocity domain with a Neumann boundary condition, we establish error estimates for the distribution function in a weighted $L^{\infty}$ norm and for the electric field in a $L^{\infty}$ norm, respectively.

This paper studies the upper and lower bounds of the lifespan for the classical solutions to the initial value problems of one dimensional wave equations with non-autonomous semilinear terms including the space-derivative of the unknown function.This is a non-trivial business comparing to the analogous results with time-derivative type semilinear terms, especially for the proof to obtain the sharp upper bound of the lifespan as we have to deal with space dependent weights among iteration procedures of the weighted functional of the solution. Also it is surprising that a part of them reaches to the same ordinary differential inequality for classical semilinear damped wave equations introduced by Li and Zhou (Discrete Contin. Dynam. Systems, 1995, 1(4): 503-520), and we show a simple proof for blow up result from this ordinary differential inequality by iteration argument and slicing method in more general situation.

In this paper, we are focusing on the proof of the blow-up result for a quadratic semilinear wave equation in two space dimensions. There is a logarithmic loss in estimating the lifespan of classical solutions if the 0th moment of the initial speed does not vanish. This result is already known with almost sharp constants. But in order to have a direct application to numerical analysis, we show a simple proof by iteration argument of point-wise estimate of the solution with the slicing technique.

We prove the existence, uniqueness, and sharp bilateral pointwise estimates for positive bounded solutions to the Lane–Emden type problem \[ \begin{cases} L u = \sum\limits_{i=1}^{m}\sigma_{i} u^{q_{i}}+\sigma_0, \quad u\geq0 & \text{in } \Omega, \liminf \limits_{x \rightarrow y} u(x) = f(y), & y \in \partial^\infty\Omega, \end{cases} \] where $0 < q_{i} < 1$. Here $Lu = - \text{div}(A \nabla u)$ is a uniformly elliptic operator with bounded coefficients, $\sigma_{i}$ is a nonnegative locally finite Borel measure on an $A$-regular domain $\Omega \subset \mathbb{R}^n$ which possesses a positive Green function associated with $L$, and $f$ is a nonnegative continuous function on the boundary $\partial^\infty\Omega$. An analogous result for positive continuous solutions to the problem is also illustrated. Our method can be adapted to address related sublinear problems with zero boundary conditions involving the fractional Laplace operator $(-\Delta)^{\alpha}$ for $0< \alpha < n/2$, in place of $L$, in $\mathbb{R}^n$ as well.

We establish a criterion for the validity of the classical (non-semiclassical) Weyl law for Schrödinger operators $ H=\Delta+V $ on complete Riemannian manifolds. In contrast to existing results, our approach does not rely on standard geometric assumptions such as bounded geometry, nor on analytic assumptions such as the doubling condition on the potential. Instead, we identify a geometric-analytic invariant that encodes the precise balance between the geometry of the manifold, the growth of $V$, and the oscillation scale of $V$. This intrinsic quantity, denoted $c_{\delta}(\lambda)$ admits effective quantitative estimates. We prove that the Weyl asymptotic holds provided $\lim_{\lambda\to\infty} c_\delta(\lambda)=0 .$ The sharpness of this criterion is demonstrated through explicit examples showing that the Weyl law can fail when the criterion is violated.

This paper presents results on the unboundedness and minimal speed of traveling wave solutions for a one-dimensional spatial reaction-diffusion equation with an asymptotically linear reaction term and a saturation parameter. By applying a Poincaré-type compactification, we reveal the full dynamics (including infinity) of the two-dimensional system of ordinary differential equations satisfied by traveling wave solutions. This yields essential information characterizing traveling wave solutions: the classification of trajectories in the phase plane, the positivity and unboundedness of front-type and sign-changing profiles, and the explicit form of the minimal speed. This paper examines a special equation with an asymptotically linear reaction term. While, our results differ from those of conventional linear determinacy. We claim that the minimal speed is derived from information at infinity within the traveling wave system.

We fix the proof of Theorem 1.17 in the quoted reference, which does not in general gives the correct coefficients in the asymptotic behavior of eigenvalues for the Dirichlet Laplacian when a small hole is removed from the domain.

We consider a system of elliptic equations, depending on a small parameter $\eps$, that models long-range segregation of populations. The diffusion is governed by the Laplacian. This system was previously investigated by Caffarelli, Patrizi, and Quitalo in \cite{CL2} as a model in population dynamics, and they established the regularity of the free boundary in two dimensions. In this paper we study the free boundary in the higher dimensional case. We extend the concept of angles and asymptotic cones to higher dimensions, and give a characterization of regular and singular points in terms of their densities and angles. We obtain a structure result of the free boundary and show that, if the angles at the singular points are away from $\frac{n\omega_n}{2}$, the regular set is open in the free boundary and locally a $C^1$ manifold of dimension $n-1$. We also show that, if the supports of the populations are convex, they are convex polytopes. A weak form of the equality of angles for the convex configuration is also derived.

We establish sharp regularity for the value function, the pressure, and the free boundary in one-dimensional first-order mean field games with power coupling and compactly supported density. Under a standard nondegeneracy assumption on the initial datum, the pressure \(p=m^\theta\) is Lipschitz continuous, the value function \(u\) is \(C^{1,1/2}\), and the two free boundary curves are smooth in time. If the initial pressure is smooth, then both \(p\) and \(u\) are smooth up to the free boundary from inside the positive phase. The proof works in Lagrangian coordinates and, through a singular change of variables, recasts the boundary degeneracy as a removable radial axis in effective dimension \(N=4+2/\theta\), allowing the application of recent estimates for even solutions to elliptic problems with degenerate weights.

The variational model for a ferroelectric nematic bears close resemblance to the well-known energy model for micromagnetics. Despite this similarity, the two models operate in fundamentally distinct parameter regimes describing different physics. In this paper we establish that the ferroelectric nematic energy functional $\Gamma$-converges to the energy of a nematic with high elastic anisotropy.

We develop a classical well-posedness and regularity theory on a finite connected weighted graph for an extended mean field game system, its associated master equation, and a Hamilton-Jacobi- Bellman equation on the probability simplex, all in the presence of an individual noise operator. The geometric structure is inherited from the logarithmic-mean activation functional of discrete optimal transport, under which the entropic Fokker-Planck equation appears as a gradient flow on the graph and the individual noise operator is a bilinear form in the probability vector and the Wasserstein gradient. A central technical step is a quantitative preservation-of-positivity estimate for the discrete continuity equation, which rules out finite-time boundary degeneracy and yields a classical solution theory for the master equation on the open simplex without imposing any boundary condition. As an application, we recover a Nash equilibrium interpretation of the discrete system in terms of Markov chains on the graph. Our setup is inspired by the computational algorithms for optimal mass transport of [10, 11] and provides a rigorous well-posedness theory for several of the equations derived in [25].

We prove that the short-period eigenfunctions of quantum cat maps constructed by Kim and the author equidistribute on $\mathbb{T}^2$ in the sense of semiclassical measures. We also show that their logarithmically large $\ell^\infty$-norm is asymptotically concentrated on a bounded number of coordinates. Thus, for this explicit family, strong coordinate localization coexists with semiclassical equidistribution. These results confirm the behavior suggested by earlier numerical evidence of Kim and the author, and contrast with the scarring phenomena for short-period eigenfunctions observed by Faure, Nonnenmacher, and De Bièvre.

We establish sharp local existence results for the Hirota-Satsuma system in $H^k(\mathbb{R}) \times H^s(\mathbb{R})$, depending on the ratio between the dispersion of the components. These theorems significantly generalize previous works, which were restricted to the diagonal case of equal regularity $s=k$. Furthermore, we extend the known global well-posedness theory to the off-diagonal regime. The argument relies on the Fourier restriction norm method coupled with the concept of integrated-by-parts strong solution - a framework that generalizes the classical notion of strong solution.

We show that intrinsic and extrinsic area density bounds are equivalent, with matching asymptotic values, for complete, connected, smooth minimal immersions $i:\Sigma^d\to\mathbb{R}^N$ of any dimension and codimension. Combining our results with a recent breakthrough by Bellettini, we extend the Schoen–Simon–Yau curvature estimates for smoothly immersed, two-sided, stable minimal hypersurfaces $i:\Sigma^n\to\mathbb{R}^{n+1}$ with bounded intrinsic area density to the missing case $n=6$, which had remained open since.

We consider the evolution and interaction of a 2-dimensional nematic liquid crystal of Ericksen-Leslie type within a 1-dimensional flexible viscoelastic structure. This is a fully macroscopic model in which the nematic liquid crystal is modelled by the simplified Ericksen-Leslie system with Ginzburg-Landau approximation. The liquid crystal is contained in a thin viscoelastic shell of arbitrary reference configuration that evolves with respect to the forces exerted by the liquid crystal. Barring any degeneracies in the shell, we construct a global weak solution for the coupled system. We then show that any family of such weak solutions that are parametrized by the Ginzburg-Landau coefficient, converges to a weak solution of the original simplified Ericksen-Leslie system without the Ginzburg-Landau term.

We investigate the formal stability of finite-amplitude non-zonal flows bifurcating from the trivial state in the unforced 2D Euler equations on the sphere. To bypass the degeneracy of the spherical Laplacian and filter out the low-frequency Fjørtoft instabilities, we restrict the functional space to the invariant subspace of the tetrahedral symmetry group. Using Arnold's Energy-Casimir method, we prove that the linearized elliptic operator derived via Liapunov-Schmidt reduction acts as the Hessian of the conserved functional. By tracking the critical eigenvalue along the bifurcating branches via the Crandall-Rabinowitz theorem, we establish a relation between the bifurcation topology and formal stability. Applying this framework to four distinct geophysical profile functions, we demonstrate that subcritical polynomial and supercritical sine-Gordon flows achieve a negative-definite second variation, that is, their formal stability. In contrast, subcritical sinh-Gordon and supercritical Liouville exponential flows generate saddle points, making them unstable. This classification identifies the specific nonlinear interactions required for the persistence of large-scale coherent waves in planetary atmospheres.

In this paper we address the existence and uniqueness of entire spacelike hypersurfaces in the Lorentz–Minkowski space $\mathbb{L}^{m+1}$ with prescribed mean curvature that are star-shaped with respect to a point and asymptotic to a light cone. We also establish a Willmore-type inequality and prove a non-existence result for spacelike radial graphs asymptotic to the light cone whose mean curvature belongs to $L^p$ for $1 \leq p\leq m$, in particular in the case of compactly supported mean curvature.

We prove a stability theorem for finite-dimensional analytic inverse problems. Let \(U\subset\R^m\) be an open parameter set, let \(F(p)\) be a boundary measurement operator, and let \(R(p)\) be the finite-dimensional quantity to be recovered. If \(F\) is real analytic and \[ F(p)=F(q)\quad\Longrightarrow\quad R(p)=R(q), \] then \(R\) satisfies a Hölder stability estimate on every compact subset of \(U\). The proof uses a Hilbert–Schmidt scalarization of the operator equation \(F(p)=F(q)\) and the Łojasiewicz distance inequality. We also prove that, after fixing countable dense families of boundary inputs and tests, finitely many scalar matrix elements of the data give the same Hölder recovery on compact parameter sets. This finite-measurement conclusion is qualitative: the proof does not give an effective measurement list, exponent, or constant. The finite-measurement statement follows from finite determinacy of real analytic zero sets. We apply the result to local Neumann-to-Dirichlet data for piecewise constant anisotropic conductivities and to localized Dirichlet-to-Neumann data for piecewise homogeneous anisotropic elasticity.

We prove small energy regularity for a parabolic boundary reaction Ginzburg-Landau problem in the full range $s\in (0,1)$, answering a question posed by Hyder, Segatti, Sire and Wang. We also obtain a similar small energy regularity result for fractional harmonic maps to spheres. Both results are uniform as $s\to 1$.

We study families $u_s$ of functions satisfying the equations $(-\Delta)^s u_s=0$, $s \in (0,1)$ in a smooth bounded open set $\Omega \subset \mathbb{R}^N$. The main purpose of this paper is twofold. First, we provide a detailed analysis of the asymptotics of these families in the zero order limit $s \to 0^+$. Second, we study the differentiability of $u_s$ as a function of $s$. Most of our results are devoted to the associated Poisson problem, where the family $u_s$ is determined by the exterior condition $u_s = g$ in $\mathbb{R}^N \setminus \Omega$ for some fixed function $g \in L^\infty(\mathbb{R}^N \setminus \Omega)$. Our results show that both the zero order asymptotics and the differentiability properties of $u_s$ can be expressed in terms of the logarithmic Laplacian of suitable extensions of $g$. This allows to deduce pointwise monotonicity properties of $u_s$ in the order parameter $s$ for a large class of functions $g$.

We consider standing waves of the nonlinear Schrödinger equation $i\partial_t u = -\Delta_\alpha u + |u|^{p-1}u$ in the defocusing case in dimensions $N=2$ and $N=3$. Here, $-\Delta_\alpha$ denotes the Laplacian with a point interaction. This operator is bounded from below by a negative constant; consequently, unlike in the free case, the associated energy functional admits non-trivial minimizers. We establish existence and uniqueness of standing waves, and prove further qualitative properties, including radial symmetry, positivity, and stability. Moreover, we build an appropriate functional space for the zero-mass case and establish sharp decay estimates in this case.

This paper studies the 1D stochastic Allen–Cahn equation on a bounded domain driven by localized white noise. We prove that the associated Markov process admits a unique invariant measure and is exponential mixing. The main challenge lies in the interaction between localized nature of the noise and non-trivial global dynamics of the system. To overcome this, our approach relies on two ingredients from PDE control theory: stabilization for the linearized system and global steady-state controllability for the nonlinear equation. The stabilization result is derived using the weak observability and Fenchel–Rockafellar duality, while the global controllability relies on quasi-static deformations combined with global dynamics.