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

从头算分子轨道法预测分子间势用于制冷剂特性计算

• 批准号: 22K03946
• 负责人: 近藤 智恵子
• 金额: $ 2.66万
• 依托单位: Nagasaki University
• 依托单位国家: 日本
• 项目类别: Grant-in-Aid for Scientific Research (C)
• 财政年份: 2022
• 资助国家: 日本
• 起止时间: 2022-04-01 至 2026-03-31
• 项目状态: 未结题
• 来源: https://kaken.nii.ac.jp/grant/KAKENHI-PROJECT-22K03946/
• 关键词: 冷媒; 温暖化係数; 分子シミュレーション; 分子軌道計算; 分子間相互作用; 熱物性

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

The next generation of Freon compounds HFO1123, HFO 1132 (E), HFO 1132a, and the molecular orbital method of non-hydrocarbon are used to optimize molecular structure and atomic charge. The binding length and binding angle of the molecule are calculated. The calculation results are used to construct the intramolecular force field. Universal force field OPLS-AA/L is compared with the combination of long and good, and the combination angle is somewhat different. MP2 is used to analyze two molecules, and intermolecular interactions are calculated. The calculation results show that the molecular interaction is very close to the Lennard-Jones relationship, and the atomic interaction is very close. The results of molecular dynamics calculations provide a reference for the latest high-precision determination of liquid density near the critical point. The results of molecular dynamics calculations for the liquid-phase equilibrium system confirm that the saturated liquid density is consistent within 1%. In addition, the saturation vapor density, saturation pressure and surface tension are calculated in a leap in accuracy. HFO1123 is calculated in a lean vapor reproduction process. The average force is calculated in a second step. The equation of state is proposed to be in good agreement with the temperature range. The equation of state for the 3rd coefficient of gravity is the equation of gravity. The third coefficient is lower than the second coefficient, and the calculation accuracy of the equation of state is lower than that of the third coefficient. The temperature of the sample is determined by the theoretical calculation of Thol.