التقلبات الحرارية

الانتشار الذري على سطح البلورة. يُعدّ اهتزاز الذرات مثالاً على التقلبات الحرارية. وبالمثل، توفر التقلبات الحرارية الطاقة اللازمة للذرات للانتقال من موقع إلى آخر مجاور. ولتبسيط الرسم، لم تُعرض التقلبات الحرارية للذرات الزرقاء.

في الميكانيكا الإحصائية ، تُعرَّف التقلبات الحرارية بأنها انحرافات عشوائية لنظام ذري عن حالته المتوسطة، تحدث في نظام في حالة اتزان. [ 1 ] وتزداد جميع التقلبات الحرارية حجماً وتواتراً مع ارتفاع درجة الحرارة، وبالمثل، تتناقص مع اقتراب درجة الحرارة من الصفر المطلق .

تُعد التقلبات الحرارية مظهرًا أساسيًا لدرجة حرارة الأنظمة: فالنظام عند درجة حرارة غير صفرية لا يبقى في حالته المجهرية المتوازنة، ولكنه بدلاً من ذلك يأخذ عينات عشوائية من جميع الحالات الممكنة، باحتمالات معطاة بواسطة توزيع بولتزمان .

تؤثر التقلبات الحرارية بشكل عام على جميع درجات حرية النظام: يمكن أن تكون هناك اهتزازات عشوائية ( فونونات )، ودورانات عشوائية ( روتونات )، وإثارات إلكترونية عشوائية، وما إلى ذلك.

تخضع المتغيرات الديناميكية الحرارية ، مثل الضغط ودرجة الحرارة والإنتروبيا ، لتقلبات حرارية أيضاً. فعلى سبيل المثال، بالنسبة لنظام ذي ضغط متوازن، يتذبذب ضغط النظام إلى حد ما حول قيمة التوازن.

إن "المتغيرات الضابطة" فقط للمجموعات الإحصائية (مثل عدد الجسيمات N والحجم V والطاقة الداخلية E في المجموعة الميكروكانونية ) لا تتقلب.

تُعدّ التقلبات الحرارية مصدرًا للتشويش في العديد من الأنظمة. وتُشكّل القوى العشوائية المُسببة لهذه التقلبات مصدرًا لكلٍّ من الانتشار والتبديد (بما في ذلك التخميد واللزوجة ) . وترتبط التأثيرات المتنافسة للانجراف العشوائي ومقاومة الانجراف بنظرية التقلب والتبديد . وتلعب التقلبات الحرارية دورًا رئيسيًا في التحولات الطورية والحركية الكيميائية .

نظرية النهاية المركزية

حجم فضاء الطورV{\displaystyle {\mathcal {V}}}، يشغلها نظام من2م{\displaystyle 2m}درجات الحرية هي ناتج حجم التكوينV{\displaystyle V}وحجم فضاء الزخم. وبما أن الطاقة هي شكل تربيعي للزخم في نظام غير نسبي، فإن نصف قطر فضاء الزخم سيكونهـ{\displaystyle {\sqrt {E}}}وبالتالي فإن حجم الكرة الفائقة سيتغير كماهـ2م{\displaystyle {\sqrt {E}}^{2m}}مما يعطي حجم طور

V=(جهـ)مΓ(م+1)،{\displaystyle {\mathcal {V}}={\frac {(C\cdot E)^{m}}{\Gamma (m+1)}},}

أينج{\displaystyle C}هو ثابت يعتمد على الخصائص المحددة للنظام و Γ{\displaystyle \Gamma }هي دالة غاما . في حالة كون هذه الكرة الفائقة ذات أبعاد عالية جدًا،2م{\displaystyle 2m}وهو الوضع المعتاد في الديناميكا الحرارية، حيث يقع كل الحجم تقريبًا بالقرب من السطح

Ω(هـ)=Vهـ=جمهـم-1Γ(م)،{\displaystyle \Omega (E)={\frac {\partial {\mathcal {V}}}{\partial E}}={\frac {C^{m}\cdot E^{m-1}}{\Gamma (m)}},}

حيث استخدمنا صيغة التكرارمΓ(م)=Γ(م+1){\displaystyle m\Gamma (m)=\Gamma (m+1)}.

مساحة السطحΩ(هـ){\displaystyle \Omega (E)}تستند هذه الكمية إلى عالمين: (أ) العالم العياني حيث تُعتبر دالةً للطاقة والمتغيرات الشاملة الأخرى، كالحجم، التي تُعتبر ثابتةً عند حساب حجم الطور، و(ب) العالم المجهري حيث تُمثل عدد التراكيب المتوافقة مع حالة عيانية معينة. هذه الكمية هي التي أشار إليها بلانك باسم الاحتمال "الديناميكي الحراري". وهي تختلف عن الاحتمال الكلاسيكي في أنها لا يمكن تطبيعها؛ أي أن تكاملها على جميع الطاقات يتباعد، ولكنه يتباعد كقوة للطاقة وليس بسرعة أكبر. ولأن تكاملها على جميع الطاقات لانهائي، يمكننا محاولة دراسة تحويل لابلاس الخاص بها.

