JOVANA
Explore Library Glossary Getting Started Three Levels Fields How it works Mission
Join the mission
All guides

反應擴散、行進波前與斑圖

在純粹的擴散上加一個反應項,這個原本只會抹平一切的方程,忽然有了自己的意志:入侵以確定的速度橫掃空間,而均勻的狀態自發碎裂成條紋與斑點。讓我們認識反應擴散方程、行進波前,以及那個驚奇——擴散竟能「失穩」而非撫平。

會反擊的擴散

你在拋物那一級學到的一切都說,擴散是個和事佬:擴散方程 u_t = k u_xx 會把任何粗糙的初始輪廓無情地抹平、推向平坦,總量守恆,凸起只會散開、淡去。現在加上一個「不」守恆任何東西的新項——一個反應,它在局部依據當地已有多少 u 而創造或摧毀 u。結果就是反應擴散方程 u_t = k u_xx + f(u):u_xx 那部分仍把物質從擁擠處搬往空曠處,但 f(u) 讓每一點都能自行繁衍或死亡。這是最簡單的半線性方程——對最高階導數而言是線性的,只透過未微分的 f(u) 才非線性——而這一個非線性項,徹底改變了整個動力學的性格。

把 u(x,t) 想成沿一條線鋪開的族群密度。擴散是個體的隨機漫遊;反應 f(u) 是出生與死亡。一個著名的選擇是邏輯斯諦增長 f(u) = r u (1 - u):當 u 很小時族群幾乎以指數增長,但當 u 趨近承載容量 1 時,增長便關閉。光有擴散,族群會逐漸消融;光有反應,族群會在各處均勻增長、毫無地點之分。把兩者合起來,你就得到本篇要談的真正嶄新的現象——一場有組織的入侵,以一個鋒利而明確的速度橫越空曠的領地推進。

一個自帶形狀的波

核心對象在此:一個行進波。它是一個解,看起來像一個固定的輪廓沿著線以等速剛性滑動,u(x,t) = U(x - c t)。形狀 U 從不改變;整個圖樣只是平移。我們在波動方程那裡已見過這個想法,那裡任何 f(x - c t) 都會行進——但在那裡速度 c 是方程直接交給我們的。反應擴散令人驚奇之處在於,方程「不」告訴你速度;波必須自己去發現 c,而且只有某些速度才是可能的。

要看清為何如此,把這個試解代進去。令 z = x - c t。則 u_t = -c U'(z)、u_xx = U''(z),於是整條偏微分方程坍縮成一個只含單一變數 z 的常微分方程:k U'' + c U' + f(U) = 0。奇蹟在於,偏微分方程變成了常微分方程——因為行進波的假設移除了時間這個獨立變數,把整幅圖景凍結在那個與波前一同移動的座標系裡。我們現在要找一個輪廓 U(z),它從一個靜止狀態爬升到另一個:當 z 趨向負無窮時 U 趨向 1(已被入侵的領地),當 z 趨向正無窮時 U 趨向 0(尚未抵達的領地)。

PDE:   u_t = k u_xx + f(u)

  ansatz  u(x,t) = U(z),   z = x - c t
        |
        |  u_t = -c U',   u_xx = U''
        v
ODE:   k U'' + c U' + f(U) = 0

  boundary conditions on the profile:
        U(-infinity) = 1   (invaded state behind the front)
        U(+infinity) = 0   (empty state ahead of the front)

  --> only special speeds c admit such a connecting profile
行進波試解把反應擴散偏微分方程化為移動座標系中的一條常微分方程;速度 c 不是給定的——它是由「輪廓必須連接兩個靜止狀態」這個要求所挑選出來的。

對於邏輯斯諦反應 f(u) = r u (1 - u),這就是著名的費雪-KPP 方程,由費雪(為一個橫掃族群的有利基因建模)以及柯爾莫哥洛夫、彼得羅夫斯基與皮斯庫諾夫各自獨立寫下。對每一個大於或等於某個最小值 c* = 2 sqrt(k r) 的速度 c,都存在一個行進波前,而你實際上看到從局部起點長出的那個波前,總是恰恰挑選那個最慢的可容許速度。所以入侵速度由僅僅兩個數字決定——個體擴散有多快、繁殖有多快——以 2 sqrt(k r) 的形式結合。這一條公式就預測了一個物種以多快速度蔓延入新棲地、火焰以多快速度吃進燃料、一個有利基因以多快速度接管全局。

雙穩波前:一場有確定贏家的拔河

費雪-KPP 有一個穩定狀態(u = 1)與一個不穩定狀態(u = 0),所以入侵是單向的:滿的狀態總會征服空的。一個更豐富、也很不一樣的故事,來自艾倫-卡恩方程 u_t = k u_xx + u (1 - u^2)。現在有「兩個」穩定的靜止狀態 u = +1 與 u = -1,中間隔著一個不穩定的狀態 u = 0。這是兩相材料的經典模型——比方兩種晶體取向,或朝上對朝下的磁化——沿一個界面相遇。行進波現在是一個波前,一側連接 -1、另一側連接 +1,它移動的方向,是讓較受偏好的相去吃掉較不受偏好的相。

