解は、いつまでも存在するとは限りません。

ある時刻までは普通の値を取っていた関数が、有限の時間で急激に大きくなり、無限大へ向かってしまうことがかし、解の爆発は、

• 化学反応の暴走
• 燃焼
• 生物や細胞の急激な増殖
• 流体の特異性
• 重力崩壊
• ランダムな外力を受ける物理現象

などを理解するために研究されている、重要な数学的テーマです。

この記事では、

1. 微分とは何か
2. 「バイバイン」のような指数的増殖
3. 常微分方程式における有限時間爆発
4. 空間を含む熱方程式の爆発
5. ランダムなノイズを加えた確率熱方程式
6. 2026年時点の研究

という順に、高校数学から最先端の研究までを一気につなげて説明します。

---

1．そもそも「微分」とは何か

自動車が2時間で100km進んだとします。

このときの平均速度は、

$$\frac{100}{2}=50$$

より、時速50kmです。

しかし、実際の自動車が2時間ずっと時速50kmで走っていたとは限りません。

信号で止まった時間もあれば、一時的に時速80kmで走った時間もあるでしょう。

そこで知りたくなるのが、

ある瞬間に、どれくらいの速さで変化しているのか

という量です。

時刻 $t$ におけるある量を $u(t)$ とします。

時刻 $t$ から少しだけ時間 $h$ が経過したとき、量は $u(t)$ から $u(t+h)$ へ変化します。

この間の平均変化率は、

$$\frac{u(t+h)-u(t)}{h}$$

です。

ここで、時間の幅 $h$ を限りなく0に近づけます。

その極限として、

$$u'(t)=\lim_{h\to0}\frac{u(t+h)-u(t)}{h}$$

を定義します。

これが関数 $u(t)$ の微分です。

微分 $u'(t)$ は、時刻 $t$ における $u$ の瞬間的な変化の速さを表します。

例えば、$u(t)$ が自動車の位置なら、$u'(t)$ はその瞬間の速度です。

$u(t)$ が人口なら、$u'(t)$ はその瞬間に人口が増えている速さです。

---

2．変化の速さから未来を予測する

ある量 $u(t)$ の増える速さが、現在の量 $u(t)$ に比例するとします。

この関係は、

$$\frac{du}{dt}=ku$$

と書けます。

ここで、$k>0$ は増殖の強さを表す定数です。

左辺の

$$\frac{du}{dt}$$

は、$u$ の時間に対する変化の速さです。

右辺の $ku$ は、現在の量 $u$ が大きいほど、増える速さも大きくなることを表しています。

このように、未知の関数とその微分の関係を表した方程式を、微分方程式といいます。

---

3．「バイバイン」を微分方程式にすると？

一定時間ごとに個数が2倍になる物体を考えます。

最初に1個だったものが、

$$1,\ 2,\ 4,\ 8,\ 16,\ 32,\ldots$$

と増えていく状況です。

これは、ドラえもんに登場する「バイバイン」のような増殖です。

個数が多いほど、次の瞬間に増える個数も多くなります。

時間を連続的に考えるなら、この現象は、

$$\frac{du}{dt}=ku$$

という微分方程式で表せます。

この方程式を解いてみましょう。

$u(t)>0$ として両辺を $u$ で割ると、

$$\frac{1}{u}\frac{du}{dt}=k$$

となります。

形式的に、$u$ に関する部分と $t$ に関する部分を分けると、

$$\frac{1}{u},du=k,dt$$

です。

この方法を変数分離法といいます。

両辺を積分すると、

$$\int\frac{1}{u},du=\int k,dt$$

です。

左辺は $\log u$ なので、

$$\log u=kt+C$$

となります。

ここで $C$ は積分定数です。

両辺の指数関数を取ると、

$$u(t)=Ae^{kt}$$

となります。

初期値を、

$$u(0)=u_0$$

とすると、

$$u_0=Ae^0=A$$

です。

したがって、

$$\boxed{u(t)=u_0e^{kt}}$$

となります。

---

4．5分ごとに2倍になる場合

5分ごとに個数が2倍になるとします。

最初の個数が $u_0$ なら、5分後には、

$$u(5)=2u_0$$

でなければなりません。

一方、

$$u(t)=u_0e^{kt}$$

なので、

$$u(5)=u_0e^{5k}$$

です。

したがって、

$$u_0e^{5k}=2u_0$$

より、

$$e^{5k}=2$$

となります。

両辺の対数を取ると、

$$5k=\log2$$

なので、

$$k=\frac{\log2}{5}$$

です。

よって解は、

$$u(t)=u_0e^{(\log2)t/5}$$

となります。

指数関数の性質から、

$$e^{(\log2)t/5}=2^{t/5}$$

なので、

$$\boxed{u(t)=u_0,2^{t/5}}$$

と書くこともできます。

---

5．指数関数は驚くほど速い

最初に1個あり、5分ごとに2倍になるとします。

1時間は60分なので、1時間後には12回倍になります。

