Blow-up of solutions to the Euler-Poisson-Darbox equation with critical power nonlinearity

1. 問題設定

対象とする方程式は、$\mathbb{R}^n$ における以下の半線形 E.P.D. 方程式のコーシー問題である：

$$\begin{cases} u_{tt} - \Delta u + \frac{\mu}{t} u_t = t^\alpha |u|^p, & t > 0, \ x \in \mathbb{R}^n, \\ u(0, x) = u_0(x), \quad u_t(0, x) = 0, & x \in \mathbb{R}^n. \end{cases}$$

ここで、$\mu > 0$ は減衰係数、$\alpha \ge 0$ は時間重み、$p > 1$ は非線形指数である。
特に、臨界ケース $p = p_S(n+\mu, \alpha)$ を扱う。ここで $p_S(n+\mu, \alpha)$ は以下の二次方程式の正の根として定義される（シフトされたストラス指数）：
$$\gamma(n, \mu, \alpha, p) := 2 + (n + \mu + 1 + 2\alpha)p - (n + \mu - 1)p^2 = 0.$$

前著 [10] では、劣臨界 ($1 < p < p_S$) な場合の有限時間爆発と存在時間の上界が確立されていたが、本論文ではその境界にあたる臨界ケースにおける爆発の証明と、より精密な存在時間推定（指数関数的な上界）を目的としている。

2. 手法と主要な技術的工夫

臨界ケースの解析は、劣臨界ケースに比べて極めて困難であり、解の存在時間が指数関数的に長くなる（爆発に近づきにくい）ため、従来の試験関数法では不十分である。本論文では、以下の 2 つの重要な改良を導入してこの困難を克服した。

(1) 改良された試験関数の構成

従来の反復法（iteration argument）において、非線形項の積分に対する下界評価を強化するために、スケーリングパラメータに関する積分を用いた試験関数を構築した。

• 修正ベッセル関数 $K_{\frac{\mu-1}{2}}(t)$ を用いて、微分方程式 $h''(t) - (\frac{\mu}{t}h(t))_t - h(t) = 0$ の解 $h(t)$ を構成する。
• これをスケーリングパラメータ $\lambda$ で積分し、新しい試験関数 $b_q(t, x)$ を定義する（式 (46)）。
  $$b_q(t, x) = \int_0^1 h(\lambda t) \phi(\lambda x) \lambda^{q-1} d\lambda$$
• この構成により、試験関数の漸近挙動（式 (49), (50)）が改善され、特異点 $t=0$ における振る舞いを制御できる。

(2) 強化された機能関数（Blow-up Functional）の構築

非線形項の時空間積分を含む機能関数 $Y(M)$ を定義し、これに対する対数増大項を含むより良い下界評価を得た。

• 機能関数 $Y(M)$ は、非線形項の積分を重み関数 $b_q$ とカットオフ関数 $\eta$ で重み付けしたものである（式 (66)）。
• 劣臨界ケースでは得られなかった「追加の対数増大（additional logarithmic growth）」を、試験関数の漸近挙動と積分順序の交換（式 (72)）を用いて導出する。
• これにより、$Y(t) \gtrsim \varepsilon^p \ln t$ という評価（式 (71)）が得られる。

3. 主要な結果

定理 1.2 において、以下の条件のもとで有限時間爆発と存在時間の上界が示された。

• 次元とパラメータの条件：
  • $n \ge 3, \alpha \ge 0$ または $n=2, \alpha > 0$
  • $0 < \mu \le \mu^*(n, \alpha)$ （ここで $\mu^*$ は式 (11) で定義される閾値）
  • 臨界指数 $p = p_S(n+\mu, \alpha)$
• 初期データ $u_0(x) \in H^1(\mathbb{R}^n)$ は正で、コンパクトな台を持つ。

結論：
十分小さな $\varepsilon > 0$ に対して、エネルギー解の存在時間 $T(\varepsilon)$ は以下の上界を満たす：
$$T(\varepsilon) \le \exp\left( C \varepsilon^{-p(p-1)} \right)$$
ここで $C$ は $\varepsilon$ に依存しない正の定数である。

4. 技術的詳細と証明の鍵

1. 特異性の除去: 方程式の減衰項 $\frac{\mu}{t}u_t$ は $t=0$ で特異である。試験関数 $b_q$ の構成において、修正ベッセル関数 $K_{\frac{\mu-1}{2}}$ の漸近挙動（式 (41)）と漸化式（式 (42)）を巧みに利用することで、$t \to 0$ における特異性が相殺され、積分可能になることを示した（式 (73), (74)）。
2. 微分不等式の導出: 機能関数 $Y(t)$ に対して、$Y'(t)$ と $Y(t)$ の関係から、対数項を含む微分不等式 $t Y'(t) \ge (\ln t)^{1-p} Y^p$ （式 (83)）を導出した。
3. 爆発の証明: この微分不等式を積分することで、$Y(t)$ が有限時間で発散すること、すなわち解が有限時間で爆発することを示した。

5. 意義と貢献

• 未解決問題への回答: D'Abbicco (2021) が提起した「半線形 E.P.D. 方程式の爆発結果の完全な理解が欠如している」というオープン問題（特に臨界ケース）に対して、部分的な回答を提供した。
• 臨界ケースの解析の進展: 臨界指数における爆発は、解の存在時間が指数関数的に長くなるため、劣臨界ケースよりも遥かに繊細な解析を要する。本論文は、試験関数のスケーリング積分と対数増大項の導出という新しい手法を確立し、この難問に挑む道筋を示した。
• 物理的・数学的応用: E.P.D. 方程式は曲面論、音波伝播、衝突する重力場などの物理分野で重要である。本結果は、これらの物理モデルにおける非線形現象の振る舞い（特に臨界点での不安定性）の理解を深めるものである。

6. 結論

本論文は、半線形オイラー・ポアソン・ダルブー方程式の臨界べき非線形性に対する有限時間爆発を証明し、存在時間の指数関数的な上界を確立した。特に、特異な減衰項を持つ問題において、改良された試験関数と対数増大項を伴う機能関数評価を用いることで、従来の手法では達成できなかった結果を導出した点に大きな学術的価値がある。