Différences
Ci-dessous, les différences entre deux révisions de la page.
Les deux révisions précédentes Révision précédente Prochaine révision | Révision précédente | ||
teaching:progappchim:matplotlib_gallery:potentiel_energy_surface [2016/05/10 13:00] – villersd | teaching: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> |
# | # | ||
# -*- coding: utf-8 -*- | # -*- coding: utf-8 -*- | ||
Ligne 69: | Ligne 69: | ||
X, Y = np.meshgrid(x_1d, | X, Y = np.meshgrid(x_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.)))/ | + | 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) | print np.min(E) | ||
Ligne 86: | Ligne 86: | ||
plt.ylabel(' | plt.ylabel(' | ||
- | fig = plt.figure(2, | + | fig = plt.figure(2, |
ax = Axes3D(fig) | ax = Axes3D(fig) | ||
ax.plot_surface(X, | ax.plot_surface(X, | ||
Ligne 93: | Ligne 93: | ||
ax.set_zlabel(' | ax.set_zlabel(' | ||
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 | ||
Ligne 117: | Ligne 117: | ||
* [[http:// | * [[http:// | ||
* [[https:// | * [[https:// | ||
+ | * W H Miller, Recent Advances in Quantum Mechanical Reactive Scattering Theory, Including Comparison of Recent Experiments with Rigorous Calculations of State-to-State Cross Sections for the H/ | ||
+ | * Peterson, Kirk A., Woon, David E., Dunning, Thom H. , Benchmark calculations with correlated molecular wave functions. IV. The classical barrier height of the H+H2→H2+H reaction, | ||
+ | |||
+ | |||
Voir aussi : | Voir aussi : | ||
* [[http:// | * [[http:// | ||
* [[http:// | * [[http:// | ||
+ | * [[http:// | ||