exercises:2014_ethz_mmm:reaction_energy
Differences
This shows you the differences between two versions of the page.
Both sides previous revisionPrevious revisionNext revision | Previous revision | ||
exercise:reaction_energy [2014/03/26 13:58] – sclelia | exercises: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 | + | In this exercise, you will calculate the reaction energy for the methane |
- | {{ : | + | \[ CH_4 + 2O_2 \rightarrow CO_2 + 2H_2O \] |
Reaction energy: | Reaction energy: | ||
- | {{ : | + | \[ \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: [[basis_sets|Basis Sets]] | ||
<code - CH4.inp > | <code - CH4.inp > | ||
+ | |||
+ | &GLOBAL | ||
+ | PROJECT CH4 | ||
+ | RUN_TYPE ENERGY | ||
+ | PRINT_LEVEL MEDIUM | ||
+ | &END GLOBAL | ||
+ | |||
+ | & | ||
+ | METHOD Quickstep | ||
+ | &DFT | ||
+ | & | ||
+ | PERIODIC NONE | ||
+ | PSOLVER | ||
+ | &END POISSON | ||
+ | & | ||
+ | METHOD GAPW ! Method: gaussian and augmented plane waves | ||
+ | &END QS | ||
+ | |||
+ | & | ||
+ | & | ||
+ | &END XC_FUNCTIONAL | ||
+ | & | ||
+ | & | ||
+ | EPS_SCHWARZ 1.0E-10 | ||
+ | &END SCREENING | ||
+ | &END HF | ||
+ | &END XC | ||
+ | &END DFT | ||
+ | |||
+ | &SUBSYS | ||
+ | &CELL | ||
+ | ABC 10 10 10 | ||
+ | PERIODIC NONE ! Non periodic calculations. That's why the POISSON section is needed | ||
+ | &END CELL | ||
+ | & | ||
+ | & | ||
+ | &END | ||
+ | &END | ||
+ | &COORD | ||
+ | C | ||
+ | H | ||
+ | H | ||
+ | H | ||
+ | H | ||
+ | &END COORD | ||
+ | &KIND H ! potential and basis for H | ||
+ | & | ||
+ | 3 | ||
+ | 1 0 0 3 1 | ||
+ | | ||
+ | 1.86870000 | ||
+ | 0.41821000 | ||
+ | 1 0 0 1 1 | ||
+ | 0.10610000 | ||
+ | 1 1 1 1 1 | ||
+ | 1.00000000 | ||
+ | & | ||
+ | POTENTIAL ALL | ||
+ | & | ||
+ | | ||
+ | | ||
+ | & | ||
+ | &END KIND | ||
+ | &KIND C ! potential and basis for C | ||
+ | &BASIS | ||
+ | 5 | ||
+ | 1 0 0 6 2 | ||
+ | | ||
+ | 188.57000000 | ||
+ | | ||
+ | | ||
+ | 3.55670000 | ||
+ | 0.54258000 | ||
+ | 1 0 0 1 1 | ||
+ | 0.16058000 | ||
+ | 1 1 1 3 1 | ||
+ | 9.14260000 | ||
+ | 1.92980000 | ||
+ | 0.52522000 | ||
+ | 1 1 1 1 1 | ||
+ | 0.13608000 | ||
+ | 1 2 2 1 1 | ||
+ | 0.80000000 | ||
+ | &END | ||
+ | POTENTIAL ALL | ||
+ | & | ||
+ | 4 2 0 | ||
+ | 0.34883045 | ||
+ | &END | ||
+ | &END KIND | ||
+ | |||
+ | &END SUBSYS | ||
+ | &END FORCE_EVAL | ||
+ | |||
</ | </ | ||
Line 48: | Line 144: | ||
==== 2.Step ==== | ==== 2.Step ==== | ||
+ | |||
Modify the input in order to perform the same calculation for: | Modify the input in order to perform the same calculation for: | ||
* H2O | * H2O | ||
* CO2 | * CO2 | ||
* O2 TRIPLET | * O2 TRIPLET | ||
- | <note tip> | + | <note tip> |
- | Remember that the O2 ground state is a triplet state, with non paired electrons | + | <note important> |
+ | MULTIPLICITY=2S+1=3. </ | ||
+ | |||
+ | For O2 triplet, the LSD and MULTIPLICITY keywords are needed in the DFT section: | ||
< | < | ||
+ | METHOD Quickstep | ||
+ | & | ||
+ | | ||
+ | | ||
+ | ... | ||
+ | | ||
</ | </ | ||
- | 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 74: | Line 180: | ||
==== Questions ==== | ==== Questions ==== | ||
* Report the total energy for all the molecules, and the overall reaction energy for CH4 combustion. | * Report the total energy for all the molecules, and the overall reaction energy for CH4 combustion. | ||
- | * (Optional) Compute the total enery for O2 singlet. | + | * (Optional) Compute the total enery for O2 singlet |
==== Additional material ==== | ==== Additional material ==== | ||
- | * O BASIS SET | + | O BASIS SET |
- | * O POTENTIAl | + | < |
+ | #O pc-1 | ||
+ | 5 | ||
+ | 1 0 0 6 2 | ||
+ | | ||
+ | 347.15000000 | ||
+ | | ||
+ | | ||
+ | 6.66460000 | ||
+ | 1.06690000 | ||
+ | 1 0 0 1 1 | ||
+ | 0.30700000 | ||
+ | 1 1 1 3 1 | ||
+ | | ||
+ | 3.68380000 | ||
+ | 0.99234000 | ||
+ | 1 1 1 1 1 | ||
+ | 0.24487000 | ||
+ | 1 2 2 1 1 | ||
+ | 1.00000000 | ||
+ | |||
+ | </ | ||
+ | |||
+ | * O potential | ||
+ | < | ||
+ | #O ALLELECTRON ALL | ||
+ | 4 4 0 | ||
+ | | ||
+ | </ | ||
* O2 coordinates | * O2 coordinates | ||
+ | < | ||
+ | O | ||
+ | O | ||
+ | </ | ||
* CO2 coordinates | * CO2 coordinates | ||
- | | + | < |
+ | C | ||
+ | O | ||
+ | O | ||
+ | </ | ||
+ | | ||
+ | < | ||
+ | O | ||
+ | H | ||
+ | H | ||
+ | </ | ||
exercises/2014_ethz_mmm/reaction_energy.1395842284.txt.gz · Last modified: 2020/08/21 10:14 (external edit)