事例2：低分子混合物（水 ＋ エタノール）

事例2：低分子混合物（水 ＋ エタノール）

この事例では、2種類の低分子（水とエタノール）を特定の個数比で混合したシステムを構築します。

対象システム

• 構成:
  • 水：300分子
  • エタノール：150分子
• 総原子数: 2,250
• 初期設定密度: 0.80 g/cm³

計算コード (run.py)

import fbtk

# 1. 各分子オブジェクトをSMILESから作成
water = fbtk.Molecule.from_smiles("O", name="water")
ethanol = fbtk.Molecule.from_smiles("CCO", name="ethanol")

# 2. Builderの初期化
builder = fbtk.Builder(density=0.80)

# 3. 各分子を指定した個数ずつ追加
builder.add_molecule(water, count=300)
builder.add_molecule(ethanol, count=150)

# 4. システムの構築と緩和
system = builder.build()
print("Starting relaxation for Water + EtOH...")
system.relax()

# 5. 結果の出力
water.to_file("water.mol")
ethanol.to_file("ethanol.mol")
system.to_file("water_etoh_cell.mol2")
print("Done.")

実行ログの確認

Atoms:     2250       | Bonds:     1800
--------------------------------------------------------------------------------
       | Fmax           | FRMS           | Total E          |
Iter   | (kcal/mol/Å)   | (kcal/mol/Å)   | (kcal/mol)       | Status
--------------------------------------------------------------------------------
     0 |     36477.1705 |      3459.2950 | 204596098964.3458 |
    10 |      1620.0353 |       134.2404 |       51220.1694 |
    20 |       813.1749 |        50.6762 |       24796.0682 |
   ...
   280 |         2.9743 |         0.2826 |        2533.5805 |
   290 |         1.5966 |         0.2099 |        2532.7015 | FRMS-Conv
--------------------------------------------------------------------------------
=== Optimization Finished ===
Reason: FRMS-Conv
Total Time: 5.678s     (Avg: 19.512ms / step)
Final Energy: 2532.7015       kcal/mol
Final Fmax:   1.5966          kcal/mol/Å
Final FRMS:   0.2099          kcal/mol/Å
Min Distance: 0.9979          Å (Atoms 1353 and 1359)
--------------------------------------------------------------------------------

結果の解説

混合状態の初期構造

緩和後の Min Distance は 0.9979 Å となっています。水素結合ネットワークが形成される前の「爆発しない」安定な初期配置が得られていることが確認できます。