したがって、

$$2^{12}=4096$$

個です。

2時間後には24回倍になるので、

$$2^{24}=16777216$$

個になります。

3時間後には、

$$2^{36}=68719476736$$

個です。

指数関数は、非常に速く増加します。

しかし、ここで重要な点があります。

どれほど急激に増加しても、有限な時刻 $t$ に対しては、

$$e^{kt}<\infty$$

です。

つまり、指数関数は有限時間では無限大になりません。

時間を無限に大きくしたときに初めて、

$$\lim_{t\to\infty}e^{kt}=\infty$$

となります。

これは数学でいう「有限時間爆発」ではありません。

---

6．数学における「解の爆発」とは何か

関数 $u(t)$ が、ある有限な時刻 $T$ までは存在しているとします。

しかし、$t$ を $T$ に近づけたとき、

$$\lim_{t\uparrow T}u(t)=\infty$$

となることがあります。

このとき、解は時刻 $T$ で爆発するといいます。

記号 $t\uparrow T$ は、$t$ が $T$ より小さい側から $T$ に近づくことを表しています。

有限時間爆発とは、単に値が大きくなることではありません。

有限な時間しか経過していないのに、解が無限大へ向かう

という現象です。

---

7．本当に有限時間で爆発する方程式

次の微分方程式を考えます。

$$\frac{du}{dt}=u^2$$

初期値を、

$$u(0)=u_0>0$$

とします。

変数分離すると、

$$\frac{1}{u^2},du=dt$$

です。

両辺を積分すると、

$$\int u^{-2},du=\int dt$$

となります。

左辺を計算すると、

$$-\frac{1}{u}=t+C$$

です。

初期値 $u(0)=u_0$ を代入すると、

$$-\frac{1}{u_0}=C$$

となります。

したがって、

$$-\frac{1}{u}=t-\frac{1}{u_0}$$

です。

両辺にマイナスを掛けると、

$$\frac{1}{u}=\frac{1}{u_0}-t$$

となります。

右辺を通分すると、

$$\frac{1}{u}=\frac{1-u_0t}{u_0}$$

です。

よって、

$$\boxed{u(t)=\frac{u_0}{1-u_0t}}$$

となります。

---

8．爆発時刻を計算する

解は、

$$u(t)=\frac{u_0}{1-u_0t}$$

でした。

分母が0になる時刻を求めると、

$$1-u_0T=0$$

です。

したがって、

$$\boxed{T=\frac{1}{u_0}}$$

となります。

$t$ が $T$ に近づくと、分母は正の値を保ったまま0へ近づきます。

そのため、

$$\lim_{t\uparrow1/u_0}\frac{u_0}{1-u_0t}=\infty$$

です。

つまり、この解は、

$$\boxed{T=\frac{1}{u_0}}$$

で爆発します。

例えば、$u_0=1$ なら、

$$u(t)=\frac{1}{1-t}$$

です。

この解は $t=1$ で爆発します。

実際、

$$u(0.9)=10$$

であり、

$$u(0.99)=100$$

さらに、

$$u(0.999)=1000$$

です。

$t$ が1へ近づくにつれて、解が急激に大きくなっていくことが分かります。

---

9．初期値が大きいほど早く爆発する

爆発時刻は、

$$T=\frac{1}{u_0}$$

でした。

初期値が $u_0=1$ なら、

$$T=1$$

です。

初期値が $u_0=10$ なら、

$$T=\frac{1}{10}=0.1$$

です。

初期値が $u_0=100$ なら、

$$T=\frac{1}{100}=0.01$$

です。

つまり、

最初の値が大きいほど、爆発する時刻が早くなる

ことが分かります。

現在の量が大きいほど増加速度 $u^2$ が大きくなり、その増加によってさらに $u^2$ が大きくなるからです。

これが自己増幅です。

---

10．一般の $p$ 乗ではいつ爆発するか

次に、より一般的な微分方程式、

$$\frac{du}{dt}=u^p$$

を考えます。

初期値は、

$$u(0)=u_0>0$$

とします。

$p\neq1$ の場合、変数分離すると、

$$u^{-p},du=dt$$

です。

両辺を積分すると、

$$\frac{u^{1-p}}{1-p}=t+C$$

となります。

初期値を代入すると、

$$\frac{u_0^{1-p}}{1-p}=C$$

です。

したがって、

$$\frac{u^{1-p}}{1-p}=t+\frac{u_0^{1-p}}{1-p}$$

となります。

両辺に $1-p$ を掛けると、

$$u^{1-p}=u_0^{1-p}+(1-p)t$$

です。

よって、

$$\boxed{u(t)=\left(u_0^{1-p}+(1-p)t\right)^{\frac{1}{1-p}}}$$

となります。

$p>1$ の場合には、

$$1-p=-(p-1)$$

なので、

$$\boxed{u(t)=\left(u_0^{1-p}-(p-1)t\right)^{-\frac{1}{p-1}}}$$

と書けます。

---

11．$p>1$ の場合

