[論文レビュー] Dynamic and static limitation in multiscale reaction networks, revisited

Alexander N. Gorban, Ovidiu Radulescu|ArXiv.org|Mar 30, 2007

Zeolite Catalysis and Synthesis被引用数 12

ひとこと要約

本稿では、離散的力学系と反応グラフ上の「サイクル修復」を用いて、線形マルチスケール反応ネットワークにおける静的および動的制限の一般理論を構築する。well-separatedな速度定数を有するネットワークに対して、固有値が実数であり、定常状態（静的）および緩和（動的）挙動の両方において、単一の制限段階によって支配される優位な遅いモードを示す。これにより、頑健で計算的に効率的なモデル簡略化が可能となる。

ABSTRACT

The concept of limiting step gives the limit simplification: the whole network behaves as a single step. However, in its simplest form this idea is applicable only to the simplest linear cycles in steady states. For such the simplest cycles the nonstationary behaviour is also limited by a single step, but not the same step that limits the stationary rate. We develop a general theory of static and dynamic limitation for all linear multiscale networks. New estimates of eigenvectors for diagonally dominant matrices are used. Multiscale ensembles of reaction networks with well separated constants are introduced and typical properties of such systems are studied. For any given ordering of reaction rate constants the explicit approximation of steady state, relaxation spectrum and related eigenvectors ("modes") is presented. The obtained multiscale approximations that we call "dominant systems" are computationally cheap and robust. These dominant systems can be used for direct computation of steady states and relaxation dynamics, especially when kinetic information is incomplete, for design of experiments and mining of experimental data, and could serve as a robust first approximation in perturbation theory or for preconditioning.

研究の動機と目的

• 単純な不可逆的サイクルを超えて、一般の線形マルチスケールネットワークへと古典的な制限段階の概念を拡張すること。
• 静的（定常状態）および動的（緩和）制限段階を区別し、複雑なネットワークでは両者が異なることがあることを示すこと。
• 不完全な動力学的データを有するマルチスケール反応ネットワークに対して、計算的に効率的で頑健な近似手法（「優位系」）を開発すること。
• well-separatedな速度定数を有する場合、固有値が実数であり、緩和スペクトルが単一の遅いモードによって支配されることを証明すること。
• スペクトル的およびグラフ論的ツールを用いて、制限段階の特定と還元モデルの構築のための明示的アルゴリズムを提供すること。

提案手法

• 有限集合上の補助的離散力学系を用いてネットワーク構造を分析し、制限段階を同定する。
• 反応グラフを変更し、優位な動力学的挙動を分離するための「サイクル修復」アルゴリズムの開発。
• 対角優勢行列の固有ベクトルに対する新しい推定式を適用し、スペクトル特性を評価する。
• 小速度定数を摂動として扱う摂動論を用いて、動力学的行列への近似を構築し、優位系を構成する。
• 固定点方程式と正規化制約を用いて、定常状態近似および緩和モードを明示的に導出する。
• 左固有ベクトルを用いて遅い変数を定義し、全動力学方程式を単一の有効な段階に還元可能にする。

実験結果

リサーチクエスチョン

• RQ1 制限段階の概念は、単純な不可逆的サイクルを超えて、任意の線形マルチスケール反応ネットワークへ一般化可能か？
• RQ2 複雑なネットワークにおいて、同じ反応段階が定常状態の速度と緩和時間の両方を制限するか？
• RQ3 マルチスケールネットワークの固有値が実数である条件は何か？これにより減衰振動が回避される。
• RQ4 どのようにして体系的に優位系を特定し、構築できるか？ これは、完全なネットワーク挙動を正確に近似する。
• RQ5 well-separatedな速度定数を有するネットワークにおける遅い多様体の数学的構造は何か？

主な発見

• well-separatedな速度定数を有するマルチスケールネットワークでは、固有値が実数であり、減衰振動が生じにくいことが示された。
• 静的（定常状態）および動的（緩和）制限段階は一般に異なる。緩和時間は別個の速度決定段階によって支配される。
• 優位系近似は明示的に構築されており、不完全な動力学的データが存在しても計算的に頑健である。
• 遅いモードは単一の有効速度定数 λ₀ ≈ −εw*₁ によって支配され、ここで ε は小さな摂動パラメータ、w*₁ は元の不変サイクルにおける定常フラックスである。
• 左固有ベクトル l⁰ は遅い変数を提供し、有効入力項を伴う還元された動力学方程式を満たす。これによりモデル簡略化が可能となる。
• 明示的な誤差境界が導出された：固有ベクトルおよび固有値近似の相対誤差は O(ε) であり、|χᵢ| < 3ε および |ςᵢ| < ε を満たす。