生成时间: 2026年01月19日 04:06:27

重点关注论文数量: 2

了解领域论文数量: 4

不相关论文数量: 8

重点关注论文（完整分析）

1. 线性化弹性中的结构变形

Structured Deformations in Linearized Elasticity

作者: Manuel Friedrich, José Matias, Elvira Zappale

类别: math.AP, math.FA, math.OC

发布日期: 2026-01-16

arXiv ID: 2601.11333v1

链接: http://arxiv.org/abs/2601.11333v1

详细分析

中文标题: 线性化弹性中的结构变形

以下是根据论文“Structured Deformations in Linearized Elasticity”提供的摘要和PDF内容进行的分析。

1. 研究对象和背景

本文的核心研究对象是线性化弹性理论中的结构变形 (Structured Deformations)。在经典弹性力学中，物体的变形通常由宏观变形场 $g$ 及其梯度 $\nabla g$ 来描述。然而，为了捕捉材料在微观尺度上的复杂行为（如微裂纹、滑移、错位等非光滑变化），Del Piero 和 Owen 引入了“结构变形”的概念。一个结构变形由一对场 $(g, G)$ 描述，其中 $g$ 是宏观位移场，$G$ 是一个附加的几何场。在这个框架下，$G$ 代表了无错位变形的宏观贡献，而 $\nabla g - G$ 则捕捉了非光滑微观变化（即“错位”，disarrangements）的宏观贡献。

本文的背景是将这一理论框架从经典的 Sobolev 空间（如 SBV 空间）扩展到有界变形 (Bounded Deformation, BD) 空间。BD 空间是线性弹性问题的自然函数空间，因为它放宽了对变形梯度的可积性要求（不需要 $\nabla u \in L^p$，只需 $Eu$ 为有界测度），这使得其能够包含更广泛的变形模式，特别是具有奇异性的变形。

数学模型与方程：

研究基于线性化弹性（小变形假设）的变分问题。论文考虑了两种类型的能量泛函：

  1. 定义在 $SBD$ 上的能量泛函：针对位移 $u \in SBD(\Omega)$（有界变形的子集，Cantor 部分为零），能量泛函 $F(u)$ 包括体能量和表面能量：

    $$ F(u) := \int_{\Omega} W\big(x, Eu(x)\big) \,dx + \int_{\Omega \cap J_u} \psi\big(x, [u](x), \nu_u(x)\big) \,dH^{N-1}(x) $$

    其中：

    • $Eu = \frac{\nabla u + (\nabla u)^T}{2}$ 是对称变形梯度。
    • $W(x, \cdot)$ 是定义在 $\mathbb{R}^{N \times N}_{\text{sym}}$ 上的体能量密度。
    • $J_u$ 是位移 $u$ 的跳跃集（断层），$[u]$ 是跳跃高度，$\nu_u$ 是法向量。
    • $\psi$ 是定义在 $\mathbb{R}^N \times S^{N-1}$ 上的表面能量密度。

  2. 定义在结构变形 $(g, G)$ 上的松弛能量 $I_p(g, G)$：为了处理非光滑变形，论文定义了松弛后的能量，即“最经济的方式”逼近给定的结构变形 $(g, G) \in StBD^p(\Omega)$：

    $$ I_p(g, G) := \inf \left\{ \liminf_{n \to \infty} F(u_n) : u_n \in SBD(\Omega), \, u_n \xrightarrow{*} (g, G) \right\} $$

    其中收敛性 $u_n \xrightarrow{*} (g, G)$ 定义为 $u_n \to g$ 强收敛于 $L^1(\Omega; \mathbb{R}^N)$，且 $Eu_n \rightharpoonup^* G$（当 $p=1$）或弱收敛（当 $p>1$）。

研究区域与函数空间：

  • 区域：$\Omega \subset \mathbb{R}^N$ 是具有 Lipschitz 边界的有界域。
  • 函数空间：
    • $BD(\Omega)$ (Bounded Deformation): 分布式对称梯度 $Eu$ 是有界向量测度的空间。
    • $SBD(\Omega)$: $BD(\Omega)$ 中满足 $Eu$ 的 Cantor 部分为零（即 $Eu = Eu^j + Esu$ 且 $Esu = 0$）的子空间。
    • $StBD^p(\Omega) = BD(\Omega) \times L^p(\Omega; \mathbb{R}^{N \times N}_{\text{sym}})$：结构变形的定义空间。

论文旨在建立结构变形理论在 $BD$ 框架下的变分公式，即找到松弛能量 $I_p$ 的积分表示（体密度和面密度），并验证其合理性。

2. 主要定理或主要结果

论文通过以下三个主要定理给出了研究结论：

定理 2.1 (逼近定理 - Approximation Theorem) 这是结构变形理论的基础。对于任意结构变形 $(g, G) \in StBD^p(\Omega)$，存在一列 $SBD(\Omega)$ 函数 $\{u_n\}$ 逼近它：

  • $u_n \to g$ 强收敛于 $L^1(\Omega; \mathbb{R}^N)$。
  • $Eu_n = G$ （在 $L^p$ 意义下，或几乎处处）。
  • $Eu_n \rightharpoonup^* Eg$ 在测度意义下弱*收敛。
  • 对 $|Eu_n|(\Omega)$ 有一致界的控制。

定理 2.8 (松弛与积分表示 - Relaxation) 这是论文的核心结果。在满足 (W1)-(W2) 和 (ψ1)-(ψ2) 等增长和正则性假设下，松弛能量 $I_p(g, G)$ 具有如下积分表示：

$$ I_p(g, G) = \int_{\Omega} f\big(x, Eg(x), G(x)\big) \,dx + \int_{\Omega \cap J_g} \Phi\big(x, [g](x), \nu_g(x)\big) \,dH^{N-1}(x) $$

