Rotational partition function

In chemistry, the rotational partition function relates the rotational degrees of freedom to the rotational part of the energy.

Definition

The total canonical partition functionZ{\displaystyle Z} of a system of N{\displaystyle N} identical, indistinguishable, noninteracting atoms or molecules can be divided into the atomic or molecular partition functions ζ{\displaystyle \zeta }:[1]Z=ζNN!{\displaystyle Z={\frac {\zeta ^{N}}{N!}}} with: ζ=jgjeEj/kBT,{\displaystyle \zeta =\sum _{j}g_{j}e^{-E_{j}/k_{\text{B}}T},} where gj{\displaystyle g_{j}} is the degeneracy of the jth quantum level of an individual particle, kB{\displaystyle k_{\text{B}}} is the Boltzmann constant, and T{\displaystyle T} is the absolute temperature of system. For molecules, under the assumption that total energy levels Ej{\displaystyle E_{j}} can be partitioned into its contributions from different degrees of freedom (weakly coupled degrees of freedom)[2]Ej=iEji=Ejtrans+Ejns+Ejrot+Ejvib+Eje{\displaystyle E_{j}=\sum _{i}E_{j}^{i}=E_{j}^{\text{trans}}+E_{j}^{\text{ns}}+E_{j}^{\text{rot}}+E_{j}^{\text{vib}}+E_{j}^{\text{e}}} and the number of degenerate states are given as products of the single contributions gj=igji=gjtransgjnsgjrotgjvibgje,{\displaystyle g_{j}=\prod _{i}g_{j}^{i}=g_{j}^{\text{trans}}g_{j}^{\text{ns}}g_{j}^{\text{rot}}g_{j}^{\text{vib}}g_{j}^{\text{e}},} where "trans", "ns", "rot", "vib" and "e" denotes translational, nuclear spin, rotational and vibrational contributions as well as electron excitation, the molecular partition functions ζ=jgjeEj/kBT{\displaystyle \zeta =\sum _{j}g_{j}e^{-E_{j}/k_{\text{B}}T}} can be written as a product itself ζ=iζi=ζtransζnsζrotζvibζe.{\displaystyle \zeta =\prod _{i}\zeta ^{i}=\zeta ^{\text{trans}}\zeta ^{\text{ns}}\zeta ^{\text{rot}}\zeta ^{\text{vib}}\zeta ^{\text{e}}.}

Linear molecules

Rotational energies are quantized. For a diatomic molecule like CO or HCl, or a linear polyatomic molecule like OCS in its ground vibrational state, the allowed rotational energies in the rigid rotor approximation are EJrot=J22I=J(J+1)22I=J(J+1)B.{\displaystyle E_{J}^{\text{rot}}={\frac {\mathbf {J} ^{2}}{2I}}={\frac {J(J+1)\hbar ^{2}}{2I}}=J(J+1)B.} J is the quantum number for total rotational angular momentum and takes all integer values starting at zero, i.e., J=0,1,2,{\displaystyle J=0,1,2,\ldots }, B=22I{\displaystyle B={\frac {\hbar ^{2}}{2I}}} is the rotational constant, and I{\displaystyle I} is the moment of inertia. Here we are using B in energy units. If it is expressed in frequency units, replace B by hB in all the expression that follow, where h is the Planck constant. If B is given in units of cm1{\displaystyle \mathrm {cm^{-1}} }, then replace B by hcB where c is the speed of light in vacuum.

For each value of J, we have rotational degeneracy, gj{\displaystyle g_{j}} = (2J+1), so the rotational partition function is therefore ζrot=J=0gjeEJ/kBT=J=0(2J+1)eJ(J+1)B/kBT.{\displaystyle \zeta ^{\text{rot}}=\sum _{J=0}^{\infty }g_{j}e^{-E_{J}/k_{\text{B}}T}=\sum _{J=0}^{\infty }(2J+1)e^{-J(J+1)B/k_{\text{B}}T}.}