Z(β)=0هـ-βهـΩ(هـ)دهـ،{\displaystyle {\mathcal {Z}}(\beta )=\int _{0}^{\infty }e^{-\beta E}\Omega (E)\,dE,}

والتي يمكن إعطاؤها تفسيراً فيزيائياً. عامل التناقص الأسي، حيثβ{\displaystyle \beta }إذا كانت قيمة المعامل موجبة، فسوف تتغلب على مساحة السطح المتزايدة بسرعة بحيث تتطور ذروة حادة للغاية عند طاقة معينةهـ{\displaystyle E^{\star }}سيأتي معظم المساهمة في التكامل من جوار مباشر حول هذه القيمة للطاقة. وهذا يُمكّن من تعريف دالة كثافة احتمالية مناسبة وفقًا لـ

و(هـ؛β)=هـ-βهـZ(β)Ω(هـ)،{\displaystyle f(E;\beta )={\frac {e^{-\beta E}}{{\mathcal {Z}}(\beta )}}\Omega (E),}

الذي يمثل تكامله على جميع الطاقات وحدةً بناءً على تعريفZ(β){\displaystyle {\mathcal {Z}}(\beta )}، والتي تُعرف بدالة التقسيم، أو الدالة المولدة. ويعود الاسم الأخير إلى حقيقة أن مشتقات لوغاريتمها تولد العزوم المركزية، وهي:

هـ=-lnZβ، (هـ-هـ)2=(Δهـ)2=2lnZβ2،{\displaystyle \langle E\rangle =-{\frac {\partial \ln {\mathcal {Z}}}{\partial \beta }},\qquad \ \langle (E-\langle E\rangle )^{2}\rangle =(\Delta E)^{2}={\frac {\partial ^{2}\ln {\mathcal {Z}}}{\partial \beta ^{2}}},}

وهكذا دواليك، حيث يمثل الحد الأول متوسط ​​الطاقة، ويمثل الحد الثاني تشتت الطاقة.

حقيقة أنΩ(هـ){\displaystyle \Omega (E)}لا يزيد معدل الزيادة عن قوة الطاقة، مما يضمن أن هذه اللحظات ستكون محدودة. [ 2 ] لذلك، يمكننا توسيع العاملهـ-βهـΩ(هـ){\displaystyle e^{-\beta E}\Omega (E)}حول القيمة المتوسطةهـ{\displaystyle \langle E\rangle }، وهو ما سيتزامن معهـ{\displaystyle E^{\star }}بالنسبة للتقلبات الغاوسية (أي أن القيم المتوسطة والأكثر احتمالاً تتطابق)، والاحتفاظ بالحدود ذات الرتبة الأدنى ينتج عنه

و(هـ؛β)=هـ-βهـZ(β)Ω(هـ)خبرة{-(هـ-هـ)2/2(Δهـ)2}2π(Δهـ)2.{\displaystyle f(E;\beta )={\frac {e^{-\beta E}}{{\mathcal {Z}}(\beta )}}\Omega (E)\approx {\frac {\exp\{-(E-\langle E\rangle )^{2}/2\langle (\Delta E)^{2}\rangle \}}{\sqrt {2\pi (\Delta E)^{2}}}}.}

هذا هو التوزيع الغاوسي، أو التوزيع الطبيعي، والذي يُعرَّف بأول عزمين له. وبشكل عام، يحتاج المرء إلى جميع العزوم لتحديد كثافة الاحتمال.و(هـ؛β){\displaystyle f(E;\beta )}، والتي يشار إليها باسم الكثافة القانونية أو الخلفية، على عكس الكثافة السابقةΩ{\displaystyle \Omega }، والتي يُشار إليها باسم دالة "البنية". [ 2 ] هذه هي نظرية النهاية المركزية كما تنطبق على الأنظمة الديناميكية الحرارية. [ 3 ]

إذا ازداد حجم الطور معهـم{\displaystyle E^{m}}، وسيتغير تحويل لابلاس الخاص به، دالة التقسيم، على النحو التاليβ-م{\displaystyle \beta ^{-m}}إعادة ترتيب التوزيع الطبيعي بحيث يصبح تعبيرًا عن دالة البنية وتقييمه عندهـ=هـ{\displaystyle E=\langle E\rangle }يعطي

Ω(هـ)=هـβ(هـ)هـZ(β(هـ))2π(Δهـ)2.{\displaystyle \Omega (\langle E\rangle )={\frac {e^{\beta (\langle E\rangle )\langle E\rangle }{\mathcal {Z}}(\beta (\langle E\rangle ))}{\sqrt {2\pi (\Delta E)^{2}}}}.}

It follows from the expression of the first moment that β(E)=m/E{\displaystyle \beta (\langle E\rangle )=m/\langle E\rangle }, while from the second central moment, (ΔE)2=E2/m{\displaystyle \langle (\Delta E)^{2}\rangle =\langle E\rangle ^{2}/m}. Introducing these two expressions into the expression of the structure function evaluated at the mean value of the energy leads to

