teaching:progappchim:matplotlib_gallery:potentiel_energy_surface

Différences

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

Lien vers cette vue comparative

Les deux révisions précédentes Révision précédente
teaching:progappchim:matplotlib_gallery:potentiel_energy_surface [2016/12/14 11:57] – [Références] villersdteaching:progappchim:matplotlib_gallery:potentiel_energy_surface [2020/12/07 17:22] (Version actuelle) – [Programme] villersd
Ligne 32: Ligne 32:
 ==== Programme ==== ==== Programme ====
  
-<sxh python; title : PES-contour-01.py>+<code python PES-contour-01.py>
 #!/usr/bin/env python #!/usr/bin/env python
 # -*- coding: utf-8 -*- # -*- coding: utf-8 -*-
Ligne 69: Ligne 69:
 X, Y = np.meshgrid(x_1d, y_1d) X, Y = np.meshgrid(x_1d, y_1d)
 print X.shape, X.dtype, X.ndim, Y.shape, Y.dtype, Y.ndim print X.shape, X.dtype, X.ndim, Y.shape, Y.dtype, Y.ndim
-E=(Q(X)+Q(Y)+Q(X+Y)-np.sqrt(2.*((a(X)-a(Y))**2.+(a(Y)-a(X+Y))**2.+(a(X+Y)-a(X))**2.)))/(1.+k)+E=(Q(X)+Q(Y)+Q(X+Y)-np.sqrt(2.*( (a(X)-a(Y) )**2.+(a(Y)-a(X+Y) )**2.+(a(X+Y)-a(X) )**2.) ))/(1.+k)
 print np.min(E)  #valeur minimale de E print np.min(E)  #valeur minimale de E
  
Ligne 86: Ligne 86:
 plt.ylabel('y (pm)') plt.ylabel('y (pm)')
  
-fig = plt.figure(2,figsize=(15, 15))+fig = plt.figure(2,figsize=(15, 15) )
 ax = Axes3D(fig) ax = Axes3D(fig)
 ax.plot_surface(X,Y,E, rstride=1,cstride=1 ,cmap=plt.cm.jet) ax.plot_surface(X,Y,E, rstride=1,cstride=1 ,cmap=plt.cm.jet)
Ligne 93: Ligne 93:
 ax.set_zlabel('E') ax.set_zlabel('E')
 plt.show() plt.show()
-</sxh>+</code>
  
 Avec les paramètres essayés, la valeur minimale de E est environ -1.603 Avec les paramètres essayés, la valeur minimale de E est environ -1.603
  • teaching/progappchim/matplotlib_gallery/potentiel_energy_surface.txt
  • Dernière modification : 2020/12/07 17:22
  • de villersd