當兩相恰好同等受偏好時,波前靜止不動——c = 0——而你能寫出輪廓的閉式解:U(z) = tanh(z / sqrt(2k)),一個從 -1 到 +1 的平滑 S 形階梯,其寬度由 sqrt(k) 決定。把平衡稍微傾斜(讓反應不對稱,譬如 f(u) = u(1-u^2) 加一個小常數),波前便開始漂移,其速度正比於兩邊能量有多不對稱。這就是雙穩波前傳播,它是相變、神經脈衝波前,以及斑圖隨時間粗化(較大的疇域緩慢吞掉較小的)背後的引擎。

圖靈的驚奇:建造斑圖的擴散

現在來到真正顛覆直覺的結果。至此為止擴散都還扮演它抹平的角色;新意是反應提供的。1952 年,艾倫·圖靈問了一個離經叛道的問題:擴散本身,這個偉大的平等主義者,會不會反而從均勻中「創造」結構?對單一種化學物,答案是不會。但取「兩種」會反應、會擴散的化學物——一個活化劑 u 既促進自身、又促進一個抑制劑 v,而 v 壓制 u——並把它們耦合起來:u_t = D_u u_xx + f(u,v)、v_t = D_v v_xx + g(u,v)。圖靈證明了,一個空間均勻的穩態,在化學物無法移動時完美穩定,卻會在你一旦讓它們擴散的那一刻被「失穩」。這就是圖靈不穩定性,也是整個這一級裡最反直覺的事實。

這個機制有個生動的名字:局部活化、遠程抑制。關鍵的要求是抑制劑擴散得遠比活化劑快,D_v 遠大於 D_u。想像一個偶然出現的小活化劑凸起。它在局部放大自己,所以凸起想長成一個尖峰。但它也產生抑制劑,而抑制劑向外奔逃的速度遠快過活化劑的擴散,於是在尖峰四周一圈寬廣的範圍內壓制活化劑。結果是一個被壓制的護城河圍住的尖峰——而這些護城河的自然間距,定下了一個偏好的波長。均勻狀態碎裂成一列規則排列的斑點或條紋,其大小由擴散常數挑選,而非由初始資料決定。

你實際上要怎麼檢驗這件事?在均勻穩態附近做線性化,看一個形如 cos(q x) 的小擾動——單一個空間波數 q。每個波數各以自己的速率增長或衰減;不穩定性出現的條件是,某一「帶」中等波數有正的增長率,而 q = 0(均勻模態)與極大的 q(細微的擺動,被擴散扼殺)兩者都衰減。那一帶中增長最快的 q,就是你將看到浮現的波長。這正是你會對任何形成斑圖的系統所做的線性穩定性加波數分析,也是斑圖形成的分析核心。

什麼為真、什麼是草圖、它棲身何處

對線性分析告訴你什麼、不告訴你什麼,要誠實。圖靈的計算只證明均勻狀態不穩定,並預測「最先」出現之斑圖的波長;它對最終的振幅或確切形狀一無所說,因為擾動一旦增長,非線性項就接管,線性化便不再有效。系統最終安頓成斑點、條紋還是六邊形,這些如何競爭與粗化,是個全然非線性的問題,通常需要電腦或細緻的弱非線性理論。線性穩定性是「斑圖從何處開始」的可靠偵測器,而非「它們最終落腳何處」的描述。

第二件誠實話:反應擴散是一個美麗而有預測力的理想化模型,但它是真實生物學與化學的一幅漫畫。別洛索夫-扎博廷斯基反應確實在培養皿裡造出行進的化學波與螺旋,這是貨真價實的勝利。但魚身上具體的條紋、或豹身上的斑點,是否真的字面上由一個兩化學物的圖靈機制所產生,仍有爭議——這些方程令人信服地捕捉了自發成圖的「邏輯」,而真正的分子角色往往未知、或更為複雜。把這個模型用作對機制的洞見,而不是「自然恰好在解這兩條方程」的主張。

最後,把這篇放回學習階梯中的位置。費雪-KPP 與艾倫-卡恩的波前是非線性的,但它們永遠保持「光滑」、行為良好——擴散讓拋物型的光滑化始終活著,沒有任何東西爆掉。那份溫和,是會飽和的反應所特有的。緊接著的下一篇要問,當反應「不」飽和會發生什麼——當 f(u) 像 u^p 那樣增長、快過擴散所能耗散的速度時——那裡的解能在有限時間內逃逸到無窮。所以這一篇是暴風雨前的寧靜:這裡的非線性是組織與建構,而本學習階梯其餘部分則展示摧毀一切的非線性。