12. よくある質問・トラブルシューティング
MOPAC, CNDO/S, GAMESS, NWChem, Gaussianに関して
Q. 系のモデリング中や、計算結果を読み込んだ際に、不要な結合が表示されます。または結合が表示されません。どのように対処したらいいですか？
Q. MOPACのログに ATOMS **, **, AND ** ARE WITHIN .**** ANGSTROMS OF A STRAIGHT LINE と表示されて異常終了します。
Q. MOPACを実行するとログに「UNRECOGNIZED KEY-WORDS: ( PM6（ハミルトニアン名）)」と出力され計算が終了してしまいます。
Q. GAMESS, Gaussian, NWChemのESP（静電ポテンシャル）の3次元表示を高速化することはできますか？
Q. GAMESSの実行時に「* ERROR: MEMORY REQUEST EXCEEDS AVAILABLE MEMORY」などとログに出力され計算が異常終了します。
Q. Windows版GAMESSをインストールした直後は問題なく利用できていたが、ある時期から計算が全て途中で異常終了するようになりました。対処方法はありますか？
Q. GAMESSのログに「 **** ERROR **** PCM SPHERE(S) MUST HAVE A POSITIVE RADIUS 」と表示され異常終了します。
Q. GAMESSの実行時に「ERROR: BAD DELOCALIZED COORDINATES GENERATED!!!」とログに出力され計算が異常終了します。
Q. NWChemの並列実行時に「Please specify an authentication passphrase for smpd: 」とログに出力され計算が流れません。
-
Q. MD計算においてSHAKE法などによる拘束は計算結果にどのような影響を与えますか？拘束方法はどのように選んだらよいですか？
Q. MD計算を実行後、アニメーションを観たり、最終構造を見ていると、分子がセルの外側に出てしまうことがあるのはなぜですか？
Q. [MD]-[自動で電荷を割り当て]メニューで「Topology file not found」というエラーが表示されます。解決方法はありますか？
Q. 界面ビルダで作成した構造からMD計算を実行しようとすると、力場の割り当てが失敗したりMD計算破綻します。解決方法はありますか？
Q. LAMMPS, Gromacsで液体（有機物）-固体（無機物）間の相互作用パラメータをどのように決めたらいいですか？
Q. Gromacsの最終構造やアニメーションを読み込んで取得した構造や、それに何かしらの編集を行ってから再びMD計算を実行すると力場の割り当てが失敗したりMD計算が破綻します。解決方法はありますか？
-
Q. Quantum ESPRESSOのSCF計算が出力ファイル（.pwoutまたは.out）に「too few bands」と表示され異常終了します。nbndの設定方法が分かりません。
Q. Quantum ESPRESSOのSCF計算が出力ファイル（.pwoutまたは.out）に「charge is wrong」と表示され異常終了します。どのように解決したらいいですか？
Q. Quantum ESPRESSOを用いて誘電関数を計算する際に、epsilon.xの出力(eps.out)に「bad band number」と表示され誘電関数を取得できません。
Q. Quantum ESPRESSOでセルサイズの構造最適化を行うと最終構造のエネルギーで不連続的な変化をしたりエネルギーが増加するのはなぜですか？
Q. 系のモデリング中や、計算結果を読み込んだ際に、不要な結合が表示されます。または結合が表示されません。どのように対処したらいいですか？
12.1. 購入に関して
12.1.1. Q. 代金の支払方法を教えてください。
A.
【法人の場合】
以下の条件での後払いとなります。
支払方法: 当社指定銀行口座への現金振込
支払期日: 納品翌月末日
【個人の場合】
PayPalにてクレジットカードでお支払いください。
12.1.2. Q. 当社から発行される書類の種類を教えてください。
A.
請求書・納品書・見積書を発行いたします。
ただしPayPalの場合のみ、PayPalから領収書を取得してください。
その他の書類を発行希望の際はご相談ください。ただし、内容によりお断りする場合もありますのでご了承ください。
12.1.3. Q. 代理店などエンドユーザ以外が注文する方法を教えてください。
A.
指定代理店をご利用ください。詳細は 価格・購入 をご確認ください。指定代理店を介したくない場合は、弊社（株式会社クロスアビリティ）からの直販のみ可能です。指定代理店を設けている理由は、Winmostarの価格とサービスを、エンドユーザに適切な形で提供するためです。
その上でなおエンドユーザが指定代理店以外の代理店との取引を希望する場合は、当該代理店から指定代理店経由での販売は可能です。ただし、その際当該代理店と弊社の間でサービス内容の保証に関する取り決めがないため、サービス内容の保証は致しかねます。
12.1.4. Q. Winmostar本体と保守の価格を分割することは可能ですか？保守のみの販売はありますか？
A.
分割することはできません。保守のみの販売はありません。教育機関向けのラインナップで保守が切れ、保守を更新したい場合は、Winmostar本体も含めて再度購入頂く必要がございます。
12.1.5. Q. 製品の保証書はありますか？
A.
保証書という形での書類はありませんが、ご購入時にご同意いただく 使用規約 と サポートサービス規約 に保守等の保証書の内容に相当する文章が書かれておりますので、使用規約とサポートサービス規約を保証書の代用としてご利用ください。
12.1.6. Q. 個人利用での割引などありますか？個人向けライセンスはありますか？
A.
特にありません。無償版については、個人でもご利用頂けます。プロフェッショナル版については、教育機関に所属し教育機関での用務に使用する目的であれば教育機関向け、それ以外であれば民間企業・官公庁向けのライセンスのご購入ください。
12.2. ライセンス・ライセンスコードに関して
12.2.1. Q. 特定ユーザライセンスの登録ユーザ変更は可能ですか？
民間企業・官公庁の場合は不可能です。教育機関の場合は、前回のユーザ変更（初回の変更の場合は購入）から1年以上経過していたら変更可能です。
12.2.2. Q. 教育機関向けライセンスを購入した時の所属機関から他の機関に移った場合、そのライセンスを使い続けることは可能ですか？
購入時の機関の所属から外れた場合、永久使用権であっても使用できません。
12.2.3. Q. 無償版、学生版、無料トライアルからプロフェッショナル版といった具合に、途中から製品を購入するなどしてエディションを変更する際に再インストールは必要ですか？
A. 不要です。これから使用したいライセンスコードを メニューの ライセンスコード に入力してください。
12.2.4. Q. 無償版とプロフェッショナル版、といった具合に異なるエディションのWinmostarを複数インストールすることは可能ですか？
A. 可能です。その場合は、Winmostarのインストーラでエディションごとに別々のインストール先を指定してください。
12.2.5. Q. すでに入力したライセンスコードを変更する方法を教えてください。
A. これから使用したいライセンスコードを メニューの ライセンスコード に入力してください。
学生の方で、無償版から学生版に切り替えたい場合は、ライセンス登録ページ で学生版にチェックを入れて再度ライセンス登録をしてください。
12.2.6. Q. MACアドレスを表示する方法を教えてください。
A. Winmostar V11の場合は、ライセンス未登録の状態で起動するとそのマシンのMACアドレスが表示されます。
Winmostar V10以下でWindows10の場合は、まず をクリックしてコマンドプロンプトを起動します。次に、コマンドプロンプトのウィンドウで ipconfig /all と入力しEnterキーを押します。様々な情報が出力されるので、その中から「物理アドレス」の行を探してください。その内容がMACアドレスです。
「物理アドレス」行が複数ある場合、Winmostarのノードロックライセンス購入時に申請するMACアドレスは、基本的にどの「物理アドレス」でも大丈夫です。
12.2.7. Q. Winmostarのライセンス料に、Gaussianのライセンス料は含まれていますか？Gaussianのライセンスが付属したWinmostarのラインナップはありますか？
A. Winmostarのライセンス料に、Gaussianのライセンス料は含まれていません。Gaussianのライセンスが付属したWinmostarのラインナップはありません。Gaussianの代理店から別途ご購入ください。
12.2.8. Q. ライセンス期限が来るとそれまで使用していたWinmostarはどうなりますか？
A. ライセンス期限が来た後Winmostarを起動すると、ライセンスコードを入力するウィンドウが立ち上がります。再びWinmostarのライセンスを更新し入力すると、ライセンス期限が来る前の状態を引き継いで使用することができます。
12.2.9. Q. 個人での学習目的に教育機関向けライセンスを購入できますか？
A.
現時点で教育機関に所属し、教育機関での用務に使用しない限り、教育機関向けライセンスをご購入頂けません。個人での利用については Q. 個人利用での割引などありますか？個人向けライセンスはありますか？ をご確認ください。
12.3. サポート・保守に関して
12.3.1. Q. 不具合サポートやパッチ提供は受けられますか?
使用規約の内容に基づき実施されます。最新の使用規約はこちら（ Winmostar 使用規約 ）です。
12.3.2. Q. 旧バージョンのサポート・保守はどこまで行われますか?
ご使用中のバージョンのWinmostarに関する、有効な使用規約に記載の内容に基づきます。また、操作方法の簡単な案内は、可能な範囲で対応します。
12.3.3. Q. Winmostarの開発元（製造元）はどこですか？製造地はどこですか？
株式会社クロスアビリティ です。製造地は日本です。
12.3.4. Q. 開発元（製造元）と直接連絡取ることはできますか?
問い合わせフォーム から連絡を取ることができますが、対応の可否は利用規約に基づきます。最新の使用規約はこちら（ Winmostar 使用規約 ）です。 有償サポート を利用することで、より進んだメールでのサポートが可能となります。
12.3.5. Q. 使用しているWinmostarのアップデート・バージョンアップ・アップグレードは可能ですか?
A. マイナーバージョン（およびリビジョン）の更新については、利用可能期間内であれば何回でも実施可能です。メジャーバージョンの更新については、永久使用権の場合はライセンスの更新が必要で、年間使用権の場合は実施可能です。
例として、「V8.039」については、「8」がメジャーバージョン、「039」がマイナーバージョンを指します。「V9.1.0」については「9」がメジャーバージョン、「1」がマイナーバージョン、「0」がリビジョンを指します。
例えば、Winmostar V9の永久使用権のライセンス取得者は、V9.1.0からV9.1.5やV9.4.4に更新することは可能ですが、V10.0.0に更新することは不可能です。
12.3.6. Q. 質問するときの注意事項はありますか？
A. 計算が上手く流れない等の質問の場合、原則として状況を再現するインプットやアウトプットファイルをお送り下さい。
12.3.7. Q. 保守のみの販売はありますか？
A.
12.4. ソフトウェアの機能・利用について
12.4.1. Q. 「Winmostar」の呼び方は何ですか？
A. 「ウインモスター」です。Wikipedia等では誤情報が掲載されることがありますが、こちらが正式な呼び方です。
12.4.2. Q. Winmostarは社内外のサーバやクラウドに接続しますか?LAN接続していない状態でも利用可能ですか?
12.4.3. Q. Winmostarの動かすPCの最小・推奨スペックを教えてください。
A. 最小・推奨スペック をご確認下さい。
12.4.4. Q. Winmostarで作成したデータを学会発表や論文に用いることは可能ですか？学会発表、論文投稿の際にどのように引用したらいいですか？
A. 使用いただいて問題ありません。発表される際には 引用について の通りに引用してください。
12.4.5. Q. Winmostarの画面を撮影した動画・画像をYouTube、SNSなどのウェブ媒体にアップロードすることは可能ですか？
A. 可能です。アップロードする際には、WinmostarのHPのURLを引用し、Winmostarを使用していること、使用したWinmostarのバージョンを明記してください。なお、ライセンスキーが表示される環境設定ウィンドウなどのアップロードは固く禁じます。
12.4.6. Q. Winmostarそのものに各種計算用ソフトがインストールされていますか？
A. MOPAC、CNDO/SのみWinmostarにインストールされています。それ以外のソフトは、ライセンスの関係上Winmostarには同梱されておらず、別途インストールする必要があります。多くのソフトは無料でインストール可能で、その手順は インストール で紹介されています。
12.4.7. Q. Winmostarはクラウド上で計算させていますか？
A. クラウド上で計算させることも可能ですが、させないことも可能です。デフォルトではクラウドを利用せず、WinmostarをインストールしたWindows PC上で計算をさせます。
12.4.8. Q. WinmostarはGPU計算に対応していますか？
A. GPU計算に対応していますが、デフォルトではGPUを使わない設定になっています。GAMESS, Gaussian, Gromacs, LAMMPS, Quantum ESPRESSOがGPUに対応していますが、動作確認および設定作業は有償での対応となります。
12.4.9. Q. IntelとAMDのどちらのCPUの方が動作に適していますか？どちらがお勧めですか？
A. 一般に、シミュレーションにおいてどちらが優れているということはありません。
12.4.10. Q. Winmostarを使って並列計算は可能ですか？
A. 可能です。詳細は、各ソルバのキーワード設定ウィンドウのページをご確認ください。
12.4.11. Q. Winmostarを使って並列計算する際、利用できるコア数の上限はありますか？利用できるコア数に応じて費用は変わりますか？
A. ユーザが用意したハードウエアの範囲内で、制限なく並列数を指定して頂けます。並列数に応じて、Winmostarのライセンス料は費用は変化しません。ローカルジョブの場合は、 Winmostar Job Manager で設定した最大コア数を上回るとジョブが流れないため、最大コア数の設定を変更してください。
12.4.12. Q. WinmostarはmacOS、Linuxで動作しますか？
A. Winmostarのアプリケーション本体はWindows OSのみサポートされています。サポートされているWindows OSの確認は 動作環境 で可能です。macOS、LinuxでWinmostarのアプリケーション本体を動かす場合は、VirtualBoxなどの仮想環境上にWindows OSをインストールした上でご使用ください。
リモートジョブを実行するコンピュータには、Linux・macOSを使用できます。
12.4.13. Q. Gaussianのインストール方法を教えてください。
Gaussianのインストール方法は、Gaussianの販売代理店より入手してください。Gaussianをインストールした後は、 において、Gaussianのプログラムパス（g03.exe, g09.exe, g16.exeなど）を選択してください。
12.4.14. Q. Winmostarで何原子あるいは何分子まで計算できますか？扱える原子数、分子数の上限はありますか？
動作速度は考えないとすると、100万原子程度までの動作確認はしています。動作速度は実行環境に強く依存するため、ご購入前に無料トライアルでご確認ください。なお、将来のバージョンではWinmostarの高速化を計画しております。
12.4.15. Q. Winmostarで粗視化モデルの計算はできますか？
純粋な量子化学計算、第一原理計算、古典分子動力学計算ではないという意味での粗視化モデルとしては、LAMMPSを用いた散逸粒子動力学（DPD）計算とKremer-Gresetモデルに対応しています（Kremer-Grestについては別途お問い合わせください）。United atomモデルやその派生の粗視化モデルは今後対応予定です。それ以外のモデルについては、個別にお問い合わせください。個別対応としている理由は、粗視化モデルを用いたほとんどの研究が、ソフトを使うだけでは有意義な結果が得られず、入念かつ慎重なコンサルティングが必要であるという認識を我々が持っているためです。Winmostarのサポート事例では、十分に検証を行った粗視化モデルのシミュレーション結果を紹介しています。
12.5. ソフトウェア動作全般に関して
12.5.1. Q. 思ったようにモデルを作成できません。計算できません。動作しません。
まず、以下の基礎的なチェックを行ってください。
インストール時の注意事項 を確認する。
使用中のWinmostarが無償版、学生版、プロフェッショナル版、プロフェッショナル版（トライアル）のいずれに該当するか確認し、問題を起こしている機能がその版で使用可能か 機能表 を見て確認する。
使用中のセキュリティ対策ソフトの活動記録を確認し、WinmostarおよびCygwinWMのインストールフォルダの下のアプリケーションの活動が妨害された記録がないか確認する。
Winmostarを最新版にアップデートし（使用中のバージョンと共存させることが可能）、 既知の不具合 、 よくある質問・トラブルシューティング に類似する状況がないか確認する。
保存するファイルやそれを含むディレクトリ（上位階層全てを含む）の名前に、日本語、全角文字などのマルチバイト文字や特殊記号（スペースも不具合の原因となります）が含まれている場合は、一部ソルバで不具合が出ることがあるため半角英数のみとなるようにする。
実行した処理で何かしらログが出力されているか作業フォルダを確認し、ログの内容を確認する。
計算が開始されたが計算結果がおかしいと感じた場合は、メインメニューで使用したソルバのメニューから「ログを表示」などをクリックし、ログの内容を確認する。
計算の不具合については、各種ソルバのバージョンが、Winmostarのインストールガイドで推奨しているバージョンと同じであるか確認する。（特にGromacs, LAMMPS, Quantum ESPRESSO）
次に、メモ帳などで以降の作業の記録を取れるようにしてください。不具合の再現方法が判明した場合、作業の記録と一緒にご報告頂くと比較的短時間で修正できることがあります。
そして、Winmostarの チュートリアル のうち、これから使いたいソルバの基礎編チュートリアルをトレースしてください。
基礎編チュートリアルのトレースに失敗する場合は、以下を試してください。
誤操作でないことを確認するため再度トレースする。
並列実行している場合は、シリアル実行（並列数1）に切り替える。
Winmostarを再起動する。
OSを再起動する。
セキュリティ対策ソフトで、Winmostar、CygwinWMのインストールフォルダ、およびソルバ（MPIを含む）が監視対象外に設定する。
CygwinWMを使用している場合は、 でCygwinWMの簡易的な診断を実行する。
Winmostar, CygwinWMおよび使用したソルバを再インストールする。
他のPCで試す。
次に、最終的に計算したいものに極力近いと思われるチュートリアルをトレースしてください。 それに成功したら、最終的に計算したいものに少しずつ寄せるように計算条件を変更し（原子数、スーパーセルのサイズ、重合度、元素の種類、相の数など）、問題発生箇所を特定したら以下を試してください。
よくある質問・トラブルシューティング に類似事例がないかご確認ください。
問題発生箇所がWinmostarが外部ソフトを呼んでいる部分の場合は、そのソフトの情報もご確認ください。
Cygwinを用いた処理で落ちている場合は、 Cygwinの一般的な不具合 をご確認ください。
12.5.2. Q. 「ERROR: I/O error 32」と表示され処理に失敗します。
A. 処理に関わるファイルがWinmostar以外のアプリケーションまたはプロセスで開かれていてロックされている場合や、削除されている可能性があります。
OSを再起動し他のアプリケーションが開いていない状況でお試しください。
12.5.3. Q. Cygwinを使う処理が異常終了します。／ 機能で ... ERROR ... と表示されます。／Cygwinの黒いウィンドウに child_info_fork::abort: ... Loaded to different address: parent ... != child ... などと表示されます。
A. 以下の手順を上から順に一つずつ実行し、その都度、エラーが起きた処理を再実施してください。
一般的な 一般的な不具合の対処 を実施する
マシンを再起動する
Windowsセキュリティ開き アプリとブラウザーコントロール から Exploit Protectionの設定 クリックする。そして、 イメージのランダム化を強制する の値を 既定でオフにする か 既定値を使用する（オフ） に変更する。
使用しているCygwinWMの
cygwin1.dll以外を検索して削除し、マシンを再起動する警告
同一マシン上にCygwinWM以外に
cygwin1.dllが存在して場合の一部のケースでこの操作が必要です。
cygwin1.dllは他にCygwinをインストールしていなくても、各種フリーウエアなどに同梱されていることがあります。
使用しているマシン上の全てのCygwinが終了している状態で、Windowsの[ファイル名を指定して実行]にて
C:\cygwin_wm\bin\ash.exe（CygwinWMをC:\cygwin_wmにインストールした場合）を実行し、 /bin/rebaseall -v というコマンドを実行しマシンを再起動する。セキュリティ対策ソフトを一時的に無効する。
CygwinのFAQ に記載されている不具合を起こしがちなソフトを無効にする。
その他、 Cygwinのfork()関連の失敗に関するFAQ に記載された方法を試す。
Cygwin公式サイト のCygwinを新規にインストールし、そこからターミナル（端末）を起動できるか確認する。
12.5.4. Q. 機能で No reference file (... filelist_cygwinwm.txt) was found... と表示されます。
A. CygwinWM診断機能ではまず最初に、CygwinWMをインストールしたフォルダの直下のfilelist_cygwinwm.txtというファイルを探しに行っており、このファイルがないというエラーが出ています。CygwinWMをインストールされた場所（デフォルトでは
C:\cygwin_wm ）でfilelist_cygwinwm.txtを探し、その直上のフォルダを[ツール]-[環境設定]-[プログラムパス]で指定してください。filelist_cygwinwm.txtが見つからない場合はCygwinWMのインストールに失敗した可能性があるので、セキュリティソフトの設定などを見直したうえでCygwinWMを再インストールしてください。12.5.5. Q. 機能で WARNING ... some files are missing と表示されます。
A. CygwinWMを再インストールしてください。
再インストールしても表示される場合は、セキュリティ対策ソフトを一時的に無効にするか、インストール先・インストーラを監視対象外に指定してください。
12.5.6. Q. ジョブマネージャに登録されたジョブが実行されません。
A. 指定したMPIの並列数がジョブマネージャのMaxCoreの設定より大きいとジョブは実行されません。
MaxCoreの初期値値は実行しているPCのコア数に設定されているはずですが、それが変更されていないか、またはMPIの並列数をそれより多く設定していないか確認してください。
ジョブマネージャを使用しないで実行したい場合は、 画面の 計算 タブの「MOPACをジョブマネージャで実行」や「その他のソルバをジョブマネージャで実行」のチェックを外します。
12.5.7. Q. ジョブの実行時に「実行できません（アクセスが拒否されました。）」と表示アラートが出現しジョブが開始されません。
A. 一般的な不具合の対処 を試してください。特に、インストールしたWinmostarおよびソルバのフォルダをセキュリティ対策ソフトの監視対象から外してください。
12.5.8. Q. Winmostarの各種機能やソルバーの実行など、黒いコンソール画面が出現する機能において、黒いコンソール画面の処理が終わらず先に進みません。
A. 黒いコンソール画面の中をたまたまクリックしてしまうと、Windowsの仕様上そこから処理がペンディングしてしまいます。
コンソール画面のウィンドウがアクティブの状態でESCキーを押すと、処理が再開されます。
12.5.9. Q. ファイルを開いたり、分子をモデリングした際に、結合が出現しなくなってしまった。または、無駄に結合が多く出現するようになってしまった。
A. まず、量子化学計算と第一原理計算（固体物理計算）では、結合の情報は計算結果に一切影響を与えないため、結合の有無はあくまで表示上の問題であることを前提にお考え下さい。Winmostarをインストールした直後の状態から結合の生成挙動が変わってしまった場合は、 の 結合判定係数 の値が適切でない可能性があります。デフォルト値に戻すか、1.15程度の値に設定してください。デフォルトの状態でも望みの状態にならない場合は、 または を使い、結合の作成または削除してください。
12.5.10. Q. Windowsのエクスプローラ上で拡張子を表示するにはどうしたらいいですか？
- Windows 7の場合:
エクスプローラを開く
Alt キーを押す
メニューの 表示 タブを開く
登録されている拡張子は表示しない のチェックが外れた状態にする
- Windows 8, 10の場合
エクスプローラを開く
表示 タブを開く
ファイル名拡張子 のチェックが付いた状態にする
12.5.11. Q. マーカー（赤丸）やグループ選択（青丸）が表示されなくなりました。原子を選択できなくなりました。
にチェックを入れてください。
12.6. ファイル入出力に関して
12.6.1. Q. Winmostar以外で生成されたファイルをWinmostarで開けません。また、Winmostarで生成したファイルを編集してからWinmostarで開こうとしても開けません。
A. 改行コードやエンコーディングが変化していないか確認してください。
12.6.2. Q. MOL形式またはSDF形式のファイルを開くと、結合長が不自然になります。水素が出現しません。
A. 次の手順で分子構造を修正してください。(1) 結合長を自動調整 (2) Z-Matrixを再生成 (3) 選択原子に付加（自動） SDFファイルの場合は SDFファイルの編集 の手順を参考に操作してください。
12.7. 分子のモデリング・系の作成に関して
12.7.1. Q. 化学結合の種類（一重、二重など）を変更する方法を教えてください。
A. 例えば以下に示す方法で変更できます。
1) またはメインウィンドウ上部の 結合を付加/変更 ボタンを複数回押すことで、結合の種類を変更できます。
2) を選択すると原子間距離から判定された結合次数で自動的に化学結合の種類が変更されます。予め により構造最適化しておくと、より妥当に自動変更されることがあります。
3) 小さい分子が一つだけしか表示されていない場合は、MOPAC計算を実行することで、Population解析結果を用いて自動的に結合次数が変更されます。
12.7.2. Q. 機能を実行すると Error : Failed to solvate. などと表示され処理に失敗します。
---------質問詳細---------
を実行した際に generate.log に下記のように出力され処理が正常終了しません。
gmx insert-molecules -try 100 -f gmx_tmp_water.gro -o gmx_tmp_water_tmp.gro -ci mol0.gro -nmol 64
...
set +v
Error : Failed to solvate.
一般的な不具合 の対処と、Cygwinの一般的な不具合 の対処に加え、分子数を減らすか、密度を減らして実行してください。
分子数が大きい場合（ケースにもよるが10,000程度）は、現在内部処理で使用している gmx solvate の処理の限界となるケースもあるので、 で分子を並べてください。
将来的には本機能で分子数が大きい場合にも対応予定です。
12.8. ローカルマシンでのMPI・並列実行に関して
12.8.1. Q. MPICHが計算途中で終了します。
---------質問詳細---------
MPICH実行中に、次のようなエラーを表示して計算が途中終了となることがあります。
op_read error on left context: Error = -1
op_read error on parent context: Error = -1
unable to read the cmd header on the left context, Error = -1
unable to read the cmd header on the parent context, Error = -1
Error posting readv, An existing connection was forcibly closed by the remote host.(10054)
connection to my parent broken, aborting.
state machine failed.
A.
このエラーはMPICHがlocalonlyでもネットワークアダプタを使うため、ネットワークアダプタが途中で切れてしまうため発生するエラーです。
しかし初めからネットワークアダプタが切れている場合、MPICHはネットワークアダプタを使用しないため、このエラーは発生しません。
MPICHを用いて長時間の計算を行う場合、ネットワークアダプタを無効にしてから計算を実行して下さい。
12.8.2. Q. LAMMPS, Quantum ESPRESSOのMPI並列実行時に Unable to open the HKEY_LOCAL_MACHINE\SOFTWARE\MPICH\SMPD\process\???? registry key, error 5, アクセスが拒否されました。 という警告が表示されます。
A. MPICHがレジストリを書き換えようとするのですが、管理者権限がないので失敗したというメッセージです。
管理者権限でWinmostarを起動すればメッセージは出なくなりますが、メッセージが出ている状態でも計算自体は正常に実行されているので、無視しても問題ありません。
12.9. リモートジョブに関して
12.9.1. Q. 社内や学内のスパコンまたはLinuxサーバにジョブを投げる方法を教えてください。
A. 接続先のコンピュータ固有の環境設定などが必要な場合も、リモートジョブ用のひな形スクリプトを作成することで可能になります。
詳しくは リモートジョブ をご確認ください。
12.9.2. Q. Test Connection での接続テストは成功するが、ジョブの投入(サブミット)に失敗します。
A. 様々な理由が考えられます。以下にいくつかの例を示します。
1. TSUBAME3.0など、SSH接続の回数制限がある場合は、 TSUBAME3.0でのSSHアクセス数制限について に記載の方法で、SSH接続を都度実行せずにつなぐ方法で回避することができます。
2. サーバ側で、秘密鍵認証だけでなく、パスワード認証もアクティブにすることで回避できる場合もあります。
3. ログインサーバの実体が複数あり、バックグラウンドで自動選択される場合は、特定のログインサーバのみを利用するか、全てのサーバがcache登録されるまで接続しておくことで回避できる場合もあります。
4. ローカルマシンからWinmostarがジョブ投入コマンド（ qsub など）を投げても、リモートサーバ上でコマンドが見つからない場合があります。 Submit Remote Job ウィンドウの Profile --> Edit Profile... の Prefix for queueing commands に、 qsub 等の実行ファイルのパスを記入することで回避できます。例えば、 qsub のフルパスが /usr/local/bin/qsub の場合は、 Prefix for queueing commands に「/usr/local/bin/」と入力してください。
12.9.3. Q. Test Connection での接続テストは成功するが、WARNINGが表示され各種の操作に失敗します。
---------質問詳細---------
TestConnectionの結果はOKにもかかわらず、各種コマンドが実行できない。
また、リモートジョブ投入画面起動時やTestConnection実施時などで以下のダイアログが表示される。
WARNING: Putty default host name was found in registry.
(\SOFTWARE\SimonTatham\PuTTY\Sessions\Default%20Settings\HostName)
This may cause errors while job submission.
Clear this setting.
A.
原因：
このWARNINGはPuttyのHostNameが設定されているときにおこります。
Puttyの設定はWindowsのレジストリに保存されるため、Winmostar同梱版以外のPuttyであってもHostNameに何らか文字列が保存されていても、この問題がおこります。
対応：
リモートジョブ投入画面の からPuttyを起動します。Default SettingsのHostName欄に文字列が設定されているか確認します。
この文字列を削除してDefault Settingsを選択した状態でSaveすると、この問題を解消できます。
(なお、Port欄の入力内容は特に影響しません。)
12.9.4. Q. リモートサーバではどの種類のMPI（MPICH、OpenMPIなど）を使用できますか？
A. 基本的にどの種類のMPIも利用可能です。MPICH、OpenMPI、MVAPICHなどで動作実績があります。
テンプレートスクリプトを編集することで、source、module、exportといったコマンドを自由に実行し、任意のMPIを実行する環境を設定できます。
使用するソルバは、使用するMPI（mpicc, mpif90）でコンパイルされている必要があります。
12.10. MOPAC, CNDO/S, GAMESS, NWChem, Gaussianに関して
12.10.1. Q. MOPAC、GAMESS、Gaussian、NWChemのキーワード設定ウィンドウやEasy Setupウィンドウで、使用したい計算手法（Hamiltonian）、基底関数（Basis Set）が選択肢の中に見つかりません。どのように設定したらいいですか？
計算手法、基底関数の設定欄に直接キーボードで入力できる場合は、直接入力することができます。分極関数（「*」や「p」「d」「f」で表現されるもの）の指定方法はソルバごとに記述が異なることがあるので、それぞれのソルバのマニュアルを確認してください。
12.10.2. Q. 系のモデリング中や、計算結果を読み込んだ際に、不要な結合が表示されます。または結合が表示されません。どのように対処したらいいですか？
A. こちらのFAQ を参照してください。
12.10.3. Q. MOPACで計算できる最大原子数はいくつですか？
A. 重原子（水素以外）70、軽原子（水素）90です。
マニュアルページ から大分子対応版MOPAC6の実行バイナリ(最大420原子)をダウンロードして使用することもできます。
WinmostarはMOPAC2016にも対応しています。
MOPAC2016は原子数の制限はなく、学位授与機関に所属する方のみ無料です。
12.10.4. Q. MOPACのログに ATOMS **, **, AND ** ARE WITHIN .**** ANGSTROMS OF A STRAIGHT LINE と表示されて異常終了します。
---------質問詳細---------
以下のように3原子が直線になったというエラーが出て止まります。
CALCULATION ABANDONED AT THIS POINT
THREE ATOMS BEING USED TO DEFINE THE
COORDINATES OF A FOURTH ATOM, WHOSE BOND-ANGLE IS
NOT ZERO OR 180 DEGREEES, ARE IN AN ALMOST STRAIGHT
LINE. THERE IS A HIGH PROBABILITY THAT THE
COORDINATES OF THE ATOM WILL BE INCORRECT.
THE FAULTY ATOM IS ATOM NUMBER 69
最後に、
ATOMS 68, 57, AND 54 ARE WITHIN .0134 ANGSTROMS OF A STRAIGHT LINE
と出ます。
A.
角度が180°近くになる角度がZ-Matrixに含まれている場合に表示されます。
メインウィンドウ右下の座標編集機能で、接続先の原子を変更し、Z-Matrixから180°に近い角度がなくなるようにしてください。
Z-Matrixに慣れていない場合は、これ以外の方法として、キーワードに"XYZ"を追加すると、このエラーを回避できることもあります。
あるいは、3原子が直線に並ぶ線上から外れた位置に、原子種XXのダミー原子を追加し、直線に並ぶ原子のZ-Matrix上の接続先として指定することで,
エラーを回避できることもあります。
12.10.5. Q. MOPACを実行するとログに「UNRECOGNIZED KEY-WORDS: ( PM6（ハミルトニアン名）)」と出力され計算が終了してしまいます。
A. MOPACキーワード設定でHamiltonian=AM1に変えると動く場合は、使しているMOPACが対応していないハミルトニアンを選択していることによるエラーが出たことになります。
WinmostarマニュアルのMOPACの各バージョンがサポートする ハミルトニアンの一覧 をご確認の上、適切なハミルトニアンを選択してください。
それでも動かない場合は 一般的な不具合 の対処を実施してください。
12.10.6. Q. Winmostarに付属のMOPACで溶媒効果（COSMO法）を使用するにはどうしたらいいでしょうか。
Winmostarに付属のMOPACは溶媒効果（COSMO法）に対応していません。MOPACより計算時間が掛かりますがより精度の高いGAMESSを用いた計算をご検討ください。
12.10.7. Q. Winmostar内蔵以外のMOPAC(MOPAC2016など)で計算すると、一部の分子軌道しか表示されません。
MOPAC SetupウィンドウのALLVECSにチェックを入れて入力ファイルを作成して、計算をしてください。
12.10.8. Q. CNDO/SでMethod=INDOを使用したときに計算できません。
A. F以降の元素は同プログラムのMethod=INDOでサポートされていません。
Method=CNDOにするか、GAMESSなどの非経験手法を使ってください。
12.10.9. Q. GAMESS, Gaussian, NWChemのESP（静電ポテンシャル）の3次元表示を高速化することはできますか？
A. Windows版Gaussianをインストールしている場合は、Cubeファイルを開いた際に出現するCubegenウインドウにおいてCubegenチェックボックスにチェックを入れると、Gaussianに付属するCubegenプログラムを使用し比較的高速に処理することが可能になります。
将来的にはWinmostar付属のcubeファイル処理プログラム（OpenCubegen）を高速化する予定です。
12.10.10. Q. GAMESSの実行時に「* ERROR: MEMORY REQUEST EXCEEDS AVAILABLE MEMORY」などとログに出力され計算が異常終了します。
A. インプットファイルで指定したメモリ量では足りていないことを意味しています。
GAMESSキーワード設定ウィンドウのAdvancedタブの$SYSTEM欄のMWORDSの数値を大きくしてください。
12.10.11. Q. Windows版GAMESSをインストールした直後は問題なく利用できていたが、ある時期から計算が全て途中で異常終了するようになりました。対処方法はありますか？
A. まずは 一般的な不具合 の対処を実施してください。それでもうまくいかない場合は、次の操作を順に実行してください。
1. GAMESSのイントーラを起動→Removeを選択→GAMESSをアンインストール
2. Windowsの設定→アプリ→アプリと機能→Microsoft MPIを選択→Microsoft MPIをアンインストール
3. Winmostarインストールガイドにある、Windows版GAMESSインストールマニュアルの内容に従ってGAMESSとMS-MPIを再インストール
特に、GAMESS実行時に、コンソールウィンドウに「mpiexec ...(中略)... server rejected credentials」などと表示される場合はこの方法が有効な場合があります。
12.10.12. Q. GAMESSにおいて、1原子だけの系の計算で「WARNING. NUMBER OF INTERNAL COORDINATES IS GREATER THAN (3N-6), BUT NO SYMMETRY COORDINATES ARE GIVEN.」と表示されて異常終了します。
A. 原子が1個だけの系においてZ-matrixを使うことによる不具合を示すメッセージになります。
この場合は直交座標を使う（COORD=UNIQUEにする）ことで解消します。
WimostarのGAMESSキーワード設定ウィンドウにおいて、COORDをUNIQUEに変更してください。
12.10.13. Q. GAMESSのログに「 **** ERROR **** PCM SPHERE(S) MUST HAVE A POSITIVE RADIUS 」と表示され異常終了します。
A. Cavity半径がGAMESSに内蔵されていない原子が含まれている可能性があります。
Cavity半径を指定するためには、$PCM行の直後に次のステートメントを追加してください。
$PCMCAV RIN(13)=1.55, RIN(15)=1.55 $END
この例では13番目と15番目の原子にCavity半径を与えます。
12.10.14. Q. GAMESSの実行時に「ERROR: BAD DELOCALIZED COORDINATES GENERATED!!!」とログに出力され計算が異常終了します。
A. WimostarのGAMESSキーワード設定ウインドウにおいて、Z-Matrixタブを選択 --> $ZMATのチェックを外してください。
12.10.15. Q. GAMESSでNMR計算を実行しようとする（RUNTYP=NMRで計算をする）とエラーが出ます。
A. GAMESSのNMR計算は閉殻のHartree-Fock法のみ対応しており、その他DFT法等では実行できません。
また、GAMESSのNMR計算の実行速度は遅いこともあるので、代わりにNWChemやGaussianを使うことをお勧めします。
どうしてもWinmostarからGAMESSで計算したい場合は、SetupウインドウのAdvancedタブの$SCF欄のDIRSCFのチェックを全て外してください。
また並列計算に対応していないので並列数は1にしてください。
エラーで止まった際のログに詳細な指示が書かれているので、そちらも参考にしてください。
12.10.16. Q. GAMESSでDiffuse関数を使うとSCF計算が収束しません。
GAMESS SetupのBasicタブの$CONTRL枠にあるOthers欄にICUT=11を追記して、2電子積分のカットオフ値を小さく(厳しく)してください。
12.10.17. Q. 一部の分子において、Fireflyのoptimaize実行後、分子軌道UVvis…を確認すると、「'*******' is not valid floating point value」というエラーがでます。
A. 基底関数に6-31+G*とdiffuse関数の+が加わっているため、基底の線形従属性が大きくなっています。
そのため、分子軌道係数の値の一部が非常に大きくなり、ログ中に****と出力されます。
解決方法としては、
1. 6-31G*基底関数を使う
2. 6-31+G*を使うのであれば、FireflyではなくGAMESSで計算する
が挙げられます。
線形従属性の処理がGAMESSには入っているため、
FireflyとGAMESSではエネルギー値が少し異なる可能性があります。
FireflyかGAMESSどちらかで統一して、一連の計算を行ってください。
12.10.18. Q. NWChemの並列実行時に「Please specify an authentication passphrase for smpd: 」とログに出力され計算が流れません。
A. MPICH2インストール時にパスフレーズ（passphrase）を省略してしまうとそのようなエラーになる場合があります。
解決方法はいくつかありますが、MPICH2を一旦アンインストールしてから、再度インストールすると解決することがあります。
その場合は、MPICH2のアンインストール前にsmpdをストップし、MPICH2の再インストール後にsmpdをインストールする必要があります。
12.10.19. Q. NWChemの実行時に「sym_geom_project: sym_center_map is inconsistent with requested accuracy」とログに出力され計算が流れません。
A. 構造がNWChemの分子対称性判断基準からわずかに外れた場合、エラーで止まります。Winmostarメインウィンドウの[ツール]-[点群解析]を選択し、[Analyze]、その後[Symmetrize]をクリックして正確な対称性を持った構造にするか、もしくは原子1つを少し動かして対称性を崩した構造にしてください。
12.10.20. Q. NWChemで溶媒効果を入れる方法を教えてください。
A. NWChem SetupウィンドウのAdvancedタブのOther Settings欄に
cosmo
solvent ○○○
end
の3行を記入してください。○○○には溶媒分子名を入れてください。対応している溶媒分子は https://nwchemgit.github.io/COSMO-Solvation-Model.html#solvents-list-solvent-keyword をご覧ください。
12.10.21. Q. Gaussianで計算が正常に実行される時とされない時があります。
A. シングルコア版のGaussianの場合は、2つ以上ジョブを実行しようとするとGaussianが強制終了します。シングルコア版のGaussianを使い、Winmostarで2つ以上のGaussianのローカルジョブを逐次実行したいときは、各ジョブの並列数（%nprocまたは%nprocshared）を1にした上で、Winmostar JMのMax Coresを1に設定してください。
12.10.22. Q. Gaussianのlogファイルを読み込んだのですが、軌道(固有)エネルギーなどが表示されません。
A. 実行したGaussianの入力ファイルにpop=fullとgfprintが抜けている場合は表示されません。
Gaussian入力もしくは出力ファイルをWinmostarで開くと、構造だけでなく計算条件も読み込まれます。Winmostar以外で作成された入力もしくはその出力ファイルでは、pop=fullとgfprintが設定されていないことがあるため、そのようなファイルを開いてキーワード設定をする場合、pop=fullとgfprintが設定されているか確認をしてください。
12.10.23. Q. Gaussianでchkファイルを読み込んだ計算を実施する方法を教えてください。
A. リモートジョブの場合はSubmitJobウィンドウで[Advance]のチェックを入れ、[Delete *.chk]のチェックを外すとchkファイルが残され、その上でchkファイルを生成した時と同じ名前でジョブを流すとchkファイルを読み込んで計算が流れます。
--Link1-- を使う方法の方が設定自体は簡便なため、こちらの使用もご検討ください。
12.10.24. Q. 数GBのGaussian fchkファイルを読み込ませて分子軌道を表示させようとするとエラーが出ます。
A. おおむね2GB程度のfchkファイルまで対応しています。上限サイズは、分子サイズ、基底関数に依存します。将来的にはfchkファイルサイズの制限を解消する予定です。
12.11. Gromacs, LAMMPSに関して
12.11.1. Q. 分子動力学計算の実行中にエラーが起こるがどうしたらいいですか。
A. まずはWinmostarのチュートリアルをトレースしてください。それで失敗する場合は環境構築に失敗している可能性があります。
次に、エラーに再現性があるか確認してください。再現性がない場合は、環境（ハード・ソフトのいずれかまたは両方）の問題か、乱数由来の問題であることが多いです。
次に、各成分1分子、電荷なしで計算してください。それでもエラーが出る場合は力場の種類も変更してください。
各成分1分子、電荷なしでエラーが出なかったら、電荷を設定したり分子数を増やして挙動を観察してください。
次に、長時間の計算でエラーが出る場合、平衡化の時間を長くしてください。
平衡化の過程でエラーが出る場合は、以下の方法を試してください。
- 並列計算をしない。並列数を1にする。
- 温度・圧力制御にBerendsen法や速度スケーリング法を使う。
- 保存量（全エネルギー、拡張ハミルトニアン）の保存を確認し、時間刻みを小さくする。
- 倍精度を使う（Gromacsの場合）
- 拘束（SHAKE法など）を使わない。その代わり時間刻みを0.1 fs程度にする。
- 計算を短いステップ数で区切り、継続を繰り返す。（特に圧力制御時に初期密度と最終密度が大きく異なる場合に有効）
12.11.2. Q. 分子動力学計算においてどのように系を平衡化したらいいですか。
A. 低分子の平衡状態の凝集系（気体ではなく液体・固体のこと）計算が目的のケースについてまず述べます。
まず初期状態の分子を並べる際には、最終的な密度に極力近い密度に設定してください。
しかし、かなり低密度でないと並べられないときはそれで構いません。
その後、ポテンシャルエネルギー、温度、密度の変化が収束するまで、エネルギー極小化、温度一定計算、温度圧力一定計算を流してください。
初期密度が低すぎた場合は、温度圧力一定計算で、目標圧力よりも高めの圧力（例えば100倍程度）で一旦圧縮してください。
最終的にアンサンブル平均の物理量に関心があり、平衡化後に目標温度・圧力に達しているならば、細かい平衡化手順の差は計算結果に大きな影響を与えることは少ないです。
高分子、ガラスの場合は、真の意味で平衡状態を得るには、現実的な計算時間では不可能な場合がほとんどのため、エネルギー、温度、密度の収束の加え、観察したい物理量に影響が大きいと思われる物理量の相関が0に到達する程度の時間平衡化計算を実施します。
気体の場合は圧力制御は不安定なため、エネルギー極小化と温度一定計算のみで平衡状態を得ます。
12.11.3. Q. MD計算においてSHAKE法などによる拘束は計算結果にどのような影響を与えますか？拘束方法はどのように選んだらよいですか？
A. SHAKE法、RATTLE法、LINCS法、SETTLE法を共有結合する原子間に適用し結合長を拘束することで、時間刻みを大きく取り、同じ計算量でもより長時間の現象をより安定して観察できるようになります。安定、というのは、ハミルトニアン（全エネルギー）の保存の観点で、になります。
拘束しない場合に共有結合を表現する関数も実現象を高精度に表現しているわけではないので、安定した計算が流れているという前提のもと、算出される各種の物性に与える影響という点では、拘束する場合・しない場合のどちらも、それぞれの事情による実現象からのずれが生じています。
分子内の振動運動自体に計算の目的がない限りは、長時間安定してハミルトニアンが保存する条件を都度選択することを基本的には推奨します。
ただし、水素原子の結合は、拘束しない場合は系内で突出して高速に運動し、ハミルトニアンのドリフトの原因になりうるので、多くの場合は水素原子の結合については拘束します。
12.11.4. Q. MD計算を実行後、アニメーションを観たり、最終構造を見ていると、分子がセルの外側に出てしまうことがあるのはなぜですか？
A. 周期境界を使用していると、分子の実体は周期境界のセルの内側に収まるべきです。
しかし、Gromacs、LAMMPSなどのソルバは、平均二乗変位などを計算するために、セルの境界を分子が跨いでも、座標を折り返さずにそのまま並進移動した値でトラジェクトリを記録しています。
どちらにしても、結果解析時には適切に考慮され同じ結果が出力されますので、結果解析への悪影響はありません。
セルの外側に分子が飛び出る様子が見た目としてよくない場合は 表示 - 周期境界条件の表現形式 の設定を調整してください。
12.11.5. Q. [MD]-[自動で電荷を割り当て]メニューで「Topology file not found」というエラーが表示されます。解決方法はありますか？
A. 電荷を割り当てようとしている分子がAcpypeに対応していない可能性があります。[ファイル名]_charge_tmpフォルダ以下のlogファイルや、sqm.outファイルを確認してください。対応していない原子がある場合は、MOPACやQMで電荷を割り当てるか、手動で入力してください。
12.11.6. Q. 力場を割り当てで処理に時間が掛かり終わりません。解決方法はありますか？
特に力場の種類にGAFF, OPLS-AA/L+GAFFを選択した場合に処理に時間が掛かることがあります。これは、GAFF, OPLS-AA/L+GAFFの割当にacpypeを使用しており、acpypeの処理が遅いために発生しています。1分子が10000原子近い高分子の場合は、5-6時間程度待てば処理が正常終了することを確認しております。環境によってはこれよりも速い場合もあり、逆に遅い場合もあります。将来のWinmostarではこれらの力場の割当処理の高速化に取り組むことを予定しております。
12.11.7. Q. 力場を割り当てでエラーが表示されます。解決方法はありますか？
A. 力場を割り当てようとしている分子がAcpypeに対応していない可能性があります。[ファイル名]_top_tmpフォルダ以下のlogファイルや、sqm.outファイルを確認してください。対応していない原子がある場合は、UFFやDreidingなどの汎用性のある力場をお試しください。
12.11.8. Q. Gromacs, LAMMPS（分子動力学計算）から誘電率を計算できますか？
A. 誘電率は外場の周波数に依存した物性であり、また周波数帯ごとにメカニズムも違うため、一概にお答えすることはできません。
WinmostarのGromacs, LAMMPSから計算される誘電率は、分子内分極が時間変化しない前提での、分子の配向に由来する成分です。
そして、その中でも、分子動力学計算のシミュレーション時間内における系全体の双極子モーメントの揺らぎから計算される、無限に遅い低周波の極限の値となります。
ポリマーのように分子量が大きく緩和が遅い物質の場合はシミュレーション時間内に観測できる範囲での情報しかわからないため注意が必要です。
WinmostarのQuantum ESPRESSOから計算される誘電関数は、原子座標が固定された状態での電子の分極に由来する高周波成分の誘電関数です。
比較対象としている誘電率の実験値の取得方法や、材料の性質、研究目的を考えたうえで、計算をプランニングする必要があります。
なお、弊社の有償サポートでプランニングのお手伝いをすることが可能です。
12.11.9. Q. 界面ビルダで作成した構造からMD計算を実行しようとすると、力場の割り当てが失敗したりMD計算破綻します。解決方法はありますか？
A. 接合箇所の間隔が短いと、両層の分子が衝突し、望まない共有結合が原子間距離から判定され生成されてしまったり、MD計算開始時に大きな力が働きMD計算が破綻することになります。平衡化に時間が掛かってしまいますが、界面ビルダの[Direction]-[Interval]の値を大きくし、長い時間を掛けて平衡化をするのが一番良い方法となります。
12.11.10. Q. 圧力制御（NPT一定またはNPH一定）を行うと、計算が途中で破綻します。解決方法はありますか？
A. まず、気相や気相中に他の相が分散しているような、分子間の相互作用が極めて弱い状況では、圧力制御は安定しにくいため、圧力制御を使わない方法も試してください。次に、破綻した計算において密度の時間変化を確認し、何が起こっているか確認してください。また、圧力制御を行う前に、密度一定で十分エネルギー・温度・圧力が平衡化している必要があります。密度一定での平衡化が終わった時点で、圧力（の平均値）は0またはマイナスである方が望ましいです。密度一定での平衡化終了時点で圧力の値が大きいと、圧力制御を開始した直後にシステムサイズが急激に変化します。密度一定での平衡化終了時点での圧力を小さくしたい場合は、初期密度を小さくしてください（最終的な密度のおおよそ50％程度）。それでもなお解決しない場合は、(1) 圧力制御をParrinello-Rahman(Nose-Hoover)法ではなくBerendsen法に切り替える、(2) 圧力制御の時定数を大きくする、(3) 圧力制御を入れた計算を短く何回かに分割する、ということで改善するかと思われます。
12.11.11. Q. LAMMPS, Gromacsで液体（有機物）-固体（無機物）間の相互作用パラメータをどのように決めたらいいですか？
A. まず、着目している液体-固体界面について分かっている実験事実や関連研究について可能な限り情報収集します。特に、表面の原子レベルの構造（ミラー指数、官能基など）や化学的性質（濡れ性や、疎水的か親水的かといった大雑把なもの）が分かることが望ましいです。
LJパラメータと電荷の値が適切な関連研究の文献に載っている場合は、その値を使うのがベストです。Winmostarでは次のように値を設定します。
- 電荷については、一度mol2形式で保存し、電荷の値を論文の値に書き換える。
- LJパラメータについては、[力場を割り当て]機能で[Exception]から各元素の値を入力します。
適切な値が載っている文献はなかったが、ある程度化学的性質が分かっている場合は、電荷、LJパラメータを何通りか変えてシミュレーションを実行し、計算結果を比較して妥当と思われるパラメータを採用します。
表面の原子レベルの構造しか分かっていない場合は、第一原理計算（Quantum ESPRESSO）を併用してパラメータを決めます。
電荷の決定には、Quantum ESPRESSOのLowdin電荷機能を流用する場合もあります。ただし、Lowdin電荷を使う場合は電荷の合計値が0とならないため値の微調整（全体的にシフトさせるなど）が必要です。また、Lowdin電荷では分極が過大評価されることもあり注意が必要です。
LJパラメータについては、既知のパラメータ（Dreiding, UFF, CLAYFFなど）を使うか、第一原理計算からForce Matchingなどのアルゴリズムを用いて算出します。
このように液体（有機物）-固体（無機物）間の相互作用パラメータが複雑となっているのは、以下の事実に由来します。
1. 古典MDでは原子位置に相互作用パラメータ（epsilon, sigma, 電荷）が依存しないという仮定を置くが、有機物-無機物界面系では有機物-有機物系に比べ、その近似が大きな誤差を生じる。
2. 現実のデバイス中の無機物表面は酸化膜などに覆われているが、実験観察も容易ではなく、原子解像度で正確にモデリングすることが難しい。
12.11.12. Q. Gromacsの最終構造やアニメーションを読み込んで取得した構造や、それに何かしらの編集を行ってから再びMD計算を実行すると力場の割り当てが失敗したりMD計算が破綻します。解決方法はありますか？
A. Gromacsの不具合により分子が分離している可能性があります。分子が分離しているか否かは、[選択]-[分子種によるグループ選択]を表示し、想定していない成分が含まれているかを確認することで判断できます。分離していた場合は Gromacsのトラジェクトリを読み込んだ構造または、最終構造において、本来繋がっているべき分子がバラバラに表示される。 をご参照ください。
12.11.13. Q. Gromacsにおいて最終構造やアニメーションを読み込むと分子がバラバラに表示されることがありますがなぜですか？
12.11.14. Q. GromacsのER法結果読み込みを実行しても結果が表示できません／エラーが出ます。
A. ER法を実行する際に指定した出力先ディレクトリに生成されるermod.outの内容を確認してください。
ermod.outの中に「 The minimum of the energy coordinate is too large; the ecdmin parameter needs to be smaller」と書かれている場合は、ER法実行ウィンドウの[Options]ボタンを押し、[For Solution System]のと[minimum value of the solute-solvent energy (ecdmin)]の値を小さくしてください。
具体的な値の設定方法など、詳しくは ERmodのwikiのFAQ を参照してください。
また、同様にermod.outの内容と ERmodのwikiのFAQ全般 の内容を照らし合わせ、ermodの設定の変更が必要な場合はER法実行ウィンドウの[Options]で設定してください。
12.11.15. Q. ある物質とある物質の間の相互作用を計算することはWinmostarで可能ですか？
A. この手の質問は大変多いのですが、「相互作用」という言葉の定義は広いため、定義によります。
まずは、着目している物質の量子化学、分子動力学、第一原理計算自体を実行できるか、という意味では、各ソルバ（GAMESS, Gaussian, LAMMPS, Gromacs, Quantum ESPRESSO）で実行可能な内容にWinmostarは依存しているので、各ソルバのマニュアルを予め調べてください。
次に、何か着目している物性や現象のメカニズムを知りたいという意味での相互作用については、その物質系固有の知識が必要となるため、計算可能かどうかの調査自体に十分な調査が必要となるため、すぐには回答できません。民間企業向けに有償で調査を行うサービスも用意しています。
最後に、Kitaura-Morokuma解析やMDにおけるCoulomb・vdW相互作用などの相互作用エネルギー解析のように、具体的な解析内容が分かっている場合は、その旨をご質問ください。
12.11.16. Q. 特定原子の軌跡をプロットする方法を教えてください。
A. 以下のように操作してください。
1) アニメーションを表示
2) 分子表示エリアで軌跡を描画したい原子をクリック
3) アニメーションウィンドウで[Custom Plot]をクリック
4) [X axis]の2番目の項目をチェックし、「Atomic position (x)」を選択（適宜適切な方向を選択）
5) [X axis]の[Apply]をクリック
6) [Y axis]の2番目の項目をチェックし、「Atomic position (y)」を選択（適宜適切な方向を選択）
7) [Y axis]の[Add]をクリック
8) [Y axis]の「Atomic position (y)」以外のチェックを外す
9) グラフの下の[X Axis]と[Y Axis]の[Autoscale]を外し[Min]と[Max]を適宜設定
なお、以上の操作を簡略化した機能を開発予定です。
12.12. Quantum ESPRESSO, OpenMXに関して
12.12.1. Q. 手順通りQuantum ESPRESSOの擬ポテンシャルファイルをインストールしたが認識されません。
A. 拡張子が.UPFの擬ポテンシャルファイル（例えば
O.pw-mt_fhi.UPF 、 Au.pbe-dn-rrkjus_psl.0.1.UPF など）をQuantum ESPRESSOのインストールフォルダの下のpseudoフォルダ（デフォルトではC:\Program Files\Quantum ESPRESSO 64-bit 5.2.1pseudoなど）にコピーしてください。ただし、旧Internet ExprolerなどのブラウザでUPFファイルをダウンロード・保存すると、拡張子が勝手に変更されたりと不具合が報告されているので、EdgeやChromeなどのブラウザも試してください。12.12.2. Q. Quantum ESPRESSOの擬ポテンシャルファイルの探し方が分かりません。どのように探したらいいですか？
A. Windows版Quantum ESPRESSOインストールマニュアル にて紹介しています。
12.12.3. Q. Quantum ESPRESSOを用いた計算が失敗します。計算結果の表示でエラーが出ます。
A. まずは 一般的な不具合 の対処を実施してください。
次に、WinmostarではQEの各モジュールをバッチ処理で連続実行しているので、Winmostarが生成したbatファイル（ローカル実行の時）またはshファイル（リモート実行の時）に記述された処理の流れを見ながら、生成された出力ファイル（pwoutまたはout）ファイルを順番に確認してください。
例えば、フォノン計算の場合はph.xの出力ログ（ph.out）を確認してください。
最初に「Error in routine ...」などのエラーが出現した箇所の対処を施し、再度ジョブを実行してください。
特定のキーワードに関するエラーは、そのキーワードの設定を 公式サイト でご確認ください。
典型的なQEのエラーの対処方法は 公式サイトのFAQ に記載されています。
12.12.4. Q. Quantum ESPRESSOを用いてPhonon計算を実行する際に、ph.xの出力(ph.out)に「third order derivatives not implemented with GGA」と表示され計算結果を取得できません。
A. GGAでない擬ポテンシャルを選択することで解消します。
12.12.5. Q. Quantum ESPRESSO, OpenMXのSCF計算または構造最適化計算が収束しません。
A. 以下の対策を順に実施してください。
必ず試すべきこと：
・第一原理計算は設定項目が多いので、適当に計算条件を変えず、きちんと記録を取りながら一連の計算を流す。
・QEの一般的な不具合 の対処を実施する。
・本当に収束しない傾向にあるかチェックする。
・QEではEstimated accuracyをSCFサイクル数に対しプロットする。両対数プロットならなおよし。
・スピン分極状態・電荷が妥当か調べる。
・Hexagonal結晶でK_POINTSのシフトを行っていたらシフトを外す。
・up/downスピンの並び方を与える。
・系全体の磁気モーメントを拘束する。
・尤もらしい初期構造を使う。
・実験や他の計算手法で得られた構造を使う。
・計算する上で配置に任意性のある原子（X線で見えない軽元素、固溶体、欠陥、非整数の組成など）がある場合は、違う配置を試す。
・固溶体・欠陥を含むようなケースでは、系内に大きなダイポールモーメントが生じないような初期構造にする。
次に試すこと：
・mixing_modeを調整する。
・smearingを使っていない場合はsmearingを使う。
・擬ポテンシャルの種類を変える。
・スピン分極の初期値を調整する。（原子単位または系全体）
・外部電場、欠陥、吸着など比較的複雑な条件を設定している場合は、それらをなくしたよりシンプルな条件で試し、その計算が収束したなら、その計算の終状態（原子配置・波動関数など）を始状態として計算を開始する。
・収束しなかった計算の途中から計算を開始する（SCFのアルゴリズムは履歴に依存するため）。
・行列計算のパラメータを調整する（収束しづらい設定のみ見直す）。
・スラブに分子が吸着するような、系内に大きなダイポールモーメントが発生してしまう場合は、ダイポールの補正を行う。
計算時間・計算精度との兼ね合いで試すこと：
・カットオフエネルギーを大きく取る。
・K点を多めにとる。
・smearingを調整する（種類・幅）。
・波動関数の更新度合（QEではmixing_beta）を小さくする。
計算精度との兼ね合いで試すこと：
・SCFの収束パラメータを緩くする。
12.12.6. Q. Quantum ESPRESSOのSCF計算が出力ファイル（.pwoutまたは.out）に「too few bands」と表示され異常終了します。nbndの設定方法が分かりません。
A. まずは QE公式のマニュアルのnbndの説明 をご確認ください。
nbndを使わずに計算を流すと、QEが自動でnbndを適当に設定して計算するので、Winmostarのキーワード設定画面で「Use nbnd」のチェックを外してください。
nbndを増やしたい場合は、nbndを使わずに実行したときにpwoutまたはoutファイルに出力される"number of Kohn-Sham states"の値よりも大きい値をnbndに設定してください。
また、Winmostarのキーワード設定画面の「Use nbnd」のところに表示される「# valence bands: 」の値も参考にしてください（詳細は 固体 ‣ Quantum ESPRESSO メニュー を参照）。
12.12.7. Q. Quantum ESPRESSOのSCF計算が出力ファイル（.pwoutまたは.out）に「fixed occupations and lsda need tot_magnetization」と表示され異常終了します。どのように解決したらいいですか？
A. occupationsにsmearingを指定するか、starting_magnetizationではなくtot_magnetizationを指定してください。
12.12.8. Q. Quantum ESPRESSOのSCF計算が出力ファイル（.pwoutまたは.out）に「charge is wrong」と表示され異常終了します。どのように解決したらいいですか？
A. まずQuantum ESPRESSOキーワード設定ウィンドウでoccupationsがsmearingになっていることを確認し、次にecutrhoを大きめ（400 Ryなど）に設定してください。
12.12.9. Q. Quantum ESPRESSOを用いて誘電関数を計算する際に、epsilon.xの出力(eps.out)に「bad band number」と表示され誘電関数を取得できません。
A. SCF計算でバンド数（nbnd）を増やすことで解消します。
12.12.10. Q. Quantum ESPRESSOを用いて誘電関数を計算する際に、epsilon.xの出力(eps.out)に「USPP are not implemented」と表示され誘電関数を取得できません。
A. SCF計算でノルム保存型の擬ポテンシャルを選択することで解消します。
12.12.11. Q. Quantum ESPRESSOを用いてPhonon計算に失敗し、計算結果を取得できません。
A. まずは、Phonon計算を実行せず、同じ擬ポテンシャルファイル、その他計算条件を使用してSCF計算が正常終了することを確認してください。次に、working directory（末尾が_qe_dataのフォルダ）の中のph.out(Phonon計算モジュールph.xの出力ファイル)を確認してください。そこに「The phonon code with US-PP and raman or elop not yet available」と書かれている場合は、ノルム保存型の擬ポテンシャルを選択することで解消します。同様に、PAWポテンシャルを利用している場合もラマン計算などがサポートされていないので、ノルム保存型の擬ポテンシャルを選択することで解消します。
12.12.12. Q. フェルミ面を出力しようとしてもそれらしきものが表示されません。
A. まず、可能なら対象の物質が金属であることを確認してください。次に、状態密度も出力し、フェルミエネルギーにおいて状態密度が0でないことを確認してください。
12.12.13. Q. Quantum ESPRESSO（バンド計算）から誘電率を計算できますか？
A. Q. Gromacs, LAMMPS（分子動力学計算）から誘電率を計算できますか？ を参照してください。
12.12.14. Q. Quantum ESPRESSOで汎関数の種類はどのように設定しますか？
A. Quantum ESPRESSOでは、汎関数ごとに擬ポテンシャルファイルが作られるので、基本的には擬ポテンシャルファイルを選んだ時点で汎関数が決定されます。一部の汎関数（HSE、vdw汎関数など）は、ベースとなる汎関数（例えばHSEの場合はPBE）で作られた擬ポテンシャルファイルを選択した上で、input_dftキーワードを使用して汎関数の設定を上書きします。
12.12.15. Q. Quantum ESPRESSOでセルサイズの構造最適化を行うと最終構造のエネルギーで不連続的な変化をしたりエネルギーが増加するのはなぜですか？
A. エネルギーの変化が最後にわずかに不連続になるのは、 Quantum ESPRESSOが、Pulay圧力に起因する問題を避けるために、最後の計算だけ利用する平面波セットを取り直しているからです。
このエネルギーの飛びは、カットオフエネルギーを大きくすることである程度軽減できます。
初期構造から最終構造が大きく異なる場合は、続けて再度構造最適化計算を実行します。
12.12.16. Q. 系のモデリング中や、計算結果を読み込んだ際に、不要な結合が表示されます。または結合が表示されません。どのように対処したらいいですか？
A. こちらのFAQ を参照してください。
12.12.17. Q. OpenMXでMPIを有効にしてローカルマシンで計算を実行すると、 tcp_peer_send_blocking: send() to socket 12 failed: Transport endpoint is not connected というエラーが表示されます。
A. CygwinのOpenMPI特有の問題で、Windowsの[設定]-[ネットワークとインターネット]-[アダプターのオプションを変更する]において使用していないネットワークアダプタを無効にしてください。また、OpenMXはローカルマシンにおいてはOpenMPで計算することを推奨します。
12.13. アドオンに関して
12.13.1. Q. 溶解度パラメータ計算モジュールを用いてポリマーのハンセン溶解度パラメータを計算する際に、ポリマーの繰り返し構造（モノマー）の取り方の違いで出力される値が変化してしまう。
A. 実装されている原子団寄与法のアルゴリズムのために発生しています。原子団を探索する際には、一番大きな原子団から探索されるようになっています。重要そうな官能基は繰り返し単位の中に入れておくことをお勧めします。
12.13.2. Q. 溶解度パラメータ計算モジュールを用いて取得したハンセン溶解度パラメータの値が、文献値と大きく異なります。
A. 溶解度パラメータ計算モジュールは、各種の文献値を学習データとしてニューラルネットワークで学習された原子団寄与法を用いてハンセン溶解度パラメータを出力しています。そのため、文献値と全く同じ値を返すわけではありません。また、文献によっては溶解度パラメータの単位が異なりますので、その点にご注意ください。