Izotermik-izobarik ansambl - Isothermal–isobaric ensemble

The izotermik-izobarik ansambl (doimiy harorat va doimiy bosim ansambli) bu a statistik mexanik ansambl doimiy haroratni saqlaydi va doimiy bosim qo'llaniladi. U shuningdek - zarralar soni bo'lgan ansambl doimiy sifatida ham saqlanadi. Ushbu ansambl kimyoda muhim rol o'ynaydi, chunki kimyoviy reaktsiyalar odatda doimiy bosim sharoitida amalga oshiriladi.[1] NPT ansambli model tizimlarining holatini tenglamasini o'lchash uchun ham foydalidir virusli kengayish chunki bosimni baholash mumkin emas yoki birinchi darajali o'zgarishlar o'tishlariga yaqin tizimlar.[2]

Asosiy xususiyatlarni ishlab chiqarish

Uchun bo'lim funktsiyasi - ansamblni tizimidan boshlab statistik mexanikadan olish mumkin a tomonidan tavsiflangan bir xil atomlar Hamiltoniyalik shaklning va bir quti ichida mavjud . Ushbu tizim. Ning bo'linish funktsiyasi bilan tavsiflanadi kanonik ansambl 3 o'lchamda:

,

qayerda , termal de Broyl to'lqin uzunligi ( va bo'ladi Boltsman doimiy ) va omil (bu zarrachalarni ajratib bo'lmaydiganligini hisobga oladi) ikkalasi ham kvaz-klassik chegarada entropiyaning normallashishini ta'minlaydi.[2] Tomonidan belgilangan yangi koordinatalar to'plamini qabul qilish qulay bo'linish funktsiyasi shunday bo'ladi

.

Agar ushbu tizim keyinchalik hajmli hammom bilan aloqa qilsa o'z ichiga olgan doimiy harorat va bosimda ideal gaz zarrachalarning umumiy soni bilan shu kabi , butun tizimning bo'linish funktsiyasi shunchaki quyi tizimlarning bo'lim funktsiyalari mahsulidir:

.
Tizim (hajmi ) doimiy haroratdagi ancha kattaroq vannaga botiriladi va zarrachalar soni sobit qolishi uchun yopiladi. Tizim vannadan erkin harakatlanadigan piston bilan ajralib turadi, shunda uning hajmi o'zgarishi mumkin.

Ning ajralmas qismi koordinatalari sodda . Bu chegarada , esa doimiy bo'lib qoladi, o'rganilayotgan tizim hajmining o'zgarishi bosimni o'zgartirmaydi butun tizim. Qabul qilish yaqinlashtirishga imkon beradi . Ideal gaz uchun, zichlik va bosim o'rtasidagi bog'liqlikni beradi. Buni qism funktsiyasi uchun yuqoridagi ifodaga almashtirish, faktorga ko'paytirish (ushbu qadamni asoslash uchun quyida ko'rib chiqing) va V hajmi bo'yicha integratsiya keyin beradi

.

Hammom uchun bo'lim vazifasi oddiygina . Ushbu atamani umumiy ifodadan ajratish uchun - ansambl:

.

Ning yuqoridagi ta'rifidan foydalanib , bo'lim funktsiyasini quyidagicha yozish mumkin

,

Bu odatda kanonik ansambl uchun ajratilgan funktsiya bo'yicha tortilgan summa sifatida yozilishi mumkin

Miqdor integralni hosil qilish uchun zarur bo'lgan teskari hajm birliklari bilan oddiygina bir necha doimiydir o'lchovsiz. Ushbu holatda, , lekin umuman olganda u bir nechta qiymatlarni qabul qilishi mumkin. Uning tanlovidagi noaniqlik hajmni hisoblash mumkin bo'lgan miqdor emasligidan kelib chiqadi (masalan, zarralar sonidan farqli o'laroq) va shuning uchun yuqoridagi lotinlashda bajarilgan yakuniy hajm integratsiyasi uchun "tabiiy o'lchov" mavjud emas.[2] Ushbu muammo turli xil mualliflar tomonidan bir necha bor hal qilingan,[3][4] bir xil teskari hajm birliklari bilan S uchun qiymatlarga olib keladi. Farqlar yo'qoladi (ya'ni tanlovi o'zboshimchalikga aylanadi) da termodinamik chegara, bu erda zarralar soni abadiylikka boradi.[5]

The -samblni Gibbs kanonik ansamblining alohida ishi sifatida ham ko'rish mumkin, unda makrostatlar tizimning harorati tashqi haroratga qarab belgilanadi va tizimga ta'sir qiluvchi tashqi kuchlar . O'z ichiga olgan bunday tizimni ko'rib chiqing zarralar. Tizimning Gamiltoniani keyin beriladi qayerda tashqi kuchlar mavjud bo'lmagan taqdirda tizimning Hamiltonianidir ular konjuge o'zgaruvchilar ning . Mikrostatlar tizimning qiymati keyin aniqlangan ehtimollik bilan sodir bo'ladi [6]