其中，体密度 $f$ 和面密度 $\Phi$ 由局部极小化问题的渐近极限给出：

  • 体密度 $f(x_0, \xi, B)$:
    $$ f(x_0, \xi, B) := \limsup_{\varepsilon \to 0^+} \frac{m(\ell_{x_0, 0, \xi, B}; Q(x_0, \varepsilon))}{\varepsilon^N} $$
    其中 $\ell_{x_0, 0, \xi, B}$ 是具有线性畸变的仿射函数，$m$ 是相关的极小化问题。
  • 面密度 $\Phi(x_0, \lambda, \nu)$:
    $$ \Phi(x_0, \lambda, \nu) := \limsup_{\varepsilon \to 0^+} \frac{m(v_{x_0, \lambda, 0, \nu}, 0; Q_\nu(x_0, \varepsilon))}{\varepsilon^{N-1}} $$
    其中 $v_{\dots}$ 是在法向量 $\nu$ 方向上发生跳跃 $\lambda$ 的函数。

定理 2.12 (含 Cantor 部分的松弛 - Relaxation with Cantor Part) 当考虑与空间变量 $x$ 无关的能量密度时，定理 2.8 的结果可以进一步扩展，不仅包含体密度和跳跃表面密度，还包含了 Cantor 部分的能量贡献：

$$ I_p(g, G) = \int_{\Omega} H_p(Eg, G) \,dx + \int_{\Omega \cap J_g} h_p([g], \nu_g) \,dH^{N-1} + \int_{\Omega} H_p^{\infty} \left( \frac{dE^c g}{d|E^c g|}, 0 \right) \,d|E^c g| $$

其中 $H_p$ 和 $h_p$ 分别对应于定理 2.8 中确定的密度，$H_p^{\infty}$ 是 $H_p$ 在无穷远处的弱 Recession 函数。

3. 研究方法、关键技术和核心工具

本论文的证明策略结合了抽象的泛函分析工具和具体的变分松弛技巧。

a) 核心分析学工具

  1. 结构变形的逼近理论 (Approximation Theory)： 论文首先建立了 $BD$ 框架下的逼近定理（Theorem 2.1）。这需要利用 Alberti 定理（涉及函数的“好”近似）和 $SBD$ 中函数的构造。在 $BD$ 空间中，由于 $Eu$ 可能是测度，构造 $SBD$ 序列逼近给定的 $(g, G)$ 需要精细的分解技术，特别是处理 $Eu$ 的奇异部分。

  2. 全局松弛方法 (Global Method for Relaxation)： 这是证明积分表示（定理 2.8）的核心工具，由 Bouchitté, Fonseca 等人引入。该方法不直接求解原始的松弛问题，而是引入一个辅助的非局部问题 $m(g, G; O)$（公式 2.7），通过证明 $m$ 与 $I_p$ 的某种对偶性或等价性来导出密度。

    • 关键步骤：利用 Brezis-Lieb 引理 和 测度的弱*收敛 性质。
    • 工具：Banach-Alaoglu 定理 用于保证在 $BD$ 空间中序列的紧性（弱收敛），虽然 $BD$ 不是自反空间，但对称梯度 $Eu$ 是有界测度，因此在弱拓扑下是紧的。

  3. 变分法与 Gamma 收敛 (Variational Methods & Gamma-Convergence)： 论文本质上研究的是 Gamma 收敛（或类似概念）下的松弛极限。为了证明局部极小化问题的解收敛到所需的密度，需要利用：

    • 周期性边界条件与缩放 (Scaling & Periodic Boundary Conditions)：通过在越来越小的立方体 $Q(x, \varepsilon)$ 上求解局部问题并取 $\varepsilon \to 0$ 的极限来获得密度。
    • 紧性原理 (Compactness Principles)：在 $BD$ 空间中，从 $SBD$ 函数序列的 $L^1$ 有界性和对称梯度的 $L^1$ 或测度有界性推导出相对紧性。

b) 技术性难题与解决技巧

  1. 非光滑分析与测度论 (Non-smooth Analysis & Measure Theory)： 难题：$BD$ 空间中的函数其对称梯度 $Eu$ 是一个矩阵值 Radon 测度，可能包含绝对连续部分、跳跃部分和 Cantor 部分。处理 Cantor 部分（非跳跃的奇异部分）非常困难，因为它不能通过传统的切平面或跳跃集来描述。 技巧：

    • 在处理一般 $BD$ 函数时，作者先集中于 $SBD$（Cantor 部分为零）的松弛结果，因为这可以通过相对简单的 $L^p$ 分析和跳跃集的 $H^{N-1}$ 测度来处理。
    • 在定理 2.12 中，专门针对 $x$-无关的能量，利用了 De Giorgi 的 Lower Semicontinuity Theorem 的推广形式。通过识别 Cantor 部分的能量密度为 $H_p^{\infty}$，作者巧妙地绕过了显式计算 Cantor 测度的复杂性，而是将其表示为体能量密度在无穷远处的渐近行为。

  2. 线性化与非线性逼近 (Linearization from Nonlinear Models)： 难题：如何从非线性弹性模型严格推导出线性化弹性中的结构变形能量？ 技巧：论文采用了 Nonsimple Materials (非简单材料) 的模型（源于 Toupin）。这类材料的能量依赖于变形的二阶梯度（$\nabla^2 y$）。

    • 二阶梯度的作用：引入 $\nabla^2 y$ 的项（如 $\int |\nabla^2 y|^2$）提供了额外的紧性（Compactness），这对于从非线性模型（$y = id + \delta u$）过渡到线性模型（$\delta \to 0$）至关重要。它确保了跳跃集的几何控制（即防止断裂在极限下过于奇异）。
    • 缩放分析：论文通过参数 $\delta$（小变形参数）对位移进行缩放，并结合松弛过程。这是一个双重极限过程（$\delta \to 0$ 同时进行结构变形的松弛），需要精细的 Gamma 收敛分析。

  3. Korn-Poincaré 不等式的适应 (Adaptation of Korn-Poincaré Inequality)： 难题：在 $BD$ 空间中，传统的 Korn 不等式（控制位移与其对称梯度的关系）形式不同，因为 $Eu$ 可能没有逆。 技巧：如附录 B 所述，论文在处理线性化非线性模型时，需要证明某些紧性估计。这依赖于针对 $BD$ 空间改进的 Poincaré 不等式，特别是处理跳跃部分 $[u]$ 和体部分 $Eu$ 之间的关系。

  4. 密度识别的具体计算： 技巧：

    • 体密度的识别（定理 2.8 (i)(iii)）利用了定义在 $SBD(Q)$ 上的局部极小化问题。对于 $p>1$，表面能量 $\psi$ 的影响在体密度 $H_p$ 的定义中是显式的；对于 $p=1$，由于 $L^1$ 弱收敛的特性，体密度中包含了表面能的渐近贡献（见 Remark 2.9，涉及 recession 函数 $W^\infty$）。
    • 面密度的识别则关注于在微观尺度上的“错位”能量，即函数在法向量 $\nu$ 方向发生跳跃 $\lambda$ 的最小能量成本。

