1 <?xml version="1.0" encoding="utf-8"?>
2 <!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Strict//EN" "http://www.w3.org/TR/xhtml1/DTD/xhtml1-strict.dtd">
3 <html xmlns="http://www.w3.org/1999/xhtml">
5 <meta http-equiv="Content-Type" content="text/html;charset=utf-8" />
7 <link rel="stylesheet" href="../etc/style.css" type="text/css" />
12 <div file="index.html">
14 <div class="centered" lang="en">
16 <h2>An Interactive Molecular Modeling Software<br />with Integrated Ruby Interpreter</h2>
17 <h3>Version 0.6.2</h3> <!-- version -->
19 <p><a href="http://molby.sourceforge.jp/index-en.html">http://molby.sourceforge.jp/index-en.html</a></p>
21 <div class="centered" lang="ja">
23 <h2>対話型分子モデリングソフトウェア<br />(Ruby インタプリタ内蔵)</h2>
24 <h3>Version 0.6.2</h3> <!-- version -->
26 <p><a href="http://molby.sourceforge.jp/index.html">http://molby.sourceforge.jp/index.html</a></p>
29 <div class="contents" lang="en">
30 <h2>Table of Contents</h2>
31 <link id="#toc_all" />
33 <div class="contents" lang="ja">
35 <link id="#toc_all" />
37 <link id="#navigation" />
40 <div file="introduction.html">
41 <link rel="Up" href="index.html" />
43 <div class="contents" lang="en">
47 <div class="contents" lang="ja">
51 <link id="#navigation" />
54 <div file="whatis.html">
55 <link rel="Up" href="introduction.html" />
57 <div class="contents" lang="en">
58 <h1>What is Molby?</h1>
60 Molby is an application to visualize and build molecular models interactively through 3D display. It provides the following capability:
63 <li>Viewing molecular models in either ball-and-stick or "line" display.</li>
64 <li>Interactive modification of molecular models by mouse operation.</li>
65 <li>Optimization of the structure by molecular mechanics calculation. Also capable is the molecular dynamics calculation.</li>
66 <li>Viewing and editing atomic coordinates, molecular mechanics atom typs, and other parameters in tables.</li>
67 <li>Importing and exporting structure/coordinate data for other program packages such as Gaussian, GAMESS, NAMD and SHELX.</li>
68 <li>Embedded Ruby interpreter for automated processing of molecular models.</li>
71 <div class="contents" lang="ja">
74 Molby は、分子を 3D 表示して、画面上で分子モデルを構築するためのアプリケーションです。以下のような機能があります。
77 <li>分子を ball-and-stick または「線」で表示できます。</li>
78 <li>分子モデルをマウス操作で編集できます。</li>
79 <li>分子力学計算による構造最適化ができます。分子動力学計算も可能です。</li>
80 <li>原子座標や分子力学の原子タイプなどのパラメータを、テーブルで編集できます。</li>
81 <li>Gaussian, GAMESS, NAMD, SHELX などの他のプログラムパッケージのファイルを読み込み/書き出しできます。</li>
82 <li>分子モデルの操作を自動化できるように Ruby インタプリタを内蔵しています。</li>
85 <link id="#navigation" />
88 <div file="installation.html">
89 <link rel="Up" href="introduction.html" />
90 <link rel="Prev" href="whatis.html" />
92 <div class="contents" lang="en">
94 <h2>1. Supported Platforms</h2>
96 Molby runs on the following platforms.
99 <li>Microsoft Windows (XP or later).</li>
100 <li>Mac OS X (10.4 or later, PowerPC or Intel).</li>
104 Installation procedures are described below for each of these platforms.
107 <h2>2. Get the Software</h2>
109 Download from the official distribution web site, <a href="http://en.sourceforge.jp/projects/molby/">http://en.sourceforge.jp/projects/molby/</a>.
112 <h2>3. Installation</h2>
114 <h3>3-1. Microsoft Windows</h3>
117 The Windows version is provided as a standard setup package (<code>SetupMolby.exe</code>). Double-click the setup package, and follow the instructions. If you are not sure, just select "OK" to go proceed.
120 After installation is finished, you will find the Molby application registered in the "Start" menu under the item "All Programs."
123 <h3>3-2. Mac OS X</h3>
125 The Mac version is provided as a disk image (<code>Molby.dmg</code>). Double-click the disk image file, and you will find a virtual disk drive mounted on the desktop. Find the Molby application inside it, and drag it to the "Applications" folder in your hard drive.
128 <span class="italic">Note:</span> The Mac version is provided as a universal binary, which runs natively on both PowerPC and Intel platforms.
131 <h2>4. Uninstallation</h2>
132 <h3>4-1. Microsoft Windows</h3>
134 Use the uninstaller in the Molby folder. You can access the Molby folder from the "Start" menu (Start -> All Programs -> Molby).
136 <h3>4-2. Mac OS X</h3>
138 Trash the Molby application in the "Applications" folder.
142 <div class="contents" lang="ja">
144 <h2>1. サポートするプラットフォーム</h2>
146 Molbyは以下のプラットフォームで動作します。
149 <li>Microsoft Windows (XPまたはそれ以降).</li>
150 <li>Mac OS X (10.4またはそれ以降, PowerPCまたはIntel).</li>
154 それぞれのプラットフォームでのインストール方法を、以下に説明します。
157 <h2>2. ソフトウェアの入手</h2>
159 公式配布サイト <a href="http://sourceforge.jp/projects/molby/">http://sourceforge.jp/projects/molby/</a> から最新版をダウンロードしてください。
164 <h3>3-1. Microsoft Windows</h3>
167 Windows版は、標準のセットアップパッケージ (<code>SetupMolby.exe</code>) で配布されています。セットアップパッケージをダブルクリックして、指示に従ってください。どうしたらよいかわからない時は、"OK"を押して進んで構いません。
170 インストールが終了したら、「スタート」メニューの「すべてのプログラム」の中に Molby が入っているはずです。
173 <h3>3-2. Mac OS X</h3>
175 Mac版は、ディスクイメージ (<code>Molby.dmg</code>) で配布されています。ディスクイメージをダブルクリックすると、仮想ディスクがデスクトップにマウントされます。Molby アプリケーションがその中にありますので、ハードディスクの「アプリケーション」フォルダにドラッグコピーしてください。
178 <span class="italic">注:</span> Mac 版は universal binary です。PowerPC, Intel の両方のマシンで動作します。
182 <h3>4-1. Microsoft Windows</h3>
184 Molbyフォルダの中にアンインストーラがあります。「スタート」メニューで「すべてのプログラム」-> Molby とたどってください。
186 <h3>4-2. Mac OS X</h3>
188 「アプリケーション」フォルダの中の Molby アプリケーションをゴミ箱に捨てます。
192 <link id="#navigation" />
195 <div file="copyright.html">
196 <link rel="Up" href="introduction.html" />
197 <link rel="Prev" href="installation.html" />
198 <link id="#header" />
199 <div class="contents" lang="en">
200 <h1>Copyright and License</h1>
201 <h2>1. Copyright Notice</h2>
203 Molby is a copyrighted product of Toshi Nagata.
206 Copyright (C) 2009-2012 Toshi Nagata <!-- copyright -->
209 Molby includes the following softwares, which are copyrighted products as described below:
213 <a href="http://ambermd.org/">AmberTools 1.3</a>: Copyright (C) Junmei Wang, Ross C. Walker, Michael F. Crowley, Scott Brozell and David A. Case
216 <a href="http://www.wxwidgets.org/">wxWidgets 2.8.9</a>: Copyright (C) 1992-2008 Julian Smart, Robert Roebling, Vadim Zeitlin and other members of the wxWidgets team. Portions (c) 1996 Artificial Intelligence Applications Institute
219 <a href="http://www.ruby-lang.org/">Ruby 1.8.7</a>: Copyright (C) 1993-2009 Yukihiro Matsumoto
222 <a href="http://www.netlib.org/clapack/">CLAPACK</a>: Copyright (C) 1992-2008 The University of Tennessee. All rights reserved.
228 Molby is distributed under the GNU General Public License (version 2).
232 Molby: An Interactive Molecular Modeling Software with Integrated Ruby Interpreter
235 Copyright (C) 2009-2012 Toshi Nagata <!-- copyright -->
238 This program is free software; you can redistribute it and/or modify
239 it under the terms of the GNU General Public License as published by
240 the Free Software Foundation; either version 2 of the License, or
241 (at your option) any later version.
244 This program is distributed in the hope that it will be useful,
245 but WITHOUT ANY WARRANTY; without even the implied warranty of
246 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
247 <a href="../etc/gpl.txt">GNU General Public License</a> for more details.
251 <div class="contents" lang="ja">
255 Molby の著作権は永田 央が保持しています。
258 Copyright (C) 2009-2012 Toshi Nagata <!-- copyright -->
261 Molby は以下のソフトウェアを含んでいます。それぞれの著作権表示は下の通りです。
265 <a href="http://ambermd.org/">AmberTools 1.3</a>: Copyright (C) Junmei Wang, Ross C. Walker, Michael F. Crowley, Scott Brozell and David A. Case
268 <a href="http://www.wxwidgets.org/">wxWidgets 2.8.9</a>: Copyright (C) 1992-2008 Julian Smart, Robert Roebling, Vadim Zeitlin and other members of the wxWidgets team. Portions (c) 1996 Artificial Intelligence Applications Institute
271 <a href="http://www.ruby-lang.org/">Ruby 1.8.7</a>: Copyright (C) 1993-2009 Yukihiro Matsumoto
277 Molby は <a href="../etc/gpl.txt">GNU General Public License (GNU 一般公衆利用許諾契約書, version 2)</a> に従って配布します。
281 Molby: 対話型分子モデルソフトウェア(Ruby インタプリタ内蔵)
284 Copyright (C) 2009-2012 Toshi Nagata <!-- copyright -->
287 このプログラムはフリーソフトウェアです。あなたはこれを、フリーソフトウェア財団によって発行された GNU 一般公衆利用許諾契約書(バージョン2か、希望によってはそれ以降のバージョンのうちどれか)の定める条件の下で再頒布または改変することができます。
290 このプログラムは有用であることを願って頒布されますが、*全くの無保証* です。商業可能性の保証や特定の目的への適合性は、言外に示されたものも含め全く存在しません。詳しくはGNU 一般公衆利用許諾契約書をご覧ください。
294 参考のため、GNU 一般公衆利用許諾契約書の<a href="../etc/gpl.ja.txt">非公式日本語訳</a>を添付します。ただし、正式な文書は英語版の方です。
297 <link id="#navigation" />
300 <div file="tutorials.html">
301 <link rel="Up" href="index.html" />
302 <link rel="Prev" href="introduction.html" />
303 <link id="#header" />
304 <div class="contents" lang="en">
308 <div class="contents" lang="ja">
312 <link id="#navigation" />
315 <div file="drawing.html">
316 <link rel="Up" href="tutorials.html" />
317 <link id="#header" />
318 <div class="contents" lang="en">
319 <h1>Step One: Draw a New Molecule</h1>
320 <h2>1. Create a Molecule</h2>
322 Start up Molby. You will see a blank window like below. This is the main window to manipulate your molecule.
324 <p><img src="../etc/tutorial_01.png" /></p>
326 Suppose you want a benzene molecule. Double-click on the black area, and you will be asked to "enter formula". Type "C6H6" and hit "OK".
328 <p><img src="../etc/tutorial_02.png" /></p>
330 Then you see a benzene molecule in the black area. The atoms and bonds are drawn in red, which means these atoms and bonds are "selected". On the left to this area, there is a table, which shows in numbers the positions of all atoms, and other useful information. We will get into this part in more details later.
332 <p><img src="../etc/tutorial_03.png" /></p>
333 <h2>2. Rotate the Molecule</h2>
335 You can rotate the molecule around by use of the three bars on the both sides and the bottom of the black area. Press the left button of the mouse on one of these bars, and try dragging.
337 <p><img src="../etc/tutorial_04.png" /></p>
339 <li>Bar 1: rotate the whole molecule around the horizontal axis.</li>
340 <li>Bar 2: rotate the whole molecule around the vertical axis.</li>
341 <li>Bar 3: rotate the selected part of the molecule. (In the present example, the whole molecule is selected, so that the bar 3 works in the same manner as the bar 2.)</li>
344 Next, look at the buttons above the black area. The left-most button "Rot" is now selected. This means the window is in the "Rotation" mode. In this case, you can drag in the black area to rotate freely the whole molecule.
347 <span class="italic">Note:</span> Be careful not to start drag from a very near point to the molecule. If you start drag on a selected atom or bond, the selected part will move along the mouse, instead of the whole molecule to rotate.
350 <h2>3. Translate and Scale</h2>
352 Press the button "Trans" above the black area. The window is now in "Translate" mode. When you drag in the black area, the whole molecule moves along with the mouse movement.
355 <span class="italic">Note:</span> the functions of the rotation bars do not change.
357 <p><img src="../etc/tutorial_05.png" /></p>
359 Next, press the button "Scale". The window is now in "Scale" mode, and you can expand and shrink the molecule by dragging in the black area.
361 <p><img src="../etc/tutorial_06.png" /></p>
363 The "Fit to Screen" menu command, avaiable under the "Show" menu, is a convenient way to fit the whole molecule to the window.
365 <p><img src="../etc/tutorial_07.png" /></p>
367 Now you have learned how to move the whole molecule around. In the next chapter, you will learn how to edit the molecule.
371 <div class="contents" lang="ja">
372 <h1>第一段階:新しい分子を描く</h1>
375 Molbyを立ち上げてください。下のような空のウィンドウが現れます。このウィンドウが、分子を取り扱うためのメインウィンドウです。
377 <p><img src="../etc/tutorial_01.png" /></p>
379 ベンゼン分子を作ってみたいとします。画面の黒いところ(編集エリア)をダブルクリックすると、"enter formula" と指示が出ます。"C6H6" とタイプして、OK を押します。
381 <p><img src="../etc/tutorial_02.png" /></p>
383 編集エリアにベンゼン分子が現れます。原子と結合は赤色で描かれていますが、これは「選択されている」という表示です。左側には表があり、原子の位置などの情報を数値で表示しています。この部分については、あとで詳しい使い方が出てきます。
385 <p><img src="../etc/tutorial_03.png" /></p>
388 編集エリアの両側と下には、分子を回転させるバーがあります。このバーでマウスボタンを押し、そのままドラッグしてみてください。
390 <p><img src="../etc/tutorial_04.png" /></p>
392 <li>バー1:水平軸の周りに分子を回転させる。</li>
393 <li>バー2:垂直軸の周りに分子を回転させる。</li>
394 <li>バー3:分子の選択部分を回転させる。(今は分子全体が選択されているため、バー2と同じ動作をします。)</li>
397 次に、編集エリアの上にあるボタンを見てください。一番左の "Rot" が選択されています。これは、このウィンドウが現在「回転 (Rotation)」モードにあることを示しています。この場合、編集エリアでマウスをドラッグすると、分子全体を自由に回転させることができます。
400 <span class="italic">注:</span> 分子に近すぎるところからドラッグを始めないように注意してください。選択している原子・結合の上でドラッグを始めると、選択部分がマウスについて動いてしまい、分子全体を動かすことができません。
402 <h2>3. Translate and scale</h2>
404 "Trans" と書かれたボタンを押してください。ウィンドウは「並進 (Translate)」モードになります。編集エリアでマウスをドラッグすると、分子全体を平行移動させることができます。
407 <span class="italic">注:</span> 回転バーの機能は変わりません。
409 <p><img src="../etc/tutorial_05.png" /></p>
411 次に、"Scale" と書かれたボタンを押してください。ウィンドウは「拡大縮小 (Scale)」モードになります。編集エリアでマウスをドラッグすると、分子全体を拡大・縮小することができます。
413 <p><img src="../etc/tutorial_06.png" /></p>
415 "Show" メニューの中に "Fit to Screen" メニューコマンドがあります。このコマンドは、分子をウィンドウサイズに合わせて表示したいときに便利です。
417 <p><img src="../etc/tutorial_07.png" /></p>
419 分子を画面上で動かす方法がわかりました。次は、編集のやり方を説明します。
422 <link id="#navigation" />
425 <div file="editing.html">
426 <link rel="Up" href="tutorials.html" />
427 <link rel="Prev" href="drawing.html" />
428 <link id="#header" />
429 <div class="contents" lang="en">
430 <h1>Step Two: Edit a Molecule</h1>
431 <h2>1. Append a Substituent</h2>
433 Suppose you have a benzene molecule now, and want to append a methyl substituent to give a toluene.
436 Press the button "Select", and click on the black area where no atom nor bond is present. The selection becomes canceled.
438 <p><img src="../etc/tutorial_08.png" /></p>
440 Select one of the hydrogen atoms. There are two ways to do this. You can click on the target atom, or you can drag to select atoms within a rectangular region.
442 <p><img src="../etc/tutorial_09.png" /> <img src="../etc/tutorial_10.png" /></p>
444 Double-click on the hydrogen atom. This time, be sure to hit the selection! Then you will see the dialog "enter formula" as shown below. Type "CH3" in the text box and hit "OK".
446 <p><img src="../etc/tutorial_11.png" /></p>
450 <p><img src="../etc/tutorial_12.png" /></p>
452 You want to save the result to a file. Select "Save as..." from the File menu, choose the file format, enter the filename, and hit "Save".
455 The "native" file format in Molby is "mbsf", i.e. <u>M</u>ol<u>b</u>y <u>S</u>tructure <u>F</u>ile. The mbsf format is a private format with no compatibility with any existing file format for chemistry. Nevertheless, it is a plain text file that is reasonably compact and easy to read. Unless you have a particular reason to choose other formats, it is recommended that you save the structues as a "mbsf" file.
457 <h2>2. Manipulate Molecular Fragments</h2>
459 After building a molecule, you may want to move some part of a molecule. Here are examples how you can do it.
462 Select the methyl group of the toluene, if you have not done so yet <span class="note">(Tip: you can add to the current selection by shift-clicking)</span>. Press the left rotation bar (marked "3" in the previous page), and drag up and down. You will see the methyl group rotates to the right and left.
464 <p><img src="../etc/tutorial_13.png" /></p>
466 The same manipulation can be achieved by selecting one bond and dragging the left rotation bar. In this case, you can also rotate the benzene ring (with the methyl group fixed) by pressing Option (Mac) or Alt (Win) key when dragging.
468 <p><img src="../etc/tutorial_14.png" /></p>
470 The selected atoms can be dragged to give a translational move. Use this feature with care, because it may result in a chemically unnatural structure.
472 <p><img src="../etc/tutorial_15.png" /></p>
473 <h2>3. Adding and Deleting Atoms</h2>
475 We have already learned one way to add atoms: double-clicking on the selected part of the molecule (or on the black area where nothing is present), and type-in the formula. Actually, this is the most convenient way to add atoms in Molby. However, there are also other ways to add atoms.
478 Look at our familiar toluene molecule. Suppose we want to convert it to indane.
480 <p><img src="../etc/indane.png" /> indane</p>
482 Rotate the molecule so that its orientation matches the chemical structure. We will start from C3, add two carbons, and then close the ring.
484 <p><img src="../etc/tutorial_16.png" /></p>
486 Press the "Erase" button above the black area. Click on the "H3" hydrogen. The hydrogen atom and the bond between C3 and H3 disappear.
488 <p><img src="../etc/tutorial_17.png" /></p>
490 Press the "Bond" button. Drag from the C3 atom to the right-bottom, and release the mouse button. A new atom and a bond to C3 are created. Drag from the new atom to the right-up, and another atom and bond are created.
492 <p><img src="../etc/tutorial_18.png" /> <img src="../etc/tutorial_19.png" /></p>
494 Press the "Erase" button again, and erase one hydrogen from the methyl group. Press the "Select" button, select the benzene-methylene bond, and rotate the methylene group as appropriate for the five-membered ring, by use of the left rotation bar.
496 <p><img src="../etc/tutorial_20.png" /> <img src="../etc/tutorial_21.png" /></p>
498 Now you can close the ring. Press the "Bond" button, and drag from the methylene carbon. When the mouse cursor comes close enough to the target carbon, the new atom will snap to the target and make a new bond. Then release the mouse button.
500 <p><img src="../etc/tutorial_22.png" /></p>
502 Finally, you want to add hydrogens to the newly created carbons. Press the "Select" button, and select these two carbons. Go to the "Edit" menu, and select the "Add Hydrogen" -> "Tetrahedral sp3" menu command.
504 <p><img src="../etc/tutorial_23.png" /></p>
508 <p><img src="../etc/tutorial_24.png" /></p>
510 A similar result can be achieved by selecting the hydrogen ortho to the methyl group, double-click it, type "CH2CH3" in, erase one hydrogen atom from each of the methyl groups, and make a bond. This is better because the newly created methylenes have reasonable bond lengths and angles as methylene groups. The above example is just for explanation of the editing features.
513 <div class="contents" lang="ja">
514 <h1>第二段階:分子を編集する</h1>
517 ベンゼン分子があり、これにメチル基をつけてトルエンを作りたいとします。
520 "Select"ボタンを押して、編集エリアの原子や結合がないところをクリックします。選択が解除され、分子の表示が赤から原子ごとの色に変わります。
522 <p><img src="../etc/tutorial_08.png" /></p>
524 水素原子の一つを選択します。やり方は二つあります。選択したい原子をクリックするか、またはドラッグで現れる四角い領域の中の原子を選択します。
526 <p><img src="../etc/tutorial_09.png" /> <img src="../etc/tutorial_10.png" /></p>
528 水素原子の上でダブルクリックします。このとき、選択部分の上でダブルクリックするように注意してください。"Enter formula"というダイアログが出てきます。"CH3"と入力し、"OK"を押します。
530 <p><img src="../etc/tutorial_11.png" /></p>
534 <p><img src="../etc/tutorial_12.png" /></p>
536 ファイルに保存しておきたい時は、通常通り File メニューから "Save as..." を選び、ファイルフォーマットを選び、ファイル名をタイプして "Save" を押します。
539 Molby の標準ファイルフォーマットは "mbsf" (<u>M</u>ol<u>b</u>y <u>S</u>tructure <u>F</u>ile) です。Mbsf フォーマットは独自のもので、既存の化学ファイルフォーマットとは互換性がありません。しかし、mbsfはテキストファイルであり、比較的コンパクトで、読みやすいものです。他のフォーマットで保存する積極的な理由がない場合は、"mbsf"フォーマットで保存しておくことをおすすめします。
542 <h2>2. 分子の一部を操作する</h2>
544 分子を組み立てたら、その一部を動かしたくなることがあります。ここでは、そのやり方の例を示します。
547 トルエンのメチル基を選択します<span class="note">(コツ:シフトキーを押しながらクリックすると、現在の選択範囲に原子を付け加えることができます)</span>。左の回転バー(前ページで③と書かれていたもの)でマウスボタンを押し、上下にドラッグしてください。メチル基が右左に回転するのがわかります。
549 <p><img src="../etc/tutorial_13.png" /></p>
551 同じ操作は、結合を1つ選択して左の回転バーをドラッグしてもできます。この場合、メチル基を固定してベンゼン環の方を回すこともできます。ドラッグの時にOptionキー(Mac)、Altキー(Win)を押してください。
553 <p><img src="../etc/tutorial_14.png" /></p>
555 選択した原子をドラッグすると平行移動できます。化学的に不自然な構造になりますので、注意が必要です。
557 <p><img src="../etc/tutorial_15.png" /></p>
558 <h2>3. 原子を追加する・削除する</h2>
560 原子を追加する方法はすでに一つ説明しました。選択部分または編集エリアの何もないところをダブルクリックして、構造式をタイプする方法です。実際、Molbyで原子を追加するにはこれが最も便利です。しかし、他の方法もありますので、それを紹介します。
563 先ほどのトルエン分子を使います。これをインダンに変えたいとします。
565 <p><img src="../etc/indane.png" /> インダン</p>
567 分子を回転させて、化学構造式と同じ向きになるようにします。C3から始めて、炭素原子を2つ追加し、環を閉じることにします。
569 <p><img src="../etc/tutorial_16.png" /></p>
571 編集エリアの上の"Erase"ボタンを押します。"H3"水素をクリックします。この水素原子と、C3-H3の結合が消えます。
573 <p><img src="../etc/tutorial_17.png" /></p>
575 "Bond"ボタンを押します。C3原子から右下にドラッグし、マウスボタンを離します。新しい原子と、その原子とC3の間の結合が新しく作られます。新しい原子から右上にドラッグすると、もう一つの原子と結合が作られます。
577 <p><img src="../etc/tutorial_18.png" /> <img src="../etc/tutorial_19.png" /></p>
579 "Erase"ボタンをもう一度押して、メチル基の原子を一つ消します。"Select"ボタンを押して、ベンゼン環とメチレン炭素(今水素原子を一つ消したところ)の間の結合を選択します。左の回転バーを使って、五員環に適切な向きになるようにメチレン基を回転させます。
581 <p><img src="../etc/tutorial_20.png" /> <img src="../etc/tutorial_21.png" /></p>
583 環を閉じます。"Bond"ボタンを押し、メチレン炭素原子から右下の炭素原子へドラッグします。マウスカーソルが目標の炭素原子に近づくと、ドラッグしてできた結合が炭素原子にくっつきます。ここでマウスボタンを離します。
585 <p><img src="../etc/tutorial_22.png" /></p>
587 最後に、新しく作った二つの炭素原子に水素原子を付加します。"Select"ボタンを押して、マウスを使って二つの炭素原子を選択します。"Edit"メニューから、"Add Hydrogen"->"Tetrahedral sp3"コマンドを実行します。
589 <p><img src="../etc/tutorial_23.png" /></p>
593 <p><img src="../etc/tutorial_24.png" /></p>
595 トルエンからインダンは次のような操作でも作ることができます。メチル基のオルト位の水素原子を選択、ダブルクリック、ダイアログが出たら"CH2CH3"と入力、二つのメチル基からそれぞれ水素原子を一つ削除、結合回転で向きを調整、結合を作成。この方法の方が、新しく作ったメチレン基の結合距離・結合角が正しくなるため、よりよい結果になります。上の例は、編集操作の説明のためのものと考えてください。
598 <link id="#navigation" />
601 <div file="cut_copy_paste.html">
602 <link rel="Up" href="tutorials.html" />
603 <link rel="Prev" href="editing.html" />
604 <link id="#header" />
605 <div class="contents" lang="en">
606 <h1>Step Three: Edit a Molecule: Cut/Copy/Paste</h1>
608 Like any other decent "editing" applications, Molby has a capability of cut/copy/paste by use of a clipboard. These functions (especially copy and paste) are also quite useful for building complex molecules.
611 We start from the toluene molecule again. Select the methyl group and do copy.
613 <p><img src="../etc/tutorial_25.png" /></p>
615 Unselect the methyl group, and do paste. Another methyl group appears, with no connection to the existing atoms. This is how "paste" works when nothing is selected beforehand.
617 <p><img src="../etc/tutorial_26.png" /></p>
619 Now we will see what happens when something is selected before pasting. As described above, start with the toluene molecule, select the methyl group, and copy it. Unselect the methyl group, and this time select the hydrogen atom that is "ortho" to the methyl group.
621 <p><img src="../etc/tutorial_27.png" /></p>
623 Do paste. The selected hydrogen atom is replaced with the pasted methyl group. Note that a new bond is created between the pasted methyl carbon and the ortho carbon, with an acceptable bond length and angles.
625 <p><img src="../etc/tutorial_28.png" /></p>
627 The dihedral angle may not be acceptable. In that case, you can rotate the pasted fragment by use of the left rotation bar until the dihedral angle looks good.
630 This "select and paste" technique is very useful for building large molecules. Suppose we want to build an oligobenzamide.
632 <p><img src="../etc/tutorial_39.png" /></p>
634 Create a "monomer", i.e. N-methyl-4-acetamidobenzamide. This is done by (1) create a benzene, (2) select H1, double-click, and enter "CONHCH3", (3) select H4, double-click, and enter "NHCOCH3". After entering each formula, you need to rotate the fragment to make the dihedral angle appropriate.
636 <p><img src="../etc/tutorial_29.png" /></p>
638 Select the whole molecule except for the COCH3 group at the left. Do copy.
640 <p><img src="../etc/tutorial_30.png" /></p>
642 Select the NHCH3 group at the right. Do paste.
644 <p><img src="../etc/tutorial_31.png" /></p>
646 The newly created amide bond has a bad dihedral angle. Rotate the pasted fragment by the left rotation bar, so that the amide bond becomes properly trans.
648 <p><img src="../etc/tutorial_32.png" /><img src="../etc/tutorial_33.png" /></p>
650 Select "Fit to Screen" from the "Show" menu, to make the whole molecule visible.
652 <p><img src="../etc/tutorial_34.png" /></p>
653 <p><img src="../etc/tutorial_35.png" /></p>
655 Repeat these procedures to make the tetramer.
657 <p><img src="../etc/tutorial_36.png" /></p>
658 <p><img src="../etc/tutorial_37.png" /></p>
659 <p><img src="../etc/tutorial_38.png" /></p>
661 Make one more iteration to make the octamer.
663 <p><img src="../etc/tutorial_39.png" /></p>
665 <div class="contents" lang="ja">
666 <h1>第三段階:分子を編集する:カット・コピー・ペースト</h1>
668 他のアプリケーションと同様に、Molby もクリップボードを使ってカット・コピー・ペーストを行う機能があります。これらの機能、特にコピー・ペーストは複雑な分子を作るときにたいへん有用です。
671 ふたたびトルエン分子から始めます。メチル基を選択して、コピーします。
673 <p><img src="../etc/tutorial_25.png" /></p>
675 メチル基の選択を解除して、ペーストを実行します。他の原子と結合していないメチル基がもう一つ現れます。原子が選択されていない状態で「ペースト」すると、このような動作になります。
677 <p><img src="../etc/tutorial_26.png" /></p>
679 今度は、ペーストする前に原子が選択されているときにどうなるかを見てみます。上と同じように、トルエンから始めて、メチル基を選択して、コピーします。メチル基の選択を解除して、今度はメチル基のオルト位の水素原子を選択します。
681 <p><img src="../etc/tutorial_27.png" /></p>
683 ペーストを実行します。選択された水素原子がメチル基に置き換わります。オルト位の炭素とメチル基の間に新しい結合ができ、その長さと結合角が適切なものになっていることがわかります。
685 <p><img src="../etc/tutorial_28.png" /></p>
687 二面角は適切とはいえないかも知れません。その場合は、ペーストされたグループを左の回転バーを使って回転させて、適切な二面角にします。
690 この「選択してペースト」のテクニックは、大きな分子を作るときにたいへん役に立ちます。たとえば、オリゴベンズアミドを作りたいとします。
692 <p><img src="../etc/tutorial_39.png" /></p>
694 まず、「モノマー」となる N-メチル-4-アセトアミドベンズアミドを作ります。次のようにします:(1) ベンゼンを作る、(2) H1 を選択してダブルクリックし、「CONHCH3」とタイプする、(3) H4 を選択してダブルクリックし、「NHCOCH3」とタイプする。化学式をタイプしたあとは、作成したフラグメントを回転させて二面角を調整しておきます。
696 <p><img src="../etc/tutorial_29.png" /></p>
698 左端の COCH3 以外全部を選択し、コピーします。
700 <p><img src="../etc/tutorial_30.png" /></p>
702 右端の NHCH3 を選択し、「ペースト」を実行します。
704 <p><img src="../etc/tutorial_31.png" /></p>
706 新しく作られたアミド結合は二面角が不適切です。左の回転バーを使ってペーストしたフラグメントを回転し、アミド結合が正しく trans になるようにします。
708 <p><img src="../etc/tutorial_32.png" /><img src="../etc/tutorial_33.png" /></p>
710 "Show" メニューから "Fit to Screen" を選び、分子全体が見えるようにします。
712 <p><img src="../etc/tutorial_34.png" /></p>
713 <p><img src="../etc/tutorial_35.png" /></p>
717 <p><img src="../etc/tutorial_36.png" /></p>
718 <p><img src="../etc/tutorial_37.png" /></p>
719 <p><img src="../etc/tutorial_38.png" /></p>
723 <p><img src="../etc/tutorial_39.png" /></p>
725 <link id="#navigation" />
728 <div file="ring_fusion.html">
729 <link rel="Up" href="tutorials.html" />
730 <link rel="Prev" href="cut_copy_paste.html" />
731 <link id="#header" />
732 <div class="contents" lang="en">
733 <h1>Step Four: Ring Fusion</h1>
734 <p>(Version 0.6.1 and later)</p>
735 <h2>1. Ring fusion by double-click and type-in</h2>
737 Molby has a capability (from version 0.6.1 and upper) to create a fused ring structure. As an example, here is presented yet another method to make indane. This time, we start from a benzene.
739 <p><img src="../etc/tutorial_4_01.png" /></p>
741 Choose "select" mode, and select C1 and C2 by clicking on the C1-C2 bond.
743 <p><img src="../etc/tutorial_4_02.png" /></p>
745 Double-click on the selection, and enter "cyclopentane" in the dialog box.
747 <p><img src="../etc/tutorial_4_03.png" /></p>
749 After pressing "OK", you will see a five-membered ring fused to the benzene ring.
751 <p><img src="../etc/tutorial_4_04.png" /></p>
753 Ring fusion also works when three or more atoms are selected in the original structure. For example, select the consecutive three carbons in indane as follows.
755 <p><img src="../etc/tutorial_4_05.png" /></p>
757 Double-click on the selection, and enter "C6H6" (or "benzene").
759 <p><img src="../etc/tutorial_4_06.png" /></p>
761 Now we have acenaphthene...almost.
763 <p><img src="../etc/tutorial_4_07.png" /></p>
765 As shown in the above figure, there is an extra, leftover hydrogen atom on one bridgehead carbon. After all, Molby is not so smart --- it just removes one hydrogen from each terminal carbon atoms in the selection, and connect a portion of the fusing fragment ("benzene" in this case) so that the newly created ring has the same number of atoms. Consequently, the stereochemistry may become strange when sp3 carbons are present either in the original structure or in the fusing fragment. You may need to remove/append hydrogen atoms and clean up structures (see <a href="mm_minimize.html">"Energy Minimization"</a> for detail).
767 <h2>2. Ring fusion by copy-and-paste</h2>
769 There is another way to make a ring fused structure, which includes copy and paste. Starting again from a benzene.
771 <p><img src="../etc/tutorial_4_01.png" /></p>
773 Open the "File" menu, and select "Open Predefined" → "Alicyclic" → "cyclopentane".
775 <p class="note">(This is for Version 0.6.1 and later)</p>
776 <p><img src="../etc/tutorial_4_08.png" /></p>
778 A new window named "*cyclopentane*" opens with one cyclopentane molecule. Select a portion containing three CH2 groups, and copy it by command-C or ctrl-C.
780 <p><img src="../etc/tutorial_4_09.png" /></p>
782 Return to the benzene molecule again, and this time select H1 and H2 (not C1 and C2).
784 <p><img src="../etc/tutorial_4_10.png" /></p>
786 Do paste by command-V or ctrl-V. You now get a indane molecule.
788 <p><img src="../etc/tutorial_4_11.png" /></p>
790 This "copy-and-paste" actually works in a similar way to previously described example (see <a href="cut-copy-paste.html">"Cut/Copy/Paste"</a>). When both the current selection and the fragment in the pasteboard have two terminals, they are connected and ring fusion takes place. This method of ring fusion is slightly more complicated than the "double-click and type-in" method, but it may be easier to understand.
793 <div class="contents" lang="ja">
794 <h1>第四段階:縮環構造の作成</h1>
795 <p>(バージョン 0.6.1 以上)</p>
796 <h2>1. ダブルクリック/キー入力による縮環</h2>
798 Molby は縮環構造を作るための機能を持っています(バージョン 0.6.1 以降)。例として、インダンを作るもう一つの方法を紹介します。今回は、ベンゼンから始めます。
800 <p><img src="../etc/tutorial_4_01.png" /></p>
802 "Select" モードにして、C1-C2 結合をクリックして C1 と C2 を選択します。
804 <p><img src="../etc/tutorial_4_02.png" /></p>
806 選択部分をダブルクリックして、ダイアログボックスに "cyclopentane" と入力します。
808 <p><img src="../etc/tutorial_4_03.png" /></p>
810 "OK" を押すと、5員環がベンゼン環に縮環します。
812 <p><img src="../etc/tutorial_4_04.png" /></p>
814 選択部分が3個以上の原子を含んでいても、縮環は可能です。例えば、インダンで図のように3つの連続した炭素原子を選択してみます。
816 <p><img src="../etc/tutorial_4_05.png" /></p>
818 選択部分をダブルクリックして、"C6H6"(または "benzene")と入力します。
820 <p><img src="../etc/tutorial_4_06.png" /></p>
822 そうすると、アセナフテンができます...だいたい。
824 <p><img src="../etc/tutorial_4_07.png" /></p>
826 上の図に示したように、橋頭位の炭素に1つ余分な水素原子が残ってしまいます。実のところ、Molby はあんまり賢くありません。選択部分の両端の原子から水素原子を1つずつ取り除いて、縮環する分子(この場合はベンゼン)の一部をかわりに結合して、生成する環の原子数が一致するようにしているだけです。このため、元の選択部分または縮環する分子が sp3 炭素を含んでいると、立体化学が変になることがあります。水素原子を除いたり付け加えたりして、さらに構造を最適化する必要があるでしょう(詳しくは「<a href="mm_minimize.html">エネルギー最小化</a>」を見てください)。
828 <h2>2. コピー/ペーストによる縮環</h2>
830 縮環構造を作るには、コピー/ペーストを使ったもう一つの方法があります。またベンゼンから始めましょう。
832 <p><img src="../etc/tutorial_4_01.png" /></p>
834 "File" メニューを開き、"Open Predefined" → "Alicyclic" → "cyclopentane" を選んでください。
836 <p class="note">(これはバージョン 0.6.1 以降の機能です。よく使われる分子構造がこのメニューから選択できます。)</p>
837 <p><img src="../etc/tutorial_4_08.png" /></p>
839 "*cyclopentane*"という名前の新しいウィンドウが開き、シクロペンタン分子が入っています。3つのメチレン基を選択して、コマンド-C またはコントロール-C でコピーしてください。
841 <p><img src="../etc/tutorial_4_09.png" /></p>
843 ベンゼン分子に戻って、今度は H1 と H2(C1 と C2 ではなく)を選択してください。
845 <p><img src="../etc/tutorial_4_10.png" /></p>
847 コマンド-V またはコントロール-V でペーストしてください。インダンができます。
849 <p><img src="../etc/tutorial_4_11.png" /></p>
851 この「コピー/ペースト」は、実は前に説明したのと (<a href="cut-copy-paste.html">カット/コピー/ペースト</a>) 同じように機能しています。現在の選択部分とクリップボードの内容がともに「2つの末端」を持っているとき、末端同士がつなぎあわされ、縮環したように見えるのです。この方法による縮環は「ダブルクリック/入力」の方法よりも少し複雑ですが、何が起きているかはこちらの方が理解しやすいかもしれません。
854 <link id="#navigation" />
857 <div file="proptable.html">
858 <link rel="Up" href="tutorials.html" />
859 <link rel="Prev" href="ring_fusion.html" />
860 <link id="#header" />
861 <div class="contents" lang="en">
862 <h1>Step Five: Edit a Molecule: Using a Property Table</h1>
863 <h2>1. The Property Table of a Molecule</h2>
865 A molecule has many properties. It has a set of atoms and bonds. An atom has a name, an atom type, an element, a weight (which can be derived from the element), a charge, coordinates, and so on. According to the custom in the biomolecular modeling, atoms are grouped in "residues." It also has a set of parameters used in molecular mechanics. Molby also has a limited support for quantum chemical calculations. In relation to that, a molecule in Molby can retain extra information such as MO coefficients.
868 Many of these information are accessible via the <span class="italic">property table</span>. The table is on the left side of the model window.
870 <p><img src="../etc/tutorial_5_01.png" /></p>
872 The table now shows the properties of each atom, namely its name, type, element, "residue" name and index, coordinates and partial charge. The property can be edited by double-clicking on the text.
874 <p><img src="../etc/tutorial_5_02.png" /></p>
876 When return (enter) or tab key is pressed, the edited value is finalized, and editing will continue on another cell. Pressing return, shift-return, tab, or shift-tab should cause editing of the bottom, above, right, or left cell, respectively. If you want to finish editing, use option- (or alt-) return key combination.
879 Several points are worth mentioning here.
882 (1) The atom index (the leftmost column) cannot be edited. If you want to change the order of the atoms, try using a "drag-and-drop" feature, as shown in the figure below.
884 <p><img src="../etc/tutorial_5_03.png" /></p>
886 (2) The residue column shows the residue name and index separated by a period. Manual editing of this cell may lead to a surprising result. For example, if you change the residue name of atom 0 from "RES" to "XXX", all atoms having residue index "1" will also have the new residue name "XXX".
888 <p><img src="../etc/tutorial_5_04.png" /></p>
890 This behavior is based on the principle that all residue names and indices should be consistent, <span class="italic">i.e.</span> the atoms having the same residue index should have the same residue name. This is sometimes convenient, but in many cases it causes confusion.
893 A more recommended way to change the residue name and index is to use the menu command "Assign Residue...", in the "Script" menu, after selecting atoms you want to assign one residue name and index.
896 <img src="../etc/tutorial_5_05.png" />
897
898 <img src="../etc/tutorial_5_06.png" />
901 <img src="../etc/tutorial_5_07.png" />
904 <h2>2. The Bond/Angle/Dihedral/Improper Table</h2>
906 The property table can also show other information. The type of information can be selected at the popup menu.
908 <p><img src="../etc/tutorial_5_08.png" /></p>
910 The bond, angle, dihedral, and improper tables show the indices, names, and types of the constituent atoms. In addition, the quantities (bond lengths, angles, etc.) and the molecular mechanics parameters (only after MM/MD calculation is performed) are also shown.
912 <p><img src="../etc/tutorial_5_09.png" /></p>
914 The bonds, angles, etc. are not editable in the property table. At present, the only way to create/delete bonds from GUI is to use the mouse operation "Bond" or "Erase".
917 <h2>3. The Parameter Table</h2>
919 The parameter table shows the molecular mechanics parameters in one table.
921 <p><img src="../etc/tutorial_5_10.png" /></p>
923 When you build a molecule from scratch, and open the parameter table before doing any MM/MD operations, this table is likely to be empty. This is because the MM parameters are looked up only when "Molecular Dynamics..." or "Minimize..." menu command is invoked, or you use MM/MD tools like Antechamber (see Step 4 for details). After doing some of these operations, you will see a set of parameters shown in this table.
926 The parameters are grouped in several classes, namely "VDWs", "Bonds", "Angles", "Dihedrals", "Impropers", and "VDW Pairs".
929 The "VDWs" parameters are for the van der Waals nonbonding interaction, described by "eps" (the energy at the potential minimum in kcal/mol), "r" (the van der Waals radius in Å, which is half of the interatomic distance at the equilibrium state), "eps14", "r14" (eps and r for the pair of atoms separated by exactly three bonds. In many applications, the same parameters with the ordinary pair are used with a specific factor multiplied), "atomNo" (the atomic number), "weight" (the atomic weight). The last two are usually overridden by the atom parameters.
932 The "Bonds" parameters consist of "k" (the force constant in kcal/mol/Å<sup>2</sup>) and "r0" (the equilibrium bond length). The "Angles" parameters consist of "k" (the force constant in kcal/mol/radian<sup>2</sup>) and "a0" (the equilibrium bond angle). The "Dihedrals" and "Impropers" parameters consist of "k" (the force constant in kcal/mol/radian<sup>2</sup>), "period" (the periodicity), and "phi0" (the equilibrium torsion angle in degree). The "VDW Pairs" parameters are (rarely) used to describe pair-specific van der Waals interactions, and consist of "eps", "r", "eps14", and "r14".
935 There is another important property for each parameter, that is, whether the parameter is taken from the "global" source (i.e. not specific to this molecule), or it is "local" (specific to this molecule), or "undefined". This property is shown by the color of the row. The global parameters are shown in white cells, the local ones in pale yellow, and the undefined ones in red.
937 <p><img src="../etc/tutorial_5_11.png" /></p>
939 You can edit the parameters (although you need to be familiar with MM parameters and to know what you are doing!), but only if they are "local" or "undefined" parameters. The "global" parameters cannot be edited because it may be also used in other molecules, where the edited parameters may not be appropriate. If you want to modify some parameters but they are "global", then you can make them "local" by copying the parameters and the pasting.
941 <p><img src="../etc/tutorial_5_12.png" /></p>
943 Another way to create a "local" parameter is to use the "Create New Parameter" menu command in the Edit menu. You can choose the parameter type in the submenu.
945 <p><img src="../etc/tutorial_5_13.png" /></p>
947 If some parameters are not to be used for calculation, you can "cut" the parameters from the table. This feature is useful when you have many duplicated parameters having different force constants. (In actual calculation, the parameter appearing later in the table will be used.)
951 <div class="contents" lang="ja">
952 <h1>第五段階:分子を編集する:属性テーブルを使う</h1>
953 <h2>1. 分子の属性テーブル</h2>
955 分子はいろいろな属性を持っています。まず原子と結合があります。原子は名前、原子タイプ、元素、原子量(元素が決まれば決まる)、電荷、座標、など。また、生化学の習慣に従い、原子は「残基」でグループ化されています。分子はまた、分子力学のパラメータを持っています。Molby はまた量子化学計算も部分的にサポートします。それと関連して、分子に分子軌道係数などの情報を持たせることができます。
958 これらの情報の多くは「属性テーブル」で見ることができます。属性テーブルは、分子モデルのウィンドウの左半分にあります。
960 <p><img src="../etc/tutorial_5_01.png" /></p>
962 図では、テーブルは原子の属性を表示しています。名前、原子タイプ、元素、残基名と残基番号、座標、部分電荷です。これらの属性はダブルクリックして編集することができます。
964 <p><img src="../etc/tutorial_5_02.png" /></p>
966 リターン(エンター)またはタブキーを押すと、編集した値は確定され、別のセルが編集状態になります。リターン、シフト+リターン、タブ、シフト+タブでそれぞれ下・上・右・左のセルに移動します。編集を終了したい時は、オプション(Windows では Alt)+リターンを使ってください。
972 (1) 原子のインデックス(一番左の列)は編集できません。原子の順序を変えたい時は、下の図にあるように「ドラッグ・アンド・ドロップ」操作を使ってください。
974 <p><img src="../etc/tutorial_5_03.png" /></p>
976 (2) 残基の列には、残基の名前と残基番号がピリオドでつないで表示されています。このセルを編集すると、妙なことが起きて驚くかもしれません。たとえば、原子0の残基名を "RES" から "XXX" に変えると、残基番号が "1" であるすべての原子の残基名が "XXX" に変わります。
978 <p><img src="../etc/tutorial_5_04.png" /></p>
980 この挙動は、すべての残基名と残基番号は統一されていなければならない、という原理によります。つまり、同じ残基番号を持つ原子の残基名は同じでなければなりません。時にはこの挙動が便利に使えますが、混乱を招くことも多くあります。
983 残基名と残基番号を変更するお勧めの方法は、変更したい原子を選択して "Script" メニューの "Assign Residue..." コマンドを使うことです。
986 <img src="../etc/tutorial_5_05.png" />
987
988 <img src="../etc/tutorial_5_06.png" />
991 <img src="../etc/tutorial_5_07.png" />
994 <h2>2. Bond/Angle/Dihedral/Improper テーブル</h2>
996 属性テーブルは他の情報も表示することができます。テーブルの種類はポップアップメニューで選択できます。
998 <p><img src="../etc/tutorial_5_08.png" /></p>
1000 Bond, angle, dihedral, improper テーブルは、それぞれ結合・結合角・二面角・improper 二面角の構成原子のインデックス、名前、原子タイプを表示します。また、現在の値(結合長、角度、など)と、分子力学パラメータ(MM/MD 計算を行った直後に限られますが)も表示されます。
1002 <p><img src="../etc/tutorial_5_09.png" /></p>
1004 このテーブルで、結合、結合角などを編集することはできません。現在、画面上のインターフェイスを使って結合を作ったり削除したりする唯一の方法は、"Bond", "Erase" のマウス操作を使うことです。
1007 <h2>3. パラメータテーブル</h2>
1009 パラメータテーブルは、分子力学パラメータを表示します。
1011 <p><img src="../etc/tutorial_5_10.png" /></p>
1013 分子を新しく作成して、MM/MD に関連する操作を何もせずにパラメータテーブルを開くと、テーブルはおそらく空のままでしょう。分子力学パラメータは、"Molecular Dynamics..." または "Minimize..." メニューコマンドを実行したときか、Antechamber のような分子力学用ツール(第四段階参照)を使った時に初めて作成されます。これらの操作をした後には、パラメータがこのテーブルに表示されます。
1016 パラメータはいくつかの種類にグループ分けされています。"VDWs", "Bonds", "Angles", "Dihedrals", "Impropers", そして "VDW Pairs" です。
1019 "VDWs" パラメータは、van der Waals 相互作用のためのもので、必要な値は "eps"(ポテンシャル極小値、kcal/mol)、"r"(van der Waals 半径、Å。ポテンシャル極小状態での原子間距離の 1/2)、"eps14", "r14"(結合3つで隔てられた2つの原子間の eps と r。多くの分子力学パッケージでは、通常の原子間と同じパラメータに一定の係数をかけて用いることが多い)、"atomNo"(原子番号)、"weight"(原子量)。最後の2つは、ふつうは各原子のパラメータで上書きされます。
1022 "Bonds" パラメータは、"k"(結合の力の定数、kcal/mol/Å<sup>2</sup>)、"r0"(平衡結合長さ、Å)から成ります。"Angles" パラメータは、"k"(結合角の力の定数、kcal/mol/radian<sup>2</sup>)、"a0"(平衡結合角、degree)から成ります。'Dihedrals", "Impropers" パラメータは、"k"(力の定数、kcal/mol/radian<sup>2</sup>)、"period"(周期)、"phi0"(平衡二面角、degree)から成ります。"VDW Pairs" パラメータは滅多に使われませんが、ある原子の組み合わせに特有の van der Waals 相互作用を記述するもので、"eps", "r", "eps14", "r14" から成ります。
1025 それぞれのパラメータにはもう一つ重要な特性があります。それは、そのパラメータが「グローバル」なのか(つまり、他の分子にも使われる共通の値なのか)、または「ローカル」なのか(この分子に特有)、または「未定義」なのか、です。この特性は、テーブルの行の色で区別できます。「グローバル」なパラメータは白、「ローカル」はうすい黄色、「未定義」は赤色のセルです。
1027 <p><img src="../etc/tutorial_5_11.png" /></p>
1029 パラメータは編集することができます(あなたが分子力学パラメータに詳しくて、何をしているかわかっていればの話だけど!)。ただし、パラメータが「ローカル」「未定義」の場合だけです。「グローバル」パラメータは、他の分子でも使われている可能性があるため、編集できません。「グローバル」なパラメータを編集したい時は、そのパラメータをコピー・ペーストすることで「ローカル」パラメータにすることができます。
1031 <p><img src="../etc/tutorial_5_12.png" /></p>
1033 「ローカル」なパラメータを作るもう1つの方法は、"Edit" メニューの "Create New Parameter" コマンドを使うことです。パラメータの種類をサブメニューから選ぶことができます。
1035 <p><img src="../etc/tutorial_5_13.png" /></p>
1037 もしあるパラメータは計算に使わないという場合には、そのパラメータをテーブルから削除することができます。この機能は、同じパラメータで異なる値を持つものが複数ある場合に有用です。(実際の計算では、テーブルの後ろの方にあるパラメータが優先されます。)
1040 <link id="#navigation" />
1043 <div file="mm_minimize.html">
1044 <link rel="Up" href="tutorials.html" />
1045 <link rel="Prev" href="proptable.html" />
1046 <link id="#header" />
1047 <div class="contents" lang="en">
1048 <h1>Step Six: Energy Minimization by Molecular Mechanics</h1>
1049 <h2>1. About Molecular Mechanics Implementation in Molby</h2>
1051 The molecular models built by hand generally include unnatural bond lengths, bond angles, van der Waals contacts, and so on. Molecular mechanics is a useful technique to remove such unnatural structures.
1054 Molby implements molecular mechanics calculation by use of basic force fields, including bond stretching, angle bending, dihedral rotation, improper torsions, van der Waals interaction, and electrostatic interaction. With the exception of the electrostatic interaction, the parameters are taken from the predefined table by looking up the "atom types" assigned to the atoms.
1057 By default, Molby uses <span class="italic">gaff</span> parameter set. Gaff, "General AMBER Force Field", is a part of <a href="http://ambermd.org/">AMBER molecular dynamics package</a>; the force field parameters are made public domain by courtesy of the AMBER developers. Molby also includes <span class="italic">parm99</span> parameter set, another AMBER force field that is widely used for biomolecules.
1060 The troublesome part in using force fields is how to assign atom types. The AMBER developers provides a useful set of programs to take care of this problem. Thanks again to the AMBER developers, these programs, named AmberTools, are now made public (open source under GPLv2). Molby includes some of these programs, and calls and imports the results.
1063 Following are the original papers published by the AMBER team. Please be sure to refer to them in your publication when you use the AMBER parameters or AmberTools.
1066 <li>Wang, J.; Wolf, R.M.; Caldwell, J.W.; Kollamn, P.A.; Case, D.A. Development and testing of a general Amber force field. <span class="italic">J. Comput. Chem.,</span> <span class="bold">2004,</span> <span class="italic">25,</span> 1157–1174.</li>
1067 <li>Wang, B.; Merz, K.M. Jr. A fast QM/MM (quantum mechanical/molecular mechanical) approach to calculate nuclear magnetic resonance chemical shifts for macromolecules. <span class="italic">J. Chem. Theory Comput.,</span> <span class="bold">2006,</span> <span class="italic">2,</span> 209–215. </li>
1069 <h2 id="minimize" name="minimize">2. Energy Minimization How-to</h2>
1071 Now we try energy minimization. We use 2,2'-dimethoxybiphenyl as an example.
1073 <p><img src="../etc/biphenyl.png" /> 2,2'-dimethoxybiphenyl</p>
1075 Build this molecule. The easiest way is, (1) double-click the empty editing area and type "C6H5C6H5", (2) select one ortho hydrogen, double-click, and type "OCH3", (3) repeat (2) for another ortho hydrogen on the other ring.
1077 <p><img src="../etc/tutorial_40.png" /></p>
1079 Open the "MM/MD" menu, and select "antechamber/parmchk..." command in the "Tools" submenu.
1081 <p><img src="../etc/tutorial_41.png" /></p>
1083 A dialog like below shows up. Turn off the first two checkboxes. The "log" directory is used by AmberTools for storing intermediate files; the default value would be acceptable, but you can change it here.
1085 <p><img src="../etc/tutorial_42.png" /></p>
1087 After pressing "OK", two dialog boxes appear in turn. They disappear so quickly that you may not recognize what they are saying; actually, the first one says "Running antechamber" and the second one "Running parmchk." These are programs included in AmberTools. In the present case, both programs complete successfully, and the following dialog appears.
1089 <p><img src="../etc/tutorial_43.png" /></p>
1091 Press "OK", and you will return to the molecule window. Do you notice what change has been made? It is the atom types that are modified. Specifically, the types of the atoms 0 and 10 are changed from "ca" to "cp".
1093 <p><img src="../etc/tutorial_44.png" /></p>
1095 If you are wondering what "ca" or "cp" mean, look at the global parameter table (MM/MD → View Global Parameters...), and find comment(s) in the "vdw" record.
1097 <p><img src="../etc/tutorial_45.png" /></p>
1099 Return to the molecule, and select MM/MD → Minimize.
1101 <p><img src="../etc/tutorial_46.png" /></p>
1103 A setting dialog opens. "Steps per frame" means the screen is updated every this number of steps. "Number of frames" means the maximum number of "frames" (i.e. screen updates) to calculate. If the minimization completes before this number of frames, the calculation will stop. The numbers 10 and 200 are reasonable choice in many cases.
1105 <p><img src="../etc/tutorial_47.png" /></p>
1107 Press "OK", and minimization starts. As you expect, the dihedral angle between the two phenyl rings becomes large. The calculation will stop after 200 frames. You can see the number "199" at the right bottom of the window, and the slider at the bottom of the window is now active. Move the slider, and you can see how the molecular structure changed during the minimization.
1109 <p><img src="../etc/tutorial_48.png" /></p>
1111 If you save this molecule at this stage, all the frames will be also saved (when you select the "mbsf" format), and the resulting file may be very large. If you do not want this, then you can remove all the frames by use of "Delete Frames..." command in the "Script" menu.
1113 <p><img src="../etc/tutorial_49.png" /></p>
1114 <p><img src="../etc/tutorial_50.png" /></p>
1115 <h2 id="electrostatic" name="electrostatic">3. Handling electrostatic interaction</h2>
1117 The above description is sufficient for initial cleanup of the molecular structure. However, we should go one further step to take electrostatic interaction into consideration. This is particularly important in molecules with polar functional groups (such as carbonyl).
1120 Continue our study on 2,2'-dimethoxybiphenyl. Open "MM/MD" → "Tools" → "Antechamber/parmchk...", and this time turn on the top checkbox. Also make sure the net molecular charge is correct.
1122 <p><img src="../etc/tutorial_51.png" /></p>
1124 Press "OK", and calculation starts. This time the calculation should take much longer than before, because semi-empirical calculation is carried out for optimizing the structure and getting the partial charges.
1126 <p><img src="../etc/tutorial_52.png" /></p>
1128 When calculation is done, the molecular structure may change, because structure optimization has been done by semi-empirical calculation. However, even more important is the "charge" values. You can see the charge values by scrolling the table to the right. By use of these "charge" values, interaction energies of the polar functional groups can be taken into account.
1130 <p><img src="../etc/tutorial_53.png" /></p>
1132 Note that the atomic charges given in the above method are derived from semi-empirical quantum chemical calculations. On the other hand, it is generally considered that the charges derived from ab initio calculations are better. Molby does not have capability to perform ab initio calculations, but it can help creating necessary input files for external quantum chemical programs. This will be described <a href="qchem.html#gamess_resp">elsewhere</a> in this User's Manual.
1136 <div class="contents" lang="ja">
1137 <h1>第六段階:分子力学計算によるエネルギー最小化</h1>
1138 <h2>1. Molby の分子力学計算について</h2>
1140 手作業で作成した分子モデルは、不自然な結合長、結合角、van der Waals 接触などを含んでいます。このような不自然な構造を修正するには、分子力学のテクニックが便利です。
1143 Molby の分子力学計算は、基本的な分子力場(結合の伸縮、結合角のたわみ、二面角の回転、improper torsion、van der Waals 相互作用、静電相互作用)に基づいています。静電相互作用以外については、パラメータのセットが用意してあり、原子に与えられた「原子タイプ」から必要なパラメータを取り出して計算に使います。
1146 デフォルトでは、Molby は <span class="italic">gaff</span> パラメータセットを使います。Gaff は "General AMBER Force Field" の略で、<a href="http://ambermd.org/">AMBER 分子動力学パッケージ</a>の一部です。AMBER 開発者の好意により、分子力場パラメータはパブリックドメインとされています。なお、Molby は <span class="italic">parm99</span> パラメータセットも同梱しています。これも AMBER の分子力場であり、生体分子の計算に広く使われています。
1149 分子力場を使うときに厄介なのは、原子タイプをどうやって指定するかです。AMBER には、このための便利なプログラム AmberTools が含まれています。これもまた AMBER 開発者の好意により、AmberTools はフリーソフトウェア(GPLv2)として公開されています。Molby は AmberTools のプログラムの一部を同梱しています。
1152 以下の論文は AMBER 開発チームによるものです。AMBER パラメータや AmberTools プログラムの機能を使った場合は、発表論文にこれらの論文を引用してください。
1155 <li>Wang, J.; Wolf, R.M.; Caldwell, J.W.; Kollamn, P.A.; Case, D.A. Development and testing of a general Amber force field. <span class="italic">J. Comput. Chem.,</span> <span class="bold">2004,</span> <span class="italic">25,</span> 1157–1174.</li>
1156 <li>Wang, B.; Merz, K.M. Jr. A fast QM/MM (quantum mechanical/molecular mechanical) approach to calculate nuclear magnetic resonance chemical shifts for macromolecules. <span class="italic">J. Chem. Theory Comput.,</span> <span class="bold">2006,</span> <span class="italic">2,</span> 209–215. </li>
1158 <h2 id="minimize" name="minimize">2. エネルギー最小化の方法</h2>
1160 さて、エネルギー最小化をやってみましょう。例として 2,2'-ジメトキシビフェニルを使うことにします。
1162 <p><img src="../etc/biphenyl.png" /> 2,2'-ジメトキシビフェニル</p>
1164 この分子のモデルを作ってください。簡単なやり方は、(1) 空の編集エリアでダブルクリックして "C6H5C6H5" とタイプする。(2) オルト位の水素原子を1つ選択し、ダブルクリックして "OCH3" とタイプする。(3) もう1つの環のオルト位の水素原子について、(2) を繰り返す。
1166 <p><img src="../etc/tutorial_40.png" /></p>
1168 "MM/MD"メニューを開き、"Tools" サブメニューから "Antechamber/parmchk..." コマンドを実行します。
1170 <p><img src="../etc/tutorial_41.png" /></p>
1172 下のようなダイアログが現れます。上の2つのチェックボックスをオフにしてください。"Log" ディレクトリは、AmberTools のプログラムが中間ファイルを保存するのに使います。デフォルトの位置で問題はないでしょうが、変更してもかまいません。
1174 <p><img src="../etc/tutorial_42.png" /></p>
1176 "OK"を押すと、2つのダイアログが順に現れます。すぐに消えてしまうので、何と書いてあるか読めないかも知れません。最初のものは "Running antechamber", 2つめのものは "Running parmchk" と書いてあります。これらは AmberTools に含まれているプログラムです。今の場合は、両方とも正しく実行され、次のダイアログが現れます。
1178 <p><img src="../etc/tutorial_43.png" /></p>
1180 "OK"を押して、分子のウィンドウに戻ってください。何が変わったかわかるでしょうか? 変わったのは原子タイプです。具体的には、原子0と原子10のタイプが "ca" から "cp" に変わっています。
1182 <p><img src="../etc/tutorial_44.png" /></p>
1184 "ca" や "cp" が何の意味か知りたければ、"MM/MD" メニューから "View Global Parameters..." を選んでみてください。Molby に同梱されているパラメータの一覧表が出てきます。"vdw" (van der Waals パラメータの意味ですが、原子タイプの定義も兼ねています)の中で "ca", "cp" を見つけて、表の右端にあるコメントを見てみてください。
1186 <p><img src="../etc/tutorial_45.png" /></p>
1188 分子のウィンドウに戻り、"MM/MD" → "Minimize" を選んでください。
1190 <p><img src="../etc/tutorial_46.png" /></p>
1192 設定ダイアログが開きます。"Steps per frame" は、構造最適化計算がこのステップ数進んだところで画面を更新する、という意味です。"Number of frames" は、この数だけ画面を更新したら計算を止める、という意味です。その前に構造最適化が終了すれば、その時点で計算は止まります。下の図にある 10, 200 を入れておけば、たいていの場合はよいでしょう。
1194 <p><img src="../etc/tutorial_47.png" /></p>
1196 "OK" を押すと、構造最適化が始まります。この分子の場合、期待される通り、2つのベンゼン環の間の二面角がだんだん大きくなります。指定通り、200回画面が更新されると計算は止まります。ウィンドウの右下に "199" という数字があり、ウィンドウの下辺にあるスライダーが有効になっています。更新された画面は、それぞれ独立したフレームとして保存されており、このスライダーを動かすと、初期状態からどのように構造が変わって行くかを再現することができます。
1198 <p><img src="../etc/tutorial_48.png" /></p>
1200 この時点で分子を "mbsf" フォーマットで保存すると、全部のフレームの情報が保存され、ファイルがとても大きくなります。これが好ましくない場合は、現在表示されている以外のフレームを削除することもできます。"Script" メニューから "Delete Frames..." コマンドを使ってください。
1202 <p><img src="../etc/tutorial_49.png" /></p>
1203 <p><img src="../etc/tutorial_50.png" /></p>
1204 <h2 id="electrostatic" name="electrostatic">3. 静電相互作用の取り扱い</h2>
1206 上の説明は、分子構造のひずみを取り除く最初の段階としては十分です。しかしながら、次の段階として、静電相互作用を考慮しなければなりません。特に、カルボニル基などの極性官能基が含まれている場合は、これは重要です。
1209 2,2'-ジメトキシビフェニルについて検討を続けましょう。"MM/MD" → "Tools" → "Antechamber/parmchk...", を開き、今度は一番上のチェックボックスをオンにします。そして、分子の総電荷が正しいことを確認してください。この場合は中性分子なので、"0" で結構です。
1211 <p><img src="../etc/tutorial_51.png" /></p>
1213 "OK" を押すと、計算が始まります。今度は、前回よりもずっと長く計算時間がかかります。半経験的分子軌道法により、構造最適化と部分電荷の計算を行うからです。
1215 <p><img src="../etc/tutorial_52.png" /></p>
1217 計算が終了すると、分子構造は少し変化することがあります。半経験的分子軌道計算で構造最適化が行われたからです。また、"charge" の値が設定されていることがわかります。テーブルを右にスクロールすると、与えられた電荷の値を見ることができます。この電荷の値を用いて静電相互作用の計算を行うことにより、極性官能基の相互作用エネルギーをより正しく見積もることができます。
1219 <p><img src="../etc/tutorial_53.png" /></p>
1221 上の方法で計算した電荷は、半経験的分子軌道計算によるものです。一方、一般的には ab initio 計算で求めた電荷の方が良いと考えられています。Molby は ab initio 計算を行う機能を持っていませんが、外部の量子化学計算プログラムへの入力ファイルを作成する機能があります。これは、このマニュアルの<a href="qchem.html#gamess_resp">別のところ</a>で解説します。
1225 <link id="#navigation" />
1228 <div file="md.html">
1229 <link rel="Up" href="tutorials.html" />
1230 <link rel="Prev" href="mm_minimize.html" />
1231 <link id="#header" />
1232 <div class="contents" lang="en">
1233 <h1>Step Seven: Molecular Dynamics (MD) Calculation</h1>
1234 <h2>1. MD Calculation within Molby</h2>
1236 Molby implements molecular dynamics (MD) calculation which uses the same force fields as the energy minimization by molecular mechanics (MM). This implementation is suitable only for preliminary calculations (to check parameters quickly, etc.); for production runs, it is strongly recommended that you use one of the established software packages.
1239 An example of preliminary MD run is presented here. We use 2,2'-dimethoxybiphenyl again. In a similar way as in the <a href="mm_minimize.html">step six</a>, create a model and assign MM parameters (including the partial charges on the atoms).
1241 <p><img src="../etc/tutorial_7_01.png" /></p>
1243 Select MM/MD → Molecular Dynamics. A setting dialog opens. Although this dialog resembles that in MM minimization, it shows different parameters that are relevant to the MD calculation.
1245 <p><img src="../etc/tutorial_7_02.png" /></p>
1247 The "timestep" parameter defines the minimum time increment in solving the equation of motions. The "target temperature" defines the temperature of the system. Before running the MD calculation, all atoms are let to have random velocities with the Boltzmann distribution corresponding to the target temperature. In addition, the velocities are modified so that the temperature is kept constant. The "steps per frame" and "number of frames" parameters have the same meaning as in the MM minimization. The screenshot shows 10 as the "steps per frame"; this is an appropriate value for MM minimization, however for MD calculations a larger value would be more appropriate (like 100).
1250 Pressing the "Advanced..." button opens another dialog with other parameters. The meanings of the parameters are described in the embedded Ruby document, in the <a href="molby_rb/MDArena.html">MDArena</a> page.
1252 <p><img src="../etc/tutorial_7_03.png" /></p>
1254 Return to the original MD dialog (by pressing "Close" in the advanced settings dialog), and press "OK". The MD run starts, and new frames are accumulated. If you find something is wrong, or want to stop before getting the specified number of frames, you can stop the MD run by pressing Command-period (Mac) or Control-period (Windows).
1256 <h2>2. Using Molby with AMBER: Creating Inputs and Importing Outputs</h2>
1258 For production runs, you can create AMBER input files from Molby. More specifically, you can create "prmtop" and "inpcrd" files for the SANDER module.
1261 <span class="italic">Note:</span> There is no guarantee that Molby creates the exactly same input as the official AMBER modeling tools, nor the generated files are valid inputs for the SANDER module. You may need to modify them by hand.
1264 To creat SANDER input files, select "Create SANDER input..." command in the "MM/MD" → "Tools" submenu.
1266 <p><img src="../etc/tutorial_7_04.png" /></p>
1268 You will be asked first for the file name for the "prmtop" file. Please be sure to add ".prmtop" extension. The other file, "inpcrd" file, is created as the same name with the ".prmtop" extension replaced by the ".inpcrd" extension.
1270 <p><img src="../etc/tutorial_7_05.png" /></p>
1272 Subsequently, you will be asked to select one of the two versions of the prmtop files. The older one, "AMBER8/NAMD", allows you to use the output file by NAMD (see below).
1274 <p><img src="../etc/tutorial_7_06.png" /></p>
1276 Now you can transfer your files to the workstation running SANDER. To perform the simulation, you still need to create the instruction file for SANDER, which you should already know if you are using AMBER.
1279 After the simulation is over, you can get the trajectory file ("mdcrd" file) back and import to Molby. Use "Import..." command in the "File" menu, select "AMBER mdcrd file (*.crd; *.mdcrd)", and specify the file.
1281 <p><img src="../etc/tutorial_7_07.png" /></p>
1282 <h2>3. Using Molby with NAMD: Creating Inputs and Importing Outputs</h2>
1284 You can also use the NAMD software package for production run. NAMD is developed by the Theoretical Biophysics Group in the University of Illinois at Urbana-Champaign, and the official information is found at their web site (<a href="http://www.ks.uiuc.edu/Research/namd/">http://www.ks.uiuc.edu/Research/namd/</a>). NAMD can use the AMBER "prmtop" as the input, by use of the instruction <code>amber yes</code>. See NAMD Users' Guide for details.
1287 You can also import the NAMD output by importing the dcd file. The file format is also listed in the "Import..." file dialog.
1289 <h2>4. Building Solvated Structures</h2>
1291 When you want to perform MD simulations in explicit solvent, you need to build a box of solvent molecules around the target molecule. Molby can help creating solvated structures.
1294 To build a solvated structure, you need to open the file containing the predefined box of the desired solvent. The following solvent box is included in the Molby package, and can be accessed from "File" → "Open Predefined" → "Solvent boxes" submenu. The tip3box was taken from the AmberTool package, and other solvent boxes were taken from <a href="http://www.pharmacy.manchester.ac.uk/bryce/amber">Amber parameter database</a>.
1296 <table border="1" cellspacing="0">
1297 <tr><th>name</th><th>solvent</th><th>reference</th></tr>
1298 <tr><td>tip3pbox</td>
1300 <td>Jorgensen, W. L.; Chandrasekhar, J.; Madura, J.; Klein, M. L.<br />
1301 J. Chem. Phys. 1983, 79, 926–935.
1303 <tr><td>chcl3box</td>
1305 <td>Cieplak, P.; Caldwell, J. W.; Kollman, P. A.<br />
1306 J. Comp. Chem. 2001, 22, 1048<br />
1308 <tr><td>dmsobox</td>
1309 <td>dimethylsulfoxide</td>
1310 <td>Fox, T.; Kollman, P. A.<br />
1311 J. Phys. Chem. B 1998, 102, 8070.
1313 <tr><td>meohbox</td>
1315 <td>Caldwell, J. W.; Kollman, P. A.<br />
1316 J. Phys. Chem. 1995, 99, 6208.
1319 <td>N-methylacetamide</td>
1320 <td>Caldwell, J. W.; Kollman, P. A.<br />
1321 J. Phys. Chem. 1995, 99, 6208.
1324 <p><img src="../etc/tutorial_7_08.png" /><img src="../etc/tutorial_7_09.png" /></p>
1326 While keeping the solvent box open, create or open the target (solute) molecule in a separate window. With this solute molecule in the front window, select "Solvate..." command from the "Script" menu. A dialog box opens.
1328 <p><img src="../etc/tutorial_7_10.png" /></p>
1330 In the popup menu "Choose solvent box:", you will find the solvent box you opened earlier. Note that this popup menu lists <span class="italic">all</span> open molecules that have the associated periodic box (or unit cell). This can lead to a confusing situation that, if you have another solvated structure also open, that structure will also be listed in this popup menu, because a solvated structure always have a periodic box. Therefore, please take care so that you choose the right solvent box.</p>
1332 If you open the solvent box from the "Open Predefined" menu, you can easily recognize it in the popup menu, because the name has asterisks at the beginning and the end, like "*CHCl3*".
1335 The "Box offset" parameters define the thickness of the solvent layer surrounding the solute molecule. More specifically, the periodic box of the solvated structure is determined as follows: the minimum cuboid that can surround the solute molecule is defined, and all the faces are offset by the "Box offset" distances to the outside direction. On the other hand, it is possible to define the size of the periodic box explicitly for any of the three (x, y, z) directions. If you wish to do this, give a negative number for that direction. (For example, if you want the periodic box to be 40 Å in the x direction, give -40 as the first "Box offset" parameter.)
1338 The size of the periodic box can be checked by selecting "MM/MD" → "Define Unit Cell" menu command.
1341 The "Exclusion limit distance" defines the minimum allowed distance between the solvent and solute molecules. The solvent molecule is removed when it has atoms with smaller distances than this parameter from the solute molecule.
1345 <div class="contents" lang="ja">
1346 <h1>第七段階:分子動力学計算</h1>
1347 <h2>1. Molby 組み込みの分子動力学 (MD) 計算</h2>
1349 Molby は分子動力学 (MD) 計算を実装しています。これは分子力学によるエネルギー最小化と同じ分子力場を使います。この機能は、パラメータのチェックなどの予備的な計算に使うことを想定しています。本格的な計算には、実績のある他のソフトウェアパッケージを使うことを強くおすすめします。
1352 予備的な MD 計算の例を示します。また 2,2'-ジメトキシビフェニルを使います。<a href="mm_minimize.html">第六段階</a>と同様に、モデルを作成して分子力場パラメータを決定してください(原子上の部分電荷も)。
1354 <p><img src="../etc/tutorial_7_01.png" /></p>
1356 MM/MD メニューの "Molecular Dynamics" を選ぶと、設定ダイアログが開きます。このダイアログは分子力場によるエネルギー最小化の時と似ていますが、MD 計算に使う他のパラメータを表示しています。
1358 <p><img src="../etc/tutorial_7_02.png" /></p>
1360 "Timestep" パラメータは、運動方程式を解くときの最小時間刻みを表します。"Target temperature" は系の温度を表します。MD 計算が始まる際に、すべての原子はこの温度に対応するボルツマン分布に従ってランダムな速度を与えられます。さらに、MD 計算中は温度が一定に保たれるように速度が調整されます。"Steps per frame" と "number of frames" パラメータは分子力学によるエネルギー最小化の時と同じ意味を持ちます。このスクリーンショットでは "steps per frame" が 10 になっています。これはエネルギー最小化では適切な値ですが、分子動力学の場合はもう少し大きな値(たとえば 100)の方がより適切です。
1363 "Advanced..." ボタンを押すと、他のパラメータを持つ別のダイアログが開きます。これらのパラメータの意味は、内蔵 Ruby インタプリタのリファレンスで <a href="molby_rb/MDArena.html">MDArena</a> のページに書かれています。
1365 <p><img src="../etc/tutorial_7_03.png" /></p>
1367 元の MD 設定ダイアログに戻り(詳細設定ダイアログで "Close" ボタンを押す)、"OK" を押します。MD 計算が始まり、新しいフレームが追加されて行きます。何か問題があったり、指定したフレーム数よりも前に計算を止めたいときは、コマンド-ピリオド (Mac) かコントロール-ピリオド (Windows) を使ってください。
1369 <h2>2. AMBER とともに使う:入力の作成と出力結果のインポート</h2>
1371 本格的な計算を行うために、Molby で AMBER への入力ファイルを作成することができます。具体的には、SANDER モジュールの入力となる "prmtop" と "inpcrd" ファイルを作ることができます。
1374 <span class="italic">注:</span> Molby が AMBER 付属のモデリングツールと同じ入力ファイルを作成する保証はありませんし、SANDER モジュールへの正しい入力になっている保証もありません。手作業で修正が必要な場合もあるかもしれません。
1377 SANDER の入力ファイルを作るには、"MM/MD" → "Tools" サブメニューから "Create SANDER input..." コマンドを選んでください。
1379 <p><img src="../etc/tutorial_7_04.png" /></p>
1381 最初に "prmtop" ファイルの名前を聞かれます。必ず ".prmtop" 拡張子をつけるようにしてください。もう一つの "inpcrd" ファイルは、".prmtop" を ".inpcrd" に置き換えた名前で保存されます。
1383 <p><img src="../etc/tutorial_7_05.png" /></p>
1385 次に、prmtop ファイルの2つのバージョンのうち1つを選ぶように促されます。古い方、"AMBER8/NAMD" を使えば、NAMD ソフトウェアパッケージ(下記参照)でも使える入力ファイルができます。
1387 <p><img src="../etc/tutorial_7_06.png" /></p>
1389 これで、SANDER が走るワークステーションにファイルを転送することができます。シミュレーションを実行するためには、SANDER の命令ファイルを作らなければなりませんが、AMBER をお使いの方ならやり方はよくご存知でしょう。
1392 計算が終了したら、トラジェクトリファイル ("mdcrd" ファイル) を取得して、Molby にインポートすることができます。"File" メニューの "Import..." コマンドを選び、ファイルタイプとして "AMBER mdcrd file (*.crd; *.mdcrd)" を選んで、ファイルを読み込んでください。
1394 <p><img src="../etc/tutorial_7_07.png" /></p>
1395 <h2>3. NAMD とともに使う:入力の作成と出力結果のインポート</h2>
1397 本格的な計算を行うために NAMD を使うこともできます。NAMD はイリノイ大学アーバナ・シャンペーン校の理論生物物理グループが開発したソフトウェアパッケージです。公式ウェブサイトは <a href="http://www.ks.uiuc.edu/Research/namd/">http://www.ks.uiuc.edu/Research/namd/</a> です。NAMD は AMBER の "prmtop" を入力として使うことができます(<code>amber yes</code> 命令を使う)。詳しくは NAMD のユーザーズガイドをご覧ください。
1400 NAMD の出力の dcd ファイルもインポートが可能です。"Import..." コマンドのファイルタイプで選択することができます。
1402 <h2>4. 溶媒和構造の作成</h2>
1404 明示的な溶媒中での MD シミュレーションを行うためには、目的分子の回りに溶媒分子の箱を作る必要があります。Molby でこのような溶媒和構造を作成することができます。
1407 溶媒和構造を作るためには、まず使用する溶媒の箱を読み込む必要があります。以下の溶媒の箱が Molby に同梱されており、"File" → "Open Predefined" → "Solvent boxes" サブメニューから選ぶことができます。Tip3box は AmberTool パッケージから得たものであり、他の溶媒箱は <a href="http://www.pharmacy.manchester.ac.uk/bryce/amber">Amber parameter database</a> で公開されているものです。
1409 <table border="1" cellspacing="0">
1410 <tr><th>name</th><th>solvent</th><th>reference</th></tr>
1411 <tr><td>tip3pbox</td>
1413 <td>Jorgensen, W. L.; Chandrasekhar, J.; Madura, J.; Klein, M. L.<br />
1414 J. Chem. Phys. 1983, 79, 926–935.
1416 <tr><td>chcl3box</td>
1418 <td>Cieplak, P.; Caldwell, J. W.; Kollman, P. A.<br />
1419 J. Comp. Chem. 2001, 22, 1048<br />
1421 <tr><td>dmsobox</td>
1422 <td>dimethylsulfoxide</td>
1423 <td>Fox, T.; Kollman, P. A.<br />
1424 J. Phys. Chem. B 1998, 102, 8070.
1426 <tr><td>meohbox</td>
1428 <td>Caldwell, J. W.; Kollman, P. A.<br />
1429 J. Phys. Chem. 1995, 99, 6208.
1432 <td>N-methylacetamide</td>
1433 <td>Caldwell, J. W.; Kollman, P. A.<br />
1434 J. Phys. Chem. 1995, 99, 6208.
1437 <p><img src="../etc/tutorial_7_08.png" /><img src="../etc/tutorial_7_09.png" /></p>
1439 溶媒箱を開いたまま、目的とする分子(溶質)を新しいウィンドウで開くか作成します。溶質分子のウィンドウを最前面にした状態で、"Script" メニューから "Solvate..." コマンドを選びます。ダイアログが開きます。
1441 <p><img src="../etc/tutorial_7_10.png" /></p>
1443 ポップアップメニュー "Choose solvent box:" には、さきほど開いた溶媒箱がリストされているはずです。注意していただきたいのは、このポップアップメニューには周期境界(または単位格子)を持つすべての分子がリストされることです。このため、紛らわしいことに、別の溶媒和構造を同時に開いていると、その構造もこのポップアップメニューに登場してしまいます(溶媒和構造は必ず周期境界を持っていますから)。ですから、正しい溶媒箱を選ぶように十分に注意してください。
1446 "Open Predefined" メニューから選んで開いた溶媒箱はすぐに区別できます。名前の前後に * がついているからです ("*CHCl3*" のように)。
1449 "Box offset" パラメータは、溶質分子の回りを囲む溶媒層の厚さを指定します。言い換えると、これから作成する溶媒和構造の周期境界は次のように決められます:まず溶質を囲む最小の直方体を算出し、それぞれの面を外側に向かって Box offsect パラメータ分だけ移動させます。一方、x, y, z のそれぞれの方向に対して、周期境界の大きさを決めたい場合もあります。このときは、その方向の Box offset パラメータに負の値を指定してください。(たとえば、x 方向の周期境界を大きさを 40 Å にしたい時は、"Box offset" の最初のパラメータを -40 にしてください。)
1452 周期境界の箱の大きさは、"MM/MD" → "Define Unit Cell" メニューを選ぶと見ることができます。</p>
1454 "Exclusion limit distance" は、溶質分子と溶媒分子の原子間距離の最小値を定めます。溶媒分子のある原子が溶質分子からこの距離以内に近づいているとき、その溶媒分子は取り除かれます。
1457 <link id="#navigation" />
1460 <div file="qchem.html">
1461 <link rel="Up" href="tutorials.html" />
1462 <link rel="Prev" href="md.html" />
1463 <link id="#header" />
1464 <div class="contents" lang="en">
1465 <h1>Step Eight: Collaboration with Other Quantum Chemistry Softwares</h1>
1467 Molby has a capability to export and import files for quantum chemistry softwares, namely <a href="http://www.gaussian.com/">Gaussian</a> and <a href="http://www.msg.ameslab.gov/GAMESS/">GAMESS</a>. At present, the capability is quite limited, and in many cases it would be much better to use other established softwares. Nevertheless, if you are already familiar with Molby, you may want to use Molby for creating input files for Gaussian and GAMESS and processing outputs. Here are instructions how to do it.
1470 Needless to say, you need to have access to Gaussian or GAMESS program packages. They can be on the same machine or on other machines (workstations) than Molby. Please learn how to use these packages before using Molby as described in this section.
1472 <h2>1. Using Gaussian</h2>
1474 The Gaussian input can be created by selecting "Export..." command in the "File" menu. The file extension is either "gjf" (as is the convention in GaussianW) or "com" (as in UNIX version of Gaussian).
1476 <p><img style="vertical-align:top;" src="../etc/tutorial_8_01.png" /><img src="../etc/tutorial_8_02.png" /></p>
1478 The output will look like this. Although Molby can create only one type of Gaussian input (optimize with PM3), it should be relatively easy to modify the generated input file by hand.
1480 <p class="code">%Chk=benzene.chk
1483 unnamed1; created by Molby at Sat Feb 11 00:30:21 +0900 2012
1486 C -0.653000 0.585000 -1.068000
1487 H -1.158000 1.039000 -1.898000
1488 C 0.729000 0.607000 -1.003000
1489 H 1.295000 1.076000 -1.783000
1490 C 1.382000 0.021000 0.069000
1491 H 2.452000 0.038000 0.119000
1492 C 0.651000 -0.586000 1.076000
1493 H 1.156000 -1.039000 1.906000
1494 C -0.732000 -0.607000 1.012000
1495 H -1.298000 -1.077000 1.792000
1496 C -1.384000 -0.021000 -0.060000
1497 H -2.455000 -0.038000 -0.110000
1501 When you do geometrical optimization, you may want to examine how the structure changes as the calculation proceeds. This can be done by importing the Gaussian output file. The extention should be either ".out" or ".log".
1503 <p><img style="vertical-align:top;" src="../etc/tutorial_8_03.png" /><img src="../etc/tutorial_8_04.png" /></p>
1504 <h2>2. Using GAMESS</h2>
1506 Creating GAMESS input can be more complicated than Gaussian, so that Molby provides a simple dialog to help creating GAMESS input. The dialog is accessible from the "Creating GAMESS input..." command in the "QChem" menu.
1508 <p><img src="../etc/tutorial_8_05.png" /></p>
1510 You can specify various settings in the dialog.
1512 <p><img src="../etc/tutorial_8_06.png" /></p>
1514 <li><b>SCF Type:</b> RHF, ROHF, or UHF.</li>
1515 <li><b>Run Type:</b> Energy, Property, or Optimize.</li>
1516 <li><b>Use internal coordinates for structure optimization:</b> Add instructions for automatic generation of internal Z-matrix. (<i>Note:</i> It cannot be used for linear molecules.)</li>
1517 <li><b>Charge:</b> The (formal) charge of the molecule.</li>
1518 <li><b>Multiplicity:</b> The spin multiplicity.</li>
1519 <li><b>Use DFT:</b> Check if you want to use DFT.</li>
1520 <li><b>DFT type:</b> At present only B3LYP is available in this dialog.</li>
1521 <li><b>Basis set:</b> One of PM3, STO-3G, 3-21G, 6-31G, 6-31G(d), 6-31G(d, p), 6-311G, 6-311G(d, p), or LanL2DZ.</li>
1522 <li><b>Use secondary basis set:</b> Check if you want to use another basis set for certain elements.</li>
1523 <li><b>Elements:</b> The elements (comma separated) to use secondary basis set.</li>
1524 <li><b>Basis set:</b> The secondary basis set.</li>
1525 <li><b>Calculate electrostatic potential:</b> This is used for RESP charge calculation for AMBER.</li>
1528 When the calculation of GAMESS is complete, you will find two output files, namely *.log and *.dat. Either format can be imported by use of the "Import..." menu command. Some informations are included in both (<i>e.g.</i> coordinates during structural optimization), but other informations are only in one of these files (<i>e.g.</i> the full description of gaussian functions is only in the *.log file, whereas the orbital coefficients with full precision are only in the *.dat file). You need to be familiar with the structure of the GAMESS output files to fully utilize the GAMESS import capability of Molby.
1530 <h2 id="gamess_resp" name="gamess_resp">3. Using GAMESS for calculation of the RESP charges</h2>
1532 We already saw <a href="mm_minimize.html#electrostatic">how to assign partial charges</a> for evaluation of electrostatic interactions in MM calculations. There we used semi-empirical calculations, although <i>ab initio</i> calculations will give better results if possible. Here are instructions how to do it using GAMESS.
1535 Select "MM/MD" → "Tools" → "GAMESS/RESP...".
1537 <p><img src="../etc/tutorial_8_07.png" /></p>
1539 The following window pops up. As the first step, press the "Create GAMESS Input..." button.
1541 <p><img src="../etc/tutorial_8_08.png" /></p>
1543 The familiar GAMESS dialog opens up. It is most important to turn on the "Calculate electrostatic potential (ESP)" checkbox (it should be turned on if you follow the steps as described here, but please double-check). Also make sure that the charge and multiplicity are correct, and select a suitable basis set (6-31G(d) is recommended).
1545 <p><img src="../etc/tutorial_8_09.png" /></p>
1547 Press OK to create the GAMESS input, and send it to GAMESS. The calculation will take time, so that you can finish Molby (after saving the molecule as a msbf file!), and work on something else at this stage.
1550 After the GAMESS calculation is complete, open the same molecule, and select "MM/MD" → "Tools" → "GAMESS/RESP..." again. This time, follow the second step by pressing the "Import GAMESS dat..." button.
1552 <p><img src="../etc/tutorial_8_10.png" /></p>
1554 Select the GAMESS dat file (which should be available when the GAMESS calculation ends successfully), and import it. When the import is complete, the "Run RESP..." button should be enabled. If it does not, the imported dat file does not have the electrostatic potential information. Start over from the step 1, and make sure that the "Calculate electrostatic potential (ESP)" checkbox is on.
1556 <p><img src="../etc/tutorial_8_11.png" /></p>
1558 Press the "Run RESP..." button, and the following dialog opens. This is almost the same as the dialog for Antechamber <a href="mm_minimize.html#minimize">described before</a>.
1560 <p><img src="../etc/tutorial_8_12.png" /></p>
1562 Press the "OK" button, and the RESP charge will be assigned to the atoms.
1565 <div class="contents" lang="ja">
1566 <h1>第八段階:他の量子化学ソフトウェアとの連携</h1>
1568 Molby は量子化学計算ソフトウェア <a href="http://www.gaussian.com/">Gaussian</a>、<a href="http://www.msg.ameslab.gov/GAMESS/">GAMESS</a> のファイルを作成/読み込みする機能を持っています。現状では、Molby には極めて制限された機能しかないため、おそらく他の実績あるソフトウェアを使った方がよいでしょう。しかし、Molby の操作に慣れているなら、Gaussian や GAMESS の入力を作成して結果を読み込むのに Molby を使いたいことがあるかもしれません。どのようにすればよいかを説明します。
1571 当然のことながら、Gaussian や GAMESS のプログラムパッケージが使えることが前提です。これらが走るのは、Molby と同じマシン上でも異なるマシン(ワークステーション)上でも構いません。この項で説明するやり方で Molby を使う前に、これらのプログラムパッケージをどのように使うかを学んでおいてください。
1573 <h2>1. Gaussian を使う</h2>
1575 Gaussian の入力は "File" → "Export..." コマンドを使えば作成することができます。ファイルの拡張子は "gjf"(GaussianW での習慣)、または "com"(UNIX 版での習慣)です。
1577 <p><img style="vertical-align:top;" src="../etc/tutorial_8_01.png" /><img src="../etc/tutorial_8_02.png" /></p>
1579 できあがったファイルは次のようになります。Molby が作成できるのは、一種類の Gaussian 入力ファイル(PM3 による構造最適化)だけですが、Gaussian の入力ファイルを手作業で修正するのは比較的簡単です。
1581 <p class="code">%Chk=benzene.chk
1584 unnamed1; created by Molby at Sat Feb 11 00:30:21 +0900 2012
1587 C -0.653000 0.585000 -1.068000
1588 H -1.158000 1.039000 -1.898000
1589 C 0.729000 0.607000 -1.003000
1590 H 1.295000 1.076000 -1.783000
1591 C 1.382000 0.021000 0.069000
1592 H 2.452000 0.038000 0.119000
1593 C 0.651000 -0.586000 1.076000
1594 H 1.156000 -1.039000 1.906000
1595 C -0.732000 -0.607000 1.012000
1596 H -1.298000 -1.077000 1.792000
1597 C -1.384000 -0.021000 -0.060000
1598 H -2.455000 -0.038000 -0.110000
1602 構造最適化を行ったあとは、構造がどのように変化するか見てみたいと思うでしょう。これは、Gaussian 出力ファイルを読み込めば実現できます。拡張子は ".out" または ".log" でなくてはなりません。
1604 <p><img style="vertical-align:top;" src="../etc/tutorial_8_03.png" /><img src="../etc/tutorial_8_04.png" /></p>
1605 <h2>2. GAMESS を使う</h2>
1607 GAMESS の入力を作成するのは Gaussian よりもずっと複雑なので、専用のダイアログが用意されています。このダイアログは、"QChem" → "Creating GAMESS input..." コマンドで開くことができます。
1609 <p><img src="../etc/tutorial_8_05.png" /></p>
1611 このダイアログでは、いろいろな設定を決めることができます。
1613 <p><img src="../etc/tutorial_8_06.png" /></p>
1615 <li><b>SCF Type:</b> RHF, ROHF, または UHF.</li>
1616 <li><b>Run Type:</b> Energy, Property, または Optimize.</li>
1617 <li><b>Use internal coordinates for structure optimization:</b> 内部的に Z-matrix を自動生成するための命令を付け加えます。(<i>注:</i> 直線状分子には使えません。)</li>
1618 <li><b>Charge:</b> 分子の形式電荷。</li>
1619 <li><b>Multiplicity:</b> スピン多重度。</li>
1620 <li><b>Use DFT:</b> DFT 計算を行うときチェックします。</li>
1621 <li><b>DFT type:</b> このダイアログでは B3LYP のみ指定できます。</li>
1622 <li><b>Basis set:</b> 以下の基底が指定できます:PM3, STO-3G, 3-21G, 6-31G, 6-31G(d), 6-31G(d, p), 6-311G, 6-311G(d, p), または LanL2DZ.</li>
1623 <li><b>Use secondary basis set:</b> 特定の元素のみ別の基底を使いたいときチェックします。</li>
1624 <li><b>Elements:</b> 別の基底を使う元素(コンマで区切って複数指定できます)。</li>
1625 <li><b>Basis set:</b> 別の基底。</li>
1626 <li><b>Calculate electrostatic potential:</b> RESP 電荷を求めるための静電ポテンシャルの計算を行う。</li>
1629 GAMESS の計算が終了すると、*.log と *.dat の2つのファイルができます。どちらも "Import..." コマンドで読み込むことができます。ある種の情報(たとえば構造最適化途中の座標)はどちらのファイルにも含まれていますが、その他の情報はどちらか一方にしか含まれません(たとえば、基底を構成する Gaussian 関数の完全な係数は *.log ファイルにしかなく、精度の高い軌道係数は *.dat にしかありません)。GAMESS 読み込み機能を十分に活用するには、GAMESS の出力が何を含んでいるかをよく理解する必要があります。
1631 <h2 id="gamess_resp" name="gamess_resp">3. GAMESS を用いて RESP 電荷を計算する</h2>
1633 分子力学計算で静電相互作用を評価するための部分電荷の計算については<a href="mm_minimize.html#electrostatic">すでに説明しました</a>。そこでは半経験的分子軌道計算を用いましたが、可能ならば ab initio 計算の方がよい結果を与えます。GAMESS を使って計算する方法を説明します。
1636 メニューより "MM/MD" → "Tools" → "GAMESS/RESP..." を選びます。
1638 <p><img src="../etc/tutorial_8_07.png" /></p>
1640 次のウィンドウが開きます。第1のステップとして、"Create GAMESS Input..." ボタンを押します。
1642 <p><img src="../etc/tutorial_8_08.png" /></p>
1644 GAMESS のダイアログが開きます。大切なのは、"Calculate electrostatic potential (ESP)" チェックボックスをオンにすることです(ここの手順に従えば、自動的にオンになっているはずですが、一応確認してください)。分子の電荷・スピン多重度が正しいことを確かめ,基底関数 (6-31G(d) がおすすめ) を指定してください。
1646 <p><img src="../etc/tutorial_8_09.png" /></p>
1648 OK ボタンを押して GAMESS 入力を作成し、GAMESS で計算を実行してください。計算には時間がかかりますから、この時点で Molby を終了して他の作業をしていただいて構いません(ただし、分子ファイルを mbsf 形式で保存することを忘れないで!)。
1651 GAMESS 計算が完了したら、同じ分子ファイルを開き、"MM/MD" → "Tools" → "GAMESS/RESP..." をもう一度選択してください。今度は、第2ステップの "Import GAMESS dat..." ボタンを押します。
1653 <p><img src="../etc/tutorial_8_10.png" /></p>
1655 GAMESS の dat ファイルを選択し(GAMESS の計算が成功すればこれができているはずです)、読み込んでください。読み込みが完了したら、"Run RESP..." ボタンが有効になっているはずです。もしそうならなかったら、読み込んだ dat ファイルに静電ポテンシャルのデータが含まれなかったということです。もう一度第1ステップからやり直し、特に "Calculate electrostatic potential (ESP)" チェックボックがオンになっていることを確かめてください。
1657 <p><img src="../etc/tutorial_8_11.png" /></p>
1659 "Run RESP..." ボタンを押すと、次のダイアログが開きます。これは、<a href="mm_minimize.html#minimize">前に説明した</a> Antechamber のダイアログとほとんど同じです。
1661 <p><img src="../etc/tutorial_8_12.png" /></p>
1663 "OK"ボタンを押すと、RESP 部分電荷が各原子に対してアサインされます。
1666 <link id="#navigation" />
1669 <div file="ruby.html">
1670 <link rel="Up" href="tutorials.html" />
1671 <link rel="Prev" href="qchem.html" />
1672 <link id="#header" />
1673 <div class="contents" lang="en">
1674 <h1>Step Nine: Using Embedded Ruby Interpreter</h1>
1676 <div class="contents" lang="ja">
1677 <h1>第九段階:組み込み Ruby インタプリタを使う</h1>
1679 <link id="#navigation" />
1682 <div file="reference.html">
1683 <link rel="Up" href="index.html" />
1684 <link rel="Prev" href="tutorials.html" />
1685 <link id="#header" />
1686 <div class="contents" lang="en">
1690 <div class="contents" lang="ja">
1694 <link id="#navigation" />
1697 <div file="ruby_ref.html">
1698 <link rel="Up" href="reference.html" />
1699 <link id="#header" />
1700 <div class="contents" lang="en">
1701 <h1>Ruby Extension Reference</h1>
1703 The following classes/module/exception are defined in the embedded Ruby interpreter in Molby. The methods defined in the module <code>Kernel</code> can be used as if they are builtin functions, like in the standard Ruby interpreter. All other classes and <code>MolbyError</code> exception are defined under the a module <code>Molby</code>. However, <code>Molby::</code> prefix is not necessary because <code>include Molby</code> is invoked on startup.
1706 On startup, <code>include Math</code> is also invoked, so that the methods and constants in the <code>Math</code> module can also be used without prefix.
1709 It is recommended for you to read the document of the most important class <a href="molby_rb/Molecule.html">Molecule</a>. Follow the links as necessary, and you will understand how other classes are defined and used.
1711 <h3 class="section-bar">Classes</h3>
1713 <li><a href="molby_rb/AtomRef.html">AtomRef</a></li>
1714 <li><a href="molby_rb/Dialog.html">Dialog</a></li>
1715 <li><a href="molby_rb/DialogItem.html">DialogItem</a></li>
1716 <li><a href="molby_rb/IntGroup.html">IntGroup</a></li>
1717 <li><a href="molby_rb/LAMatrix.html">LAMatrix</a></li>
1718 <li><a href="molby_rb/MDArena.html">MDArena</a></li>
1719 <li><a href="molby_rb/MolEnumerable.html">MolEnumerable</a></li>
1720 <li><a href="molby_rb/Molecule.html">Molecule</a></li>
1721 <li><a href="molby_rb/Parameter.html">Parameter</a></li>
1722 <li><a href="molby_rb/ParameterRef.html">ParameterRef</a></li>
1723 <li><a href="molby_rb/ParEnumerable.html">ParEnumerable</a></li>
1724 <li><a href="molby_rb/Transform.html">Transform</a></li>
1725 <li><a href="molby_rb/Vector3D.html">Vector3D</a></li>
1727 <h3 class="section-bar">Module</h3>
1729 <li><a href="molby_rb/Kernel.html">Kernel</a></li>
1730 <li><a href="molby_rb/MolbyModule.html">Molby</a></li>
1732 <h3 class="section-bar">Exception</h3>
1734 <li><a href="molby_rb/MolbyError.html">MolbyError</a></li>
1737 <div class="contents" lang="ja">
1738 <h1>Ruby 拡張リファレンス</h1>
1740 以下のクラス/モジュール/例外が Molby 内蔵の Ruby インタプリタに組み込まれています。<code>Kernel</code> モジュールのメソッドは、通常の Ruby インタプリタと同様に組み込み関数のように使うことができます。その他のクラスと <code>MolbyError</code> 例外は <code>Molby</code> モジュールの下に定義されています。ただし、起動時に <code>include Molby</code> が呼び出されるため、<code>Molby::</code> 接頭辞は必要ありません。
1743 起動時には <code>include Math</code> も呼び出されるため、<code>Math</code> モジュールのメソッドや定数も接頭辞なしで使うことができます。
1746 最も重要なクラスは <a href="molby_rb/Molecule.html">Molecule</a> なので、このドキュメントを最初に読むことをお勧めします。ここから必要に応じてリンクをたどっていけば、他のクラスがどのように定義され使われるかが理解しやすくなります。
1748 <h3 class="section-bar">Classes</h3>
1750 <li><a href="molby_rb/AtomRef.html">AtomRef</a></li>
1751 <li><a href="molby_rb/Dialog.html">Dialog</a></li>
1752 <li><a href="molby_rb/DialogItem.html">DialogItem</a></li>
1753 <li><a href="molby_rb/IntGroup.html">IntGroup</a></li>
1754 <li><a href="molby_rb/LAMatrix.html">LAMatrix</a></li>
1755 <li><a href="molby_rb/MDArena.html">MDArena</a></li>
1756 <li><a href="molby_rb/MolEnumerable.html">MolEnumerable</a></li>
1757 <li><a href="molby_rb/Molecule.html">Molecule</a></li>
1758 <li><a href="molby_rb/Parameter.html">Parameter</a></li>
1759 <li><a href="molby_rb/ParameterRef.html">ParameterRef</a></li>
1760 <li><a href="molby_rb/ParEnumerable.html">ParEnumerable</a></li>
1761 <li><a href="molby_rb/Transform.html">Transform</a></li>
1762 <li><a href="molby_rb/Vector3D.html">Vector3D</a></li>
1764 <h3 class="section-bar">Module</h3>
1766 <li><a href="molby_rb/Kernel.html">Kernel</a></li>
1767 <li><a href="molby_rb/MolbyModule.html">Molby</a></li>
1769 <h3 class="section-bar">Exception</h3>
1771 <li><a href="molby_rb/MolbyError.html">MolbyError</a></li>
1774 <link id="#navigation" />