CP2K Open Source Molecular Dynamics

User Tools

Site Tools

Trace:
exercises:2014_ethz_mmm:reaction_energy

Differences

This shows you the differences between two versions of the page.

Link to this comparison view

Both sides previous revisionPrevious revision
Next revision		Previous revision
exercise:reaction_energy [2014/03/26 16:04] – scleliaexercises:2014_ethz_mmm:reaction_energy [2020/08/21 10:15] (current) – external edit 127.0.0.1
Line 1: Line 1:
 ====== Reaction Energy ====== ====== Reaction Energy ======
-In this exercise, you will calculate the reaction energy for the mehane combustion reaction: +In this exercise, you will calculate the reaction energy for the methane combustion reaction: 
-{{ :exercise:f1.jpg?300 |}}+\[ CH_4 + 2O_2 \rightarrow CO_2 + 2H_2O \]
  
 Reaction energy: Reaction energy:
-{{ :exercise:f2.jpg?500 |}}+\[ \sum E_\text{products} -  \sum E_\text{rectants} = \left (2\cdot E_{H_2O} + E_{CO_2} \right) - \left(E_{CH_4} + 2\cdot E_{O_2}\right) \] 
  
 Ground state oxygen, O2, is a triplet diradical, a property which can explain why liquid oxygen is paramagnetic and attracted to the poles of a magnet.  Ground state oxygen, O2, is a triplet diradical, a property which can explain why liquid oxygen is paramagnetic and attracted to the poles of a magnet. 
-{{ :exercise:o3.png?600 |+{{ o3.png?600 |
 }} }}
  
Line 15: Line 16:
  
 Run a single point calculation for CH4, using the given input file.  Run a single point calculation for CH4, using the given input file. 
-Note that the file contains explicit basis sets and potential for all-electron calculations. An explanation of the basis set formats is given here: [[exercise:basis_sets|Basis Sets]]+Note that the file contains explicit basis sets and potential for all-electron calculations. An explanation of the basis set formats is given here: [[basis_sets|Basis Sets]]
  
 <code - CH4.inp > <code - CH4.inp >
  
 &GLOBAL &GLOBAL
-  PROJECT CH4_simple4+  PROJECT CH4
   RUN_TYPE ENERGY   RUN_TYPE ENERGY
   PRINT_LEVEL MEDIUM   PRINT_LEVEL MEDIUM
Line 28: Line 29:
   METHOD Quickstep              ! Electronic structure method (DFT,...)   METHOD Quickstep              ! Electronic structure method (DFT,...)
   &DFT   &DFT
-  BASIS_SET_FILE_NAME EMSL_BASIS_SETS 
-  POTENTIAL_FILE_NAME POTENTIAL 
     &POISSON                    ! Solver requested for non periodic calculations     &POISSON                    ! Solver requested for non periodic calculations
       PERIODIC NONE       PERIODIC NONE
Line 52: Line 51:
     &CELL     &CELL
       ABC 10 10 10       ABC 10 10 10
-      PERIODIC NONE              ! Non periodic calculations. That's why the POISSON scetion is needed +      PERIODIC NONE              ! Non periodic calculations. That's why the POISSON section is needed 
     &END CELL     &END CELL
     &TOPOLOGY                    ! Section used to center the atomic coordinates in the given box. Useful for big molecules     &TOPOLOGY                    ! Section used to center the atomic coordinates in the given box. Useful for big molecules
Line 59: Line 58:
     &END     &END
     &COORD     &COORD
 +  C         4.6425962273        5.0574874650        5.2069537560 
 +  H         5.7240587065        5.0555482951        5.2189766147 
 +  H         4.2766068912        5.8773176685        5.8100567767 
 +  H         4.2759350196        4.1226994019        5.6087492584 
 +  H         4.2938562590        5.1744089096        4.1899119266
     &END COORD     &END COORD
-    &KIND H +    &KIND H                                      ! potential and basis for H  
-    POTENTIAL ALL +      &BASIS    
-    BASIS_SET pc-1+        3 
 +  1  0  0  3  1 
 +         12.25200000          0.02282200 
 +          1.86870000          0.15564000 
 +          0.41821000          0.48898000 
 +  1  0  0  1  1 
 +          0.10610000          1.00000000 
 +  1  1  1  1  1 
 +          1.00000000          1.00000000 
 +      &END                       
 +      POTENTIAL ALL 
 +      &POTENTIAL 
 +       1    0    0 
 +       0.20000000    0 
 +      &END 
     &END KIND     &END KIND