For all but the lightest molecules or the very lowest temperatures we have BkBT{\displaystyle B\ll k_{\text{B}}T}. This suggests we can approximate the sum by replacing the sum over J by an integral of J treated as a continuous variable. ζrot0(2J+1)eJ(J+1)B/kBTdJ=kBTB.\displaystyle \zeta ^{\text{rot}}\approx \int _{0}^{\infty }(2J+1)e^{-J(J+1)B/k_{\text{B}}T}dJ={\frac {k_{\text{B}}T}{B}}.}

This approximation is known as the high temperature limit. It is also called the classical approximation as this is the result for the canonical partition function for a classical rigid rod.

باستخدام صيغة أويلر-ماكلورين، يمكن إيجاد تقدير محسّن [ 3 ]ζتعفن=كبتيب+13+115(بكبتي)+4315(بكبتي)2+1315(بكبتي)3+.{\displaystyle \zeta ^{\text{rot}}={\frac {k_{\text{B}}T}{B}}+{\frac {1}{3}}+{\frac {1}{15}}\left({\frac {B}{k_{\text{B}}T}}\right)+{\frac {4}{315}}\left({\frac {B}{k_{\text{B}}T}}\right)^{2}+{\frac {1}{315}}\left({\frac {B}{k_{\text{B}}T}}\right)^{3}+\cdots .}

بالنسبة لجزيء أول أكسيد الكربون عندتي=300 ك{\displaystyle T=\mathrm {300~K} }المساهمة (بدون وحدة)ζتعفن{\displaystyle \zeta ^{\text{rot}}}لζ{\displaystyle \zeta }تبين أن نطاقها102{\displaystyle 10^{2}}.

يمكن الآن حساب متوسط ​​الطاقة الدورانية الحرارية لكل جزيء عن طريق اشتقاقζتعفن{\displaystyle \zeta ^{\text{rot}}}فيما يتعلق بدرجة الحرارةتي{\displaystyle T}في تقريب حد درجة الحرارة العالية، تكون الطاقة الدورانية الحرارية المتوسطة لدوار صلب خطي هيكبتي{\displaystyle k_{\text{B}}T}.

تأثيرات التناظر الكمي

بالنسبة لجزيء ثنائي الذرة ذي مركز تناظر، مثلح2،شمال2،جيا2،{\displaystyle {\rm {H_{2},N_{2},CO_{2},}}}أوح2ج2{\displaystyle \mathrm {H_{2}C_{2}} }(أيدح{\displaystyle D_{\infty h}}المجموعة النقطية )، دوران الجزيء بواسطةπ{\displaystyle \pi }سيؤدي الدوران بزاوية π/π حول محور عمودي على محور الجزيء ويمر بمركز كتلته إلى تبادل أزواج من الذرات المتكافئة. تنص نظرية الإحصاء الدوراني في ميكانيكا الكم على أن الدالة الموجية الجزيئية الكلية إما متناظرة أو غير متناظرة بالنسبة لهذا الدوران، وذلك اعتمادًا على ما إذا كان عدد أزواج النوى الفرميونية المتبادلة زوجيًا أم فرديًا. ستكون الدالة الموجية الإلكترونية والاهتزازية المعطاة إما متناظرة أو غير متناظرة بالنسبة لهذا الدوران. ستشهد الدالة الموجية الدورانية ذات العدد الكمي J تغيرًا في الإشارة.(-1)ج{\displaystyle (-1)^{J}}يمكن تصنيف حالات اللف المغزلي النووي إلى حالات متناظرة أو غير متناظرة بالنسبة للتباديل النووية الناتجة عن الدوران. في حالة جزيء ثنائي الذرة متناظر ذي عدد كمي لللف المغزلي النووي I لكل نواة، يوجد(أنا+1)(2أنا+1){\displaystyle (I+1)(2I+1)}دوال الدوران المتناظرة وأنا(2أنا+1){\displaystyle I(2I+1)}هي دوال مضادة للتناظر لعدد إجمالي من الدوال النوويةزns=(2أنا+1)2{\displaystyle g^{\text{ns}}=(2I+1)^{2}}. Nuclei with an even nuclear mass number are bosons and have integer nuclear spin quantum number, I. Nuclei with odd mass number are fermions and had half integer I. For the case of H2, rotation exchanges a single pair of fermions and so the overall wavefunction must be antisymmetric under the half rotation. The vibration-electronic function is symmetric and so the rotation-vibration-electronic will be even or odd depending upon whether J is an even or odd integer. Since the total wavefunction must be odd, the even J levels can only use the antisymmetric functions (only one for I = 1/2) while the odd J levels can use the symmetric functions ( three for I = 1/2). For D2, I = 1 and thus there are six symmetric functions, which go with the even J levels to produce an overall symmetric wavefunction, and three antisymmetric functions that must go with odd J rotational levels to produce an overall even function. The number of nuclear spin functions that are compatible with a given rotation-vibration-electronic state is called the nuclear spin statistical weight of the level, often represented as gJ{\displaystyle g_{J}}. Averaging over both even and odd J levels, the mean statistical weight is (1/2)(2I+1)2{\displaystyle (1/2)(2I+1)^{2}}, which is one half the value of gns{\displaystyle g^{\text{ns}}} expected ignoring the quantum statistical restrictions. In the high temperature limit, it is traditional to correct for the missing nuclear spin states by dividing the rotational partition function by a factor σ=2{\displaystyle \sigma =2} with σ{\displaystyle \sigma } known as the rotational symmetry number which is 2 for linear molecules with a center of symmetry and 1 for linear molecules without.