$p>1$ のとき、解は、

$$u(t)=\left(u_0^{1-p}-(p-1)t\right)^{-\frac{1}{p-1}}$$

です。

括弧の中は時間とともに小さくなり、やがて0になります。

0になる時刻を $T$ とすると、

$$u_0^{1-p}-(p-1)T=0$$

です。

したがって、

$$\boxed{T=\frac{u_0^{1-p}}{p-1}}$$

となります。

この時刻に近づくと、

$$u(t)\to\infty$$

です。

よって、

$$\boxed{p>1\text{ なら有限時間で爆発する}}$$

ことが分かります。

---

12．$p=1$ の場合

$p=1$ のとき、方程式は、

$$\frac{du}{dt}=u$$

となります。

解は、

$$u(t)=u_0e^t$$

です。

この解は急激に増加しますが、すべての有限時刻において有限です。

したがって、有限時間では爆発しません。

---

13．

13．$0<p<1$ の場合

<p<1$ の場合

$0<p<1$ のとき、

$$1-p>0$$

です。

解は、

$$u(t)=\left(u_0^{1-p}+(1-p)t\right)^{\frac{1}{1-p}}$$

となります。

括弧の中は時間とともに増加しますが、有限な時刻では有限です。

したがって、この場合も有限時間爆発は起こりません。

以上をまとめると、

$$\frac{du}{dt}=u^p,\qquad u(0)=u_0>0$$

に対して、

$$\boxed{p>1\text{ のとき、そしてそのときに限って有限時間爆発する}}$$

となります。

---

14．爆発するかどうかを判定する一般公式

さらに一般的な微分方程式、

$$\frac{du}{dt}=b(u)$$

を考えます。

$b(u)>0$ とすると、$u$ は時間とともに増加します。

変数分離すると、

$$\frac{du}{b(u)}=dt$$

です。

初期値を $u(0)=u_0$ とし、$u_0$ から $u(t)$ まで積分すると、

$$t=\int_{u_0}^{u(t)}\frac{ds}{b(s)}$$

となります。

解が無限大になるまでに必要な時間は、形式的には、

$$T=\int_{u_0}^{\infty}\frac{ds}{b(s)}$$

です。

この積分が有限なら、無限大に到達するまでに有限時間しか必要ありません。

したがって、

$$\boxed{\int_{u_0}^{\infty}\frac{ds}{b(s)}<\infty}$$

なら、有限時間爆発が起こります。

反対に、

$$\boxed{\int_{u_0}^{\infty}\frac{ds}{b(s)}=\infty}$$

なら、無限大に到達するためには無限の時間が必要です。

この判定条件は、一般にオズグッド条件と呼ばれます。

例えば、

$$b(u)=u^p$$

なら、

$$\int_{u_0}^{\infty}\frac{ds}{s^p}$$

を考えます。

この広義積分が有限になるのは、

$$p>1$$

の場合です。

先ほど求めた爆発条件と一致しています。

---

ここから空間を取り入れる

15．一点の量から、場所ごとの量へ

これまでは、

$$u(t)$$

という関数だけを考えてきました。

これは、時間によって変化する1つの量です。

しかし、実際の温度や化学物質の濃度は、場所によって異なります。

そこで、

$$u(t,x)$$

という関数を考えます。

ここで、

• $t$ は時間
• $x$ は空間上の位置
• $u(t,x)$ は時刻 $t$、場所 $x$ における温度や濃度

を表します。

時間と空間の両方に依存する関数を扱うとき、通常の微分ではなく偏微分を使います。

---

16．熱方程式とは何か

一次元の熱方程式は、

$$\frac{\partial u}{\partial t}=\kappa\frac{\partial^2u}{\partial x^2}$$

です。

簡略化して、

$$\boxed{u_t=\kappa u_{xx}}$$

とも書きます。

ここで、$\kappa>0$ は熱の広がりやすさを表す定数です。

$u_t$ は、ある場所の温度が時間とともにどのくらい変化しているかを表します。

一方、

$$u_{xx}$$

は、温度分布のグラフがその場所でどのように曲がっているかを表します。

---

17．熱方程式は「でこぼこを平らにする」

ある場所の温度が周囲より高く、温度分布の山になっているとします。

山の頂上では、

$$u_{xx}<0$$

となります。

熱方程式より、

$$u_t=\kappa u_{xx}<0$$

なので、その場所の温度は時間とともに下がります。

反対に、ある場所の温度が周囲より低く、谷になっているとします。

谷の底では、

$$u_{xx}>0$$

となります。

したがって、

$$u_t=\kappa u_{xx}>0$$

なので、その場所の温度は上昇します。

つまり熱方程式は、

$$\boxed{\text{山を低くし、谷を高くする}}$$

働きを持っています。

その結果、温度分布のでこぼこは時間とともになだらかになります。

---

18．多次元の熱方程式

空間が $n$ 次元の場合には、ラプラシアンと呼ばれる微分作用素、

