Inter molecular potential prediction by ab initio molecular orbital method for refrigerant property calculation

次世代冷媒として期待されているフッ化エチレン系物質HFO1123，HFO1132（E），およびHFO1132aについて，非経験的分子軌道法の一つであるカプルドクラスター法を用いて分子構造最適化および原子電荷の割り当てを行った．分子内の結合長・結合角について変形エネルギーを計算し，計算結果にフィッティングを行ったパラメータを用いて分子内力場の構築を行った．汎用力場OPLS-AA/Lと比較したところ，結合長は良く一致するが，結合角については幾分乖離が見られた．次いで摂動法MP2を用いて二量体解析を行い，分子間相互作用エネルギーの増減を計算した．この計算結果から，分子間相互作用エネルギーを近似するLennard-Jones関数のパラメータについて，パラメータ調整を行うべき原子を選択した．この微調整には，分子動力学計算結果が，臨界点近傍で提供されている最新の気液密度の高精度測定値を参照した．気液平衡系の分子動力学計算の結果，それぞれの物質で，飽和液密度がおおむね1％以内で一致することを確認した．さらに，飽和蒸気密度，飽和圧力および表面張力についても計算精度が飛躍的に向上することを確認した．HFO1123では，稀薄蒸気についても再現を行い，平均力ポテンシャルから第2ビリアル係数を算出し，提案されている状態方程式が与える値と比較して広い温度範囲にわたって良く一致していることを確認した．しかしながら，第３ビリアル係数については状態方程式から大きく乖離した．第３ビリアル係数は，第2ビリアル係数に比して状態方程式の計算精度が著しく低下することが要因として考えられる．さらにボイル温度についても検証を行い，Tholらの理論計算値と概ね一致することを確認した．

次世代冷媒として期待されているフッ化エチレン系物質HFO1123，HFO1132（E），およびHFO1132aについて，非経験的分子軌道法の一つであるカプルドクラスター法を用いて分子構造最適化および原子電荷の割り当てを行った．分子内の結合長・結合角について変形エネルギーを計算し，計算結果にフィッティングを行ったパラメータを用いて分子内力場の構築を行った．汎用力場OPLS-AA/Lと比較したところ，結合長は良く一致するが，結合角については幾分乖離が見られた．次いで摂動法MP2を用いて二量体解析を行い，分子間相互作用エネルギーの増減を計算した．この計算結果から，分子間相互作用エネルギーを近似するLennard-Jones関数のパラメータについて，パラメータ調整を行うべき原子を選択した．この微調整には，分子動力学計算結果が，臨界点近傍で提供されている最新の気液密度の高精度測定値を参照した．気液平衡系の分子動力学計算の結果，それぞれの物質で，飽和液密度がおおむね1％以内で一致することを確認した．さらに，飽和蒸気密度，飽和圧力および表面張力についても計算精度が飛躍的に向上することを確認した．HFO1123では，稀薄蒸気についても再現を行い，平均力ポテンシャルから第2ビリアル係数を算出し，提案されている状態方程式が与える値と比較して広い温度範囲にわたって良く一致していることを確認した．しかしながら，第３ビリアル係数については状態方程式から大きく乖離した．第３ビリアル係数は，第2ビリアル係数に比して状態方程式の計算精度が著しく低下することが要因として考えられる．さらにボイル温度についても検証を行い，Tholらの理論計算値と概ね一致することを確認した．