...

Page 1 京都大学 京都大学学術情報リポジトリ 紅

by user

on
Category: Documents
6

views

Report

Comments

Transcript

Page 1 京都大学 京都大学学術情報リポジトリ 紅
Title
生体モデル膜のカオス的自励発振(カオスとその周辺,研
究会報告)
Author(s)
渕上, 信子; 沢島, 信介; 神原, 武志; 八木沢, 亨一; 内藤, 正
美; 佐々木, 直幸
Citation
Issue Date
URL
物性研究 (1990), 53(5): 551-558
1990-02-20
http://hdl.handle.net/2433/93974
Right
Type
Textversion
Departmental Bulletin Paper
publisher
Kyoto University
「カオスとその周辺」
生体 モデ ル膜 の カオス的 自励 発 振
都立大理 測上信子,沢島信介,電通大 神原武志,八木沢亨-,
日立基礎研 内藤正美,日歯大 佐々木直幸
1. は じめ に
'
,K'
,
C
l等のイオン濃度差 によっ
生体系 において,神経興菅現象が膜の内側 と外側 の Na
て駆動 されていることはよく知 られている 【
1
】. 一方,生体膜の二分子構造 を摸 したモデル膜
において も,浪度勾 配や直流電流印加 な ど種 々の条件下で,膜電位 が振動することが最近多 く
の実験 によ り示 されている 【
2日 31. ことに,浪度勾配下での 自励発振 は,特定のチ ャネル黄白
を持 たない比較的単純 な合成膜が,イオン濃度差のみによる電気 的興奮 とい う神経類似機能を
持つ ことを意味する点で興味深 い.
自励発振の振幅や周波数 ,波形は実験条件 によ り様々で,周期 的な ものか ら概周期的.また
4】
. 我々はこの ような自励発振現象 を再現する簡単
明 らかに*.
オス的な場合 も報告 されている【
なモデルシステムを提案 し,計算機 シミュレーションを行 なった. 系 は倍周期分岐 カオスや間
欠性 カオスへの転移 を含 む多彩 な分岐構造 を示す. ここではお もにカオス的運動 に焦点 をあて
て報告する.
2. モデル
濃度の異 なる二つ?溶液層 とこれ らを仕切
る膜 (脂質二分子膜 を想定 )か ら成 る図 1の よ
うな 1次元系 を考え ,x十 と Y とい う二種類の
イオンの運動 に注 目す る. 陰イオン Y をあら
高浪度溶液
佐渡度溶液
(X,
=Y
,
)
(XI=YI)
ヘルムホルツ
ヘルムI
i
5ルツ
わに含めることによ り,イオ ン分布 と電場をポ
電気二重層
電気二重層
ワソン方程式 を通 じて整合的に関係づけること
図1
. モデル .
がで きる. この点がこのモデルの特徴 である.
溶液内および膜内のイオ ンの流 れの駆動力 としては,浪度勾 配のみならず,電場 による ドリフ
トの影響 も考慮する. また,膜の表面でのイオ ンの出入 りについて非線形性 を仮定す る.
イオン濃度 は場所 によって連続的 に変化 しているので,本来ならば無限自由度系 であるが,
a
.
b
,
Bにおけるイオ ン濃度
ここでは近似的 に,膜の状態 を表面近傍 の薄い面状の領域 A.
xp/A.yp/A(F
L=A,a.b
.
B;Aは両の厚 さ )のみで記述する・ 膜 に比べて電解質溶液の電
気抵抗 は十分小 さいので,膜表面の電荷 は溶液中のイオンによって遮蔽 され,膜から十分遠方
.
Rを考 え,浴
では至 るところ電気的に中性であるとしてよい. そこで,膜の両側 に二つの面 L
液中における電荷分布の不均一の影響 を,これ らの面上の電荷 (面密度 eoL.eon)のみで表
.R]内の全電荷 は運動のあいだ じゅうゼロである:
わす ことにする. 区間 【L.
qL+qA+qd+qb
+qB+qR=0
qv=Xp-Yp,
(F
L=A,
B);
qy=Xp-Yp-qn,
(
(F
L=a,b)
1)
(2)
而電荷 qL とq^(および qR と qB)は左側 (およ甲右側 )のヘルムホルツ電気二重層 【
5]に対応
する. 固定電荷 屯 は,溶液中では膜の極性基か ら陽イオンが解放 して膜が負 に帯電すること
を考慮 した ものである. イオンの連動の各 プロセスは次のように表わされる:
-
5
5
1-
研究会報告
十
Aa
ky
o
一
年 ーJ
l
+
dy2
xb #
k
x
b
▲b
、
xB
+
kyb
溶液 お よび膜 内の速度定数 dを,電位差 4
'による ドリフ トの効 果 を含 め て次 の形で表わす:
土
(p=x,Y;Z
'
=1,2,3)
d
l
d-dl
d(l±eP中.
・
/2),
¢2= ¢
O
-¢b
-e6(
qL+JA+qd)/en
^-eLqL/CW,
¢3≡QB- Q,
--eLq斤/E"
¢1= ¢L
-¢
6:膜 の厚 さ;
A:面 の厚 さ
L:ヘル ムホル ツ電気二重層 の幅 ;
P
≡1/k
B T
6
.
,ら :水 お よび膜 の誘電率 ; e:素電荷 ;
電位差 Q^一九 (-A)お よび 中.-QB(-A)は A<<6,
tに よ り無視 で きる. 膜 表面での
イオ ンの出入 りについて は次の よ うな非線 形性 を仮 定す る:
(p-a
,b)
p=k完.ki
,2(Y
y)2,
k;y=k
;
.
・k;
2(
Xp)2,
A;
(5)
第 2項 は 3次の 自己触媒的過程
+
X^+2Xa
kx
d
子
-二二三
k
x
a
3Xd
に対 応 してお り,この形 は膜 の 自励発振 につい ての Ka
wa
kub
oのモデルで も用 い られ た 【
5
1.
x
,
(
膜 か ら十分 遠方 の溶液 中のイオ ン浪度
-Y,
)
,
x,
(-Y
,) を与えた とき.イオ ン分布 お
よび電荷分布 は次の時間発展方定式 (連続の式 )に従 う:
慧
I
x,
慧
r Ix.
・
・
・l・
慧
ニーI l
・lr・ (
+
x ・
1
Io
-′Y,
・
・
-Ir,-
),
(i-1・
- ,4)
慧 -lx.5-Ir.5(-)
0,
(
xl
,
X2
,.・,
X9,
X1
0
)≡(
XA,..
1
,
XB,Y^,‥ ,YB,qL ,qR).
ここで た とえば
+
l
x.
1≡d
x
I
Xl
-d
x
I
XA ,
十
J
x.
2
≡kx
q
XA- kx
d
Xq,
e
t
c
.
(7)
(eI
。は外部電流 .
) ただ し中性条件 (
I
)式 か ら,独 立 な方程式は 9個 であ り,これ らは (
7)
,(
3)
,
(
4)
,
(
2)
,
(
5
)式 を代 入 して次の形 にま とめられ る:
-
552-
「カオス とその周辺」
2 -,
( )
I
(8)
x = (xl・・・・
X9)
I
正負 イオ ンの連動 は βを通 じて結合 している. β=0 の場合 はこれ らは分離 し,4変数
(x ^,..,
xB)についての運動方程 式 は,Kawakuboのモデルで P=0 と した場合に-敦する 【
5
]
.
なお, 方程 式の形 か らパ ラメ タ一 についてい くつかの端尺別 が成 り立 ち,た とえばあるバ ラメ
,x,
,
qh,
β,a,k)について解 x(t)が得 られた とす れば ,a を定数 として,これらを
タ-の観 (x ,
.
,/α,βう aP.k
'
x
2う d k'
x
2,
x,
.
,一 X,
k'
r
2う a2k'
r
2,(残 りの速度定数 d.kは不
変 )の ように変換 した ときの解 は ∫(り /α
で与 え られ る. 以下の数値計算 では,次 の値
l
o
o
(6
,I
,
A)i
nA =
(
5
0, ,1),
(Xl
/A,X,
/A,
qJ A)
i
nmM凡=
(
3.
5,
0.
1
,
1
.
0).
cp-0
.
48
..
を用 いた:
3. 結 果
(
3
1
)
陽イオ ンについてのみ非線形性 を仮定
した とき,バ ラメ タ-の適 当な領域で ,系 は
周期振動 をす る. 振動す るの は主 と して陽
イオ ン濃度 であるが ,これ に引 きず られて陰
イオ ンが弱 い振動 をす る場合 もある.
(
32
) 直流 電流 J
oを印加す る と,振動 の振幅
は増大 し,周期 はわずかに短 くなる.
O3
)
P.
陽イオ ンと陰 イオ ンの両方 に非線形性
がある場合 ,二種類 の振動 は結合す る. 以下
ではことに断 らない限 り.陰イオ ンの膜 内取
P,
P.C
図2
. 周期倍化 カオスを含む分岐構造 .
横軸:0.
1
270<k'
,
。<0.
1
3
02;縦軸:1
.
4<Y.._ <2.
9.
(
a
)拡大臥 胡 軸:0
.
1
28
88くA◆r
.く0.
1
2938;
縦軸:1
.
58<Y._ <1
.
68.
'
r
. をコン トロール 'パ
り込みの速度 定数 k
ラメ タ- と し,残 りのバ ラメ タ-は次 の値 に
固定 した:(d
n,
d
x
2
.
d
x
,
.
d,I.
a,
2
.
d,
,)i
J
IS-1=
(0.
2,1
.
1
,2.
0,
0.
3,1
.
5
,2.
8),(k'
x
.
,
Ax
.
,
k'
r
.)
i
ns-1-(0.
1
3,
2.
0,2.
8),(k'
x
2
,
k'
x
2
,
k'
,
2
,
A
n)i
n S-I(mMAP.
)
-2=
(
2.
6,1
.
6,
3.
7,
2.
2).
≡T.q =0.
5(T.
=298K).
倍周期分 岐 カオス
図 2は区間 0.
1
27く
β
k
1.<0.
1
302 における分岐 図である. 縦軸
はk
'
,
.
の値 を決め た ときに得 られる時系列
Y.( I) の極小値 Y
^.
血 を表 わす. k
1.∼
0.
1
291
5 で周期倍化 カスケー ドに よるカオス
への転移が見 られる. 周期振動 Pl(k
'
Y
.=
0.
1
2
7),
P2 (k
'
Y
.=0.
1
28),
P.(k
'
Y
.=0.
1
29)
お よびカオス振動 C (k'
Y
.=0.
1
29
2
4)におけ
るイオ ン浪度
x
.,Y.の時間変化 と,対応す
る位相空 間 (
xA,
Y.)での軌道 を図 3
,
4に示
x^
-Y^
す. これ らはいず れ も初期状態 (
-5
5
3-
図 3. X. お よび Y. の振動波形 .
横 軸:o<L<50 ;縦軸 :I
.
0<(X.
.
Y
.)<3
j.
'
'
o
P':
A FO・
1
27.P暮:
上 ◆,. EO
.
1
28.P/:k●
.
.LO.
1
29,
C:
上●,. Z
10・
1
292
4. 周期の長い方が Y.(I).
研究会報告
図 5. ㌔ (り のパ ワー ・スペ ク トル .
図 4. 振動状態の位相空P
I
J
での軌道 .
,縦軸 ・
.l<Y.<3・
収 軸:1<X.<3・
X,
XB =Y
B
=X,,Xa=Xb-6 Y.-Yb=0)
l
.
.
H
か ら出発 して十分時間が経 過 してか らの もの
である.図 5に 211個 の デ-タか ら計算 した
l
うP2
う P
パ ワー・
スペ ク トル を示す . P
.
と
周期倍化す るにつれて、次 々に低調波成分 が
現われてい る.時系列
㌦ (り
の 〝番 目と
(
a+1)番 目の極小値 Yi.. ,Y^.A.Iを横 臥
縦軸 に とってプロッ トした リター ン ●マ ップ
(ポア ンカ レ断面 )を Pl
,P2,P.,P8(k'r
。=
0.
1
291
2),
PIG(k
'
Y
。=0.
1
291
3),
C について
図 6に示す. C の リター ン ●マ ップか ら,
この カオス ・ア トラクターの次元が 2に近い
ことが伺 える.
カオス ・ア トラクターの相関次元
カ レ断面上 のサ ンプル点 を
がァ ン
x l,
X2
,..,
X〟
とす る とき,積分相 関関数
C (r)= t
i
m
Nサ
Jち∑.0(r一 一XL-XjJ)
仙N .
I
,
・
J
♂(r)=1,(r>0); ♂(r)≡0,(rく0)
と相 関次元
Vは次の関係 にある 【
7
】
:
図 6. Y
. の リター ン ・マ ップ (d
!7 ンカレ断面 ):
.
.
.
.の プロ ッ ト・
極小値 Y.一 一 Y
l
i
m C
(r)牝r
y
pl:
A◆
,
o
岩
;0.
1
291
2.PIt:A●
,
o
=0・12913・
(p.
.pl
.
P
.
.
Cのバ ラメ タ-掛 ま図 3を劉 臥 )
∫
◆0
ここでは x と して 9次 元 の ベ
クト ル (
X^,
-
5
5
4
-
「カオスとその周辺J
C
P
i
図 7. ポア ンカ レ断面の相 関次元 V.
図8
. 閃欠性 カオス を含 tl分岐構 造 .
横軸 :0.
1
2
9
20<k
'
,
。く0.
1
92
42;縦軸 :0・
9<y<l
・
3・
破軸 :0.
1
3
5<A
●
T
。く0.
1
3
9;縦 軸:1
.
5
5<Y.,∼.< 1
.
7
5.
p.(A'
,
。=0.
1
2
92
8)は 6周J
y
J
点なので
y=0・
‥,
XB,YA,‥ .YB,
6,
.) を用い,N=2
000について
C(r)を計算 L l
ogr に対する l
o
gC(r) の プロ
003くr<0.
03 における傾 きか ら
ッ トの,区間 0.
V を求めた. アル ゴリズムは p.
,p2
,P.,P8 につ
いて V-0 となることでチェック した. カオス
領域 におけるい くつかの点での
Vの億 を図 7に示
す. ア トラクターの相 関次元 は β =(y+1) で与
えられる.
1
35<k
'
r
。<0.
1
39 に
間欠性 カオス
区間 0.
おける分岐図 を図 8に示す. k'r。-0.
1
35
9 と k1.
∼0.
1
377 のあた りで転移がみ られるが,こ れはサ
ドル ■ノー ド分岐であ りヒステ リシス はない. (す
なわち,バ ラメタ-を逆向 きに変化 させで も同 じ
図 9. 朕花位 l.…中.+Ql+中,の振動波形 ・
横 軸:0<t
<0
0
0;縦 軸.
I0.
2<eP
丸 く-0
・
l・
L
1
p・:
k◆
.
o=0.
1
36.C.:よ
◆
,
〇三0・
1
358・
分岐図が得 られる.
)周期 振動 p(k
'
r
D
=0.136)と
カオス振動 C'(k'r。
=0.
1
358)
における膜電位 礼
…¢l+¢2+¢,の時間変化 を図 9に,Clにおける
Q
mの極小値の リター ン ・
マ ップを図 1
0に示す.
C-は間欠性 カオスである.
ヒステ リシス
区間 0.
1
03<k'r。<0.
1
27 にお
ける分 岐図を図 11に示す. k'
r
。,
vO.
11
0 のあた り
で見 られる二つの周期振動 間の不 連続転移 p
.H
Pb はヒステ リシス を伴 う. 図 11は k'r. を左側 か
ら右側 に増加 させて得 られたものであるが,逆 に右
側か ら減少 させ たときの分岐図 はこれ と完全 には一
致 しない. すなわち,区間 0.
1
085<k
'
r
o<0.
11
00
では PAと pb が共 に安定軌道 と して共存 している.
k
'
,
o
=0.109 で異 なる初期条件 か ら出発 した二つの
周期振動 を図 1
2に示す. なお ,k
'
r
。∼0
.
1
23 の 3
-
5
55-
\
、
ヽ\
.
・
図1
0・ 問欠性 カオス C・における +_の リター ン ・
マ ップ :穣小催 九.一 九.
.
.0
)プロット
研究会報告
図 11
. 不達扶転移 を含む分岐梢造 .
胡 軸:o・
1
03<上◆
,
。<0.
1
271
. 縦軸:1
.
0<Y.
.
山.
く3.
4.
(
a
)拡大図. 磯 軸:0.
)
23
4<k●
,
.<0.
1
2
44・
.
縦 軸:2.
4<Y.
.
.
h
hく2.
9. (
b)拡大 臥
r
^
h<2.
0.
胡 軸:0.
1
234<l◆
,
。<0.
1
244; 縦軸:1
.
5<
,
.
固I
S. 声一
・O で秩周JVI振馴 こ判 る分岐構 造 ・
城 軸:0.
1<声<0.
9; 縦軸:0<
h<4・
0・
C点は図 2のそれ と同 じ.
r
^
,
A
図 1
2. 共存す る異なる剛 V
I
振動 ・
p..p
.非 に 上
◆
,
o
=0・
1
09・
城 軸:0<l< 1
0
0 ;縦軸 :11
0<Y.<2・
5・
-5
5
6-
「カオスとその周辺」
1(
a
)
,
(
b
))
,この先 ,〃周期
周期振動 か らカオスへ到 る分 岐 は,周期倍化 に よる ものであ り(図 1
振動の領域 とカオスの領域 が P.う Cj
ph.1
→
C→ の ように交互 に現 われ てい る. これは
一種 の p
er
io
d
i
c
c
h
a
o
d
c的n
s
i
L
i
o
ns
e
q
u
e
n
∝ と考 え られる 【
81.
概周期 振動 か らの カオス
これまで は k
'
,
。を コン トロール '
バ ラメタ一 に選 んだ. 図 1
3はこ
'
r
。=0
.
1
2
9
2
4 に固定 し,
・
Pを区間 0.
1<β< 0
.
9 で変化 させ た ものである.
れを k
声
=Oの と
き x の運動 と Y の運動 は独 立であ り系全体 で は概周期 振動 になっている. 声?小 さい領域
を拡大 す る と,初 めの うちは概周期 的 であるが ,ある ところで周期 振動が現 われ (周 波数 固定 )
,
やがてカオス に到 る様子 が見 て取 れる.
4.考 察
図 11や図 1
3に見 られ る ように.方程 式 (
8
)の分 岐構造 は 3変数や 4変数 な どの場合 に比
べ てか な り複雑 である・ これは自由度 の多 い ことを反映 している と思 われる・ だが基本 的 に阜
系 は陽 イオ ン系 と陰 イオ ン系 とい う 2種類 の部分系 が電気的 に結合 した もの と考 えて よい.
/dE
=J:
・
(x) におい て f.
I
(X) が少 な くと
振動が安定 に存在す るため には,時間発展方程 式 dx.
も 3次 以上の非線形項すなわち x
.
I
)に比例す る項 を含 む ことが必 要 である 【9]. この非線形
項 は,膜 表両 でのイオ ンの出入 りの速度定数が-定借 でな く,イオ ン法度 に非線 形 に依存 す る
とい う仮定 によって取 り入 れ られた. この点 は直接 実験 によ り証明 されてはいない. だが,
To
k
o等 は,イオ ン濃度 の変化 によ り膜 のキ ャパ シタンスが急激 に変化 す る機構 を,同種 の仮 定
によって概念 的 に説明 している 【
1
0】.
振動の仕組 み は定牲的 には次 の ように理解 される :表面のイオ ン濃度 が低 い間は膜 内取 り
込みの速度 が遅 く,表面 にはイオ ンが どんどんたまる. イオ ンがある程度 た まって くると,膜
内取 り込みの速度 が急 に速 くなってイオ ンの流 れ も急 に増大す る. す る と表面 でのイオ ン強度
は一気 に減少す る. 振動 はこの ようなプロセスの放 り返 しである. ただ しこれだけでは高々周
期 的 な振動 しか得 られない. カオスや概周期振動 は今 の場合 ,2組 の系 の振動が結合す ること
によって実現 されている. ここでは簡単の ため陽イオ ン X と陰 イオ ン Y を-種類ずつ考 えた.
モデル をよ り現実 に近づ けるには (
Na
'と K'の両方 を想定 して )さらにもうー種類 の陽 イオン
Ⅹ●
を含 め ることが考 え られ る. モデル をこの ように拡張す ることは容易 であ り、その場合 も β
→ O の極 限で X,Ⅹ,Y の運動 は完全 に分離す る. この とき 二種類以上のイオ ンが振動 してい
れば ,系全体 の運動 は一般 に概周期的 であ り,したが って
移が可能 となる.
5. む す び
β を有限 に した とき,カオスへの転
孜 々は,自励発振す る膜のモデルについて計算機 シミュ レー シ ョンを行 ない,カオス を含 む
種 々の振動状態 を得 た. 生体系 の振動 を記述す るモデル方程 式の代 表例 としては, フアン ・デ
ル ●ポル振動子 【
11
】やホジキ ン ーハ クス レイ方程式 【
1
2
】な どがあるが ,それ らは周期 的 に振動
す る外力 の もとでカオス となる 【
1
3日 1
4】. 一方我 々が扱 ったの は,一定の法度勾 配のみ で駆動
され る自律系 の カオスである. ここで扱 った ような単純 なモデル シス テムの研 究が,脳波 な ど
複雑 な生物 の カオス的振舞 いについて,その メカニズ ムの基本 的 な部分 の解明 に役 立 つ ことを
期待 してい る.
【
文献】
t
l
】 松本元 : 「
神経興奮の現象 と実体 」 上 ,第 4章 (丸善 ,1
981年 ).
[
2
] K.
To
k
o
,
K Ry
u
,S.
Ez
a
k
ia
n
dK.
Ya
ma
f
u
j
i:
I
.
Ph
y
s
.
S
o
c
.
J
p
n
.
51
,3
3
9
8(
1
9
8
2
)
.
,
61
4(
1
98
9
)
.
【
3
】林健 司 ,都 甲潔 ,山藤者 :物性研 究 51
【
41 植村寿公 、石井淑夫 :日本化学会 第 5
6春季年会予稿集 p.
4
6
0(
1
9
8
8年 )
.
-
557 -
研究 会報告
[
5
] 花井哲也 : 「
痕 とイオ ン」,第 6章 (化学 同人 ,1
978年 )
.
【
6】 T.
Ka
wa
ku
b
o:Bi
o
ph
ys
.
Ch
e
m.2
3,2
2
9(
1
986
)
.
.
Pf
O
C
a
C
C
i
a‥
Phy
s
.
R占
V.
Le
l
t
.
50,
3
46(
1
983
)
.
[
7] P.
Gms
s
b
c
r
g
e
randI
L.Swi
n
n
e
y:Phy
s
.L
et
t
.85A,
9(
1
981
)
.
[
8J I
.S.
Tu
r
n
e
r
,
I
.
C.
Rou
x,W.D.Mc
Co
r
mi
c
ha
ndH.
dI
.
Pr
i
g
o
in
g
e:
Se
l
f
Or
gani
z
at
i
ot
li
nNone
qu
L
'
L
i
b
r
L
'
u
mS
ys
t
e
ms
,
c
h
a
p.3(Wi
l
e
y,
Ne
w
[
9] G.
Ni
c
ol
i
san
Yor
k,1
97
7
)
.
[
1
0] K.Tok
o.
J
.Ni
L
t
aa
n
dY.Ya
ma
f
u
j
i:
J
.
Phy
s
.So
c
.
J
p
n.
50,1
3
43(
1
981
)
.
l
.
Ma
g.
6,
7
63(
1
9
2
8
)
.
l
ll
】 B.v
a
nd
e
rPola
ndJ
.
v
a
r
L
d
e
rMa
r
k:
Ph
i
[
1
2] A.
L Hod
g
ki
na
n
dA.
F.Hu
x
l
e
y:J
.
Ph
ys
i
ol
.1
1
7,5
(
氾(
1
952
)
.
[
1
3] F.
A.
Ro
be
r
g
e,
R.
A.
Na
de
a
ua
n
dT.
N.
J
a
me
s'
.Ca
r
di
ov
a
s
c
.
Re
s
.
2,1
9(
1
96
8
)
.
[
1
4】 K.Ai
h
a
r
a
,
G.
Ma
t
s
u
mol
oa
n
dY.
I
k
e
ga
y
a:J
.
l
h
q
o
r
.
Bi
o
l
.1
09,
2
49(
1
98
4)
.
-5
5
8-
Fly UP