4. 与之前工作的比较

本文在多个方面推进了现有的结构变形理论和线性化弹性理论。

a) 函数空间的扩展：从 $SBV$ 到 $BD$

  • 先前工作：经典的结构变形理论（如 Choksi & Fonseca, Owen & Paroni 的工作）主要建立在 $SBV$ (Special Bounded Variation) 或 $SBD$ (Special Bounded Deformation) 空间上。这些空间要求变形梯度属于 $L^p$，这在物理上对应于变形梯度在某种意义上是“正则”的，排除了 Cantor 型奇异。
  • 本文的突破：本文将理论推广到了完整的 $BD$ (Bounded Deformation) 空间。$BD$ 空间包含 Cantor 奇异测度，这是处理某些物理现象（如非局部断裂、深层微观结构）所必需的。
    • 技术改进：在 $BD$ 中建立逼近定理（Theorem 2.1）比在 $SBV$ 中更难，因为不能简单地通过 Sobolev 嵌入来控制。作者使用了 Alberti 定理和全局松弛方法，证明了即使在 $Eu$ 为测度的情况下，也可以通过 $SBD$ 序列逼近结构变形。
    • 更一般的收敛性：之前的工作通常要求 $p > 1$ 的 $L^p$ 收敛，而本文在 $BD$ 框架下处理了 $p=1$ 的情况（弱*测度收敛），这对于线性弹性模型是自然的（能量通常是 1 次或 2 次的，1 次对应于 $L^1$ 控制）。

b) 线性化结果的严格推导

  • 先前工作：虽然线性化弹性是标准理论，但将其作为非线性弹性的严格 Gamma 极限（在结构变形的背景下）来推导的工作很少。已有的工作（如 Dal Maso 等人）处理了标准的线性化，但未涉及结构变形中的“错位”部分。
  • 本文的突破：
    • 双重极限过程：本文同时处理了两个极限：一是从非线性（有限变形）到线性（无穷小变形）的极限；二是从微观非光滑变形到宏观结构变形的松弛极限。
    • 非简单材料 (Nonsimple Materials) 的引入：这是相对于标准线性化理论的一个显著改进。标准的线性化往往面临紧性不足的问题，特别是在处理断裂或奇异表面能时。通过引入依赖于二阶梯度的能量项（$\nabla^2 y$），本文建立的模型（Theorem 2.15 的前导模型）具有更好的紧性，从而允许从非线性模型严格推导出包含体密度和面密度的线性化能量。这比仅依赖一阶梯度的线性化模型更稳健，特别是对于涉及“Cantor 部分”贡献的情况。

c) 密度识别的完整性与物理意义

  • 先前工作：在 $SBV$ 框架下的结构变形理论通常只识别了体密度和跳跃表面密度（对应于断裂面）。对于 Cantor 部分的能量贡献，往往是隐含的或未被显式识别。
  • 本文的突破：
    • Cantor 部分的显式化：定理 2.12 显式地给出了 Cantor 部分的能量密度 $H_p^{\infty}$。这揭示了在小变形极限下，即使没有宏观断裂（跳跃集 $J_g = \emptyset$），微观的 Cantor 型错位（如微观晶格错位）也会对宏观能量产生贡献。这在物理上对应于材料的“织构”或微观缺陷演化。
    • 面密度的解耦与重耦：对于 $p>1$，本文证明了表面密度 $h_p$ 可以从体能量 $W$ 中解耦定义（定理 2.8 (iii)），这意味着表面能主要由跳跃 $\psi$ 决定。而对于 $p=1$（线性弹性中最常见的生长条件），体密度和表面密度通过 recession 函数 $W^\infty$ 相互作用（Remark 2.9），这比之前的理论更精确地捕捉了线性化极限下的能量行为。

d) 技术方法的推广

  • 全局松弛方法在 $BD$ 的应用：此前全局松弛方法多用于 $BV$ 或 $SBV$。本文证明了该方法在 $BD$ 空间（非自反空间，具有测度值导数）同样有效，并且可以处理结构变形这种双重场 $(g, G)$ 的情形。这为处理更广泛的非光滑分析问题提供了通用框架。
  • Korn-Poincaré 不等式的适用性：论文在附录中处理了 $BD$ 空间的 Poincaré 不等式，这对于证明线性化极限的紧性至关重要。相比于传统弹性力学中标准的 Korn 不等式，这里需要适应包含跳跃项和 Cantor 项的测度形式，这是对经典不等式的重要推广。

综上所述，本文通过将结构变形理论扩展到 $BD$ 空间，并结合非简单材料的模型，严格推导了线性化弹性中的松弛能量。其主要贡献在于提供了更广泛的函数空间框架、显式识别了 Cantor 部分的能量密度，以及建立了非线性到线性化在结构变形背景下的严格联系。

2. 低正则性折射率反谱问题

Inverse Spectral Problem With Low Regularity Refractive Index

作者: Kewen Bu, Youjun Deng, Yan Jiang, Kai Zhang

类别: math.AP

发布日期: 2026-01-16

arXiv ID: 2601.11146v1

链接: http://arxiv.org/abs/2601.11146v1

详细分析

中文标题: 低正则性折射率反谱问题

研究对象和背景

本文研究的是一个反谱问题，其具体背景是内传输问题（Interior Transmission Problem）。该问题源于电磁波或声波在介质中的散射，目标是通过频谱数据（传输本征值）来唯一地重构或确定介质内部的折射率（refractive index）函数 $n$。

