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 refrigerant is expected to be HFO1123, HFO1132 (E), HFO1132a, and HFO1132a, a futurized chemical material.ついて, non-経験 molecular orbital method の一つであるカプルドクラスターmethod をUsing いて molecular structure optimization およびatomic charge の Cut り てを行った. Intra-molecular binding length and binding angle are calculated, and the calculation results areを行った． The general force field OPLS-AA/L is a comparative one, a long one is a good one and a consistent one, and a combined angle one is a little bit eccentric. The sub-action method MP2 uses the two-dimensional body analysis and calculation of intermolecular interactions.このCalculation resultsから, intermolecular interaction エネルギーをapproximationするLennard-Jo nes off number のパラメータについて, パラメータ adjustment を行うべきatomic を选択した． We provide fine adjustment, molecular dynamics calculation results, and the latest high-precision measurement of liquid density near the critical point. The results of the molecular dynamics calculation of the liquid equilibrium system are consistent with each other and the saturated liquid density is within 1%, and the results are confirmed.さらに, saturated vapor density, saturated pressure force および surface tension についても calculation accuracy が に upward す る こ と を confirmation し た. HFO1123では, thin steamed water reappears, average force HFO1123 second hand coefficient is calculated, ProposalされているState equationがComparison withえる値として広いTemperature range囲にわたって好くConsistencyしていることをConfirmationした．しかしながら, the 3rd ビリアル coefficient については state equation から大きくdeviationした． The calculation accuracy of the 3rd and 2nd ビリアル coefficients is lower than that of the してstate equation.さらにボイル Temperature についても検证を行い, TholらのTheoretical calculation 〤とapproximately consistent することをconfirmed した.