$$\Delta u=\frac{\partial^2u}{\partial x_1^2}+\frac{\partial^2u}{\partial x_2^2}+\cdots+\frac{\partial^2u}{\partial x_n^2}$$

を使います。

多次元の熱方程式は、

$$\boxed{u_t=\kappa\Delta u}$$

です。

$\Delta u$ は、ある場所の値と周囲の値の違いを測り、空間的なでこぼこを平らにする方向に働きます。

---

19．拡散と増殖を同時に起こす

熱方程式に、先ほどの増殖項 $u^p$ を加えます。

定数を簡単にするため、ここからは $\kappa=1$ とします。

すると、

$$\boxed{u_t=\Delta u+u^p}$$

という方程式が得られます。

これは半線形熱方程式の代表例です。

また、拡散と反応を同時に表すため、反応拡散方程式とも呼ばれます。

右辺には、2つの反対方向の働きがあります。

$$\underbrace{\Delta u}{\text{周囲へ拡散させる}}+\underbrace{u^p}{\text{その場所で増殖させる}}$$

$\Delta u$ は、温度や濃度を周囲へ広げ、極端に大きな値を小さくしようとします。

一方、$u^p$ は、値が大きくなるほど、さらに強く値を増加させます。

したがって、この方程式では、

拡散が増殖を抑え込むのか
増殖が拡散に勝って爆発するのか

という競争が起こります。

---

20．空間がある場合の「解の爆発」

常微分方程式では、

$$u(t)\to\infty$$

となることを爆発と呼びました。

空間を含む場合には、場所によって値が異なります。

そこで通常は、

$$|u(t,\cdot)|_{\infty}=\sup_x|u(t,x)|$$

という量を考えます。

これは、時刻 $t$ における $|u(t,x)|$ の最大値です。

有限な時刻 $T$ に対して、

$$\lim_{t\uparrow T}|u(t,\cdot)|_{\infty}=\infty$$

となるとき、解は時刻 $T$ で爆発するといいます。

これは、少なくともどこかの場所で、温度や濃度が限りなく大きくなることを意味します。

---

21．拡散があるのに、なぜ爆発するのか

熱方程式だけなら、温度分布は平らになっていきます。

ところが、

$$u_t=\Delta u+u^p$$

では、値が大きい場所ほど $u^p$ が大きくなります。

例えば $p=2$ なら、

$$u=10\quad\Longrightarrow\quad u^2=100$$

ですが、

$$u=100\quad\Longrightarrow\quad u^2=10000$$

です。

値が10倍になると、増殖速度は100倍になります。

一度ある場所の値が十分大きくなると、反応項 $u^p$ による増殖が拡散よりも速くなります。

その結果、拡散が値を周囲へ逃がすよりも速く、その場所の値が増加し、有限時間爆発が起こることがあります。

---

22．空間全体では「次元」が運命を決める

空間全体 $\mathbb{R}^n$ で、

$$u_t=\Delta u+u^p$$

を考えます。

初期条件は、

$$u(0,x)=u_0(x)\geq0$$

とします。

さらに、$u_0$ は恒等的に0ではないとします。

この方程式には、

$$\boxed{p_{\mathrm F}=1+\frac{2}{n}}$$

という特別な指数が現れます。

これを藤田指数といいます。

---

23．藤田指数より小さい場合

指数 $p$ が、

$$1<p\leq1+\frac{2}{n}$$

を満たすとします。

このとき、0ではない非負の初期値から出発した解は、有限時間で爆発します。

重要なのは、初期値をどれほど小さくしても爆発することです。

つまり、

$$\boxed{1<p\leq1+\frac{2}{n}}$$

の範囲では、拡散だけでは増殖を抑えきれません。

---

24．藤田指数より大きい場合

一方、

$$p>1+\frac{2}{n}$$

の場合には、状況が変わります。

初期値が十分小さければ、解がすべての時刻で存在する場合があります。

このように、すべての時刻で存在する解を大域解といいます。

ただし、初期値が十分大きければ、同じ指数 $p$ でも有限時間爆発が起こります。

したがって、藤田指数は、

どんなに小さな初期値でも爆発する領域

と、

初期値が十分小さければ爆発を避けられる領域

を分ける境界です。

---

25．空間次元ごとの藤田指数

藤田指数は、

$$p_{\mathrm F}=1+\frac{2}{n}$$

です。

空間次元が $n=1$ のとき、

$$p_{\mathrm F}=1+\frac{2}{1}=3$$

です。

$n=2$ のとき、

$$p_{\mathrm F}=1+\frac{2}{2}=2$$

です。

$n=3$ のとき、

$$p_{\mathrm F}=1+\frac{2}{3}=\frac{5}{3}$$

です。

$n=4$ のとき、

$$p_{\mathrm F}=1+\frac{2}{4}=\frac{3}{2}$$

です。

まとめると、