在标准的波散射模型中，设 $D$ 是 $\mathbb{R}^d$ ($d=2,3$) 中的一个有界单连通区域，边界 $\partial D$ 分段光滑。设 $n(x)$ 为折射率，满足 $n(x)=1$ 当 $x \notin D$。$k$ 为波数（与频率相关），$u$ 为总波场，$u_s$ 为散射波场。平面波入射散射问题可表示为：

$$ \begin{cases} \Delta u + k^2 n(x) u = 0, & x \in \mathbb{R}^d, \\ u(x) = e^{ikx \cdot l} + u_s(x), \\ \lim_{r \to \infty} r (\frac{\partial u_s}{\partial r} - i k u_s) = 0, \end{cases} $$

其中 $l$ 为单位入射向量。

内传输问题与上述散射问题的远场算子的核密切相关。在采样方法（如线性采样法）中，若入射波位于远场算子的核中，则不产生散射波，因此无法直接通过远场模式重构介质。内传输问题为研究这一情况提供了有力工具。

内传输问题的具体数学模型（在本文中针对 $D = B(0,1)$ 单位球）如下：寻找非平凡的 $u, v$ 满足：

$$ \begin{cases} \Delta u + k^2 n(x) u = 0, & x \in D, \\ \Delta v + k^2 v = 0, & x \in D, \\ u = v, & \text{on } \partial D, \\ \frac{\partial u}{\partial \nu} = \frac{\partial v}{\partial \nu}, & \text{on } \partial D. \end{cases} $$

使得存在非零解的 $k$ 值被称为传输本征值（Transmission Eigenvalues）。传输本征值是介质的指纹，理论上可以通过测量获得。

本文研究的区域 $D$ 是单位球 $B(0,1)$，且假设折射率 $n$ 是径向函数，即 $n(x) = n(r)$，其中 $r = |x|$。在此径向对称假设下，我们可以寻找径向对称的本征函数 $u$ 和 $v$，形式为 $u = y(r)/r$ 和 $v = y_0(r)/r$。代入方程组，得到两个常微分方程（ODE）系统：

$$ \begin{cases} y'' + k^2 n(r) y = 0, & r \in (0,1), \\ y_0'' + k^2 y_0 = 0, & r \in (0,1), \\ y(1) = y_0(1), \\ y'(1) = y_0'(1). \end{cases} $$

这也就是文中提到的方程组 (2.2)。

函数空间与正则性假设： 本文的核心贡献之一在于放宽了对折射率 $n$ 的正则性要求。传统反谱问题通常假设 $n$ 是 $C^2$ 或 $L^\infty$ 的。本文考虑了以下两类折射率：

  1. 低正则性：$n \in C^2[0,1]$（二阶连续可微）或 $C^{1,1}[0,1]$（一阶导数 Lipschitz 连续），且满足 $0 < n_* < n(r) < n^*$。
  2. 分段结构：允许 $n$ 具有分段定义的结构，记为 $n \in M$。集合 $M$ 定义为分段连续函数，且对于任意 $f_1, f_2 \in M$，存在有限分割 $0 = x_0 < x_1 < \dots < x_S = 1$，使得在每个区间 $(x_s, x_{s+1})$ 内，Heaviside 阶跃函数 $H(f_1 - f_2)$ 保持常数。这意味着两个 $M$ 类函数之差的符号改变次数是有限的。

主要定理或主要结果

本文的主要结果由两个定理构成，分别处理折射率 $n$ 在 $C^2$ 和 $C^{1,1}$ 正则性下的唯一性确定问题。证明唯一性需要结合特殊的传输本征值（对应径向对称本征函数的本征值）以及部分先验信息（即 $n$ 在某个子区间上的值）。

定理 2.1 (Theorem 2.1): 假设折射率 $n$ 具有 $L$ 层结构（见式 (2.4)），满足 $n \in C^2[0,1]$ 且 $n \in M$。设 $0 < n_* < n(r) < n^*$， $|n'(r)| < \tilde{n}^*$， $|n''(r)| < \tilde{n}^{**}$。 已知参数 $\{\hat{\delta}_l\}_{l=1}^L$（定义见原文，与层厚度和折射率积分有关）以及 $n$ 在区间 $[\alpha, 1]$ 上的值（其中 $\alpha$ 由 $\int_\alpha^1 \sqrt{n(t)} dt$ 的特定方程定义，或者 $\alpha=1$），且已知 $n(1) \neq 1$。 如果满足以下条件之一：

  1. $\tilde{\beta}_1$ 与 $\tilde{\beta}_{2L}$ 在有理数域上线性独立，且 $\varepsilon_0$ 足够小（满足特定不等式）；
  2. $\tilde{\beta}_1$ 与 $\tilde{\beta}_{2L}$ 在有理数域上线性相关，且满足特定的 $\delta_L$（总光程）与 $n_L(1)$ 的关系及 $\varepsilon$ 的界限（具体条件见原文 (a)(b)）；

那么，特殊的传输本征值集合唯一地确定了折射率 $n$。

定理 2.2 (Theorem 2.2): 假设折射率 $n$ 满足 $n \in C^{1,1}[0,1] \cap M$，且满足 $0 < n_* < n(r) < n^*$， $|n'(r)| < n^*$。 已知 $n$ 在区间 $[\alpha, 1]$ 上的轮廓，且满足边界条件 $n(1) \neq 1$ 和 $n'(1) = 0$。 在此条件下，特殊的传输本征值唯一地确定了 $n$。

推论/反例 (Remark 2.3): 对于仅分段连续（如分段常数）的折射率 $n$，仅凭特殊的传输本征值无法唯一确定 $n$。这说明了正则性的降低会增加确定 $n$ 所需的信息量。

研究方法、关键技术和核心工具

本文的证明融合了常微分方程理论、谱理论、李维尔变换（Liouville Transform）以及数论中的丢番图逼近（Diophantine Approximation）。

