teaching:progappchim:lennard-jones

Différences

Ci-dessous, les différences entre deux révisions de la page.

Lien vers cette vue comparative

Prochaine révision
Révision précédente
Dernière révisionLes deux révisions suivantes
teaching:progappchim:lennard-jones [2014/01/29 12:35] – créée villersdteaching:progappchim:lennard-jones [2015/03/05 11:08] villersd
Ligne 4: Ligne 4:
  
 Exemple du [[http://fr.wikipedia.org/wiki/Potentiel_de_Lennard-Jones|potentiel de Lennard-Jones]] de l'argon : Exemple du [[http://fr.wikipedia.org/wiki/Potentiel_de_Lennard-Jones|potentiel de Lennard-Jones]] de l'argon :
 +
 +$V_{LJ} = 4\varepsilon \left[ \left(\frac{\sigma}{r}\right)^{12} - \left(\frac{\sigma}{r}\right)^{6} \right] = \varepsilon \left[ \left(\frac{r_{m}}{r}\right)^{12} - 2\left(\frac{r_{m}}{r}\right)^{6} \right]$
 +
 +où σ est la distance à laquelle le potentiel entre les particules s'annule et ε est l'énergie du puits de potentiel d'interaction. La distance r<sub>m</sub> à laquelle le potentiel a cette valeur minimale (pour une interaction entre seulement deux atomes) est reliée à σ par la relation suivante : $r_{m} = 2^{1/6} \sigma$ 
  
 <sxh python; title : Lennard-Jones-01.py> <sxh python; title : Lennard-Jones-01.py>
Ligne 33: Ligne 37:
 show() show()
 </sxh> </sxh>
 +
 +<note tip>Suggestion : récrire ce programme en utilisant des directives d'importation standard des librairies Matplotlib/NumPy</note>
 +
 +===== Application : forces de cohésion dans les cristaux de gaz rares =====
 +
 +À basse température, les gaz rares peuvent donner des cristaux de structure cubique à face centrée ou hexagonale compacte. La phase cubique centrée peut également être investiguée. Dans chaque de ces structures, chaque atome possède alors d'autres atomes dans son voisinage immédiat (la coordination 12 ou 8), mais aussi au delà. Pour caractériser l'énergie de cohésion et relier le paramètre du réseau cristallin à σ, il faut considérer la somme de toutes ces interactions, ou du moins de celles correspondant à un voisinage suffisant pour converger. En sommant pour N atomes, on obtient :
 +
 +$U_{tot} = \frac{1}{2} N (4\epsilon) \left[ \sum_j \left( \frac{\sigma} {p_j R}\right)^{12} - \sum_j \left( \frac{\sigma} {p_j r}\right)^{6} \right]$
 +
 +où R est la distance de séparation entre plus proches voisins et $p_j$ représente les distances non dimensionnelles (réduites par rapport à R) entre les paires possibles d'atomes par rapport à l'atome de référence.
 +
 +La distance d'équilibre $R_0$ pourra être obtenue facilement en minimisant $U_tot$. Cela nécessite de calculer pour le réseau envisagé les valeurs de $\sum_j p_j^{-12}$ et $\sum_j p_j^{-6}$.
 +
 +==== Références ====
 +  * Kittel, physique de l'état solide, chapitre 3
 +  * [[http://www.cmmp.ucl.ac.uk/~ikr/3225/Section%203.pdf]]
 +  * [[http://th.fhi-berlin.mpg.de/th/lectures/materialscience-2004/vorlesung_2004/fest-k6-2004.pdf]]
 +  * Données sur les gaz, σ du potentiel de Lennard-Jones ([[http://www.hindawi.com/journals/jther/2013/828620/]]) :
 +    * He : 262.8 pm
 +    * Ne : 277.5 pm
 +    * Ar : 340.1 pm
 +    * Kr : 360.1 pm
 +    * Xe : 405.5 pm
 +  * Données cristaux (cubic close packed ou face centered, [[http://www.webelements.com/]]
 +    * He : 424.2 pm
 +    * Ne : 442.9 pm
 +    * Ar : 525.6 pm
 +    * Kr : 570.6 pm
 +    * Xe : 620.23 pm
 +
 +
 +
  • teaching/progappchim/lennard-jones.txt
  • Dernière modification : 2015/03/05 12:14
  • de villersd