• $n=1$ では $p_{\mathrm F}=3$
• $n=2$ では $p_{\mathrm F}=2$
• $n=3$ では $p_{\mathrm F}=\dfrac{5}{3}$
• $n=4$ では $p_{\mathrm F}=\dfrac{3}{2}$

となります。

次元が高くなるほど、熱や物質は多くの方向へ広がることができます。

そのため、拡散の影響が相対的に強くなり、すべての非自明な解が爆発する指数の範囲は狭くなります。

---

26．有限の箱の中では事情が変わる

今度は空間全体ではなく、有限の長さの棒を考えます。

例えば、

$$0<x<1$$

とします。

棒の両端の温度を常に0に保つなら、境界条件は、

$$u(t,0)=u(t,1)=0$$

です。

これをディリクレ境界条件といいます。

方程式は、

$$u_t=u_{xx}+u^p$$

です。

有限区間では、任意の $p>1$ に対して、

• 初期値が十分小さければ大域解が存在する
• 初期値が十分大きければ有限時間爆発する

という状況が起こります。

境界で $u=0$ に固定されているため、熱や濃度が境界から外へ逃げられるからです。

したがって、爆発するかどうかは指数 $p$ だけでは決まりません。

• 初期値の大きさ
• 空間次元
• 領域の大きさや形
• 境界条件
• 反応項の形

などにも左右されます。

---

27．この方程式は何を表しているのか

方程式、

$$u_t=\Delta u+u^p$$

は非常に単純化されたモデルですが、多くの現象の基本構造を表しています。

発熱を伴う化学反応

$u(t,x)$ を温度と考えます。

温度が高くなるほど化学反応が速くなり、反応によってさらに熱が発生するとします。

このとき、

• $\Delta u$ は熱が周囲へ逃げる効果
• $u^p$ は化学反応によって熱が発生する効果

を表します。

熱が逃げる速さよりも発熱の速さが大きければ、温度は急激に上昇します。

数学的な爆発は、このような熱暴走を理想化して表していると考えられます。

生物の個体数や細胞密度

$u(t,x)$ を生物や細胞の密度と考えることもできます。

この場合、

• $\Delta u$ は個体や細胞が周囲へ移動する効果
• $u^p$ は密度が高い場所で急激に増殖する効果

を表します。

ただし、現実には食料や空間に限界があります。

そのため、生物モデルでは通常、増殖を抑える項も加えます。

燃焼

$u(t,x)$ を温度や反応物の濃度と考えます。

高温になるほど燃焼反応が進み、さらに熱が発生する場合、反応項と拡散項の競争が起こります。

反応が拡散を上回ると、急激な温度上昇が発生します。

---

さらにランダムな揺らぎを加える

28．現実の世界にはノイズがある

現実の現象は完全に決定論的ではありません。

温度には微小な揺らぎがあります。

化学反応では、分子がランダムに衝突します。

生物の誕生や死亡にも偶然性があります。

外部から加わる力も、正確には予測できない場合があります。

そこで、熱方程式にランダムなノイズを加えます。

代表的な方程式は、

$$\boxed{\frac{\partial u}{\partial t}=\Delta u+b(u)+\sigma\dot W}$$

です。

ここで、

• $b(u)$ は反応や増殖を表す項
• $\sigma>0$ はノイズの大きさ
• $\dot W$ は時空間ホワイトノイズ

を表します。

---

29．時空間ホワイトノイズとは何か

時空間ホワイトノイズ $\dot W$ は、時間と空間の両方において非常に細かく変動するランダムなノイズです。

直感的には、テレビ画面の砂嵐のようなものを時間と空間の両方に加えるイメージです。

ただし、数学的な $\dot W$ は通常の関数ではありません。

非常に不規則であり、各点で値を持つ普通の関数としては扱えません。

そのため、

$$\frac{\partial u}{\partial t}=\Delta u+b(u)+\sigma\dot W$$

という式も、通常の微分方程式とまったく同じ意味では解釈できません。

厳密には、熱核を使って積分したマイルド解として定義します。

---

30．確率熱方程式のマイルド解

熱方程式の解を進める作用を $S(t)$ と書きます。

形式的には、

$$S(t)=e^{t\Delta}$$

です。

初期値を $u_0$ とすると、確率熱方程式のマイルド解は、おおよそ、

$$u(t)=S(t)u_0+\int_0^tS(t-s)b(u(s)),ds+\sigma\int_0^tS(t-s),dW(s)$$

と表されます。

第1項、

$$S(t)u_0$$

は、初期値が熱方程式によって拡散した部分です。

第2項、

$$\int_0^tS(t-s)b(u(s)),ds$$

は、反応項による増殖を表します。

第3項、

$$\sigma\int_0^tS(t-s),dW(s)$$

は、ランダムなノイズの影響を表します。

高校数学の範囲を超える式ですが、

拡散、増殖、ランダムな揺らぎを全部足したもの

だと考えればよいでしょう。

---

31．小さなノイズでも爆発を引き起こすのか

有限区間における決定論的な方程式、

$$u_t=u_{xx}+u^p$$

