ヌルクライン

ヌルクライン（英: nullcline）またはゼロ等傾線（ゼロとうけいせん）とは、常微分方程式系、特に平面自励系の相平面解析において用いられる幾何学的な概念であり、ある変数の時間微分が0となる相空間上の集合（曲線や曲面）のことである。

ヌルクラインは、解析的に解くことが困難な非線形微分方程式の解の挙動（定性理論）を理解するために不可欠なツールである。これを描くことで、解の軌道（ベクトル場）の概略的な方向、平衡点の存在と安定性、リミットサイクルの有無などを視覚的に把握することが可能となる。

定義

一般に、 $n$ 次元の自励系（力学系）：

{\frac {dx_{i}}{dt}}=f_{i}(x_{1},x_{2},\dots ,x_{n})\quad (i=1,\dots ,n)

において、第 $i$ 成分のヌルクライン $N_{i}$ は、幾何学的には以下の条件を満たす集合として定義される^[1]。

N_{i}=\{(x_{1},\dots ,x_{n})\in \mathbb {R} ^{n}\mid {\frac {dx_{i}}{dt}}=f_{i}(x_{1},\dots ,x_{n})=0\}

2次元系における性質

ヌルクラインが最も頻繁に用いられるのは、2次元の自励系（平面力学系）である。以下の系を考える。

{\begin{cases}{\frac {dx}{dt}}=f(x,y)\\{\frac {dy}{dt}}=g(x,y)\end{cases}}

このとき、 $x$ -ヌルクラインと $y$ -ヌルクラインは以下のように定義される。

$x$ -ヌルクライン: 曲線 $f(x,y)=0$ のこと。この線上では $x$ の変化率が0であるため、ベクトル場は「垂直」方向（ $y$ 軸に平行）となる^[1]。
$y$ -ヌルクライン: 曲線 $g(x,y)=0$ のこと。この線上では $y$ の変化率が0であるため、ベクトル場は「水平」方向（ $x$ 軸に平行）となる^[1]。

領域の分割と符号

ヌルクラインは相平面をいくつかの基本領域に分割する。関数 $f(x,y)$ や $g(x,y)$ が連続である場合、ヌルクラインで区切られた各領域内では、微分の符号（ $dx/dt$ の正負、 $dy/dt$ の正負）は変化しない^[1]。

これにより、各領域における解軌道の概略的な方向を以下のように決定できる。

$dx/dt>0,dy/dt>0$ の領域：軌道は「北東」へ進む。
$dx/dt<0,dy/dt>0$ の領域：軌道は「北西」へ進む。
$dx/dt<0,dy/dt<0$ の領域：軌道は「南西」へ進む。
$dx/dt>0,dy/dt<0$ の領域：軌道は「南東」へ進む。

平衡点とその安定性

すべてのヌルクラインの交点は、すべての成分の微分が0となる点であるため、系の不動点（平衡点）となる。

P\in N_{x}\cap N_{y}\iff f(P)=g(P)=0

2次元系においては、 $x$ -ヌルクラインと $y$ -ヌルクラインの交点が平衡点に対応する。

線形化不能な場合の判定: 平衡点近傍でのヌルクラインの幾何学的配置と、その間の領域におけるベクトル場の向きを確認することで、線形化解析を行わずに平衡点が吸引的（安定）か反発的（不安定）かを判定できる場合がある^[2]。
傾きによる判定: 特定の生物モデル（競争モデルなど）では、交点における2つのヌルクラインの傾きの大小関係が、その平衡点の安定性を決定する条件となる^[3]。

解析への応用

ヌルクラインは、微分方程式の解のグローバルな挙動を証明・説明するために用いられる。

トラッピング領域とリミットサイクル

ポアンカレ・ベンディクソンの定理を用いてリミットサイクル（閉軌道）の存在を示す際、ヌルクラインは「トラッピング領域（Trapping Region）」を構築するための境界線として利用される。例えば、解糖系のモデルであるセルコフ (Sel'kov) モデルにおいては、ヌルクラインで囲まれた領域を設定し、ベクトル場がその領域から外へ出ないことを示すことで、振動解の存在を証明できる^[2]。

弛緩振動

ファン・デル・ポール振動子や神経パルスのモデルにおいて、パラメータが極端な値をとる場合（特異摂動）、解の軌道はヌルクラインに非常に近い場所をゆっくりと移動（クロール）し、ヌルクラインの折り返し点などで急速に別の場所へ移動（ジャンプ）するという挙動を示す。このような弛緩振動 (Relaxation Oscillation) のメカニズムは、ヌルクラインの形状（S字型など）によって幾何学的に説明される^[2]。

分岐現象

パラメータの変化によって2つのヌルクラインの配置が変化し、交点（平衡点）の数が増減したり、交点の位置関係が変わったりすることがある。これは力学系における分岐（サドルノード分岐、トランスクリティカル分岐など）に対応する^[1]^[3]。

具体的なモデル例

数理生物学や神経科学などの分野では、ヌルクラインを用いた解析が標準的な手法となっている。

ロトカ・ヴォルテラの方程式: 捕食者と被食者の個体数変動モデル。自明なヌルクライン（座標軸）と非自明なヌルクライン（直線）を持ち、その交点として共存平衡点が現れる。ヌルクライン解析により、個体数が周期的変動をする閉軌道の存在が示唆される^[3]。

競争モデル: 2種の生物が資源を争うモデル。それぞれの種のヌルクライン（通常は直線）が相平面上でどのように交差するか（あるいは交差しないか）によって、「競争排除則（一方が絶滅する）」が成立するか、あるいは「共存」が可能かが決定される^[3]。

フィッツフュー-南雲モデル: 神経細胞の膜電位変化を記述するモデル。膜電位 $v$ のヌルクラインが3次関数（N字型または逆N字型）の形状をしており、回復変数 $w$ の直線ヌルクラインとの交差位置によって、系が「静止状態」にあるか、外部刺激によって「発火（一回だけのパルス）」するか、あるいは自発的に「反復発火（リミットサイクル）」するかが決まる^[3]^[2]。

化学反応モデル: 酵素反応における基質阻害モデルや、ベロウソフ・ジャボチンスキー反応（BZ反応）のモデルなどでは、ヌルクラインが複雑な曲線（ベル型など）を描き、それが振動現象やスイッチング現象（複数の安定状態の切り替わり）を引き起こす要因となる^[3]^[2]。