システム準備¶
GROMACS を使用してシミュレーションシステムを準備する方法はいくつかあります。これらの方法は、検討している科学的な質問や、関わる物理モデルによって異なります。例えば、タンパク質-リガンドの原子レベルの自由エネルギーシミュレーションには、マルチステートのトポロジーが必要となる場合があります。一方、粗粒モデルのシミュレーションでは、より密なシステムに適したデフォルトの設定を管理する必要がある場合があります。
考慮すべきステップ¶
以下の一般的なガイダンスは、成功したシミュレーションの計画に役立ちます。一部のシミュレーションの種類では、一部の段階はオプションです。
シミュレーションを通じて、研究対象となるプロパティまたは現象を明確に特定してください。 これが明確になるまで、シミュレーションを実行せずに進めないでください。 シミュレーションを実行した後で、その結果を使って仮説を検証する方法を探すのではなく、シミュレーションが適切でない場合や、必要な情報が保存されていない場合に備えてください。
適切なツールを選択することで、シミュレーションを実行し、関心のある特性や現象を観察することができます。類似のシステムに関する他の研究者の出版物を読み、内容を理解することが重要です。 選択できるツールには以下のようなものがあります。
シミュレーションを実行するためのソフトウェア（力場を考慮する必要がある）
フォースフィールドは、システム内の粒子が互いにどのように相互作用するかを記述するものです。 調査対象のシステムと、関心のある特性または現象に適したものを選択してください。 これは非常に重要で、容易ではありません！
シミュレーションデータを分析し、観察を行うための手順。
システム内に配置する各分子に対して、初期の座標ファイルを入手または生成します。多くの異なるソフトウェアパッケージは、分子構造を構築し、それらを適切な構成に組み立てることができます。
システム用の基本的な構造を生成するには、分子を適切な座標ファイルに配置します。分子は、特定の場所に配置することも、ランダムに配置することもできます。いくつかの非|Gromacs|ツールが役立ちます。|Gromacs|では、gmx solvate、gmx insert-molecules、および:ref:`gmx genconf`が、よくある問題を解決します。
システム用のトポロジーファイルを取得または生成します。例えば、次のツールを使用できます：gmx pdb2gmx、gmx x2top、
SwissParam、`Automated Topology Builder <https://atb.uq.edu.au/>`（GROMOS96 53A6用）またはお好みのテキストエディタと、|Gromacs|の`Reference Manual`の第5章:ref:`Chapter 5 <ff>`を組み合わせて使用します。AMBER力場を使用する場合は、`antechamber <https://ambermd.org/antechamber/antechamber.html>`または`acpype <https://github.com/alanwilter/acpype>`が適切です。最終的に希望する密度に適したサイズを持つシミュレーションボックス（例：:ref:`gmx editconf`を使用）を作成し、溶媒（例：:ref:`gmx solvate`を使用）で満たします。また、システムを中和するために必要なカウンターイオン（例：:ref:`gmx grompp`と:ref:`gmx insert-molecules`を使用）を追加します。これらの手順では、座標ファイルと一致するようにトポロジーファイルを編集する必要がある場合があります。
システム上でエネルギー最小化を実行します（:ref:`gmx grompp`と:ref:`gmx mdrun`を使用）。これは、生成中に発生した不適切な初期構造を解消するために必要です。これにより、シミュレーションがクラッシュするのを防ぐことができます。また、溶媒分子（またはリピッドバイレイヤーなど）を導入する前に、溶媒中でソリューション構造を最小化する必要がある場合があります。柔軟な水モデルを使用し、結合制約や凍結グループを使用しないことを検討してください。位置制約と/または距離制約の使用は、慎重に評価する必要があります。
適切なシミュレーションパラメータを選択してください（
:ref:`mdp`ファイルで定義）。シミュレーションパラメータは、フォースフィールドがどのように導出されたかに基づいて選択する必要があります。NVTで位置制約を適用して溶媒と/または分子の温度をほぼ調整した後、密度を固定するためにNPTに移行する必要がある場合があります（c-rescaleバロスタットの使用を推奨）。必要に応じて、さらにシミュレーションを進めて、目的の生産シミュレーションのエンsembles（例：NVT、NVE）に到達します。システムが:ref:`blowing-up`の問題を引き起こしている場合は、そのページにある提案を検討してください。例えば、分子に対する位置制約の使用、結合制約の使用を避ける、またはより小さな積分ステップサイズを使用する、またはより穏やかな加熱ステージを複数使用することを検討してください。十分な時間をかけて平衡シミュレーションを実行し、システムがターゲットのエンsemblesに十分に緩和されるようにし、生産シミュレーションを開始できるようにします（gmx grompp および gmx mdrun、その後 gmx エネルギー および 可視化ソフトウェア を使用）。
適切なシミュレーションパラメータを、生産シミュレーション（:ref:`mdp`ファイルで定義）に合わせて選択してください。特に、速度を再生成しないように注意してください。フォースフィールドがどのように導出され、および、関心のあるプロパティまたは現象を測定する方法について一貫性を保つ必要があります。
ヒントとコツ¶
データベースファイル¶
GROMACS のインストールに含まれる share/top ディレクトリには、.dat という拡張子を持つ多数のプレーンテキストのヘルパーファイルが含まれています。これらのファイルは、一部のコマンドラインツール（コマンドラインリファレンス を参照）で使用されており、各ツールは使用されているファイルと、それらの使用方法をドキュメントに記載しています。
もし、これらのファイルを修正する必要がある場合（たとえば、VDW半径を含む新しい原子タイプを vdwradii.dat に導入する場合）、ファイルはインストールディレクトリからワークディレクトリにコピーできます。その場合、GROMACS のツールは、ワークディレクトリにあるコピーを読み込みます。標準の定義をすべて無効にするには、ワークディレクトリに空のファイルを作成します。