では、初期値が十分小さければ爆発を避けられる場合がありました。

ところが、一定の強さの加法ノイズを加えた、

$$u_t=u_{xx}+b(u)+\sigma\dot W,\qquad \sigma>0$$

では、状況が変化します。

適切な仮定のもとでは、解が有限時間で爆発するかどうかが、

$$\int^\infty\frac{ds}{b(s)}$$

の有限性によって決まることが示されています。

つまり、

$$\boxed{\int^\infty\frac{ds}{b(s)}<\infty}$$

であることが、確率偏微分方程式でも爆発条件として現れます。

これは、常微分方程式で登場したオズグッド条件と同じ形です。

---

32．反応項が $u^p$ のように増加する場合

大きな正の $u$ に対して、

$$b(u)\approx u^p$$

とします。

$p>1$ なら、

$$\int^\infty\frac{ds}{s^p}<\infty$$

です。

そのため、一定の技術的条件のもとで、加法ノイズを含む方程式にも有限時間爆発が生じます。

ノイズの係数 $\sigma$ が非常に小さくても、$\sigma>0$ である限り、ノイズは解を一時的に大きな正の値へ押し上げる可能性があります。

一度解が十分大きくなると、反応項 $b(u)$ の自己増幅が強くなります。

その後は、

$$u’=u^p$$

の場合と同じように、増加が増加を呼ぶ状態になります。

やがて拡散だけでは抑えられなくなり、爆発へ向かうことがあります。

---

33．ノイズを加えた瞬間に爆発するわけではない

「小さなノイズでも爆発を引き起こす」といっても、

ノイズを加えた瞬間、必ず解が無限大になる

という意味ではありません。

爆発時刻はランダムです。

ノイズが小さい場合には、解が危険な大きさに達するまで長い時間がかかる可能性があります。

反対に、偶然大きな正の揺らぎが発生すれば、比較的早く爆発することもあります。

決定論的な方程式では、同じ初期値なら未来は1つに決まります。

一方、確率偏微分方程式では、ノイズの実現ごとに異なる未来が生じます。

そのため、爆発時刻そのものが確率変数になります。

---

34．ノイズは爆発を起こすだけではない

ノイズは、必ず解を不安定にするとは限りません。

次のような方程式を考えます。

$$\boxed{u_t=u_{xx}+u^\beta+u^\gamma\dot W}$$

ここでは、ノイズの大きさ自体が $u^\gamma$ によって変化します。

このようなノイズを乗法ノイズといいます。

$u$ が大きいほど、ノイズも強くなります。

直感的には、ノイズが強くなれば、さらに爆発しやすくなるように思えます。

しかし、確率偏微分方程式では、適切な強さの乗法ノイズが、反応項による爆発を抑えることがあります。

---

35．ノイズが爆発を防ぐ仕組み

乗法ノイズは、解を上向きにも下向きにも揺らします。

解が大きくなったときにノイズも強くなると、解を大きく下向きへ押し戻す揺らぎも強くなります。

その結果、決定論的な反応項だけなら爆発する状況でも、確率的な揺らぎによって爆発を回避できることがあります。

2024年に公表された研究では、一次元の確率熱方程式について、反応項が超線形に増加する場合でも、乗法ノイズが反応項に比べて十分に強く、なおかつ速く増加しすぎない場合には、解が爆発せず大域的に存在することが証明されました。

大まかには、

適度に強いノイズは爆発を抑える可能性がある

ということです。

この現象は、確率安定化と呼ばれる考え方につながります。

---

36．ノイズが強すぎると再び爆発する

一方、乗法ノイズが解の大きさに対して速く増加しすぎると、ノイズそのものが爆発の原因になります。

一次元の時空間ホワイトノイズを含む確率熱方程式では、

$$\sigma(u)\approx u^\gamma$$

という場合に、指数、

$$\gamma=\frac{3}{2}$$

が重要な境界として現れます。

大まかには、

• $\gamma\leq\dfrac{3}{2}$ の範囲では非爆発を示せる場合がある
• $\gamma>\dfrac{3}{2}$ ではノイズ自身が爆発を引き起こし得る

という構造があります。

したがって、ノイズには少なくとも3つの顔があります。

1. 加法ノイズが反応による爆発のきっかけを作る
2. 適切な乗法ノイズが反応による爆発を防ぐ
3. 強すぎる乗法ノイズが、それ自体で爆発を引き起こす

爆発するかどうかは、単にノイズの有無だけでは決まりません。

反応項とノイズ項の増加速度のバランスが重要です。

---

爆発問題の最先端

37．研究者は何を調べているのか

「爆発するか、しないか」が分かった後にも、多くの問題が残ります。

例えば、

• どの場所で爆発するのか
• 最初に爆発する場所は1点なのか複数点なのか
• 爆発時刻はどのくらいか
• 解はどの速度で無限大へ向かうのか
• 爆発直前のグラフはどのような形か
• 初期値を少し変えても爆発は起こるのか
• 爆発する初期値と爆発しない初期値の境界は何か
• ノイズが爆発時刻をどのように変えるのか
• 爆発後に解を何らかの意味で延長できるのか

