データファイルからのデータの読み込み

	データファイルからのデータの読み込み
	ファイル

CHDメニュ   CHDコマンド   CHDダイアローグ   作図例

通常、CHDで使用する熱力学データはMALTデータベースからMALTダイレクト によって提供されます。
更に、CHDでは特殊な形式で保存されている熱力学データを読むことができます。現在、次のデータファイルを正確に読み込むことができます

1. MALT 化合物データファイル
   
2. MALT 第１版データファイル 拡張子".sc"を持つ     例(2)
   
3. CHDユーザデータファイル 拡張子".dat"を持つ     例(3) (理想溶液だけを混合相として使うことができます。).
   
4. テストデータファイル 拡張子 ".abs"を持つ　    例(4) （一般化された任意の空間に適用することができます。）
   例(4)-2（例としてLiFePO系のDFT計算データを使っています。

(2) MALT 第１版データファイル　"*.sc"

例(2)   トップへ
典型的な内容が以下に示されています。データはすべてＡＳＣＩＩテキストとして与えられ、カンマ","が区切り記号として用いられます。

1. 最初の行はファイルのタイトルを含むコメント行です。この行はデータを正しく読むために必要ですが、現在のＣＨＤで使われることはありません。
   
2. 各化合物に関する熱力学データが次に続きます。ひとつの化合物のデータは次の４つの異なる行で表されます。
   
   "*"から始まるコメント行;
   
   1. 化合物タイトル行 (一行) ６項目含まれます
      
      例: "malt1.0,La,La,g,g, 1"
      
      1. データソース名
         
      2. 化学式
         
      3. 印刷用化学式
         
      4. 298 Kにおける相
         
      5. 温度領域全体での相
         
      6. 高温相の数
         
      
      
   2. 化学量論数(7行) ２項目含まれます。
      
      例: "La, 1"
      
      1. 元素記号
         
      2. その化学量論数
         
      
      
   3. 高温データNo. 1
      
      例: "138.905 , 0 , 0 , 56.9 , 0 , 1"
      
      1. 熱容量式
         
         C_p / J K^-1mol^-1 = a + 10^-3 b T/K + 10⁵ c (K/T)² + 10^-6 d (T/K)² + 10⁸ e (K/T)³
         
         1. 熱容量式の係数　  a
            
         2. 熱容量式の係数　  b
            
         3. 熱容量式の係数　  c
            
         4. 熱容量式の係数　  d
            
         5. 熱容量式の係数　  e
            
         6. 文献情報
            
      
      
   4. 高温データNo. 2
      
      例: "$,c1, 298 , 550 , .364 ,tp,"
      
      1. 高温相の名前
         
      2. 高温相の相記号
         
      3. 有効温度領域の下限
         
      4. 有効温度領域の上限
         
      5. 相転移エンタルピー変化（上限で相転移があれば）
         
      6. 転移記号（上限が相転移であれば）
         
      
      
   データ行 (3) および (4) は各高温相毎に繰り返されます。
   

(2a) 例（MALT 第１版データファイル）

トップへ

* 87/08/12 [87/08/12]
malt1.0,La,La,g,g, 1
La, 1
0, 0
0, 0
0, 0
0, 0
0, 0
0, 0
138.905 , 431 , 393.59 , 182.268 , 22.753 , 1
$,g, 298 , 2500 , 0 ,$,
32.372 , .282 ,-28.21 , 0 , 5.869 ,-15.8248 ,-.193013 , 11,
malt1.0,La,La,cs,sl, 4
La, 1
0, 0
0, 0
0, 0
0, 0
0, 0
0, 0
138.905 , 0 , 0 , 56.9 , 0 , 1
$,c1, 298 , 550 , .364 ,tp,
27.07 , 1.34 ,-.33 , 0 , 0 ,-8.24116 ,-.154819 , 5,
$,c2, 550 , 1134 , 3.121 ,tp,
17.66 , 15.02 , 3.93 , 0 , 0 ,-3.99622 ,-.101601 , 5,
$,c3, 1134 , 1193 , 6.197 ,mp,
39.54 , 0 , 0 , 0 , 0 ,-16.3762 ,-.235862 , 5,
$,l, 1193 , 3730 , 0 ,bp,
34.31 , 0 , 0 , 0 , 0 ,-3.93978 ,-.193617 , 5,

