root / examples / Gaussian / SN2_zmat / SN2_zmat.path @ 1
Historique | Voir | Annoter | Télécharger (1,89 ko)
| 1 | 1 | equemene | &path |
|---|---|---|---|
| 2 | 1 | equemene | nat=6, ! Number of atoms |
| 3 | 1 | equemene | ngeomi=2, ! Number of initial geometries |
| 4 | 1 | equemene | ngeomf=18, !Number of geometries along the path |
| 5 | 1 | equemene | OptReac=.F., ! Do you want to optimize the reactants ? |
| 6 | 1 | equemene | OptProd=F, ! Don't optimize the products |
| 7 | 1 | equemene | Align=T, ! Do not align the molecules when interpolating geometries |
| 8 | 1 | equemene | EReac=-0.21076866, ! Energy of the reactant, because it is not calculated |
| 9 | 1 | equemene | EProd=-0.21076866, ! As we don't optimize products, we can give their energy |
| 10 | 1 | equemene | FFrozen=F, ! Some atoms are frozen. See FrozenList. |
| 11 | 1 | equemene | PathOnly=.F., ! do NOT stop after generating the first path |
| 12 | 1 | equemene | AutoCart=T, ! Determine automatically atom described in cartesian. |
| 13 | 1 | equemene | HesUpd="MS", ! Use Murtagh-Sargent update for the Hessian matrix |
| 14 | 1 | equemene | PathName="Test_Gaussian", ! Geometries along the path store in 'Test_Gaussian.XXX' |
| 15 | 1 | equemene | prog="gaussian", ! Use Gaussian to get E and forces. |
| 16 | 1 | equemene | ProgExe="g09", ! Use g09 for Gaussian executable |
| 17 | 1 | equemene | coord='zmat', ! we use mixed coordinates |
| 18 | 1 | equemene | maxcyc=51, ! Launch 50 iterations |
| 19 | 1 | equemene | IReparam=5, ! re-distribution of points along the path every 5 iterations |
| 20 | 1 | equemene | ISpline=500, ! Start using spline interpolation at iteration 500 (ie no spline !) |
| 21 | 1 | equemene | SMax=0.02 ! Max displacement will be 0.02 a.u. |
| 22 | 1 | equemene | / |
| 23 | 1 | equemene | 6 |
| 24 | 1 | equemene | Cl- + H3CCl |
| 25 | 1 | equemene | C 0.0 0.0 0.0 |
| 26 | 1 | equemene | Cl 1.791689 0.0 0.0 |
| 27 | 1 | equemene | H -0.358375 1.037571 0.0 |
| 28 | 1 | equemene | H -0.356525 -0.518921 -0.898137 |
| 29 | 1 | equemene | H -0.360734 -0.518109 0.897921 |
| 30 | 1 | equemene | Cl -6.0 0.001 0. |
| 31 | 1 | equemene | 6 |
| 32 | 1 | equemene | ClCH3 + Cl- |
| 33 | 1 | equemene | C 0.0 0.0 0.0 |
| 34 | 1 | equemene | Cl 6.0 0.001 0. |
| 35 | 1 | equemene | H 0.358375 1.037571 0.0 |
| 36 | 1 | equemene | H 0.356525 -0.518921 -0.898137 |
| 37 | 1 | equemene | H 0.360734 -0.518109 0.897921 |
| 38 | 1 | equemene | Cl -1.791689 0.0 0.0 |
| 39 | 1 | equemene | %chk=Test |
| 40 | 1 | equemene | #P HF/3-21G FORCE |
| 41 | 1 | equemene | |
| 42 | 1 | equemene | SN2 Reaction |
| 43 | 1 | equemene | |
| 44 | 1 | equemene | -1,1 |
| 45 | 1 | equemene | C 0.0 0.0 0.0 |
| 46 | 1 | equemene | Cl 1.791689 0.0 0.0 |
| 47 | 1 | equemene | H -0.358375 1.037571 0.0 |
| 48 | 1 | equemene | H -0.356525 -0.518921 -0.898137 |
| 49 | 1 | equemene | H -0.360734 -0.518109 0.897921 |
| 50 | 1 | equemene | Cl -6.0 0. 0. |
| 51 | 1 | equemene |