a) 核心分析学工具

  1. 李维尔变换 (Liouville Transform): 这是处理径向 Schrödinger 型方程 $y'' + k^2 n(r) y = 0$ 的标准工具。通过变量代换 $\xi(r) = \int_0^r \sqrt{n(t)} dt$ 和函数缩放 $z(\xi) = [n(r)]^{1/4} y(r)$，将变系数方程转化为具有势函数的方程：

    $$ z'' + (k^2 - p(\xi)) z = 0 $$

    其中势函数 $p(\xi)$ 由 $n$ 的导数构成（见式 (3.7)）。这一变换消除了方程的主要部分中的变系数，使得后续的渐近分析和估计更为标准。

  2. 常微分方程的渐近分析 (Asymptotic Analysis of ODEs): 为了利用本征值谱确定 $n$，需要分析特征函数 $y(1)$ 和 $y'(1)$ 以及特征行列式 $d(k)$ 的渐近行为（当 $k \to \infty$ 时）。文章通过李维尔变换后的方程结构，结合引理 3.9，得到了 $z$ 的渐近展开式（正弦/余弦项加上高阶小量），进而推导出 $y(1)$ 和 $y'(1)$ 的主导项（Lemma 4.2）。这允许作者将特征行列式 $d(k)$ 表示为一系列振荡项的和，其频率由光程 $\delta \hat{}$ 决定。

  3. Sturm-Liouville 理论与解的先验估计: 文章首先建立了低正则性下解的 $L^2$ 估计（Lemma 4.1）。利用归纳法和分层介质的匹配条件，证明了解 $y(r)$ 及其导数 $y'(r)$ 的增长受 $|k|$ 和积分 $\int_0^r \sqrt{n(t)} dt$ 的控制。这些估计对于控制后续的积分余项至关重要。

  4. 丢番图逼近 (Diophantine Approximation): 这是处理反演问题中“唯一性”与“非唯一性”边界的关键数论工具。利用 Dirichlet 逼近定理 的某种形式（文中引用的 Lemma 2.1），作者处理了特征频率（光程 $\beta$）的有理线性依赖性。 在定理 2.1 的条件中，区分了 $\tilde{\beta}_1$ 和 $\tilde{\beta}_{2L}$ 是否有理线性相关。如果它们是无理相关（独立），则可以通过本征值序列的分布唯一确定相位；如果存在有理关系，则需要更强的先验信息（如 $\varepsilon$ 的界限）来消除多解性。这实际上是利用频率的非共振性质来锁定物理参数。

b) 解决技术性难题的主要技巧

  1. 处理低正则性 ($C^{1,1}$) 的技巧: 对于 $C^{1,1}$ 的折射率，二阶导数 $n''$ 不存在（仅在分布意义下存在），因此标准的 Liouville 变换公式（包含 $n''$ 项）不再适用。作者采用的技巧是：

    • 逼近论证： 在证明 Theorem 2.2 时，利用 $C^{1,1}$ 函数可以用 $C^\infty$ 函数序列一致逼近的性质（通过卷积或磨光算子）。
    • 极限过程： 构造一列光滑折射率 $n_\epsilon$，它们的本征值收敛于 $C^{1,1}$ 折射率的本征值。通过在 $C^2$ 框架下证明 $n_\epsilon$ 的唯一性，然后取极限得到 $C^{1,1}$ 情形的唯一性。这种“先光滑后逼近”的方法避免了直接处理 $n''$ 不存在带来的技术困难。

  2. 构造非唯一性反例 (Counterexample): 为了证明仅凭特殊传输本征值不足以确定低正则性 $n$（Section 3.1），作者构造了一个精巧的分段常数反例（$n^{(1)}$ 和 $n^{(2)}$）。

    • 技巧： 选取两组不同的折射率参数（如 $n=4,16$ 和 $n=16,4$，切换层的位置），通过计算 $d(k)$ 的显式表达式（式 3.3），证明尽管 $n^{(1)} \neq n^{(2)}$，但它们的特征行列式 $d(k)$ 仅相差一个常数因子，因此本征值零点完全相同。
    • 意义： 这个反例严格界定了正则性要求的必要性：当 $n$ 的正则性降低到分段常数时，频谱信息丢失了层位置的信息（仅保留了某种加权平均）。

  3. 积分恒等式与 Wronskian: 在证明唯一性时，核心技巧是利用两个不同折射率 $n^{(1)}, n^{(2)}$ 对应的解 $y_1, y_2$ 建立积分恒等式。

    • 通过将两个方程分别乘以 $y_2, y_1$ 并相减，得到 Wronskian 的导数表达式。
    • 积分得到（见 Section 3.3）：
      $$ \int_0^\alpha k^2 (n^{(2)} - n^{(1)}) y_1 y_2 dr = [y_1'y_2 - y_1 y_2']|_{r=1} $$
    • 由于在 $[\alpha, 1]$ 上已知 $n^{(1)} = n^{(2)}$（先验信息），积分上限截断在 $\alpha$。右边的边界项由 $d(k)$ 决定。通过分析该积分在 $k \to \infty$ 时的行为，结合 $d(k)$ 的渐近展开，可以推导出 $n^{(2)} - n^{(1)}$ 必须为零。

与之前工作的比较

本文在反谱问题的唯一性理论方面做出了显著的改进，主要体现在降低了正则性要求和引入了新的先验信息模式。

a) 放宽假设条件：从 $C^2$ 到 $C^{1,1}$ 及分段 $C^2$

  1. 正则性的降低:

    • 先前工作： 经典结果（如 Colton, Kress, McLaughlin, Polyakov 等人的工作）通常要求折射率 $n$ 具有较高的光滑性，例如 $n \in C^\infty$ 或至少 $C^2[0,1]$。Gintides (2010) 的工作虽然允许了分段 $C^2$（两层介质），但仍然依赖于 $C^2$ 的局部性质。
    • 本文改进： 本文证明了对于 $C^{1,1}$（一阶导数 Lipschitz 连续，二阶导数可能不存在或仅为 $L^\infty$）的折射率，唯一性依然成立（Theorem 2.2）。这在物理上更有意义，因为实际介质的折射率变化往往是有限平滑的，未必处处二阶可导。文中通过将 $C^{1,1}$ 视为 $C^2$ 的逼近极限，成功将唯一性理论推广到了弱正则性情形。

  2. 分层结构的推广:

    • 先前工作： 针对低正则性的研究（如 Gintides）主要局限于两层介质（Two-layer medium），这使得分界面的数量和位置受到严格限制。
    • 本文改进： 本文引入了集合 $M$，允许任意有限层（$L$ 层）的分段连续结构（Theorem 2.1）。这使得模型能描述更复杂的多层涂层或层状地质结构。虽然仍需要“小跳跃”条件（即 $\varepsilon$ 较小，层间性质接近），但这相比于之前仅处理单一分界面的工作是一个巨大的推广。