(3) CHD ユーザーデータファイル

例(3)   トップへ
このデータファイルは次のフォーマットのユーザーデータファイルを読むためのものです。

1. データ全て、適切な順番で各行に対して決められた形式で与えられていなければなりません。
   
2. コメント行; (*) で始まる行はコメント行と定義され読み飛ばされます。
   
3. 何の文字情報を含まない行（空行）は読み飛ばされます。
   
4. 全てのデータ行は空白を区切り記号としていくつかのデータに分離されます。
   
系に関するデータは次のようになっています。

システム情報

例　  " 2 1 5 3"
１）構成要素の数(N)
２）混合相の数(MX)
３）から３+MX) I番目混合相の化学種の数
MX+4) 化学量論化合物の数

構成要素（元素）名前

例   "O Co "
構成要素の名前が順番にN要素分並んでいる。

構成要素（元素）分子量 : この情報はCHDでは使いません。

この情報はCHDでは使いません。
現在、CHDでは原子量は使いません。この行は読み飛ばされます。

ギブズエネルギー表現情報

例   " 6 1 2 3 4 5 6"
通常のギブズエネルギー表現
G = A +
     B * T +
     C * ln(T/K)*T +
     D * T² +
     E * T³ +
     F / T.
1) : 項の数
2) から 2+n) : 上の表現式における順序を用いて表現した用いる項

追加情報 この情報はCHDでは使いません。

例   " 1 1"
ギブズエネルギー表現としての追加情報です。
"1 1 " と入力して下さい。

主たるデータに関する行は次のようになっています。

混合相/化学量論化合物の名前

例   "GAS"
混合相の名前で、次の候補から選ぶ必要があります。

GAS
AQUEOUS
SOLID
LIQUID

混合相の種類

例   "IDMX"
次の中から選びます。

IDMX

各化合物の熱力学は次のようになっています。

化合物の名前

例   "O"
化合物の名前を入力します。

データの種類, 相の数, 化学量論数

例   "4 1 1.0000 0.0000"
次の情報を提供します。

1. データの種類 : 4 を入力します。
   
2. 高温相の数
   
3. 各元素毎の化学量論数
   
ギブズエネルギー No. 1..

例   "2.5000000E+03 2.4328149E+05 -2.1885825E+01 -2.0845000E+01 8.5000000E-06"
ギブズエネルギーのNo. 1です。

最初の項 : 有効温度領域の上限です。
データNO.1の残りの項とNo.2のデータとで、ギブズエネルギーを温度の関数として表すＮ＋３項を与えます。
最初の N 項 (N はこの前の入力で与えられます) + 追加の３項

ギブズエネルギーデータ No. 2

例   "0.0000000E+00 -4.8550000E+04 0.0000000E+00 0.0000000E+00 0.0000000E+00"
データNO.1の残りの項とNo.2のデータとで、ギブズエネルギーを温度の関数として表すＮ＋３項を与えます。
最初の N 項 (N はこの前の入力で与えられます) + 追加の３項

ギブズエネルギーデータ No. 3

例   "-2.0 0.0 0.0 "
-2は追加の項がT_-2（Ｔは温度です）であることを意味しています。

3a) 例 （CHDユーザーデータファイル）

