化学反応速度論シミュレーター 戻る
化学工学・反応工学

化学反応速度論シミュレーター

分子が衝突して反応する様子と、濃度-時間プロファイルをリアルタイムにアニメーション。反応速度定数・半減期・平衡定数・転化率を即時表示。

反応系の選択
d[A]/dt = −k[A]
速度定数
k (正方向)
s⁻¹
初期濃度 (mol/L)
[A]₀
M
[B]₀
M
シミュレーション時間 (s)
s
再生コントロール
プリセット
温度依存(アレニウス)

アレニウス式

温度 T (K)
K
活性化エネルギー Ea (kJ/mol)
kJ
k(T) = —
計算結果(ライブ)
0.00
経過時間 t (s)
反応速度 −d[A]/dt
半減期 t½ (s)
平衡定数 K
転化率 X (%)
収率 Y (%)
反応容器(分子アニメーション)
濃度-時間プロファイル(進行点つき)
理論・主要公式

速度則: $-\dfrac{d[A]}{dt}=k[A]^n$ / 一次半減期 $t_{1/2}=\dfrac{\ln 2}{k}$

可逆平衡: $[A]_{eq}=\dfrac{k_r}{k_f+k_r}[A]_0,\quad K=\dfrac{k_f}{k_r}$

アレニウス: $k(T)=k_0\exp\!\left[-\dfrac{E_a}{R}\left(\dfrac1T-\dfrac1{T_0}\right)\right]$

RK4: $k_1=f(t,y)$, $k_2=f(t+h/2,\,y+hk_1/2)$

$k_3=f(t+h/2,\,y+hk_2/2)$, $k_4=f(t+h,\,y+hk_3)$

$y_{n+1}=y_n + \dfrac{h}{6}(k_1+2k_2+2k_3+k_4)$

化学反応速度論シミュレーターとは

🙋
反応速度定数「k」って何ですか?シミュレーターで変えると何が起こるんですか?
🎓
大まかに言うと、反応の「速さ」を決める数字だよ。上の容器アニメーションを見て。「k」を大きくすると、青い分子(反応物A)が赤い分子(生成物B)にどんどん置き換わって、下の濃度曲線でも[A]が急降下する。実務では、触媒を加えるとこのkが大きくなって、効率よく製品を作れるんだ。
🙋
え、そうなんですか!じゃあ「半減期」って表示は、kと関係あるんですか?
🎓
その通り!一次反応A→Bなら、半減期は $t_{1/2}= \ln(2)/k$ で計算できる。曲線グラフ上の点線「t½」が、[A]が初期値の半分になる時刻だよ。kが2倍になれば半減期は半分になる。スライダーを動かして点線が左右に動くのを確かめてみて。
🙋
「可逆反応」のタブで、正方向と逆方向のkを両方設定できますね。これはどういう場面で使うんですか?
🎓
例えば、エステル化反応(酸+アルコール⇄エステル+水)みたいに、行ったり来たりする反応だね。正方向のkを大きく、逆方向のkを小さく設定すると、曲線が途中で水平になって平衡濃度 $[A]_{eq}=k_r/(k_f+k_r)\,[A]_0$ に落ち着く。容器の中でも青と赤の比率がある一定値で止まるのが見えるよ。現場では、このバランスをシミュレーションして最適な反応条件を探すんだ。

よくある質問

各反応物の初期濃度(mol/L)と反応速度定数k(一次反応なら1/s、二次反応ならL/(mol·s))を数値入力してください。可逆反応の場合は平衡定数も設定可能です。数値は半減期の自動計算にも反映されます。
一次反応(A→B)、二次反応(A+B→C)、可逆反応(A⇌B)、連続反応(A→B→C)、自触媒反応(A+B→2B)に対応しています。各反応の濃度-時間プロファイルをリアルタイムでアニメーション表示し、速度定数や半減期も自動計算します。
はい、シミュレーション時間(例:0~100秒)を設定できます。また、容器内の分子は反応物・生成物ごとに色分けされ、下の濃度曲線には現在時刻を示す進行点と半減期マーカーが重ねて描かれます。
本シミュレーターは4次のルンゲ=クッタ法で数値的に解いており、一般的な条件で高精度です。ただし、反応速度定数が極端に大きい場合や時間刻みが粗いと誤差が生じるため、適切な時間刻み(例:0.01秒)を推奨します。

実世界での応用

医薬品開発:薬の体内での代謝(分解)速度を一次反応モデルで予測し、適切な投与間隔や持続時間を設計します。シミュレーターで半減期を推定することで、臨床試験前のスクリーニングが効率化されます。

化学プラント設計:連続反応(A→B→C)をシミュレーションし、目的製品Bの収率が最大になる反応器の滞留時間や温度を決定します。これにより、原料の無駄を減らしコストを削減できます。

環境工学:排水処理における有害物質の分解反応を二次反応モデルで解析します。処理剤の最適な投入量と反応時間を求めることで、環境基準を満たす効率的な処理プロセスを構築します。

食品化学:食品の保存中に起こる栄養素の分解や風味の変化を可逆反応モデルで追跡します。温度や包装条件(平衡の変化)が品質に与える影響を予測し、賞味期限を科学的に設定します。

よくある誤解と注意点

まず、「反応速度定数kは物質固有の定数ではない」という点を押さえよう。kは温度や圧力、触媒の有無で大きく変わるんだ。例えば、同じエステル化反応でも、温度を50℃から80℃に上げると、アレニウスの式に従ってkは数倍になる。シミュレーターで「k」をいじるのは、実は「温度や触媒を変えたらどうなるか」を仮想的に試しているのと同じことなんだよ。

次に、「初期濃度を変えると、反応の『見かけの速さ』が変わる」という現象に注意。二次反応A+B→Cで、Aの初期濃度だけを10倍にすると、反応開始直後の速度は確かに速くなる。でも、速度定数k自体は変わっていない。グラフの見た目だけで「反応が速くなった!」と判断せず、速度則そのものを確認するクセをつけよう。

最後に、実務でシミュレーション結果を使う時の落とし穴。「理想的なバッチ反応器」を想定した計算であることを忘れないで。実際のプラントでは、撹拌ムラや熱移動の制限、副反応などで計算通りにならないことが多い。例えば、シミュレーターで最適と出た滞留時間10分に対して、現場では安全マージンをみて15分に設定する、といった調整が必要なんだ。

使い方ガイド

  1. 反応タイプ(一次・二次・可逆・連続・自触媒)を選択し、速度定数k(単位:s⁻¹またはL/(mol·s))をスライダーまたは数値入力で設定
  2. 初期濃度A₀(mol/L)を指定し、可逆反応の場合は逆反応速度定数kᵣを入力。容器アニメーションと濃度曲線が自動で再生される
  3. 「一時停止」「最初から」で再生を制御。半減期t½、平衡定数K、転化率X、収率Y、ライブの反応速度がリアルタイムに更新される

具体的な計算例

酢酸エステル加水分解(一次反応):初期濃度A₀=0.1mol/L、速度定数k=0.05s⁻¹の場合、半減期t½は約14.0秒です。可逆反応(A⇌B)でk=0.08s⁻¹、kᵣ=0.02s⁻¹、A₀=1.0mol/L、B₀=0のとき、平衡定数K=4.0、時間t=100秒での転化率Xは約80.0%です。二次反応は本ツールではA+B→Cとして扱い、k=0.02L/(mol·s)、A₀=B₀=2.0mol/Lならt½は約25.2秒、t=100秒でC=1.6mol/L(収率80.0%)となります。

実務での注意点