Contact set of solutions to gradient flow of Landau-de Gennes energy with a singular entropy potential

技术摘要：具有奇异熵势的 Landau-de Gennes 能量梯度流的接触集解

问题陈述
本文研究了向列型液晶 Landau-de Gennes 能量梯度流解的正则性与几何结构。该模型使用矩阵值序参数 $Q$（Q-张量）在物理定义域 $\mathbb{Q}_{phy} \subset \mathbb{R}^{3\times3}$ 内演化，该定义域由迹为零且特征值严格限制在 $(-1/3, 2/3)$ 之间的对称矩阵组成。能量泛函包含各向异性弹性能量 $G(\nabla Q)$ 和由奇异熵势 $f(Q)$ 驱动的体能量 $B(Q)$。该势函数 $f(Q)$ 由前人工作引入，在物理定义域内是有限的，但当 $Q$ 趋于物理边界 $\partial \mathbb{Q}_{phy}$ 时会趋于 $+\infty$，从而强制执行物理约束。

主要目标是分析“接触集” $C(Q) = \{x \mid f(Q(x)) = \infty\}$，它代表了解违反物理约束（即特征值达到奇异极限）的空间位置。虽然之前的研究已经建立了全局解的存在性，但在动态设定下，使奇异势发生爆炸的集合的大小和结构仍是一个悬而未决的问题，特别是关于解是否在几乎处处都严格保持在物理定义域内。

方法论
作者结合凸分析、容量理论和球面平均技术来解决这一问题。

1. 对偶形式与强制性：
   论文首先为奇异势 $f(Q)$ 建立了严谨的对偶形式。通过使用拉格朗日乘子和偏移特征值对问题进行重参数化，作者推导出了 $f(Q)$ 在物理边界附近的下界。至关重要的是，他们提供了一个简化的“类自共轭”（self-concordant-like）强制性不等式的证明：
   $$A : D^2_Q f(Q) A \geq C (D_Q f(Q) : A)^2$$
   对于所有 $Q \in \mathbb{Q}_{phy}$ 和测试矩阵 $A$ 均成立。该不等式此前是通过复杂的积分计算建立的，本文则利用配分函数的 Legendre 变换性质证明了该不等式，展示了势函数的 Hessian 能够控制势函数梯度的平方。

2. 存在性与正则性（定理 1.1）：
   为了处理奇异性，作者结合 Moreau-Yosida 正则化方案与 Galerkin 近似构造了全局弱解。他们推导出了此前无法获得的更高阶能量估计。具体而言，他们证明了对于具有有限能量和变分的初始数据，解 $Q(t, x)$ 满足：

   • $Q \in L^\infty(0, T; H^2(\mathbb{T}^n))$ 且 $\partial_t Q \in L^\infty(0, T; L^2) \cap L^2(0, T; H^1)$。
   • 关键正则性： 奇异势本身具有 $H^1$ 正则性，即 $f(Q) \in L^2(0, T; H^1(\mathbb{T}^n))$。这是通过结合强制性不等式与高阶能量估计得出的，证明了 $f$ 的 Hessian 对 $Q$ 和 $\partial_t Q$ 梯度的积分是有限的。

3. 接触集分析（定理 1.2）：
   本文的核心贡献在于仅基于 $f(Q)$ 的 $H^1$ 正则性和 $Q$ 的 $H^2$ 正则性来表征接触集 $C(Q)$。作者利用球面平均和容量理论来分析 $f(Q)$ 在奇异点附近的爆炸速率。

   • 2D 情况： 他们表明，如果存在接触点，则 $f(Q)$ 的对数爆炸速率（由势函数的下界导出）将与 $f(Q)$ 梯度的 $H^1$ 界在二维空间中不相容。具体而言，$H^1$ 能量界限限制了 $f(Q)$ 球面平均值的增长为 $O(\sqrt{|\log r|})$，而奇异性要求为 $O(|\log r|)$。这种矛盾意味着接触集为空集。
   • 3D 情况： 在三维空间中，$H^1$ 嵌入并不能完全排除奇异性的存在。然而，通过定义 $Q$ 和 $\nabla Q$ 的“良好行为”勒贝格点集，并应用关于梯度分数密度的引理 4.1，作者证明了接触集的 Hausdorff 维数至多为 1。

主要结果

• 定理 1.1： 确立了具有 $H^2$ 空间正则性的梯度流的唯一全局解的存在性，并且关键性地确定了奇异势 $f(Q)$ 的 $H^1$ 正则性。
• 定理 1.2： 对于任何满足 $f(Q) \in H^1$ 的 $H^2$ 函数 $Q$，证明了：
  • 在 2D 中，接触集 $C(Q)$ 为空 ($C(Q) = \emptyset$)。
  • 在 3D 中，接触集的 Hausdorff 维数至多为 1 ($\mathcal{H}^{1+\epsilon}(C(Q)) = 0$，对于任何 $\epsilon > 0$)。
• 推论 1.1： 将定理 1.2 应用于梯度流解，确认了对于几乎所有时间 $t$，二维中的解严格满足物理约束，而在三维中，违反约束的集合是维度 $\leq 1$ 的“薄”集。

意义与主张
本文声称在两个主要方面改进了以往的研究（特别是 [13] 以及 Evans–Kneuss–Hung [6] 的静态变分工作）：

1. 动态扩展： 虽然 [6] 建立了静态能量极小值的局部正则性，但本研究将接触集估计扩展到了动态演化系统，在动态系统中，静态变分技术并不直接适用。
2. 物理真实性： 作者处理了特定的对数奇异势 (1.7) 情况，该势函数表现出比早期一般框架中所假设的代数爆炸更弱的发散性。通过证明该特定物理势的强制性不等式并利用由此产生的 $H^1$ 正则性，他们获得了更精细的结果：在 2D 中接触集完全为空（这改进了原本暗示 2D 中为 0 维集合的一般性 $H^{n-2+\epsilon}$ 界限，因为如果没有特定的爆炸分析，该界限并不能保证集合为空）。

作者强调，3D 中的结果（维度 $\leq 1$）是由于物理势特有的对数爆炸速率与自身 $H^1$ 正则性（通过新的对偶形式和球面平均论证得出）之间相互作用的直接结果。这项工作为带有奇异熵势的 Landau-de Gennes 梯度流在定义域内几乎处处自然地保持物理约束提供了严格的理论依据。