といった問題です。

---

38．爆発速度

常微分方程式、

$$u’=u^p$$

の爆発時刻を $T$ とします。

解は、

$$u(t)=\left((p-1)(T-t)\right)^{-\frac{1}{p-1}}$$

という形で表せます。

したがって、$t$ が $T$ に近づくとき、

$$u(t)\sim C(T-t)^{-\frac{1}{p-1}}$$

という速さで大きくなります。

ここで $C>0$ は定数です。

半線形熱方程式、

$$u_t=\Delta u+|u|^{p-1}u$$

でも、常微分方程式と同じ程度の速度で爆発する場合があります。

この基準的な爆発を、一般にタイプI爆発と呼びます。

---

39．タイプI爆発とタイプII爆発

半線形熱方程式の解が時刻 $T$ で爆発するとします。

タイプI爆発では、解の大きさが、

$$|u(t,\cdot)|_\infty\leq C(T-t)^{-\frac{1}{p-1}}$$

のように、常微分方程式から予想される速度の範囲内で増大します。

一方、この基準を超えて、より特殊な速度や形で爆発するものをタイプII爆発と呼びます。

偏微分方程式には空間構造があります。

そのため、

• 解が一点へ集中する
• 特別な定常解や自己相似解へ近づく
• 空間的な形を変えながら爆発する

など、常微分方程式にはない複雑な現象が起こります。

---

40．2025年に報告された爆発速度の研究

次の半線形熱方程式を考えます。

$$u_t=\Delta u+|u|^{p-1}u$$

2025年に公表された研究では、凸とは限らない領域や非有界領域において、符号が変化する解も含めた爆発速度が研究されました。

エネルギー劣臨界と呼ばれる範囲、

$$(n-2)p<n+2$$

において、タイプII爆発が起こらないことが証明されています。

つまり、この範囲では、解が爆発するとしても、その速度は常微分方程式から予測されるタイプIの範囲に収まります。

これは1980年代から残されていた問題の解決として報告されています。

爆発研究では、

解が爆発するかどうか

だけでなく、

爆発するときに、どれほど速く大きくなれるのか

まで分類しようとしているのです。

---

41．爆発か大域存在かを分ける条件

藤田指数より大きい、

$$p>1+\frac{2}{n}$$

の領域では、すべての解が爆発するわけではありません。

小さな初期値では大域解が存在し、大きな初期値では爆発することがあります。

そこで重要になるのが、

与えられた初期値が、爆発側にあるのか、大域存在側にあるのか

を判定する問題です。

---

42．2026年に提案された新しい判定方法

2026年5月に公表された研究では、

$$u_t-\Delta u=|u|^{p-1}u$$

という半線形熱方程式について、相似変換と重み付きソボレフ空間を用いた新しい判定条件が提案されました。

対象となる指数の範囲は、

$$p_{\mathrm F}<p<p_{\mathrm S}$$

です。

ここで、

$$p_{\mathrm F}=1+\frac{2}{n}$$

は藤田指数です。

また、$n\geq3$ のとき、

$$p_{\mathrm S}=\frac{n+2}{n-2}$$

をソボレフ臨界指数といいます。

この研究では、スケール不変性のため通常のポテンシャル井戸法を直接使いにくい問題を、相似変換によって別の方程式へ移し、重み付き空間上で解析しています。

高校数学から始まった、

$$u’=u^p$$

という単純な方程式が、最先端では、

• 相似変換
• ソボレフ空間
• スケール不変性
• エネルギー法
• ポテンシャル井戸法

といった高度な解析へつながっています。

---

43．反応と乗法ノイズの境界

確率熱方程式、

$$u_t=u_{xx}+u^\beta+\sigma(u)\dot W$$

を考えます。

大きな $u$ に対して、

$$\sigma(u)\approx u^\gamma$$

とします。

この方程式では、反応項の指数 $\beta$ と、ノイズ項の指数 $\gamma$ の組み合わせによって、

• 解が大域的に存在する領域
• 正の確率で爆発する領域
• まだ十分に理解されていない領域

に分かれます。

---

44．2026年に報告された確率的爆発の範囲

2026年5月に公表された研究では、一次元の周期境界条件を持つ確率熱方程式について、次の範囲で解が正の確率で爆発することが示されました。

まず、

$$\beta>1,\qquad 0<\gamma\leq\frac{\beta}{2}$$

という範囲です。

さらに、

$$1<\beta<3$$

かつ、

$$\frac{\beta}{2}<\gamma<\frac{\beta+3}{4}$$

という範囲でも、適切な初期値から出発した解が正の確率で爆発することが示されています。

ここで「正の確率」とは、必ず爆発するという意味ではありません。

ノイズの実現によっては爆発しない可能性もありますが、爆発する確率が0ではないという意味です。

---

45．まだ未解決の領域が残っている

