2025-04-02 | | Total: 23
We show that mean field optimal controls satisfy a first order optimality condition (at a.e. time) without any a priori requirement on their spatial regularity. This principle is obtained by a careful limit procedure of the Pontryagin maximum principle for finite particle systems. In particular, our result applies to the case of mean field selective optimal control problems for multipopulation and replicator dynamics.
We construct a time-independent, incompressible, and Lipschitz-continuous velocity field in R3 that generates a fast kinematic dynamo - an instability characterized by exponential growth of magnetic energy, independent of diffusivity. Specifically, we show that the associated vector transport-diffusion equation admits solutions that grow exponentially fast, uniformly in the vanishing diffusivity limit ε→0. Our construction is based on a periodic velocity field U on T3, such as an Arnold-Beltrami-Childress flow, which satisfies a generic spectral instability property called alpha-instability, established via perturbation theory. This provides a rigorous mathematical framework for the alpha-effect, a mechanism conjectured in the late 1960s to drive large-scale magnetic field generation. By rescaling with respect to ε and employing a Bloch-type theorem, we extend the solution to the whole space. Finally, through a gluing procedure that spatially localizes the instability, we construct a globally defined velocity field u in R3 that drives the dynamo instability.
Korn-Maxwell-Sobolev (KMS) inequalities represent a tool for estimating differential expressions and have gained particular importance in recent years, especially concerning elliptic operators. In my Master's thesis, together with Peter Lewintan (University of Duisburg-Essen), we extended this concept to also apply to operators of constant rank. This makes it possible to cover more complex structures such as the curl or divergence of vector fields. A key difference from the elliptic theory is that in the constant rank case, a correction term ΠB is necessary on the left-hand side of the inequality. Results were also obtained for the limiting case p=1, although additional assumptions are required here. This article provides an illustrative introduction to KMS inequalities and demonstrates their application in both the elliptic and constant rank cases.
This paper concerns the eigenvalues of the Neumann-Poincaré operator, a boundary integral operator associated with the harmonic double-layer potential. Specifically, we examine how the eigenvalues depend on the support of integration and prove that the map associating the support's shape to the eigenvalues is real-analytic. We then compute its first derivative and present applications of the resulting formula. The proposed method allows for handling infinite-dimensional perturbation parameters for multiple eigenvalues and perturbations that are not necessarily in the normal direction.
We consider the existence of stationary wave solutions with prescribed mass to a supercritical nonlinear Schrödinger equation on a noncompact connected metric graph without a small mass assumption.
We consider the convected Helmholtz equation with a generalized Myers boundary condition (a boundary condition of the second-order) and characterize the set of physical parameters for which the problem is weakly well-posed. The model comes from industrial applications to absorb acoustic noise in jet engines filled with absorbing liners (porous material). The problem is set on a 3D cylinder filled with a d-upper regular boundary measure, with a real 1 < d ≤ 2. This setup leads to a parametric shape optimization problem, for which we prove the existence of at least one optimal distribution for any fixed volume fraction of the absorbing liner on the boundary that minimizes the total acoustic energy on any bounded wavenumber range.
The first-order approach to boundary value problems for second-order elliptic equations in divergence form with transversally independent complex coefficients in the upper half-space rewrites the equation algebraically as a first-order system, much like how harmonic functions in the plane relate to the Cauchy-Riemann system in complex analysis. It hinges on global Lp -bounds for some p > 2 for the resolvent of a perturbed Dirac-type operator acting on the boundary. At the same time, gradients of local weak solutions to such equations exhibit higher integrability for some p > 2, expressed in terms of weak reverse Hölder estimates. We show that the optimal exponents for both properties coincide. Our proof relies on a simple but seemingly overlooked connection with operator-valued Fourier multipliers in the tangential direction.
This paper aims to develop a theory for linear-quadratic Nash systems and Master equations in possibly infinite-dimensional Hilbert spaces. As a first step and motivated by the recent results in [31], we study a more general model in the linear quadratic case where the dependence on the distribution enters just in the objective functional through the mean. This property enables the Nash systems and the Master equation to be reduced to two systems of coupled Riccati equations and backward abstract evolution equations. We show that solutions for such systems exist and are unique for all time horizons, a result that is completely new in the literature in our setting. Finally, we apply the results to a vintage capital model, where capital depends on time and age, and the production function depends on the mean of the vintage capital.
We prove optimal Hölder boundary regularity for a non-local operator with a singular, symmetric kernel that depends on the distance to the boundary of the underlying domain. Additionally, we prove higher boundary regularity of solutions.
In this paper, we study the parabolic equations of the form {Lu(y,t)=f,(y,t)∈Q,u(y,t)=0,(y,t)∈∂Q,u(y,t)is uniformly bounded from below,(y,t)∈Q, where Q=Ω×R⊂Rn+1 and Ω⊂Rn is a bounded Lipschitz domain with 0∈Ω. Here L is a general second order uniformly parabolic differential operator in non-divergence form or divergence form. For f=0, we establish the structure of the solution space, which is one dimensional and the solutions in this space grow exponentially at one end and decay exponentially at the other. For f≠0, we show that all solutions can be presented by the solutions corresponding to the homogenous equations(f=0) and a bounded special solution of the inhomogeneous equations. Our method is based on maximum principle in Q and the Harnack type inequalities.
We consider the free boundary problems of degenerate elliptic equations that describe the level set formulation of the interface motion evolved by anisotropic forced mean curvature flows. The type of free boundary problems in this paper was initially studied as the first-order Hamilton-Jacobi-Isaacs equations arising in pursuit-evasion differential games and applied to the models of first-order front propagation in Soravia (1994). In this paper, we consider an extension of these free boundary problems to the second-order equations and give a deterministic game representation based on a discrete approximation scheme in Kohn and Serfaty (2006). Furthermore, we prove the comparison principle for our free boundary problems by using the framework of time-discrete games.
We show the existence of local and global in time weak martingale solutions for a stochastic version of the Othmer-Dunbar-Alt kinetic model of chemotaxis under suitable assumptions on the turning kernel and stochastic drift coefficients, using dispersion and stochastic Strichartz estimates. The analysis is based on new Strichartz estimates for stochastic kinetic transport. The derivation of these estimates involves a local in time dispersion analysis using properties of stochastic flows, and a time-splitting argument to extend the local in time results to arbitrary time intervals.
The sine-Gordon equation is a fundamental nonlinear partial differential equation that governs soliton dynamics and phase evolution in a variety of physical systems, including Josephson junctions and superconducting circuits. In this study, we analyze the effects of external perturbations such as damping and driving forces on the stability of soliton solutions. Using a rigorous Hilbert space framework, we establish well-posedness and derive regularity results for the perturbed equation. In particular, we provide sufficient conditions for the boundedness of the perturbation function, which plays a crucial role in determining the persistence of the soliton structure. Furthermore, we propose a neural network-based approach for solving the inverse problem of estimating unknown initial conditions. By training a data-driven model on simulated PDE solutions, we demonstrate that the network can accurately recover the initial states from limited or noisy observations. Numerical simulations validate the theoretical results and highlight the potential of combining mathematical analysis with machine learning techniques to study nonlinear wave phenomena. This approach offers valuable insights into soliton behavior and has potential applications in the design of quantum computing systems based on Josephson junctions.
The generator L of the linearized evolution equation of adiabatic oscillations of a gaseous star, ELASO, is a second order integro-differential operator and is realized as a self-adjoint operator in the Hilbert space of square integrable unknown functions with weight, which is the density distribution of the compactly supported background. Eigenvalues and eigenfunctions of the operator L have been investigated in practical point of view of eigenmode expansion of oscillations. But it should be examined whether continuous spectra are absent in the spectrum of L or not. In order to discuss this question, the existence of essential spectra in a closely related evolution problem is established.
This paper is devoted to the analysis of the incompressible Euler equation in a time-dependent fluid domain, whose interface evolution is governed by the law of linear elasticity. Our main result asserts that the Cauchy problem is globally well-posed in time for any irrotational initial data in the energy space, without any smallness assumption. We also prove continuity with respect to the initial data and the propagation of regularity. The main novelty is that no dissipative effect is assumed in the system. In the absence of parabolic regularization, the key observation is that the system can be transformed into a nonlinear Schrödinger-type equation, to which dispersive estimates are applied. This allows us to construct solutions that are very rough from the point of view of fluid dynamics-the initial fluid velocity has merely one-half derivative in L2. The main difficulty is that the problem is critical in the energy space with respect to several key inequalities from harmonic analysis. The proof incorporates new estimates for the Dirichlet-to-Neumann operator in the low-regularity regime, including refinements of paralinearization formulas and shape derivative formulas, which played a key role in the analysis of water waves.
This paper investigates the Dirichlet problem for a non-divergence form elliptic operator L in a bounded domain of R2. Assuming that the principal coefficients satisfy the Dini mean oscillation condition, we establish the equivalence between regular points for L and those for the Laplace operator. This result bridges a gap in the authors' recent work on higher-dimensional cases. Furthermore, we construct the Green's function for L in regular two-dimensional domains, extending a result by Dong and Kim.
In this work, we propose an exponentially convergent numerical method for the Caputo fractional propagator Sα(t) and the associated mild solution of the Cauchy problem with time-independent sectorial operator coefficient A and Caputo fractional derivative of order α∈(0,2) in time. The proposed methods are constructed by generalizing the earlier developed approximation of Sα(t) with help of the subordination principle. Such technique permits us to eliminate the dependence of the main part of error estimate on α, while preserving other computationally relevant properties of the original approximation: native support for multilevel parallelism, the ability to handle initial data with minimal spatial smoothness, and stable exponential convergence for all t∈[0,T]. Ultimately, the use of subordination leads to a significant improvement of the method's convergence behavior, particularly for small α<0.5, and opens up further opportunities for efficient data reuse. To validate theoretical results, we consider applications of the developed methods to the direct problem of solution approximation, as well as to the inverse problem of fractional order identification.
Throughout physics Brownian dynamics are used to describe the behaviour of molecular systems. When the Brownian particle is confined to a bounded domain, a particularly important question arises around determining how long it takes the particle to encounter certain regions of the boundary from which it can escape. Termed the first passage time, it sets the natural timescale of the chemical, biological, and physical processes that are described by the stochastic differential equation. Probabilistic information about the first passage time can be studied using spectral properties of the deterministic generator of the stochastic process. In this work we introduce a framework for bounding the leading eigenvalue of the generator which determines the exponential rate of escape of the particle from the domain. The method employs sum-of-squares programming to produce nearly sharp numerical upper and lower bounds on the leading eigenvalue, while also giving good numerical approximations of the associated leading eigenfunction, the quasi-stationary distribution of the process. To demonstrate utility, the method is applied to prototypical low-dimensional problems from the literature.
We describe the transformation of a polynomial planar dynamical system into a second order differential equation by means of a polynomial change of variables. We then, by means of the Krylov-Bogoliubov-Mitropolsky averaging method, identify sufficient conditions involving said change of variables so that a limit cycle exists.
This article considers the problem of reconstructing unknown driving forces based on incomplete knowledge of the system and its state. This is studied in both a linear and nonlinear setting that is paradigmatic in geophysical fluid dynamics and various applications. Two algorithms are proposed to address this problem: one that iteratively reconstructs forcing and another that provides a continuous-time reconstruction. Convergence is shown to be guaranteed provided that observational resolution is sufficiently high and algorithmic parameters are properly tuned according to the prior information; these conditions are quantified precisely. The class of reconstructable forces identified here include those which are time-dependent and potentially inject energy at all length scales. This significantly expands upon the class of forces in previous studies, which could only accommodate those with band-limited spectra. The second algorithm moreover provides a conceptually streamlined approach that allows for a more straightforward analysis and simplified practical implementation.
We investigate in this chapter the mathematical models for electromagnetic wave propagation in dispersive isotropic passive linear media for which the dielectric permittivity ε and magnetic permeability μ depend on the frequency. We emphasize the link between physical requirements and mathematical properties of the models. A particular attention is devoted to the notions of causality and passivity and its connection to the existence of Herglotz functions that determine the dispersion of the material. We consider successively the cases of the general passive media and the so-called local media for which ε and μ are rational functions of the frequency. This leads us to analyse the important class of non dissipative and dissipative generalized Lorentz models. In particular, we discuss the connection between mathematical and physical properties of models through the notions of stability, energy conservation, dispersion and modal analyses, group and phase velocities and energy decay in dissipative systems.
In this second chapter, we analyse transmission problems between a dielectric and a dispersive negative material. In the first part, we consider a transmission problem between two half-spaces, filled respectively by the vacuum and a Drude material, and separated by a planar interface. In this setting, we answer to the following question: does this medium satisfy a limiting amplitude principle? This principle defines the stationary regime as the large time asymptotic behavior of a system subject to a periodic excitation. In the second part, we consider the transmission problem of an infinite strip of Drude material embedded in the vacuum and analyse the existence and dispersive properties of guided waves. In both problems, our spectral analysis enlighten new and unusual physical phenomena for the considered transmission problems due to the presence of the dispersive negative material. In particular, we prove the existence of an interface resonance in the first part and the existence of slow light phenomena for guiding waves in the second part.
We investigate the zero-noise limit for SDE's driven by Brownian motion with a divergence-free drift singular at the initial time and prove that a unique probability measure concentrated on the integral curves of the drift is selected. More precisely, we prove uniqueness of the zero-noise limit for divergence-free drifts in L1loc((0,T];BV(Td;Rd))∩Lq((0,T);Lp(Td;Rd)) where p and q satisfy a Prodi-Serrin condition. The vector field constructed by Depauw [C. R. Acad. Sci. Paris, 2003] lies in this class and we show that for almost every intial datum, the zero-noise limit selects a probability measure concentrated on several distinct integral curves of this vector field.