Štejfa Vojtěch, Fulem Michal, Růžička Květoslav
Department of Physical Chemistry, University of Chemistry and Technology, Prague, Technická 5, CZ-166 28 Prague 6, Czech Republic.
J Chem Phys. 2019 Jun 14;150(22):224101. doi: 10.1063/1.5093767.
First-principles calculations, coupled with statistical thermodynamics, can provide ideal-gas thermodynamic properties but get complicated and less reliable with an increasing number of conformers. An approach designed for calculation of ideal-gas thermodynamic properties of long-chain molecules, R1SM, and its simplified version, sR1SM, is tested in this work by calculation of ideal-gas heat capacities and entropies for a homologous series of n-alkanes up to n-tetradecane. The R1SM approach incorporates the rigid rotor-harmonic oscillator approximation in combination with a correction for internal rotations of methyl tops using the one-dimensional hindered rotor scheme and the mixing model accounting for the population of conformers based on the Boltzmann distribution. The R1SM approach is applicable for compounds with up to hundreds of conformers, while the simplified sR1SM approach can be used for molecules with up to 10 conformers when coupled with rules for enumeration of stable conformers and estimation scheme for their energies. The obtained results for n-alkanes are compared with experimental values and previously employed computational schemes. As the conformational behavior and conformer energies are inherent parts of the proposed approaches, a thorough conformational study of n-alkanes is performed and compared with experiments and the Tasi rules for enumeration of n-alkane conformers. Finally, the standard uncertainty of the R1SM-calculated ideal-gas thermodynamic properties is estimated based on the error propagation from the used input quantities and approximations as well as on comparison to experimental values and amounts to less than 1% for both ideal-gas heat capacity and standard ideal-gas entropy.
第一性原理计算与统计热力学相结合,可以提供理想气体的热力学性质,但随着构象异构体数量的增加,计算会变得复杂且可靠性降低。本文通过计算直至十四烷的一系列正构烷烃的理想气体热容和熵,对一种专为计算长链分子理想气体热力学性质而设计的方法R1SM及其简化版本sR1SM进行了测试。R1SM方法结合了刚性转子-谐振子近似,并使用一维受阻转子方案对甲基顶的内旋转进行校正,以及基于玻尔兹曼分布的混合模型来考虑构象异构体的分布。R1SM方法适用于具有多达数百个构象异构体的化合物,而简化的sR1SM方法在与稳定构象异构体的枚举规则及其能量估计方案相结合时,可用于具有多达10个构象异构体的分子。将正构烷烃的计算结果与实验值和先前使用的计算方案进行了比较。由于构象行为和构象异构体能量是所提出方法的固有组成部分,因此对正构烷烃进行了全面的构象研究,并与实验以及正构烷烃构象异构体枚举的塔西规则进行了比较。最后,基于所用输入量和近似值的误差传播以及与实验值的比较,估计了R1SM计算的理想气体热力学性质的标准不确定度,对于理想气体热容和标准理想气体熵,该不确定度均小于1%。