b) 信息需求的优化：频谱完备性 vs. 部分先验信息

  1. 从完全谱到部分谱:

    • 先前工作： 大多数唯一性结果（如 Cakoni 2006, Gintides 2010）要求知道所有传输本征值（或无穷子集）。在实际中，获取完整的谱是非常困难的。
    • 本文改进： 本文明确指出，仅凭特殊的传输本征值（对应径向对称模式的子集）加上部分先验信息（$n$ 在 $[\alpha, 1]$ 上的值），即可唯一确定 $n$。这是一种“互补信息”的策略：如果谱信息较弱（只有特殊本征值），则需要少量的先验信息；反之，如果先验信息少，则需要更完备的谱。

  2. 先验信息的具体形式:

    • 先前工作： 传统的唯一性证明往往依赖较强的全局正则性假设（如 $C^2$）来通过 Liouville 变换的反演公式直接重构 $n$。如果正则性降低，Liouville 变换的反演公式失效（因为 $n''$ 项发散）。
    • 本文改进： 本文通过引入 $[\alpha, 1]$ 区间上的已知值作为锚点，利用 Wronskian 积分恒等式。这个积分恒等式不需要对 $n''$ 进行显式反演，只需利用 $n$ 的有界性和连续性。这表明，在低正则性情况下，局部的先验信息可以弥补全局光滑性的不足。

  3. 折射率范围的灵活性:

    • 先前工作： 很多文献假设 $0 < n < 1$（低于背景）或 $n > 1$（高于背景）。
    • 本文改进： 本文利用丢番图逼近，仅需假设 $0 < n_* < n(r) < n^*$，允许 $n$ 在背景值 1 上下波动（只要波动范围有界），这更符合实际介质的复杂性。

总结： 这篇文章通过结合微分方程的渐近分析、数论逼近和逼近论方法，显著降低了传输本征值反问题中对折射率正则性的要求，并阐明了在低正则性下，局部先验信息与频谱信息之间的精确互补关系。

了解领域论文（摘要翻译）

1. 基于三维Koiter线性弹性的均匀化中等皱褶壳理论

Homogenized moderately wrinkled shell theory from 3D Koiter’s linear elasticity

作者: Pedro Hernández-Llanos, Rajesh Mahadevan, Ravi Prakash

类别: math.AP

发布日期: 2026-01-16

arXiv ID: 2601.11384v1

链接: http://arxiv.org/abs/2601.11384v1

摘要翻译

中文标题: 基于三维Koiter线性弹性的均匀化中等皱褶壳理论

摘要翻译: 本文通过双尺度收敛方法，从三维线性弹性理论出发，推导出周期性皱褶壳模型。根据小参数 $\varepsilon>0$ 和 $p>1$ 的行为不同，会出现不同的理论体系。我们假设壳的中面由 $\displaystyle ψ(x_1,x_2)+\varepsilon^pθ\left(\frac{x_1}{\varepsilon},\frac{x_2}{\varepsilon}\right)\vect{a}_{3}(x_1,x_2)$ 给出，其中 $θ$ 是 $[0,1)^2$ 周期函数，且 $p=2$。同时假设壳的应变能服从Koiter模型。

2. 一种用于模拟原生细胞生长的活性液滴Mullins-Sekerka模型：稳定性分析与数值计算

On a Mullins-Sekerka model for the growth of active droplets modelling protocells: Stability analysis and numerical computations

作者: Harald Garcke, Kei Fong Lam, Robert Nürnberg, Andrea Signori

类别: math.AP

发布日期: 2026-01-16

arXiv ID: 2601.11155v1

链接: http://arxiv.org/abs/2601.11155v1

摘要翻译

中文标题: 一种用于模拟原生细胞生长的活性液滴Mullins-Sekerka模型：稳定性分析与数值计算

摘要翻译: 具有化学反应的Mullins-Sekerka模型可能导致液滴生长、变得不稳定、分裂、生长并经历进一步分裂的情景。这些生长和分裂循环已被提出作为原生细胞的模型，并被认为在活体系统中扮演着基础性角色，通过提供对活体系统组织至关重要的化学隔室。本文分析了化学活性的Mullins-Sekerka模型。证明了径向对称解的存在性，并给出了径向和平面情况下的详细稳定性分析。特别是，我们还分析了导致壳层型情况的多层解。最后，我们介绍了一种基于参数有限元方法的数值方法，该方法明确考虑了拓扑变化，从而允许液滴分裂和合并。多项数值模拟验证了理论稳定性分析的发现，并展示了复杂的动力学行为，包括多重不稳定性、液滴分裂和壳层型解的出现。

3. 一种基于低秩逼近与Neumann矩阵级数求解随机系数非稳态扩散型偏微分方程的高效求解器

An efficient solver based on low-rank approximation and Neumann matrix series for unsteady diffusion-type partial differential equations with random coefficients

作者: Yujun Zhu, Min Li, Yulan Ning, Ju Ming

类别: math.NA, math.AP

发布日期: 2026-01-16

arXiv ID: 2601.11152v1

链接: http://arxiv.org/abs/2601.11152v1

摘要翻译

中文标题: 一种基于低秩逼近与Neumann矩阵级数求解随机系数非稳态扩散型偏微分方程的高效求解器

摘要翻译: 本文针对具有随机系数的非稳态扩散型偏微分方程，开发了一种高效数值求解器。此类问题的主要计算挑战在于，在不确定性条件下，需反复处理由空间和时域离散化产生的大规模线性系统。为解决这一问题，我们提出了一种新颖的广义低秩矩阵逼近方法，用于表示随机刚度矩阵，并利用Neumann矩阵级数展开逼近其逆矩阵。该方法将高维矩阵求逆转化为一系列低维矩阵乘法。因此，该求解器在保持高数值精度的同时，显著降低了计算成本和存储需求。论文还提供了所提求解器的误差分析。最后，我们将该方法应用于两个经典的不确定性量化问题：非稳态随机扩散方程及相关分布式最优控制问题。数值结果证明了所提求解器的可行性与有效性。

