3.8.4 IROTSYS

Next: 3.9 PPの指定と原子初期位置の入力 Up: 3.8 データ形式などの入力 Previous: 3.8.3 IUNIT

3.8.4 IROTSYS

セルに沿った直交座標系に配置したクラスターを座標軸に回転操作を施すことにより、クラスターを傾いた配置で計算するときに使用する。クラスターなどをセルの[111]方向などに合わせて配置するときに便利である。IROTSYS=1のとき、原子配置座標をinit1.fで読み込んだ後、ファイル(29)よりオイラー角( $\theta$ , $\phi$ , $\psi$ )を読み込んで、次に示す回転操作で座標変換を行う。回転後の座標は、

$\displaystyle {\bf r}$	$\textstyle =$	$\displaystyle x {\bf e}_{1} + y {\bf e}_{2} + z {\bf e}_{3}$	(19)
	$\textstyle =$	$\displaystyle ( {\bf e}_{1},{\bf e}_{2},{\bf e}_{3}) \left( \begin{array}{c} x y z \end{array}\right)$	(20)
	$\textstyle =$	$\displaystyle ( {\bf e}_{x},{\bf e}_{y},{\bf e}_{z})\ R_{z}(\phi) R_{x}(\theta) R_{z}(\psi)\ \left( \begin{array}{c} x y z \end{array}\right)$	(21)

$\displaystyle R_{z}(\phi)$	$\textstyle =$	$\displaystyle \left( \begin{array}{ccc} \cos(\phi) & -\sin(\phi) & 0 \\ \sin(\phi) & \cos(\phi) & 0 \\ 0 & 0 & 1 \end{array}\right)$	(22)
$\displaystyle R_{x}(\theta)$	$\textstyle =$	$\displaystyle \left( \begin{array}{ccc} 1 & 0 & 0 \\ \cos(\theta) & -\sin(\theta) & 0 \\ \sin(\theta) & \cos(\theta) & 0 \end{array}\right)$	(23)

下記のオイラー角は、z軸を[111]軸に向ける場合とy軸を[111]軸に向ける場合のための入力である。

cat > ${tmpdir}/eangle.dat << end
  -54.73561 -45.0 0.0         THETA, PHI, PSI (z-axis --> [111]axis)
   35.26439 -45.0 0.0         THETA, PHI, PSI (y-axis --> [111]axis)
end

初期座標をxyz形式の原子座標をファイル(16)に出力しているので、そのファイルで初期座標を確認することができる。

Next: 3.9 PPの指定と原子初期位置の入力 Up: 3.8 データ形式などの入力 Previous: 3.8.3 IUNIT