ホタルのきらきらランダムウォーク
草桶翠

1(指導教員：郡宏)

はじめにシンクロとは，リズムの集団が協調性を持ち，大きなリズムを作り出すことである．
ホタルの集団は空をふわふわ飛びながら，別のホタルを見ると光るタイミングをシンクロさせる．また，飛ぶ方向も合わせているように見えることがある．この様子をシミュレーションで再現し，どのような条件下でどのような集団秩序が生まれるのかを調べる．

2モデル

2.1

ホタルの発光リズム

ホタルはそれぞれ固有のリズムを持っており，それにしたがって周期的に光るとする．それを次のようなモデルで表す．
(φ)φ̇i = ωi + µ(φ) ξi

図 1: シミュレーションのスナップショット
ここで，~ri (t) = (xi (t), yi (t)) は個体 i の位置，θi は個体の進行方向，vi は進行速度である．
本研究では，飛ぶ方向を合わせる現象をふわふわシンクロと呼ぶ．ホタルは光っている別のホタルを認識すると，飛ぶ方向を変化させるとする．ふわふわシンクロを導く相互作用を以下のようにする．

(1)

ここで，φi (t) はホタルの発光リズムの位相で，ωi は固有振動数，右辺第 2 項は発光周期のゆらぎを表し，(φ)µ(φ) は発光リズムに対するノイズ強度，ξi は白色ノイズである．個体数は N とする (i = 1, 2, . . . , N )．この位相を用いて，ホタル i の発光強度を次の関数で表現する．
cos φi −1p(φi ) = e τ 2(2)N∑dθi(θ)= µ(θ) ξi + H
G(Rij )p(φj ) sin(θj − θi ) (6)dt
j=1

ここで，µ(θ) は飛ぶ方向に関するノイズ強度，ξ (θ) は白色ノイズである．H はホタル同士の飛ぶ方向の結合強度，項 sin(θj − θi ) はホタル i がホタル j に方向を揃える効果を表している．

τ は光る時間の長さである．
2.2

きらきらシンクロ

本研究では，光るタイミングを合わせる現象をきらきらシンクロと呼ぶ．ホタルは別のホタルを認識すると，光るタイミングを変化させるとする．まず，ホタル i に到達する光の量はQi =N∑G(Rij )p(φj )以上のモデルを使い，数値シミュレーションを行った．個体数 N = 30，固有振動数 ωi = 1，進行速度
vi = 5.0，発光リズムに対するノイズ強度 µ(φ) = 0.1，飛ぶ方向に対するノイズ強度 µ(θ) = 0.5，光る時間の長さ τ = π/6，発光リズムの結合強度 K = 1.0，飛ぶ方向の結合強度 H = 0.1，初期条件は，xi と yi については範囲 [0, L] の一様乱数で与え，φi と θi については範囲 [0, 2π] の一様乱数で与えた．processing 言語でホタルの動きを可視化した結果が，図 1 である．

(3)

j=1である．ここで，Rij = |~ri − ~rj | はホタル i とホタルj の距離，G(Rij ) はホタル j からホタル i に到達する光の強度である．ここでは，2 次元空間に光が放射する状況，つまり G(Rij ) = R12 とする．きらきらシン3.1ij

クロを導く相互作用を以下のようにする．

dφi
Qi(φ)
= ωi + µ(φ) ξi − Ksin(φi )dt
Qi + 1

シンクロ率

きらきらシンクロとふわふわシンクロが，それぞれどの程度のものなのかを計算する．その計算式を以下のようにする．
vuN
Nu 1 ∑1 ∑Rφ (t) = t(cos φj )2 + (sin φj )2 (7)NN(4)

ここで，K は発光リズムの結合強度である．ホタルが光を感じる量に限界があるため，ホタルが感じる量をQi
Qi +1 としている．

2.3

シミュレーション

3j=1j=1vuNNu 1 ∑1 ∑Rθ (t) = t(cos θj )2 + (sin θj )2N j=1N j=1

ふわふわシンクロ空間を２次元の周期境界条件 (L × L) で考える．ホタルのランダムウォークモデルを次のモデルで表現する．
()
d~ri
cos θi= vi(5)dtsin θi(8)

