双安定な系のポテンシャルエネルギー グラフ。極小値は
x
1
{\displaystyle x_{1}}
と
x
2
{\displaystyle x_{2}}
の二つある。このように二つの「谷」を持つ曲面は双安定系になりうるもので、ボールは谷の位置でのみ安定に静止できる(図の①、②)。間にある極大値
x
3
{\displaystyle x_{3}}
の位置に置かれたボール③は平衡状態にあるが不安定で、わずかにでも擾乱を受けると安定点のどちらかに移ってしまう。
力学系 における双安定性 (そうあんていせい、英 : bistability )とはその系が二つの安定平衡状態を持つことを言う[ 1] 。双安定 (そうあんてい)な系は二つの状態のいずれかを取り続けることができる。双安定的なふるまいは機械的なリンク機構 、電子回路 、非線形光学系 、化学反応 、生理学的ないし生物学的システムなどで見られる。双安定な機械装置の例には照明スイッチ (英語版 ) がある。スイッチのレバーは「オン」位置か「オフ」位置のどちらかで静止するよう設計されており、中間では止まらない。
概要
保存力 場における双安定性は、安定平衡点であるポテンシャルエネルギー 極小点 が二つあることに由来する[ 2] 。二つの極小値は異なっていても構わない。粒子が静止し続けられるのはエネルギーの局所最低である極小位置のいずれかのみである。数学的な議論によると、二つの極小点の間には必ず不安定平衡点である極大点 が存在する。極大は二つの平衡位置を隔てる障壁と見られる。障壁を乗り越えるのに十分な励起エネルギーを与えられると、系は一方のエネルギー極小状態からもう一方に遷移することができる(化学系については活性化エネルギー とアレニウスの式 を参照せよ)。障壁位置に達した系は緩和時間 と呼ばれる時間を経てもう一方の極小状態へ緩和していく。
双安定系は本質的にヒステリシス を伴う。すなわち、系の出力はその時点での入力の強さだけでなく、過去の履歴によって系がどちらの状態を取っているかに依存する[ 3] 。
双安定性はバイナリ データを記憶するデジタル回路素子 で広く利用されており、コンピュータやある種の半導体メモリ の基本構成要素であるフリップフロップ 回路の本質的な特性でもある。双安定デバイスの一つの状態に「0」、もう一つの状態に「1」を割り当てることで1ビット のバイナリデータを格納できる。弛張発振器 (英語版 ) やマルチバイブレーター 、シュミットトリガ にも双安定性が利用されている。光双安定性 (英語版 ) はある種の光デバイス が持つ特性で、入力に応じて2つの共振伝送状態が安定になるというものである。生化学システムでも双安定性が発現することがあり、構成化学物質の濃度と活量によって0か1かのスイッチ的な出力が得られる。
数理モデル化
以下は力学系理論 で扱われるもっとも単純な双安定系の一つを数学的に表現したものである。
d
y
d
t
=
y
(
1
−
y
2
)
{\displaystyle {\frac {dy}{dt}}=y(1-y^{2})}
この系は
y
4
4
−
y
2
2
{\displaystyle {\frac {y^{4}}{4}}-{\frac {y^{2}}{2}}}
の形を持つ曲線上を転がり落ちるボールを表しており、3つの平衡点
y
=
1
{\displaystyle y=1}
、
y
=
0
{\displaystyle y=0}
、
y
=
−
1
{\displaystyle y=-1}
を持つ。真ん中の平衡点
y
=
0
{\displaystyle y=0}
は不安定、ほかの二点は安定である。時間とともに
y
(
t
)
{\displaystyle y(t)}
がどのように増減するかは初期条件
y
(
0
)
{\displaystyle y(0)}
に依存する。初期条件が正の場合(
y
(
0
)
>
0
{\displaystyle y(0)>0}
)解
y
(
t
)
{\displaystyle y(t)}
は時間とともに1に近づくが、初期条件が負の場合(
y
(
0
)
<
0
{\displaystyle y(0)<0}
)は−1に近づく。この意味で系のダイナミクスは双安定である。終状態は初期条件によって
y
=
1
{\displaystyle y=1}
か
y
=
−
1
{\displaystyle y=-1}
のどちらかになる[ 4] 。この双安定領域の出現については、分岐パラメータ
r
{\displaystyle r}
の値によって超臨界ピッチフォーク分岐 が起きるモデル系
d
y
d
t
=
y
(
r
−
y
2
)
{\displaystyle {\frac {dy}{dt}}=y(r-y^{2})}
を通して理解できる。
パラメータの特定範囲でのみ双安定性が発現する生化学システムも存在し、そのパラメータはフィードバックの強さと解釈されることが多い。典型的ないくつかの例では、パラメータの値が小さいときには一つの安定不動点 しか存在しない。パラメータが臨界値を超えるとサドルノード分岐 が起きて新たな不動点の対(一つは安定、もう一つは不安定)が生まれる。さらに増加すると別のサドルノード分岐により不安定解が最初の安定解と結合して消滅し、後で生まれた安定解のみが残ることがある。パラメータがそれらの臨界値の間にあるなら系は2つの安定解を持つことになる。そのような性質を持つ力学系の例には
d
x
d
t
=
r
+
x
5
1
+
x
5
−
x
{\displaystyle {\frac {\mathrm {d} x}{\mathrm {d} t}}=r+{\frac {x^{5}}{1+x^{5}}}-x}
がある。ここでは
x
{\displaystyle x}
が出力であり、パラメータの
r
{\displaystyle r}
が入力としてはたらく[ 5] 。
非線形な結合振動子にノイズを加えた系では、二つの安定なリミットサイクル の間を前後に飛び移る「モードホッピング」と呼ばれる不安定性が現れることがあり、そのポアンカレ断面 上では通常の双安定性と同様のふるまいが見られる[ 6] 。
生物学的・化学的システム
双安定な化学系は緩和速度論や非平衡熱力学 、確率共鳴 、気候変動 の分野とのかかわりで広く研究対象とされてきた[ 7] 。細胞周期 進行・細胞分化 [ 8] ・アポトーシス における意思決定プロセスという細胞機能の基本現象を理解する上でも双安定性は重要である。癌 発生やプリオン 病の初期段階にともなう細胞恒常性の喪失や、新しい種の発生(種分化 )ともかかわりがある[ 7] 。キイロショウジョウバエ の胚発生 を例に取ると、前後軸[ 9] や背腹軸[ 10] [ 11] の形成、ならびに眼の発生に双安定性が関与していることが報告されている[ 12] 。
生物学的・化学的なシステムが双安定性を持つには三つの必要条件を満たさなければならない。正フィードバック の存在、弱い刺激をフィルターするメカニズム、無限の出力増加を防ぐメカニズムである[ 7] 。空間的な広がりを持つ双安定なシステムでは局所相関の発生や進行波の出現が研究されている[ 13] [ 14] 。
「XがYを活性化させ、YがXを活性化させる」という単純な正フィードバックモチーフ (英語版 ) でも双安定性を作り出せるが、実際の細胞シグナル伝達では複数のフィードバックループが組み合わされてスイッチを構成し、重要な調節ステップの役割を担っている[ 15] 。これまでの研究により、Xenopus (ツメガエル属)の卵母細胞 の成熟[ 16] 、哺乳類のカルシウムシグナル 伝達、出芽酵母 の極性形成など多くの生物学的システムに、正のテンポラル・フィードバック (英語版 ) (速いループと遅いループの組み合わせ)、もしくは別タイミングで発動する複数のフィードバックループが組み込まれていることが分かっている。速さの異なるフィードバックの組み合わせには活性化時間と不活性化時間を別々に調節したりノイズへの過敏な反応を抑えたりといった利点がある[ 15] 。
反応系が反応活性因子と阻害因子の両者にフィードバックを行うと、反応物濃度の大きな変化に耐えるロバストな双安定スイッチを実現できる。細胞生物学においては、細胞周期を有糸分裂 の段階に進ませる役割を持つサイクリン依存性キナーゼ1 (CDK1) が活性化されると自身の活性因子 Cdc25 を活性化し、同時に不活性因子 Wee1 を不活性化する例がある。この二重フィードバックがなくても系は双安定だが、それほど広い範囲の濃度には耐えられないと考えられている[ 17] 。
分泌シグナル伝達分子の一つソニック・ヘッジホッグ (Shh) のシグナル伝達ネットワークは鋭敏な感受性を持つ正負のフィードバックループの組が双安定性を生み出している好例である。Shhシグナル伝達ネットワークは肢芽組織の分化のパターニングなど多くの発生プロセスで双安定スイッチとしてはたらいており、細胞は決められたShh濃度で敏速に状態を切り替える。このシグナルによって転写 が活性化された gli1 および gli2 の遺伝子産物は、gli1 , gli2 自身の発現をいっそう活発にするとともに、Shhシグナル伝達経路の下流にある標的遺伝子の転写活性化因子としてもはたらく[ 18] 。この正フィードバックと並行して、Gli 転写因子が転写抑制因子の一つ (Ptc) の転写を促進するという負のフィードバックループがShhシグナル伝達ネットワークを調節している。
機械的システム
ラチェットの動作。歯車のそれぞれの歯と溝が一つの安定状態に当たる。
機械的システムの設計に組み込まれた双安定性は一般に「オーバーセンター (英 : over-center )」と呼ばれることがある。機構を動かして部品がピーク位置を越えると(オーバーセンターすると)次の安定位置に移行し、変位が足りなければ最初の安定位置に戻る。このような機構はトグルスイッチとして機能する。身近な例として、壁に取り付ける一般的な電気スイッチがある。スイッチレバーは「オン」か「オフ」の位置にのみ静止するよう設計されており、その中間を越えて動かすと切り替わる。
ばね は双安定機構を実現する手段として一般的である。二つの安定位置を持つラチェット 構造にばねを組み込んだボタン 機構はノック式ボールペン の多くに採用されている。回転ラチェット (歯車と歯止め)は精巧な多安定機構であり、不可逆運動を作るために使われる。歯車を順方向に回していくにつれて歯止めは一つの歯を乗り越えて次の安定状態に移る。
柔軟な梁 が複数の安定な座屈 状態を持つことを利用した双安定コンプライアント機構 (英語版 ) は、がたつきや摩擦 がなく、部品構成が単純で小型化に向いていることからマイクロサイズの機構 設計で広く用いられている[ 19] [ 20] 。
関連項目
脚注
^ Morris, Christopher G. (1992). Academic Press Dictionary of Science and Technology . Gulf Professional publishing. pp. 267. ISBN 978-0122004001 . https://books.google.com/books?id=nauWlPTBcjIC&pg=PA267&dq=bistable+bistability
^ Nazarov, Yuli V.; Danon, Jeroen (2013). Advanced Quantum Mechanics: A Practical Guide . Cambridge University Press. pp. 291. ISBN 978-1139619028 . https://books.google.com/books?id=w20gAwAAQBAJ&pg=PA291&dq=bistability+minimum
^ Tao Cui, Yong Wang, and Daren Yu (2014). “Bistability and Hysteresis in a Nonlinear Dynamic Model of Shock Motion”. Journal of Aircraft 51 (5): 1373-1379. doi :10.2514/1.C032175 .
^ Ket Hing Chong; Sandhya Samarasinghe; Don Kulasiri; Jie Zheng (2015). “Computational techniques in mathematical modelling of biological switches”. MODSIM2015 : 578–584. For detailed techniques of mathematical modelling of bistability, see the tutorial by Chong et al. (2015) http://www.mssanz.org.au/modsim2015/C2/chong.pdf The tutorial provides a simple example illustration of bistability using a synthetic toggle switch proposed in Collins, James J.; Gardner, Timothy S.; Cantor, Charles R. (2000). “Construction of a genetic toggle switch in Escherichia coli”. Nature 403 (6767): 339–42. Bibcode : 2000Natur.403..339G . doi :10.1038/35002131 . PMID 10659857 . . The tutorial also uses the dynamical system software XPPAUT http://www.math.pitt.edu/~bard/xpp/xpp.html to show practically how to see bistability captured by a saddle-node bifurcation diagram and the hysteresis behaviours when the bifurcation parameter is increased or decreased slowly over the tipping points and a protein gets turned 'On' or turned 'Off'.
^ Angeli, David; Ferrell, JE; Sontag, Eduardo D (2003). “Detection of multistability, bifurcations, and hysteresis in a large calss of biological positive-feedback systems” . PNAS 101 (7): 1822–7. Bibcode : 2004PNAS..101.1822A . doi :10.1073/pnas.0308265100 . PMC 357011 . PMID 14766974 . https://www.ncbi.nlm.nih.gov/pmc/articles/PMC357011/ .
^ Mathias L. Heltberg, Sandeep Krishna, and Mogens H. Jensen (2017). “Time Correlations in Mode Hopping of Coupled Oscillators”. Journal of Statistical Physics 167 : 792-805. doi :10.1007/s10955-017-1750-x .
^ a b c Wilhelm, T (2009). “The smallest chemical reaction system with bistability” . BMC Systems Biology 3 : 90. doi :10.1186/1752-0509-3-90 . PMC 2749052 . PMID 19737387 . https://www.ncbi.nlm.nih.gov/pmc/articles/PMC2749052/ .
^ “Multistable switches and their role in cellular differentiation networks” . BMC Bioinformatics 15 : S7+. (2014). doi :10.1186/1471-2105-15-s7-s7 . PMC 4110729 . PMID 25078021 . https://www.ncbi.nlm.nih.gov/pmc/articles/PMC4110729/ .
^ Lopes, Francisco J. P.; Vieira, Fernando M. C.; Holloway, David M.; Bisch, Paulo M.; Spirov, Alexander V.; Ohler, Uwe (26 September 2008). “Spatial Bistability Generates hunchback Expression Sharpness in the Drosophila Embryo” . PLOS Computational Biology 4 (9): e1000184. Bibcode : 2008PLSCB...4E0184L . doi :10.1371/journal.pcbi.1000184 . PMC 2527687 . PMID 18818726 . https://www.ncbi.nlm.nih.gov/pmc/articles/PMC2527687/ .
^ Wang, Yu-Chiun; Ferguson, Edwin L. (10 March 2005). “Spatial bistability of Dpp–receptor interactions during Drosophila dorsal–ventral patterning”. Nature 434 (7030): 229–234. Bibcode : 2005Natur.434..229W . doi :10.1038/nature03318 . PMID 15759004 .
^ Umulis, D. M.; Mihaela Serpe; Michael B. O’Connor; Hans G. Othmer (1 August 2006). “Robust, bistable patterning of the dorsal surface of the Drosophila embryo” . Proceedings of the National Academy of Sciences 103 (31): 11613–11618. Bibcode : 2006PNAS..10311613U . doi :10.1073/pnas.0510398103 . PMC 1544218 . PMID 16864795 . https://www.ncbi.nlm.nih.gov/pmc/articles/PMC1544218/ .
^ Graham, T. G. W.; Tabei, S. M. A.; Dinner, A. R.; Rebay, I. (22 June 2010). “Modeling bistable cell-fate choices in the Drosophila eye: qualitative and quantitative perspectives” . Development 137 (14): 2265–2278. doi :10.1242/dev.044826 . PMC 2889600 . PMID 20570936 . https://www.ncbi.nlm.nih.gov/pmc/articles/PMC2889600/ .
^ Elf, J.; Ehrenberg, M. (2004). “Spontaneous separation of bi-stable biochemical systems into spatial domains of opposite phases”. Systems Biology 1 (2): 230–236. doi :10.1049/sb:20045021 . PMID 17051695 .
^ Kochanczyk, M.; Jaruszewicz, J.; Lipniacki, T. (Jul 2013). “Stochastic transitions in a bistable reaction system on the membrane” . Journal of the Royal Society Interface 10 (84): 20130151. doi :10.1098/rsif.2013.0151 . PMC 3673150 . PMID 23635492 . https://www.ncbi.nlm.nih.gov/pmc/articles/PMC3673150/ .
^ a b Onn Brandman, James E Ferrell Jr, Rong Li, Tobias Meyer (2005). “Interlinked fast and slow positive feedback loops drive reliable cell decisions”. Science 310 (5747): 496–498. doi :10.1126/science.1113834 .
^ J. E. Ferrell Jr.; Machleder EM (1998). “The biochemical basis of an all-or-none cell fate switch in Xenopus oocytes.”. Science 280 (5365): 895–8. Bibcode : 1998Sci...280..895F . doi :10.1126/science.280.5365.895 . PMID 9572732 .
^ Ferrell JE Jr. (2008). “Feedback regulation of opposing enzymes generates robust, all-or-none bistable responses” . Current Biology 18 (6): R244–R245. doi :10.1016/j.cub.2008.02.035 . PMC 2832910 . PMID 18364225 . https://www.ncbi.nlm.nih.gov/pmc/articles/PMC2832910/ .
^ Lai, K., M.J. Robertson, and D.V. Schaffer, The sonic hedgehog signaling system as a bistable genetic switch. Biophys J, 2004. 86(5): p. 2748-57.
^ A.G. Dunning, N. Tolou, P.J. Pluimers, L.F. Kluit, and J.L. Herder (2012). “Bi-stable Compliant Mechanisms: Correction for Finite Element Modeling, Tuning the Stiffness and Preloading Incorporation”. Journal of Mechanical Design 134 (8): 084502. doi :10.1115/1.4006961 .
^ 小塚裕明他「コンプライアントメカニズムを用いた円弧ばね関節を有する ばね―パラレルメカニズムによる精密位置決め装置の開発」『日本機械学会論文集(C編)』第78巻第793号、2012年、3216-3226頁、doi :10.1299/kikaic.78.3216 。
外部リンク