4. 发展方程框架下非局部积分算子的数值处理

Numerical Treatment of Non-local Integral Operators in the Framework of Evolutionary Equations

作者: Sebastian Franz, Sascha Trostorff

类别: math.NA, math.AP

发布日期: 2026-01-16

arXiv ID: 2601.11132v1

链接: http://arxiv.org/abs/2601.11132v1

摘要翻译

中文标题: 发展方程框架下非局部积分算子的数值处理

摘要翻译: 本文运用发展方程理论，研究包含非局部积分算子的抽象微分方程。在给出该方程适定性条件后，我们提出了一种求解其近似解的数值方法。在积分算子的核及解满足特定条件的情况下，我们给出了收敛性证明，并通过若干模拟结果完成本文论述。

不相关论文（基本信息）

1. 基于共振硬夹杂的弹性Calderón问题：N维映射的线性化与密度重建

Elastic Calderón Problem via Resonant Hard Inclusions: Linearisation of the N-D Map and Density Reconstruction

作者: Huaian Diao, Mourad Sini, Ruixiang Tang

类别: math.AP, math-ph

发布日期: 2026-01-16

arXiv ID: 2601.11356v1

链接: http://arxiv.org/abs/2601.11356v1

摘要: We study an elastic Calderon-type inverse problem: recover the mass density $ρ(x)$ in a bounded domain $Ω\subset\mathbb{R}^3$ from the Neumann-to-Dirichlet map associated with the isotropic Lamé system $\mathcal{L}_{λ,μ}u+ω^2ρ(x)u=0$. We introduce a constructive strategy that embeds a subwavelength periodic array of resonant high-density (hard) inclusions to create an effective medium with a uniform negative density shift. Specifically, we place a periodic cluster of inclusions of size $a$ and density $ρ_1\asymp a^{-2}$ strictly inside $Ω$. For frequencies $ω$ tuned to an eigenvalue of the elastic Newton (Kelvin) operator of a single inclusion, we show that as $a\to0$ and the number of inclusions $M\to\infty$, the Neumann-to-Dirichlet map $Λ_D$ converges to an effective map $Λ_{\mathcal{P}}$ corresponding to a background density shift $-\mathcal{P}^2$, with the operator norm estimate $\|Λ_D-Λ_{\mathcal{P}}\|\le Ca^α\mathcal{P}^6$ for some $α>0$ determined by the geometric scaling. Around this negative background we derive a first-order linearization of $Λ_{\mathcal{P}}$ in terms of $ρ$ and the Newton volume potential for the shifted Lamé operator. Testing the linearized relation with complex geometric optics solutions yields an explicit reconstruction formula for the Fourier transform of $ρ$, and hence a global density recovery scheme. The results provide a metamaterial-inspired analytic framework for inverse coefficient problems in linear elasticity and a concrete paradigm for leveraging nanoscale resonators in reconstruction algorithms.

2. 有界区域上分数阶Hamilton-Jacobi方程的遍历对：大解

Ergodic pairs for fractional Hamilton-Jacobi equations on bounded domains: large solutions

作者: Alexander Quaas, Erwin Topp

类别: math.AP

发布日期: 2026-01-16

arXiv ID: 2601.11241v1

链接: http://arxiv.org/abs/2601.11241v1

摘要: In this article, we study the ergodic problem associated to viscous Hamilton-Jacobi equation where the diffusion is governed by the censored fractional Laplacian, a nonlocal elliptic operator restricted to a bounded domain $Ω\subset \mathbb{R}^N$. We restrict ourselves to the case in which the nonlinear gradient term has a scaling less or equal than the fractional order of the diffusion. In similarity to its second-order counterpart, we provide existence of ergodic pairs involving solutions that blow-up on $\partial Ω$. We use the celebrated vanishing discount method, where the analysis of the approximated solutions have its own interest, leading to qualitative properties for the ergodic problem such as precise blow-up rates for the solution and characterization of the ergodic constant. The main difficulties arise from the state-dependency of the operator, from which the arguments of the local case based on well-known invariance properties of the Laplacian are not longer at disposal.

3. 完备非紧黎曼流形上的最优Trudinger-Moser不等式：从局部不等式到全局不等式的论证再探

Optimal Trudinger-Moser inequalities on complete noncompact Riemannian manifolds: Revisit of the argument from the local inequalities to global ones

作者: Jungang Li, Guozhen Lu

类别: math.AP

发布日期: 2026-01-16

arXiv ID: 2601.10996v1

链接: http://arxiv.org/abs/2601.10996v1

摘要: The main purpose of this short note, on the one hand, to is clarify some part of the proof of Theorem 1.3 in [8] in a simple way, and on the other hand, to give an alternative argument from local inequalities to global ones.

4. 不定变分问题中 Nehari-Pankov 型归一化解

Normalized solutions of Nehari-Pankov type to indefinite variational problems

作者: Damien Galant, Tobias Weth

类别: math.AP

发布日期: 2026-01-16

arXiv ID: 2601.10941v1

链接: http://arxiv.org/abs/2601.10941v1

摘要: We consider abstract nonlinear equations of the form $A u = λu + I'(u)$, where $A$ is a self-adjoint operator with compact resolvent on a Hilbert space $H$, $λ\in \mathbf{R}$ is a parameter, and $u \mapsto I'(u)$ is a superlinear term of variational nature. In this abstract setting, we develop a new approach to existence and multiplicity of solutions with prescribed norm in $H$. We then consider various applications of this approach. First, we obtain, under fairly general assumptions including the mass-supercritical case, the existence of infinitely many solutions to a class of nonlinear (time-independent) Schrödinger equations on a compact graph $\mathcal{G}$ with prescribed (arbitrarily large) mass. Moreover, we derive a similar existence and multiplicity result for a biharmonic semilinear equation in the $2$-torus. For a larger class of second order and higher order equations in a bounded domain with Dirichlet boundary conditions, we also show the existence of multiple solutions with prescribed small mass. The solutions we obtain are detected as ground states of Nehari-Pankov type for the associated $λ$-dependent action functional, where $λ$ varies in a spectral gap between sufficiently large eigenvalues of $A$. The key observation in this abstract framework is the fact that the $H$-norms of these $λ$-dependent solution families form connected sets. To estimate the size of these connected sets in specific settings, we use Weyl type estimates for the length of spectral gaps, variational characterizations of eigenvalues, bounds for associated eigenfunctions and a classical bound from analytic number theory.

