反応速度式
出典: フリー百科事典『ウィキペディア(Wikipedia)』 (2023/01/12 20:44 UTC 版)
化学反応 の反応速度式 (はんのうそくどしき、英語 : rate equation)あるいは速度式 (rate law)[1] とは、反応速度 と反応物 の濃度 または圧力 および定数パラメーター(主に反応速度定数 と反応次数 )の関係式である[2] 。多くの反応では、反応速度rは次のような指数関数 で与えられる。
r
=
k
[
A
]
x
[
B
]
y
{\displaystyle r\;=\;k[\mathrm {A} ]^{x}[\mathrm {B} ]^{y}}
AとBの濃度の、平衡に達するまでの時間変化(A
0 = 0.25 mol/l)。正反応の速度定数k
f = 2 min
−1 、逆反応の速度定数k
r = 1 min
−1
簡単な例
以下のように、2つの化学種 の間に平衡が成立しているとする。
A
↽
−
−
⇀
B
{\displaystyle {\ce {A <=> B}}}
反応はt =0でAの初期濃度が
[
A
]
0
{\displaystyle {\ce {[A]_0}}}
、Bの初期濃度が0の状態から始まる。
このとき、平衡定数Kは以下のように書ける。
K
=
d
e
f
k
f
k
b
=
[
B
]
e
[
A
]
e
{\displaystyle K\ {\stackrel {\mathrm {def} }{=}}\ {\frac {k_{f}}{k_{b}}}={\frac {\left[{\ce {B}}\right]_{e}}{\left[{\ce {A}}\right]_{e}}}}
ここで、
[
A
]
e
{\displaystyle [{\ce {A}}]_{e}}
と
[
B
]
e
{\displaystyle [{\ce {B}}]_{e}}
は平衡状態でのAとBの濃度である。
時刻tにおけるAの濃度を
[
A
]
t
{\displaystyle [{\ce {A}}]_{t}}
、Bの濃度を
[
B
]
t
{\displaystyle [{\ce {B}}]_{t}}
とすると、両者は次の平衡反応の等式を満たす。
[
A
]
t
=
[
A
]
0
−
[
B
]
t
{\displaystyle {\ce {[A]_{\mathit {t}}=[A]_{0}-[B]_{\mathit {t}}}}}
ここで、
[
B
]
0
{\displaystyle {\ce {[B]_0}}}
は0であることに注意する。
これは、時刻が無限大となり、平衡に達した状態でも成立する。
[
A
]
e
=
[
A
]
0
−
[
B
]
e
{\displaystyle {\ce {[A]_{\mathit {e}}=[A]_{0}-[B]_{\mathit {e}}}}}
これは平衡定数Kの定義より、
[
B
]
e
=
x
=
k
f
k
f
+
k
b
[
A
]
0
{\displaystyle [{\ce {B}}]_{e}=x={\frac {k_{f}}{k_{f}+k_{b}}}[{\ce {A}}]_{0}}
ゆえに、
[
A
]
e
=
[
A
]
0
−
x
=
k
b
k
f
+
k
b
[
A
]
0
{\displaystyle \ [{\ce {A}}]_{e}=[{\ce {A}}]_{0}-x={\frac {k_{b}}{k_{f}+k_{b}}}[{\ce {A}}]_{0}}
これらの等式により、微分方程式 を解かずともAの濃度を求めることができる。
反応速度式は以下のように与えられる。
r
=
k
1
[
A
]
s
[
B
]
t
−
k
2
[
X
]
u
[
Y
]
v
{\displaystyle r={k_{1}[{\ce {A}}]^{s}[{\ce {B}}]^{t}}-{k_{2}[{\ce {X}}]^{u}[{\ce {Y}}]^{v}}\,}
−
d
[
A
]
d
t
=
k
f
[
A
]
t
−
k
b
[
B
]
t
{\displaystyle -{\frac {d[{\ce {A}}]}{dt}}={k_{f}[{\ce {A}}]_{t}}-{k_{b}[{\ce {B}}]_{t}}\,}
微分係数が負なのは正反応がAからBに変わる反応なので、Aの濃度は減少しているからである。簡略化するため、時刻tでのAの濃度
[
A
]
t
{\displaystyle [{\ce {A}}]_{t}}
をxとおく。また平衡時のAの濃度を
x
e
{\displaystyle x_{e}}
とする。このとき、
−
d
[
A
]
d
t
=
k
f
[
A
]
t
−
k
b
[
B
]
t
−
d
x
d
t
=
k
f
x
−
k
b
[
B
]
t
=
k
f
x
−
k
b
(
[
A
]
0
−
x
)
=
(
k
f
+
k
b
)
x
−
k
b
[
A
]
0
{\displaystyle {\begin{aligned}-{\frac {d[{\ce {A}}]}{dt}}&={k_{f}[{\ce {A}}]_{t}}-{k_{b}[{\ce {B}}]_{t}}\\-{\frac {dx}{dt}}&={k_{f}x}-{k_{b}[{\ce {B}}]_{t}}\\&={k_{f}x}-{k_{b}([{\ce {A}}]_{0}-x)}\\&={(k_{f}+k_{b})x}-{k_{b}[{\ce {A}}]_{0}}\end{aligned}}}
だから
k
f
+
k
b
=
k
b
[
A
]
0
x
e
{\displaystyle k_{f}+k_{b}={k_{b}{\frac {{\ce {[A]_0}}}{x_{e}}}}}
よって反応速度 は
d
x
d
t
=
k
b
[
A
]
0
x
e
(
x
e
−
x
)
{\displaystyle \ {\frac {dx}{dt}}={\frac {k_{b}[{\ce {A}}]_{0}}{x_{e}}}(x_{e}-x)}
つまり、
ln
(
[
A
]
0
−
[
A
]
e
[
A
]
t
−
[
A
]
e
)
=
(
k
f
+
k
b
)
t
{\displaystyle \ln \left({\frac {[{\ce {A}}]_{0}-[{\ce {A}}]_{e}}{[{\ce {A}}]_{t}-[{\ce {A}}]_{e}}}\right)=(k_{f}+k_{b})t}
という結果になる[7] 。
例の一般化
t=0での濃度が上と異なる場合、上式のような簡略化は使えず、微分方程式を解くことが必要になる。しかし、その微分方程式は解くことができ、その解は以下のように一般化したものとなる。
[
A
]
=
[
A
]
0
1
k
f
+
k
b
(
k
b
+
k
f
e
−
(
k
f
+
k
b
)
t
)
+
[
B
]
0
k
b
k
f
+
k
b
(
1
−
e
−
(
k
f
+
k
b
)
t
)
{\displaystyle \left[{\ce {A}}\right]=\left[{\ce {A}}\right]_{0}{\frac {1}{k_{f}+k_{b}}}\left(k_{b}+k_{f}e^{-\left(k_{f}+k_{b}\right)t}\right)+\left[{\ce {B}}\right]_{0}{\frac {k_{b}}{k_{f}+k_{b}}}\left(1-e^{-\left(k_{f}+k_{b}\right)t}\right)}
[
B
]
=
[
A
]
0
k
f
k
f
+
k
b
(
1
−
e
−
(
k
f
+
k
b
)
t
)
+
[
B
]
0
1
k
f
+
k
b
(
k
f
+
k
b
e
−
(
k
f
+
k
b
)
t
)
{\displaystyle \left[{\ce {B}}\right]=\left[{\ce {A}}\right]_{0}{\frac {k_{f}}{k_{f}+k_{b}}}\left(1-e^{-\left(k_{f}+k_{b}\right)t}\right)+\left[{\ce {B}}\right]_{0}{\frac {1}{k_{f}+k_{b}}}\left(k_{f}+k_{b}e^{-\left(k_{f}+k_{b}\right)t}\right)}
平衡定数が温度によらず一定に近く、反応速度がとても速い場合、例えば分子の立体配座異性体同士の平衡 の分析では、反応速度を求めるのには別の方法が必要になる。それは例えば、核磁気共鳴分光法 などである。
連鎖反応
反応
A
⟶
B
⟶
C
{\displaystyle {\ce {A -> B -> C}}}
についてそれぞれの反応の速度定数が
k
1
{\displaystyle k_{1}}
と
k
2
{\displaystyle k_{2}}
であるとき、それぞれの物質の時間当たりの変化量は以下のようになる。
反応物A:
d
[
A
]
d
t
=
−
k
1
[
A
]
{\displaystyle {\frac {d[{\ce {A}}]}{dt}}=-k_{1}[{\ce {A}}]}
反応物B:
d
[
B
]
d
t
=
k
1
[
A
]
−
k
2
[
B
]
{\displaystyle {\frac {d[{\ce {B}}]}{dt}}=k_{1}[{\ce {A}}]-k_{2}[{\ce {B}}]}
反応物C:
d
[
C
]
d
t
=
k
2
[
B
]
{\displaystyle {\frac {d[{\ce {C}}]}{dt}}=k_{2}[{\ce {B}}]}
それぞれの濃度が反応物全体の物質量で測られる場合、これらのような線形微分方程式はマスター方程式 として計算される。その微分方程式は解析的に解くことができ、解は以下のようになる。
[
A
]
=
[
A
]
0
e
−
k
1
t
{\displaystyle [{\ce {A}}]=[{\ce {A}}]_{0}e^{-k_{1}t}}
[
B
]
=
{
[
A
]
0
k
1
k
2
−
k
1
(
e
−
k
1
t
−
e
−
k
2
t
)
+
[
B
]
0
e
−
k
2
t
k
1
≠
k
2
[
A
]
0
k
1
t
e
−
k
1
t
+
[
B
]
0
e
−
k
1
t
otherwise
{\displaystyle \left[{\ce {B}}\right]={\begin{cases}\left[{\ce {A}}\right]_{0}{\dfrac {k_{1}}{k_{2}-k_{1}}}\left(e^{-k_{1}t}-e^{-k_{2}t}\right)+\left[{\ce {B}}\right]_{0}e^{-k_{2}t}&k_{1}\neq k_{2}\\\left[{\ce {A}}\right]_{0}k_{1}te^{-k_{1}t}+\left[{\ce {B}}\right]_{0}e^{-k_{1}t}&{\text{otherwise}}\\\end{cases}}}
[
C
]
=
{
[
A
]
0
(
1
+
k
1
e
−
k
2
t
−
k
2
e
−
k
1
t
k
2
−
k
1
)
+
[
B
]
0
(
1
−
e
−
k
2
t
)
+
[
C
]
0
k
1
≠
k
2
[
A
]
0
(
1
−
e
−
k
1
t
−
k
1
t
e
−
k
1
t
)
+
[
B
]
0
(
1
−
e
−
k
1
t
)
+
[
C
]
0
otherwise
{\displaystyle \left[{\ce {C}}\right]={\begin{cases}\left[{\ce {A}}\right]_{0}\left(1+{\dfrac {k_{1}e^{-k_{2}t}-k_{2}e^{-k_{1}t}}{k_{2}-k_{1}}}\right)+\left[{\ce {B}}\right]_{0}\left(1-e^{-k_{2}t}\right)+\left[{\ce {C}}\right]_{0}&k_{1}\neq k_{2}\\\left[{\ce {A}}\right]_{0}\left(1-e^{-k_{1}t}-k_{1}te^{-k_{1}t}\right)+\left[{\ce {B}}\right]_{0}\left(1-e^{-k_{1}t}\right)+\left[{\ce {C}}\right]_{0}&{\text{otherwise}}\\\end{cases}}}
この方程式は定常状態近似によって簡単に解けるようになっている。
並行反応・競合反応
1つの物質から2種類の生成物が生まれる場合、並行反応または競合反応が起こっている。
2つの一次反応が起こっている場合
反応
A
⟶
B
{\displaystyle {\ce {A -> B}}}
と
A
⟶
C
{\displaystyle {\ce {A -> C}}}
の速度定数がそれぞれ
k
1
{\displaystyle k_{1}}
と
k
2
{\displaystyle k_{2}}
であるとする。この時それぞれの濃度の時間変化の式は
−
d
[
A
]
d
t
=
(
k
1
+
k
2
)
[
A
]
{\displaystyle -{\frac {d[{\ce {A}}]}{dt}}=(k_{1}+k_{2})[{\ce {A}}]}
d
[
B
]
d
t
=
k
1
[
A
]
{\displaystyle {\frac {d[{\ce {B}}]}{dt}}=k_{1}[{\ce {A}}]}
d
[
C
]
d
t
=
k
2
[
A
]
{\displaystyle {\frac {d[{\ce {C}}]}{dt}}=k_{2}[{\ce {A}}]}
と表される。
したがって、積分形の反応速度式は
[
A
]
=
[
A
]
0
e
−
(
k
1
+
k
2
)
t
{\displaystyle \ [{\ce {A}}]=[A]_{0}e^{-(k_{1}+k_{2})t}}
[
B
]
=
k
1
k
1
+
k
2
[
A
]
0
(
1
−
e
−
(
k
1
+
k
2
)
t
)
{\displaystyle [{\ce {B}}]={\frac {k_{1}}{k_{1}+k_{2}}}[{\ce {A}}]_{0}(1-e^{-(k_{1}+k_{2})t})}
[
C
]
=
k
2
k
1
+
k
2
[
A
]
0
(
1
−
e
−
(
k
1
+
k
2
)
t
)
{\displaystyle [{\ce {C}}]={\frac {k_{2}}{k_{1}+k_{2}}}[{\ce {A}}]_{0}(1-e^{-(k_{1}+k_{2})t})}
と表される。
この場合、
[
B
]
[
C
]
=
k
1
k
2
{\displaystyle {\frac {{\ce {[B]}}}{{\ce {[C]}}}}={\frac {k_{1}}{k_{2}}}}
が重要な関係式となる。
一次反応と二次反応が1つずつ起こっている場合
[8]
これは2分子による反応と、擬一次反応とみなせる加水分解が同時に起こっている場合に適用できる。並行反応によって反応物が一部消費されるため、加水分解の反応速度を調べるのは難しい。例えば、AとRが反応してCが生成するが、同時に加水分解が進行してAがBに変わると言った場合である。反応式 で表せば、
A
+
H
2
O
⟶
B
{\displaystyle {\ce {{A}+ H2O -> B}}}
と
A
+
R
⟶
C
{\displaystyle {\ce {{A}+ R -> C}}}
となる。
反応速度式は以下のとおりになる。
d
[
B
]
d
t
=
k
1
[
A
]
[
H
2
O
]
=
k
1
′
[
A
]
{\displaystyle {\frac {d[{\ce {B}}]}{dt}}=k_{1}{\ce {[A][H2O]}}=k_{1}'[{\ce {A}}]}
d
[
C
]
d
t
=
k
2
[
A
]
[
R
]
{\displaystyle {\frac {d[{\ce {C}}]}{dt}}=k_{2}{\ce {[A][R]}}}
ただし
k
1
′
{\displaystyle k_{1}'}
は擬一次速度定数である。
主な生成物Cの濃度について積分すると、以下の式が得られる。
[
C
]
=
[
R
]
0
[
1
−
e
−
k
2
k
1
′
[
A
]
0
(
1
−
e
−
k
1
′
t
)
]
{\displaystyle [C]=[R]_{0}\left[1-e^{-{\frac {k_{2}}{k_{1}'}}[{\ce {A}}]_{0}(1-e^{-k_{1}'t})}\right]}
これは
l
n
[
R
]
0
[
R
]
0
−
[
C
]
=
k
2
[
A
]
0
k
1
′
(
1
−
e
−
k
1
′
t
)
{\displaystyle ln{\frac {\ce {[R]_{0}}}{\ce {[R]_{0}-[C]}}}={\frac {k_{2}[{\ce {A}}]_{0}}{k_{1}'}}(1-e^{-k_{1}'t})}
と等価である。
[B]と[C]の濃度の関係は次のようになっている。
[
B
]
=
−
k
1
′
k
2
l
n
(
1
−
[
C
]
v
[
R
]
0
)
{\displaystyle [{\ce {B}}]=-{\frac {k_{1}'}{k_{2}}}ln\left(1-{\frac {\ce {[C]}}{v}}{[R]_{0}}\right)}
これは解析的に得られた解であるが、次の近似が用いられている。
[
A
]
0
−
[
C
]
≈
[
A
]
0
{\displaystyle {\ce {{[A]_{0}}-[C]\approx [A]_{0}}}}
そのため、前の式における[C]は[C]が[A]0 に比べ非常に小さい時のみ使うことができる。
反応の量的ネットワーク
化学反応ネットワーク理論(英語版 ) の最も一般的な考え方は、
R
{\displaystyle R}
個の反応に関わる異なる化学種の数
N
{\displaystyle N}
を考えることである[9]
[10] 。一般に、
j
{\displaystyle j}
番目の反応について次のように記述できる。
s
1
j
X
1
+
s
2
j
X
2
…
+
s
N
j
X
N
→
k
j
r
1
j
X
1
+
r
2
j
X
2
+
…
+
r
N
j
X
N
,
{\displaystyle s_{1j}{\ce {X}}_{1}+s_{2j}{\ce {X}}_{2}\ldots +s_{Nj}{\ce {X}}_{N}{\ce {->[k_{j}]}}\ r_{1j}{\ce {X}}_{1}+\ r_{2j}{\ce {X}}_{2}+\ldots +r_{Nj}{\ce {X}}_{N},}
これは、上式と同値な下式で表されることも多い。
∑
i
=
1
N
s
i
j
X
i
→
k
j
∑
i
=
1
N
r
i
j
X
i
.
{\displaystyle \sum _{i=1}^{N}s_{ij}{\ce {X}}_{i}{\ce {->[k_{j}]}}\sum _{i=1}^{N}\ r_{ij}{\ce {X}}_{i}.}
ここで
j
{\displaystyle j}
は1から
R
{\displaystyle R}
までの反応の番号。
X
i
{\displaystyle {\ce {X}}_{i}}
は
i
{\displaystyle i}
番目の化学種。
k
j
{\displaystyle k_{j}}
は
j
{\displaystyle j}
番目の反応の速度定数
s
i
j
{\displaystyle s_{ij}}
と
r
i
j
{\displaystyle r_{ij}}
は反応式における反応物および生成物の係数
この反応の反応速度は化学平衡の法則(英語版 ) から推測される。
f
j
(
[
X
→
]
)
=
k
j
∏
z
=
1
N
[
X
z
]
s
z
j
{\displaystyle f_{j}([{\vec {\ce {X}}}])=k_{j}\prod _{z=1}^{N}[{\ce {X}}_{z}]^{s_{zj}}}
これは単位時間・単位体積あたりの物質の変化量で表される。ここで、
[
X
→
]
=
(
[
X
1
]
,
[
X
2
]
,
⋅
⋅
⋅
,
[
X
N
]
)
{\displaystyle {\ce {[{\vec {X}}]=([X1],[X2],...,[X_{\mathit {N}}])}}}
は濃度のベクトル [要曖昧さ回避 ] である。ここで、この式が定義される反応は素反応である事に注意する。
零次反応
s
z
j
=
0
{\displaystyle s_{zj}=0}
が全ての
z
{\displaystyle z}
について成り立つ。
一次反応
s
z
j
=
1
{\displaystyle s_{zj}=1}
がある1つの
z
{\displaystyle z}
について成り立つ。
二次反応
2分子の反応では2つの
z
{\displaystyle z}
について
s
z
j
=
1
{\displaystyle s_{zj}=1}
が成り立つ。また二量化 では
s
z
j
=
2
{\displaystyle s_{zj}=2}
がある1つの
z
{\displaystyle z}
について成り立つ。
それぞれについて次のように議論される。この時、それぞれの反応について反応の量的関係についての行列(stoichiometric matrix )を定義することができる。
S
i
j
=
r
i
j
−
s
i
j
,
{\displaystyle S_{ij}=r_{ij}-s_{ij},}
これは
j
{\displaystyle j}
番目の反応について存在する正味の
i
{\displaystyle i}
の物質量を表す。この時、反応速度式は以下のようなより一般的な形に書き直すことができる。
d
[
X
i
]
d
t
=
∑
j
=
1
R
S
i
j
f
j
(
[
X
→
]
)
.
{\displaystyle {\frac {d[{\ce {X}}_{i}]}{dt}}=\sum _{j=1}^{R}S_{ij}f_{j}([{\vec {\ce {X}}}]).}
ここで、これは反応の量的関係を表す行列と反応速度の関数の積である事に注意する。
系内で起こっている反応が可逆反応のみであり、反応が平衡状態にある場合、この方程式には簡単な解が存在する。(
d
[
X
i
]
d
t
=
0
{\displaystyle {\frac {d[{\ce {X}}_{i}]}{dt}}=0}
)この場合性反応と逆反応の反応速度は等しいので、詳細釣り合い が成り立っている。ただし、詳細釣り合いは反応の量的行列
S
i
j
{\displaystyle S_{ij}}
のみについて成り立つ性質であり、反応速度関数
f
j
{\displaystyle f_{j}}
には依存しない。詳細釣り合いが成り立たない場合については、代謝経路 を理解するために開発された流速均衡解析(英語版 ) によって研究されている。[11] [12]
一般的な1分子の変換反応についての動力学
一般的に、1分子が
N
{\displaystyle N}
種類の化学種に変換される反応について、時刻
t
{\displaystyle t}
での化学種1 - Nの濃度を
X
1
(
t
)
{\displaystyle X_{1}(t)}
through
X
N
(
t
)
{\displaystyle X_{N}(t)}
とおくと、各時刻ごとの各化学種の濃度が分かる。ここで、
X
i
{\displaystyle X_{i}}
から
X
j
{\displaystyle X_{j}}
に変わる反応の速度定数を
k
i
j
{\displaystyle k_{ij}}
とおき、
k
i
j
{\displaystyle k_{ij}}
などのそれぞれの反応の速度定数を成分とする行列
K
{\displaystyle K}
を作る。
また、時間の関数として濃度のベクトル
X
(
t
)
=
(
X
1
(
t
)
,
X
2
(
t
)
,
.
.
.
,
X
N
(
t
)
)
T
{\displaystyle X(t)=(X_{1}(t),X_{2}(t),...,X_{N}(t))^{T}}
をおく。
そして、ベクトル
J
=
(
1
,
1
,
1
,
.
.
.
,
1
)
T
{\displaystyle J=(1,1,1,...,1)^{T}}
をおく。
さらに、
I
{\displaystyle I}
をN次の単位行列 とする。
また、
D
i
a
g
{\displaystyle Diag}
を関数とする。ただしこの関数は対角行列 を作り、その対角線上の成分があるベクトルの成分となっているものとする。
そして、
L
−
1
{\displaystyle \displaystyle {\mathcal {L}}^{-1}}
は
s
{\displaystyle s}
から
t
{\displaystyle t}
への逆ラプラス変換 (英語版 )とする。
この時
X
(
t
)
{\displaystyle X(t)}
は
X
(
t
)
=
L
−
1
[
(
s
I
+
D
i
a
g
(
K
J
)
−
K
T
)
−
1
X
(
0
)
]
{\displaystyle X(t)=\displaystyle {\mathcal {L}}^{-1}[(sI+Diag(KJ)-K^{T})^{-1}X(0)]}
,
となる。
このようにして、初期状態と時刻
t
{\displaystyle t}
での状態の関係が示される。
脚注
^ Peter Atkins 、Loretta Jones、Leroy Laverman 『アトキンス一般化学(下)』渡辺正 訳、東京化学同人 、2011年。ISBN 9784807908554 。
^ IUPAC Gold Book definition of rate law IUPAC のCompendium of Chemical Terminology(英語版 ) も参照
^ Peter Atkins 、Julio de Paula、Ronald Friedman 『アトキンス基礎物理化学(下)―分子論的アプローチ―』千原 秀昭・稲葉 章 訳、東京化学同人 、2011年。ISBN 9784807907519 。
^ a b
『アトキンス基礎物理化学(下)―分子論的アプローチ』p.661
^ Walsh R, Martin E, Darvesh S. A method to describe enzyme-catalyzed reactions by combining steady state and time course enzyme kinetic parameters... Biochim Biophys Acta.(英語版 ) 2010 Jan;1800:1-5
^ a b c NDRL Radiation Chemistry Data Center [リンク切れ ]
。次のページも参照のこと。Christos Capellos and Bennon H. Bielski "Kinetic systems: mathematical description of chemical kinetics in solution" 1972, Wiley-Interscience (New York) Archived 2013年4月14日, at Archive.is
^ これを解いた例としては次のような文献がある。Determination of the Rotational Barrier for Kinetically Stable Conformational Isomers via NMR and 2D TLC An Introductory Organic Chemistry Experiment Gregory T. Rushton, William G. Burns, Judi M. Lavin, Yong S. Chong, Perry Pellechia, and Ken D. Shimizu J. Chem. Educ. 2007 , 84, 1499. Abstract
^ José A. Manso et al."A Kinetic Approach to the Alkylating Potential of Carcinogenic Lactones" Chem. Res. Toxicol.(英語版 ) 2005, 18, (7) 1161-1166
^ Heinrich, R. and Schuster, S. (1996) The regulation of cellular systems.チャップマン・アンド・ホール(英語版 ) 、ニューヨーク
^ Chen, L. and Wang, R. and Li, C. and Aihara, K. (2010) Modeling biomolecular networks in cells: structures and dynamics.シュプリンガー・サイエンス・アンド・ビジネス・メディア
^ Szallasi, Z. and Stelling, J. and Periwal, V. (2006) System modeling in cell biology: from concepts to nuts and bolts. MIT Press Cambridge.
^ Iglesias, P.A. and Ingalls, B.P. (2010) Control theory and systems biology. MIT Press Cambridge.
関連項目