反応項 $u^\beta$ と乗法ノイズ $u^\gamma\dot W$ のどちらが強いかによって、解の運命は大きく変化します。

しかし、すべての $\beta$ と $\gamma$ の組み合わせについて、爆発するか大域存在するかが完全に分類されているわけではありません。

反応による自己増幅と、乗法ノイズによる不安定化または安定化が拮抗する範囲には、現在も未解決の領域が残されています。

つまり、2026年時点でも、

ノイズはどの強さまで爆発を防ぎ、どの強さから爆発を促進するのか

という基本的な問題が完全には解決していません。

---

46．数学的な爆発は現実の無限大なのか

現実の温度や個体数が、本当に無限大になるとは限りません。

実際の現象では、極端な状態へ近づく前に、それまで無視していた別の効果が現れます。

例えば、

• 燃料がなくなる
• 化学反応の法則が変化する
• 物質の状態が変化する
• 圧力や密度の影響が無視できなくなる
• 生物の食料や空間が不足する
• 相対論的効果が必要になる
• 量子力学的効果が必要になる
• 方程式そのものの適用範囲を超える

といったことが起こります。

したがって、数学的な爆発は、

$$\boxed{\text{現在の数理モデルが有限時間で極端な状態へ達する}}$$

ことを表しています。

爆発時刻は、その方程式だけでは現象を説明できなくなる限界時刻とも解釈できます。

---

47．解が爆発することは「数学の失敗」ではない

解が有限時間で存在しなくなると聞くと、方程式や数学に欠陥があるように感じるかもしれません。

しかし、爆発は計算の失敗ではありません。

むしろ、

どの仮定が極端な増加を生み出しているのか
どの時刻までモデルが有効なのか
どの効果を新たに加える必要があるのか

を教えてくれる重要な情報です。

爆発の研究によって、モデルの限界や、現象を安定化させる仕組みが見えてきます。

---

まとめ

微分とは、ある量がその瞬間にどれくらい変化しているかを測る操作です。

増加速度が現在の量に比例する方程式、

$$u’=u$$

の解は、

$$u(t)=u_0e^t$$

です。

指数関数は急速に増加しますが、有限時間では爆発しません。

一方、

$$u’=u^p$$

では、$p>1$ のとき、

$$u(t)=\left(u_0^{1-p}-(p-1)t\right)^{-\frac{1}{p-1}}$$

となり、

$$T=\frac{u_0^{1-p}}{p-1}$$

で有限時間爆発します。

空間変数を加えた半線形熱方程式、

$$u_t=\Delta u+u^p$$

では、拡散と増殖の競争が起こります。

空間全体 $\mathbb{R}^n$ では、

$$p_{\mathrm F}=1+\frac{2}{n}$$

という藤田指数が、すべての非自明な非負解が爆発する領域と、小さな初期値なら大域解が存在する領域を分けます。

さらにノイズを加えた確率熱方程式、

$$u_t=\Delta u+b(u)+\sigma(u)\dot W$$

では、

• 加法ノイズが爆発のきっかけを作る
• 適切な乗法ノイズが爆発を防ぐ
• 強すぎる乗法ノイズが爆発を引き起こす

という複雑な現象が現れます。

研究の最先端では、

• 爆発時刻
• 爆発場所
• 爆発速度
• 爆発直前の形
• 初期値による爆発と大域存在の境界
• 反応項とノイズ項の臨界的なバランス

が調べられています。

高校数学で解ける単純な方程式、

$$u’=u^p$$

から始まった問題が、偏微分方程式、関数解析、確率解析、ソボレフ空間、相似変換といった最先端の数学へつながっているのです。

解の爆発は、

解はいつまで存在できるのか
無限大になるなら、いつ、どこで、どの速さで爆発するのか
拡散やランダムな揺らぎは、爆発を防ぐのか、それとも促進するのか

を問う研究です。

これは高校数学から現代数学の最前線までを一本につなぐ、非常に奥深いテーマなのです。

---

参考文献

M. Salins and Y. Zhang, “Positive Probability of Explosion for Stochastic Heat Equation with Superlinear Accretive Reaction Term and Polynomially Growing Multiplicative Noise”, arXiv:2605.11319.

H. Fujita, “On the Blowing Up of Solutions of the Cauchy Problem for $u_t=\Delta u+u^{1+\alpha}$”, 1966.

M. Foondun and E. Nualart, “The Osgood Condition for Stochastic Partial Differential Equations”, arXiv:1907.12096.

M. Salins, “Global Solutions to the Stochastic Heat Equation with Superlinear Accretive Reaction Term and Polynomially Growing Multiplicative Noise”, arXiv:2409.15527.

H. Miura, J. Takahashi and E. Zhanpeisov, “Blow-up Rate for the Subcritical Semilinear Heat Equation in Non-convex Domains”, arXiv:2510.17229.

K. Zhang and Z. Li, “On the Existence and Nonexistence of Global Solutions of the Semilinear Heat Equation”, arXiv:2605.11933.