-    &KIND C +    &KIND C                                      ! potential and basis for C 
-    POTENTIAL ALL +      &BASIS 
-    BASIS_SET pc-1+         5 
 +  1  0  0  6  2 
 +       1252.60000000          0.00557360          0.00000000 
 +        188.57000000          0.04149600         -0.00027440 
 +         42.83900000          0.18263000         -0.00255830 
 +         11.81800000          0.46129000         -0.03337500 
 +          3.55670000          0.44931000         -0.08730500 
 +          0.54258000          0.00000000          0.53415000 
 +  1  0  0  1  1 
 +          0.16058000          1.00000000 
 +  1  1  1  3  1 
 +          9.14260000          0.04449900 
 +          1.92980000          0.23108000 
 +          0.52522000          0.51227000 
 +  1  1  1  1  1 
 +          0.13608000          1.00000000 
 +  1  2  2  1  1 
 +          0.80000000          1.00000000 
 +      &END 
 +      POTENTIAL ALL 
 +      &POTENTIAL 
 +      4    2    0 
 +      0.34883045    0 
 +      &END
     &END KIND     &END KIND
  
Line 109: Line 149:
   * CO2   * CO2
   * O2 TRIPLET    * O2 TRIPLET 
-<note tip> Atomic coordinates for all the moelcules, POTENTIAL and BASIS SET for KIND O are given at the end of the exercise. </note>+<note tip> Atomic coordinates for all the molecules, POTENTIAL and BASIS SET for KIND O are given at the end of the exercise. </note>
 <note important>Remember that the O2 ground state is a triplet state, with non paired electrons. <note important>Remember that the O2 ground state is a triplet state, with non paired electrons.
  MULTIPLICITY=2S+1=3. </note>  MULTIPLICITY=2S+1=3. </note>
Line 125: Line 165:
 </code> </code>
  
-Another example can be found here [[exercise:basis_sets|Basis Sets]]+Another example can be found here [[basis_sets|Basis Sets]]
  
 ==== 3.Step ==== ==== 3.Step ====
Line 145: Line 185:
 O BASIS SET O BASIS SET
 <code> <code>
-#O  pc-2 +#O  pc-1 
-  9 +  5 
-  1  0  0  8  2 +  1  0  0  6  2 
-      14782.00000000          0.00053519          0.00000000 +       2306.70000000          0.00539400          0.00000000 
-       2217.30000000          0.00413750         -0.00000193 +        347.15000000          0.04024800         -0.00031692 
-        504.74000000          0.02124500         -0.00005796 +         78.89000000          0.17921000         -0.00259440 
-        142.87000000          0.08245300         -0.00079494 +         21.87600000          0.45978000         -0.03624100 
-         46.30000000          0.23671000         -0.00731250 +          6.66460000          0.45234000         -0.08779000 
-         16.33700000          0.44039000         -0.04057400 +          1.06690000          0.00000000          0.53320000
-          5.98280000          0.36465000         -0.09159400 +
-          1.67180000          0.00000000          0.20940000+
   1  0  0  1  1   1  0  0  1  1
-          0.64662000          1.00000000 +          0.30700000          1.00000000 
-  1  0  0  1  1 +  1  1  1  3  1 
-          0.21669000          1.00000000 +         17.02200000          0.04891900 
-  1  1  1  4  1 +          3.68380000          0.24962000 
-         60.42400000          0.00689490 +          0.99234000          0.51347000
-         13.93500000          0.04900500 +
-          4.15310000          0.18255000 +
-          1.41580000          0.37633000+
   1  1  1  1  1   1  1  1  1  1
-          0.47549000          1.00000000 +          0.24487000          1.00000000
-  1  1  1  1  1 +
-          0.14529000          1.00000000 +
-  1  2  2  1  1 +
-          2.20000000          1.00000000+
   1  2  2  1  1   1  2  2  1  1
-          0.65000000          1.00000000 +          1.00000000          1.00000000
-  1  3  3  1  1 +
-          1.10000000          1.00000000+
  
 </code> </code>
Line 186: Line 215:
  
   * O2 coordinates   * O2 coordinates
-<code> +<code>   
-  O         4.4705353672        4.7577727756        4.9999951398 +O         4.4720538104        4.7584649515        4.9999999998 
-  O         5.5294646329        5.2422272244        5.0000048599+O         5.5279461896        5.2415350485        4.9999999995
 </code> </code>
   * CO2 coordinates   * CO2 coordinates
 +<code>  
 +  C         4.9999776408        4.9999662056        4.9999894728
 +  O         5.6486993295        5.9339540261        5.0004691016
 +  O         4.3512530072        4.0659797648        4.9995464311
 +</code>
   * H2O coordinates   * H2O coordinates
 +<code>  
 +  O         4.6926974603        4.7525411835        4.6307067609
 +  H         5.6350172910        4.8022721035        4.7052454388
 +  H         4.3528571397        5.2445222023        5.3644975249
 +</code>
  
exercises/2014_ethz_mmm/reaction_energy.1395849851.txt.gz · Last modified: 2020/08/21 10:14 (external edit)