5. 二维 Ambrosio-Tortorelli 泛函在相场能量收敛条件下的临界点研究

Critical points of the two-dimensional Ambrosio-Tortorelli functional with convergence of the phase-field energy

作者: Jean-François Babadjian, Martin Rakovsky, Rémy Rodiac

类别: math.AP

发布日期: 2026-01-15

arXiv ID: 2601.10875v1

链接: http://arxiv.org/abs/2601.10875v1

摘要: We consider a family $\{(u_\varepsilon, v_\varepsilon)\}_{\varepsilon>0}$ of critical points of the Ambrosio-Tortorelli functional. Assuming a uniform energy bound, the sequence $\{(u_\varepsilon, v_\varepsilon)\}_{\varepsilon>0}$ converges in $L^2(Ω)$ to a limit $(u, 1)$ as $\varepsilon \to 0$, where $u$ is in $SBV^2(Ω)$. It was previously shown that if the full Ambrosio-Tortorelli energy associated to $(u_\varepsilon,v_\varepsilon)$ converges to the Mumford-Shah energy of $u$, then the first inner variation converges as well. In particular, $u$ is a critical point of the Mumford-Shah functional in the sense of inner variations. In this work, focusing on the two-dimensional setting, we extend this result under the sole convergence of the phase-field energy to the length energy term in the Mumford-Shah functional.

6. 基于 Robin-to-Dirichlet 数据的 EIT 内部 Robin 界面定性重构方法

Qualitative reconstruction methods for imaging interior Robin interfaces in EIT from Robin-to-Dirichlet data

作者: Rafael Ceja Ayala, Malena I. Español, Govanni Granados

类别: math.NA, math.AP

发布日期: 2026-01-15

arXiv ID: 2601.10839v1

链接: http://arxiv.org/abs/2601.10839v1

摘要: We consider an inverse shape problem arising in electrical impedance tomography (EIT) for nondestructive testing, in which interior defects are modeled through Robin transmission conditions. Unlike classical formulations, we impose Robin boundary conditions on both the exterior measurement surface and the interior interface, and use the Robin-to-Dirichlet (RtD) map as the available data. Within this setting, we develop qualitative (non-iterative) reconstruction methods based on the Linear Sampling Method (LSM) and the Regularized Factorization Method (RFM), and derive new analytical characterizations that enable these methods to identify interior regions. We further propose a numerical implementation that incorporates regularization strategies and demonstrate, through experiments, that the methods reliably reconstruct interior regions of interest.

7. 无界符号、热流与Toeplitz算子

Unbounded symbols, heat flow, and Toeplitz operators

作者: Sam Looi

类别: math.FA, math.AP, math.CV

发布日期: 2026-01-15

arXiv ID: 2601.10711v1

链接: http://arxiv.org/abs/2601.10711v1

摘要: We disprove the natural domain extension of the Berger–Coburn heat-flow conjecture for Toeplitz operators on the Bargmann space and identify the failure mechanism as a gap between pointwise and uniform control of a Gaussian averaging of the squared modulus of the symbol, a gap that is invisible to the linear form $T_g$. We establish that the form-defined operator $T_g$ and the natural-domain operator $U_g$ decouple in the unbounded symbols regime: while $T_g$ is governed by linear averaging, $U_g$ is controlled by the quadratic intensity of $|g|^2$. We construct a smooth, nonnegative radial symbol $g$ satisfying the coherent-state admissibility hypothesis with bounded heat transforms for all time $t>0$; for this symbol, $T_g$ is bounded, yet $U_g$ is unbounded. This is a strictly global phenomenon: under the coherent-state hypothesis, local singularities are insufficient to cause unboundedness, leaving the ``geometry at infinity’’ as the sole obstruction. Boundedness of $U_g$ is equivalent to the condition that $|g|^2 dμ$ is a Fock–Carleson measure, a condition strictly stronger than the linear average $g dμ$ governing $T_g$. Finally, regarding the gap between the known sub-critical sufficiency condition and the critical heat time, we prove that heat-flow regularity is irreversible in this context and show that bootstrapping strategies cannot resolve the gap between sufficiency and critical time.

8. 通过多线性Kakeya型界给出的接近面积的各向异性能量的Michael-Simon不等式

Michael-Simon inequality for anisotropic energies close to the area via multilinear Kakeya-type bounds

作者: Guido De Philippis, Alessandro Pigati

类别: math.AP, math.DG

发布日期: 2026-01-15

arXiv ID: 2601.10647v1

链接: http://arxiv.org/abs/2601.10647v1

摘要: Given an anisotropic integrand $F:\text{Gr}_k(\mathbb R^n)\to(0,\infty)$, we can generalize the classical isotropic area by looking at the functional

$$ \mathcal{F}(Σ^k):=\int_ΣF(T_xΣ)\,d\mathcal{H}^k. $$

While a monotonicity formula is not available for critical points, when $k=2$ and $n=3$ we show that the Michael-Simon inequality holds if $F$ is convex and close to $1$ (in $C^1$), meaning that $\mathcal{F}$ is close to the usual area. Our argument is partly based on some key ideas of Almgren, who proved this result in an unpublished manuscript, but we largely simplify his original proof by showing a new functional inequality for vector fields on the plane, which can be seen as a quantitative version of Alberti’s rank-one theorem. As another byproduct, we also show Michael-Simon for another class of integrands which includes the $\ell^p$ norms for $p\in(1,\infty)$. For a general $F$ satisfying the atomic condition, we also show that the validity of Michael-Simon is equivalent to compactness of rectifiable varifolds.