bu erda normalizatsiya omili bilan belgilanadi

.

The - ansamblni qabul qilish orqali topish mumkin va . Keyin normalizatsiya omili bo'ladi

,

bu erda Hamiltonian zarralar momentumi bo'yicha yozilgan va lavozimlar . Ushbu summani ikkalasining ham integraliga etkazish mumkin va mikrostatlar . Oxirgi integral uchun o'lchov standart o'lchovdir fazaviy bo'shliq bir xil zarralar uchun: .[6] Integral tugadi muddatli a Gauss integrali, va aniq tarzda baholanishi mumkin

.

Ushbu natijani kiritish tanish ifoda beradi:

.[6]

Bu deyarli uchun funktsiya - ansambl, lekin u hajm birliklariga ega, bu yuqoridagi summani ajralmas qismga aylantirishning muqarrar natijasidir. Doimiylikni tiklash uchun tegishli natijani beradi .

Oldingi tahlillardan ko'rinib turibdiki, ushbu ansamblning o'ziga xos holat vazifasi Gibbs bepul energiya,

Bu termodinamik potentsial bilan bog'liq Helmholtsning erkin energiyasi (kanonik bo'lim funktsiyasining logarifmi), , quyidagi tarzda:[1]

Ilovalar

  • Doimiy bosim simulyatsiyalari aniqlash uchun foydalidir davlat tenglamasi sof tizim. Monte-Karlo simulyatsiyalari - ansambl, ayniqsa, taxminan 1 atm bosimdagi suyuqlik holatining tenglamasini aniqlash uchun foydalidir, bu erda ular boshqa ansambllarga qaraganda ancha kam hisoblash vaqti bilan aniq natijalarga erishishlari mumkin.[2]
  • Nolinchi bosim - ansambl simulyatsiyalari aralash fazali tizimlarda bug 'va suyuqlikning birgalikda yashash egri chiziqlarini baholashning tezkor usulini taqdim etadi.[2]
  • Monte-Karlo ansambli simulyatsiyasi qo'llanilgan ortiqcha xususiyatlar [7] va holat tenglamalari [8] suyuqlik aralashmalarining turli xil modellari.
  • The ansambl ham foydalidir molekulyar dinamikasi simulyatsiyalar, masalan. atrof-muhit sharoitida suvning xatti-harakatlarini modellashtirish.[9]

Adabiyotlar

  1. ^ a b Dereotu, Ken A .; Bromberg, Sarina; Stigter, Dirk (2003). Molekulyar haydash kuchlari. Nyu York: Garland fani.
  2. ^ a b v d e Frenkel, Daan.; Smit, Berend (2002). Molekulyar simulyatsiya haqida tushuncha. Nyu York: Akademik matbuot.
  3. ^ Attard, Fil (1995). "Izobarik ansambldagi tovush holatlarining zichligi to'g'risida". Kimyoviy fizika jurnali. 103 (24): 9884–9885. doi:10.1063/1.469956.
  4. ^ Koper, Ger J. M.; Reys, Xovard (1996). "Doimiy bosim ansambli uchun uzunlik ko'lami: kichik tizimlarga qo'llanilishi va Eynshteynning tebranish nazariyasiga aloqasi". Jismoniy kimyo jurnali. 100 (1): 422–432. doi:10.1021 / jp951819f.
  5. ^ Tepalik, Terrens (1987). Statistik mexanika: tamoyillar va tanlangan qo'llanmalar. Nyu York: Dover.
  6. ^ a b v Kardar, Mehran (2007). Zarrachalarning statistik fizikasi. Nyu York: Kembrij universiteti matbuoti.
  7. ^ McDonald, I. R. (1972). "- ikkilik suyuqlik aralashmalari uchun Monte-Karlo hisob-kitoblarini yig'ing ". Molekulyar fizika. 23 (1): 41–58. doi:10.1080/00268977200100031.
  8. ^ Wood, W. W. (1970). "- Qattiq disk suyuqligi uchun Monte-Karlo hisob-kitoblari. " Kimyoviy fizika jurnali. 52 (2): 729–741. doi:10.1063/1.1673047.
  9. ^ Shmidt, Yoxen; VandeVondele, Joost; Kuo, I. F. Uilyam; Sebastiani, Doniyor; Siepmann, J. Ilja; Xutter, Yurg; Muni, Kristofer J. (2009). "Zichlik funktsional nazariyasidan foydalangan holda izobarik-izotermik molekulyar dinamikaning simulyatsiyalari: atrof-muhit sharoitida suvning tuzilishi va zichligini baholash". Jismoniy kimyo jurnali B. 113 (35): 11959–11964. doi:10.1021 / jp901990u.