トップへ
　 * ユーザが生成するデータファイル
2 1 5 3
O Co
15.99940 58.93320
6 1 2 3 4 5 6
1 1
GAS
IDMX
O
4 1 1.0000 0.0000
2.5000000E+03 2.4328149E+05 -2.1885825E+01 -2.0845000E+01 8.5000000E-06
0.0000000E+00 -4.8550000E+04 0.0000000E+00 0.0000000E+00 0.0000000E+00
-2.0 0.0 0.0
O2
41 2.0000 0.0000
2.5000000E+03 -1.4231317E+04 3.0877261E+01 -3.4350000E+01 -9.6000000E-04
0.0000000E+00 9.2250000E+05 -6.7666667E+07 0.0000000E+00 0.0000000E+00
-2.0 0.0 0.0
O3
41 3.0000 0.0000
2.5000000E+03 1.1552677E+05 1.6258143E+02 -5.7780000E+01 -3.5500000E-04
0.0000000E+00 2.1245000E+06 -1.2850000E+08 0.0000000E+00 0.0000000E+00
-2.0 0.0 0.0
Co
41 0.0000 1.0000
2.5000000E+03 4.1556568E+05 2.1115903E+00 -2.6819000E+01 3.1500000E-05
0.0000000E+00 1.7010000E+05 0.0000000E+00 0.0000000E+00 0.0000000E+00
-2.0 0.0 0.0
CoO
4 1 1.0000 1.0000
4.0000000E+03 2.8028824E+05 1.4209312E+01 -3.6490000E+01 -2.6000000E-04
0.0000000E+00 2.6970000E+05 0.0000000E+00 0.0000000E+00 0.0000000E+00
-2.0 0.0 0.0
Co
4 4 0.0000 1.0000
7.0000000E+02 -6.7380728E+03 1.1042648E+02 -2.0280000E+01 -7.7800000E-03
0.0000000E+00 0.0000000E+00 0.0000000E+00 0.0000000E+00 0.0000000E+00
-2.0 0.0 0.0
1.3940000E+03 -1.2804053E+04 2.6397436E+02 -4.5270000E+01 2.3620000E-02
-6.4300000E-06 0.0000000E+00 0.0000000E+00 0.0000000E+00 0.0000000E+00
-2.0 0.0 0.0
1.7680000E+03 2.6707448E+05 -1.0849677E+03 1.2803000E+02 -3.2700000E-02
0.0000000E+00 -7.8695000E+07 0.0000000E+00 0.0000000E+00 0.0000000E+00
-2.0 0.0 0.0
3.2010000E+03 -1.5016034E+03 2.4415313E+02 -4.0500000E+01 0.0000000E+00
0.0000000E+00 0.0000000E+00 0.0000000E+00 0.0000000E+00 0.0000000E+00
-2.0 0.0 0.0
CoO
4 2 1.0000 1.0000
2.0780000E+03 -2.5477621E+05 3.1589943E+02 -5.5100000E+01 2.3000000E-03
-9.9000000E-07 8.3500000E+04 0.0000000E+00 0.0000000E+00 0.0000000E+00
-2.0 0.0 0.0
2.5000000E+03 -2.0404540E+05 3.3456010E+02 -6.0670000E+01 0.0000000E+00
0.0000000E+00 0.0000000E+00 0.0000000E+00 0.0000000E+00 0.0000000E+00
-2.0 0.0 0.0
Co3O4
4 1 4.0000 3.0000
1.2200000E+03 -9.4106306E+05 8.3592218E+02 -1.3665000E+02 -1.3765000E-02
-7.9216667E-06 1.1445000E+06 0.0000000E+00 0.0000000E+00 0.0000000E+00
-2.0 0.0 0.0

(4) テストデータファイル

例   （ トップへ）

通常の（熱力学空間ではない）空間で多面体を構築するためのファイルです。
空間内の各面は次の式で表されます。

la + mb + nc = const

ここで、次元数は３です。変数 a, b ,c, が空間を構成します。係数, l, m, n, が平面の勾配をきめます。Constが平面の位置を決めます。

データファイルの最初の行が次元数を与えます。
各面のデータが次の３行から成ります。

1. 面の名前
   
2. 係数データ
   
3. 定数値のデータ
   

(4a) 例1

トップへ
（テストデータファイル） （テストデータファイル例２）

3
A
1, 0, 0
0
B
0, 1, 0
0
C
0, 0, 1
0
AB
1, 1, 0
-10
AC
1, 0, 1
-20


例2: 作図例9ｂで LiFePO系のDFTデータを lifepo.abs に記述し、3次元図を作成する例を提示している。 この例では、最初の行に次元数を入力するとともに、元素名も入力している。

トップへ
（テストデータファイル） （テストデータファイル例）