Rφ (t) と Rθ (t) をそれぞれ，きらきら度，ふわふわ度と呼ぶ．これらの量は，0 がばらばらな状態，1 に近づくほどシンクロしている状態を表す．
63

11

Kirakira
Fuwafuwa

0.8

0.8

R

0.6

R

0.6
0.4

0.4

0.2

0.2

000400800
time12001600

KirakiraFuwafuwa
00.02

0.04

0.06

0.080.1

H

図 2: シンクロ率の時系列

図 4: 飛ぶ方向の結合強度 H とシンクロ率

0.8

0.8

0.6

0.6
R1R

10.4

0.4

0.2
00.2

Kirakira
Fuwafuwa
00.2

0.4

0.6

0.8

01

Kirakira
Fuwafuwa
020

K406080100N

図 3: 発光リズムの結合強度 K とシンクロ率

図 5: 個体数 N とシンクロ率

これらのシンクロ率を長時間観察した結果が図 2 である．右上がりのグラフがきらきら度，大きく振動しているグラフがふわふわ度である．
シミュレーションを注意深く観察し，以下のことがわかった．きらきらシンクロは，まず複数のグループごとにシンクロし，光るタイミングを揃える．その後，グループごとに近くのグループとシンクロし始める．
その結果，グループが段々大きなものになり，最終的には１つの集団となって７割前後のシンクロ率を保つ．
その一方で，ふわふわシンクロは時間が経過していってもシンクロ率が安定しない．きらきらシンクロと同様に，複数のグループごとに飛ぶ方向を揃えるが，それがより大きなグループになろうとする時に，グループ内での飛ぶ方向が定まらず，また，どこか別のグループの飛んでいる方向に影響されたりするので，お互いのメンバーが入れ替わったりグループの規模が変化したりはするが，なかなか大きなグループではシンクロしない．たまたま大きなグループができたとしても，他の個体の影響と，ホタル自身の持つ飛行のランダム性により，そのシンクロ率は持続しない．これは，発光リズムに対するノイズ強度よりも，移動方向に対するノイズ強度を大きく設定しているためだと考えられる．

ふわ度においては，発光リズムの結合強度の変化に依存しない．しかしきらきら度においては，飛ぶ方向の結合強度にも依存してシンクロ率を上げている．これは，ふわふわ度が上がることでホタルが認識する光が増えるため，きらきら度も同時に上がると考えられる．

3.2

3.3

密度とシンクロ率

個体数 N を変化させたものが図 5 である．きらきら度，ふわふわ度ともに，個体数を増やす (密度を上げる) とシンクロ率が上がる．これは限られた空間内でのシミュレーションなので，密度が上がるとホタル同士の距離が近くなり，認識し合うホタルも増えてお互いがシンクロし合うので，よりシンクロ率が上がると考えられる．特にふわふわシンクロについては，最初のシミュレーションではグループがなかなか大きくならなかったが，密度が上がるにつれてホタルのグループ自体が大きくなり，高いシンクロ率を保つことができるのではないかと推測できる．

4

まとめと今後の課題

本研究では，結合強度や個体数 (密度) とシンクロ率の相互関係について，シミュレーションを通して調べた．その結果，結合強度や個体数が上がるとシンクロ率も上がることがわかった．今後は，ノイズや周囲の明るさなどの環境が与える作用によって，シンクロ率がどのように変化していくか，ホタルの光り方に変化があるかを調べ，より現実に近いものにしたい．

結合強度とシンクロ率パラメーターの値を変化させて，シンクロ率がどのように変化するかを調べる．変化させるパラメーター
以外の値は，上のシミュレーションと同じ条件にする．
10000 単位時間経過したのち，10000 単位時間にわたって平均したシンクロ率を図 3〜5 にプロットした．
発光リズムの結合強度 K ，飛ぶ方向の結合強度 H を変化させたものがそれぞれ図 3，図 4 である．きらきら度，ふわふわ度とも，それぞれの結合強度を上げるに伴い，そのシンクロ率も上がることがわかる．ふわ参考文献[1] 郡宏，森田善久．生物リズムと力学系．共立出
版，2011．
[2] Tamas Vicsek，Anna Zafeiris．Physics Reports：Collective motion，pp.71-140，2012．

64