19/3/2020 Dinamica molecolare - Lennard-Jones
Résumé de section
-
https://meet.google.com/oth-hqpz-asd
Drive:
https://drive.google.com/open?id=1nTSH-d8P4cidmIVaSC5EEK8iEE5PLz0
scaricamento:
https://drive.google.com/open?id=1W-dpZbJ7Cif4tuoBZmrq9XeG6YA7Wxa7
Avrete visto che quando si inizializza il sistema a volte l'energia esplode. il problema è legato a thre e deltat. Prendiamo thre =10^-3. Al limite può siccedere che due pasrticelle siano così vicine. Quindi la forza che è circa 48*r^-13 per r=thre=10^-3 diventa 48*10^33. Dto che la massa è 1 questa è anche l'accelerazione. Supponiamo che la velocità sia zero. Nel tempo deltat =10^-2 la velocità diventa 48^10^31, e in un tempo deltat può benissimo arrivare molto più vicino ad un'altra particella di thre, dopodiché prende una botta ancora peggiore. Mi sa che thre dev'essere al più 0.1.