Ω(E)=Em1m2πmmmem{\displaystyle \Omega (\langle E\rangle )={\frac {\langle E\rangle ^{m-1}m}{{\sqrt {2\pi m}}m^{m}e^{-m}}}}.

The denominator is exactly Stirling's approximation for m!=Γ(m+1){\displaystyle m!=\Gamma (m+1)}, and if the structure function retains the same functional dependency for all values of the energy, the canonical probability density,

f(E;β)=β(βE)m1Γ(m)eβE{\displaystyle f(E;\beta )=\beta {\frac {(\beta E)^{m-1}}{\Gamma (m)}}e^{-\beta E}}

will belong to the family of exponential distributions known as gamma densities. Consequently, the canonical probability density falls under the jurisdiction of the local law of large numbers which asserts that a sequence of independent and identically distributed random variables tends to the normal law as the sequence increases without limit.

Distribution about equilibrium

The expressions given below are for systems that are close to equilibrium and have negligible quantum effects.[4]

Single variable

Suppose x{\displaystyle x} is a thermodynamic variable. The probability distribution w(x)dx{\displaystyle w(x)dx} for x{\displaystyle x} is determined by the entropy S{\displaystyle S}:

w(x)exp(S(x)).{\displaystyle w(x)\propto \exp \left(S(x)\right).}

If the entropy is Taylor expanded about its maximum (corresponding to the equilibrium state), the lowest order term is a Gaussian distribution:

w(x)=12πx2exp(x22x2).{\displaystyle w(x)={\frac {1}{\sqrt {2\pi \langle x^{2}\rangle }}}\exp \left(-{\frac {x^{2}}{2\langle x^{2}\rangle }}\right).}

The quantity x2{\displaystyle \langle x^{2}\rangle } is the mean square fluctuation.[4]

Multiple variables

The above expression has a straightforward generalization to the probability distribution w(x1,x2,,xn)dx1dx2dxn{\displaystyle w(x_{1},x_{2},\ldots ,x_{n})dx_{1}dx_{2}\ldots dx_{n}}:

w=i,j=1n1(2π)n/2xixjexp(xixj2xixj),{\displaystyle w=\prod _{i,j=1\ldots n}{\frac {1}{\left(2\pi \right)^{n/2}{\sqrt {\langle x_{i}x_{j}\rangle }}}}\exp \left(-{\frac {x_{i}x_{j}}{2\langle x_{i}x_{j}\rangle }}\right),}

where xixj{\displaystyle \langle x_{i}x_{j}\rangle } is the mean value of xixj{\displaystyle x_{i}x_{j}}.[4]

Fluctuations of the fundamental thermodynamic quantities

In the table below are given the mean square fluctuations of the thermodynamic variables T,V,P{\displaystyle T,V,P} and S{\displaystyle S} in any small part of a body. The small part must still be large enough, however, to have negligible quantum effects.

Averages xixj{\displaystyle \langle x_{i}x_{j}\rangle } of thermodynamic fluctuations. CP{\displaystyle C_{P}} is the heat capacity at constant pressure; CV{\displaystyle C_{V}} is the heat capacity at constant volume.[4]
ΔT{\displaystyle \Delta T}ΔV{\displaystyle \Delta V}ΔS{\displaystyle \Delta S}ΔP{\displaystyle \Delta P}
ΔT{\displaystyle \Delta T}kBCVT2{\displaystyle {\frac {k_{\rm {B}}}{C_{V}}}T^{2}}0{\displaystyle 0}kBT{\displaystyle k_{\rm {B}}T}kBCVT2(PT)V{\displaystyle {\frac {k_{\rm {B}}}{C_{V}}}T^{2}\left({\frac {\partial P}{\partial T}}\right)_{V}}
ΔV{\displaystyle \Delta V}0{\displaystyle 0}kBT(VP)T{\displaystyle -k_{\rm {B}}T\left({\frac {\partial V}{\partial P}}\right)_{T}}kBT(VT)P{\displaystyle k_{\rm {B}}T\left({\frac {\partial V}{\partial T}}\right)_{P}}kBT{\displaystyle -k_{\rm {B}}T}
ΔS{\displaystyle \Delta S}kBT{\displaystyle k_{\rm {B}}T}kBT(VT)P{\displaystyle k_{\rm {B}}T\left({\frac {\partial V}{\partial T}}\right)_{P}}kBCP{\displaystyle k_{\rm {B}}C_{P}}0{\displaystyle 0}
ΔP{\displaystyle \Delta P}kBCVT2(PT)V{\displaystyle {\frac {k_{\rm {B}}}{C_{V}}}T^{2}\left({\frac {\partial P}{\partial T}}\right)_{V}}kBT{\displaystyle -k_{\rm {B}}T}0{\displaystyle 0}kBT(PV)S{\displaystyle -k_{\rm {B}}T\left({\frac {\partial P}{\partial V}}\right)_{S}}

See also

Notes

  1. In statistical mechanics they are often simply referred to as fluctuations.
  2. 12Khinchin 1949
  3. Lavenda 1991
  4. 1234Landau & Lifshitz 1985

References