◎１．水素の原子軌道関数を見てみましょう。
DV-Xα計算環境（コマンドプロンプト）を開きます。
cd calc[Enter]
・・・C:\dvxa から C:\dvxa\calc に移動（チェンジディレクトリ）します。
md H[Enter]
・・・C:\dvxa\calc に“H”というディレクトリを作成（メイクディレクトリ）します。
cd H[Enter]
・・・C:\dvxa\calc から C:\dvxa\calc\H に移動（チェンジディレクトリ）します。
atomn[Enter]
・・・プログラム“atomn”を起動します。
Input Atomic Number(Z) (from 1 to 94), Z=?
と聞いてきますので、水素の原子番号 = 1 を入力してください。
ここまでの作業で、C:\dvxa\calc\Hには、
・F01
・F09
・F39
・c04d
・H1s--00.sca 水素の1s軌道(n = 1,
(l = 0, m = 0))の原子軌道関数（波動関数）が出来ています。
※n = 主量子数
※
(l = 方位量子数, m = 磁気量子数) は、実数型の球面調和関数です。詳しくは、教科書・量子材料学の初歩のp.14〜25を参照してください。
次に、VENUS\Programs\VEND.exeをダブルクリックして、VENDを起動します。
VEND: Menu
Dialog
Graphics
VEND
の４つのウィンドウが立ち上がります。
VEND: Menu の 【Open】をクリックし、Format: のプルダウンメニューで、
上から５行目の【SCAT (*.scat, *.sca)】を選択、【Browse...】をクリックし、
ファイルの場所(I): にc:\dvxa\calc\H を選んで、H1s--00.sca を選択し、
【開く(O)】をクリックします。
Dialogウィンドウの【Properties + 】をクリック、【Isosurfaces】をクリックし、
Isosurface level: の値を、例えば0.02a0-3/2など、
小さい値にして、【Apply】をクリックします。
この値を小さくすればするほど、原子軌道関数の拡がりは大きく表示されます。
この値を大きくすればするほど、原子軌道関数の拡がりは小さく表示されます。
さらに、Smoothingの【More】をクリックすれば、原子軌道関数の等値表面が
スムージングされて、滑らかになります。
原子軌道関数の等値表面の内部を見たい場合は、
Opacity (%): 100
となっている箇所が右下の欄にありますので、これを例えば、
Opacity (%): 50
として、【Apply】をクリックします。
そうしますと、50 %ぐらい透明に透けて見えるようになります。
Opacityは不透明度という意味です。
100 %が普通の状態で、内部は全く見えませんが、
50 %だと内部がかなり透けて見えるようになります。
【OK】をクリックしてIsosurfacesウィンドウを閉じます。
Graphicsウィンドウをマウスの左ドラッグで好きな位置に回転させてやり、
VEND: Menuの【Export】をクリックし、Export Imageウィンドウで適当な名前をつけて
ファイルの種類(T):
BMP (*.bmp)
EPS (*.eps)
JPEG (*.jpg)
JPEG 2000 (*.jp2)
PPM (*.ppm)
RAW (*.raw)
RGB (*.rgb)
TGA (*.tga)
TIFF (*.tif)
EPS (vector data) (*.eps)
の中から、適当なファイル書式を選択して保存してやれば、各種お絵描きソフト、ワープロ、
プレゼンテーションソフト等に持っていけるかと思います。
以下に、以上の操作で描いた、水素の原子軌道関数を以下に示します。
| 図 タ イ ト ル |
フ ァ イ ル 名 |
iso sur face level *1 |
n 主 量 子 数 |
実 数 型 球 面 調 和 関 数 |
l 方 位 量 子 数 |
m 磁 気 量 子 数 |
Smo oth ing |
con trd # me sh |
V I C S |
Opa city (%) |
Animated GIF 画像 画像の下の ファイルサ イズの数字 (** KB)を クリックし てください |
JPEG 静止画像 [1]〜[5] をクリックし てください いずれも 40 KB 前後です |
[eV] | 電 子 数 |
| H 1s |
H1s --00 .sca |
+0.02 -0.02 |
1 | |
0 | 0 | ○ | 61 x 61 x 61 |
− | 100 | H1s--00 62 KB 270 KB 545 KB |
H1s--00 [1] [2] [3] [4] [5] |
-1.19 | 1 |
*1) isosurface levelの単位は、(bohr)-3/2 です。
ここでbohrは、長さの原子単位のことで、
1 bohr = a0 = 5.291772108(18)×10-11 m = 0.5291772108(18) Åです。
a0は基礎物理定数の１つで、ボーア半径のことです。a0については、こちらを参照してください。
→ もくじに戻る
岡山理科大学 理学部 化学科 坂根弦太
ご意見・ご質問は下記までｇｓａｋａｎｅ＠ｃｈｅｍ．ｏｕｓ．ａｃ．ｊｐ