Discussion :

[flxpez]

Bonjour,
nous avons un TP à l'école où l'on doit calculer l'énergie d'un gaz (atomique) par la méthode de Montecarlo.
Voici le concept "en gros":
-On travail dans un espace canonique
-On place aléatoirement des atomes de gaz dans la boite
-On calcul l'énergie du système E1
-Aléatoirement on sélectionne un atome et il bouge
-On calcul l'énergie nouvellement définie du système E2
-Si E2-E1<0 on accepte le nouveau système
-Si E2-E1>0 on accepte le nouveau système avec la probabilité exp(-(E2-E1)/kbT)
-On répète ça x fois, le but étant d'avoir une énergie moyenne qui converge

L'énergie du système est pour nous la somme des potentiels de Lennard-Jones.

Notre problème est que l'on obtient des résultats incohérents. Si on met un pas trop élevé, il y a des sauts d'énergie qui devraient être refusés. On a codé deux programmes l'un deux utilisant deux tableaux pour les coordonnées x et y, et l'énergie baisse très bas puis remonte ce qui ne devrait pas être le cas. De plus l'énergie finale obtenue est négative alors qu'elle devrait être positive.

On essaye de débuguer le programme depuis hier déjà et on est a court d'idées.

Merci !

[ternel]

Bonjour,

Quelle est la formule mathématique du potentiel?
attention, dans selection, rand() ne doit pas être utilisé avec un modulo, la sélection risque d'être biaisée.
il faut faire comme dans move_atom.
Essaie donc:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
 
int selection() {
    return ( ((double)rand())/RAND_MAX)*ATOM_MAX;
}

Dans move_atom, la validation de la boite est mal faite pour les cas x ou y < 0
0 qui devrait être XMIN et YMIN

[flxpez]

Bonjour,

merci pour les indications, en fait le prof s'était trompé et la valeur attendue était bien négative.
On a cependant pu améliorer quelques petites choses, dont le MIN.

à bientôt,

[ternel]

Super.

Comme ca, tu te souviendras: si un résultat est faux, ca peut venir du calcul, des données, ou de la manière de juger le résultat.

N'hésite pas à passer la discussion comme résolue