Nonlinear molecules

A rigid, nonlinear molecule has rotational energy levels determined by three rotational constants, conventionally written A,B,{\displaystyle A,B,} and C{\displaystyle C}, which can often be determined by rotational spectroscopy. In terms of these constants, the rotational partition function can be written in the high temperature limit as [4]ζrotπσ(kBT)3ABC{\displaystyle \zeta ^{\text{rot}}\approx {\frac {\sqrt {\pi }}{\sigma }}{\sqrt {\frac {(k_{\text{B}}T)^{3}}{ABC}}}} with σ{\displaystyle \sigma } again known as the rotational symmetry number [5] which in general equals the number ways a molecule can be rotated to overlap itself in an indistinguishable way, i.e. that at most interchanges identical atoms. Like in the case of the diatomic treated explicitly above, this factor corrects for the fact that only a fraction of the nuclear spin functions can be used for any given molecular level to construct wavefunctions that overall obey the required exchange symmetries. Another convenient expression for the rotational partition function for symmetric and asymmetric tops is provided by Gordy and Cook: ζrot5.34×106σT3ABC{\displaystyle \zeta ^{\text{rot}}\approx {\frac {5.34\times 10^{6}}{\sigma }}{\sqrt {\frac {T^{3}}{ABC}}}} where the prefactor comes from (πkB)3h3=5.34×106{\displaystyle {\sqrt {\frac {(\pi k_{\text{B}})^{3}}{h^{3}}}}=5.34\times 10^{6}} عندما يتم التعبير عن A و B و C بوحدات ميغاهرتز. [ 6 ]

التعبيرات لـζتعفن{\displaystyle \zeta ^{\text{rot}}}يعمل مع الدوارات العلوية غير المتماثلة والمتماثلة والكروية.

مراجع

  1. دونالد أ. ماكواري، الميكانيكا الإحصائية ، هاربر آند رو، 1973
  2. دونالد أ. ماكواري، المرجع نفسه
  3. جي. هيرزبرغ، أطياف الأشعة تحت الحمراء وأطياف رامان ، فان نوستراند رينهولد، 1945، المعادلة (V,21)
  4. ج. هيرزبرغ، المرجع نفسه ، المعادلة (V,29)
  5. ج. هيرزبرغ، المرجع نفسه ؛ انظر الجدول 140 للاطلاع على قيم مجموعات النقاط الجزيئية الشائعة
  6. كوك، روبرت ل.؛ جوردي، والتر (1970). أطياف الجزيئات الميكروية (  الطبعة الثانية). دار إنترساينس للنشر. الصفحات 56-57 . ISBN  0-471-08681-9.

انظر أيضاً