Potențialul Lennard-Jones ( potențialul 6-12 ) este un model simplu de interacțiune perechi a moleculelor nepolare care descrie dependența energiei de interacțiune a două particule de distanța dintre ele. Acest model transmite în mod destul de realist proprietățile interacțiunii reale a moleculelor sferice nepolare și, prin urmare, este utilizat pe scară largă în calcule și simulări pe computer. Acest tip de potențial a fost propus pentru prima dată de Lennard-Jones în 1924. [unu]
Potențialul Lennard-Jones este scris după cum urmează:
unde este distanța dintre centrele particulelor, este adâncimea puțului de potențial, este distanța la care energia de interacțiune devine egală cu zero. Parametrii și sunt caracteristicile atomilor substanței corespunzătoare. Forma caracteristică a potențialului este prezentată în figură, minimul său se află în punctul .
La mari molecule sunt atrase, ceea ce corespunde unui termen din formulă. Această dependență poate fi fundamentată teoretic și se datorează forțelor van der Waals (interacțiunea indusă dipol-dipol).
La distanțe scurte, moleculele se resping unele pe altele datorită interacțiunii de schimb (atunci când norii de electroni se suprapun, moleculele încep să se respingă puternic), ceea ce corespunde termenului . Acest tip particular de potențial respingător, în contrast cu tipul de potențial de atracție, nu are o justificare teoretică. Mai rezonabilă este dependența exponențială . Cu toate acestea, potențialul de repulsie Lennard-Jones este mai convenabil în calcule, deoarece , ceea ce justifică utilizarea acestuia.
Potențialul Lennard-Jones este adesea scris în următoarea formă, cea mai simplă:
Unde
Există și această formă de scriere:
unde = este punctul minim al potențialului.
Modelul Lennard-Jones poate fi folosit pentru a descrie fazele gazoase, lichide și solide ale materiei . Cea mai mică valoare a energiei libere pentru o substanță condiționată, pentru care modelul Lennard-Jones este valabil, este obținută cu împachetare strânsă hexagonală . Pe măsură ce temperatura crește, structura cu cea mai scăzută energie liberă se schimbă într-o împachetare strânsă centrată pe față cubică și apoi se observă o tranziție la un lichid. Sub presiune, pentru structura cu cea mai scăzută energie, există o tranziție de la împachetarea cubică la împachetare hexagonală. [2]
Punctul critic pentru substanța condiționată considerată în variabile adimensionale ( ) se află la următoarele valori ale temperaturii și concentrației: [3]
Fliegenhart și Lekkerkerker au propus următoarea expresie pentru relația dintre punctul critic și al doilea coeficient virial: [4]
Poziția punctului triplu a fost stabilită de Mastney și Pablo: [5]
(lichid); (solid).Pentru a accelera calculele, potențialul Lennard-Jones este adesea întrerupt la distanță . Alegerea se datorează faptului că la această distanță valoarea energiei de interacțiune este doar ≈0,0163 din adâncimea sondei .
Cu toate acestea, reducerea potențialului în acest fel nu este întotdeauna convenabilă. Și anume, o astfel de rupere înseamnă că atunci când o moleculă traversează o sferă de rază , energia sistemului se modifică brusc sau, ceea ce este același, o forță infinit de mare acționează asupra moleculei. Pentru a evita această situație non-fizică, atunci când potențialul se sparge, acesta este și deplasat, astfel încât :
este valoarea potențialului Lennard-Jones neîntrerupt la distanță .
O altă modalitate de a accelera calculele este utilizarea spline . În acest caz, potențialul de interacțiune este împărțit în mai multe secțiuni, pe fiecare dintre acestea fiind aproximat printr-o funcție simplă. Următoarea aproximare este adesea folosită: [6]
aici .
Uneori, potențialul Lennard-Jones înseamnă forma sa mai generală, și anume:
unde . Coeficientul este ales astfel încât valoarea minimă a energiei potențiale să fie [ 7] .
Parametrii pot fi găsiți în două moduri.
O metodă se bazează pe compararea valorilor experimentale ale celui de-al doilea coeficient virial al funcției stărilor unui gaz neideal [8] pentru energia potențială a configurațiilor cu valorile obținute cu ajutorul parametrilor , [9] . Valorile , găsite în acest fel sunt în bună concordanță cu valorile găsite prin coeficientul Joule-Thomson .
Conform unei alte metode, parametrii , pot fi aflați prin compararea valorii experimentale a coeficientului de vâscozitate cu valoarea obținută din formula energiei potențiale [9] .
Tabelele cu valorile parametrilor pentru unele substanțe sunt în carte [10] .
Deși potențialul Lennard-Jones este utilizat în modelarea lichidelor și solidelor, strict vorbind, interacțiunea moleculelor la densități mari nu mai este o interacțiune de pereche. În mediile condensate, perechea de molecule luate în considerare este influențată de moleculele mediului. Așa că s-a constatat că pentru argonul solid , contribuția la energie din interacțiuni triple poate ajunge la 10%. [11] Cu toate acestea, contabilizarea interacțiunilor triple este prea costisitoare din punct de vedere computațional, așa că de obicei se mulțumește cu un potențial efectiv de pereche, unde parametrii și diferă de cei pentru gazele rarefiate.