CP2K Open Source Molecular Dynamics

User Tools

Site Tools

Trace:
exercises:2015_ethz_mmm:basis_sets

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
exercises:2015_ethz_mmm:basis_sets [2015/04/16 14:09] – oschuettexercises:2015_ethz_mmm:basis_sets [2020/08/21 10:15] (current) – external edit 127.0.0.1
Line 1: Line 1:
  
 ====== Basis Sets ====== ====== Basis Sets ======
-In this exercise you will compare different basis sets and use them for computing the binding energy of an H2 molecule.+In this exercise you will compare different basis sets and use them for computing the binding energy of an H$_2$ molecule.
  
 The cp2k basis set format is described in detail [[:basis_sets|here]]. The cp2k basis set format is described in detail [[:basis_sets|here]].
  
-===== Part I: Different basis sets for H and H2 =====+===== Part I: Different basis sets for H and H$_2$ =====
 ==== 1.Step ==== ==== 1.Step ====
  
Line 77: Line 77:
  
  
-^ Basis set                     ^ Energy  H (Eh)         ^      +^ Basis set                     ^ Energy  H ($E_h$)      ^      
 | mybasis (from given input)    | ....                   | | mybasis (from given input)    | ....                   |
 | basis try 1                   | ....                   |  | basis try 1                   | ....                   | 
Line 89: Line 89:
  
 ==== 3.Step ==== ==== 3.Step ====
-Repeat the procedure for H2. For this you will have to add a second H atom to the coordinate section and run a geometry optimization to determine the equilibrium distance. Not that this equilibrium distance might also depend on your basis-set. +Repeat the procedure for H$_2$. For this you will have to add a second H atom to the coordinate section and run a **geometry optimization** to determine the equilibrium distance. Howto run a geometry optimization was covered in a [[geometry_optimization|previous exercise]]. Note that the equilibrium distance will depend on your basis set.
  
-<note important> The H2 molecule does not have unpaired electrons. Remember to take out the LSD and MULTIPLICITY keywords.</note>+<note important> The H$_2$ molecule does not have unpaired electrons. Remember to take out the LSD and MULTIPLICITY keywords.</note>
  
-===== Part II: Estimate the binding energy of H2 =====+===== Part II: Estimate the binding energy of H$_2$ =====
  
-Binding energy: \\ +Based on the formula for the //binding energy//, you can now update your table.
  
 \[ \sum E_\text{products} -  \sum E_\text{reactants} = E(H_2) - 2 \cdot E(H) \] \[ \sum E_\text{products} -  \sum E_\text{reactants} = E(H_2) - 2 \cdot E(H) \]
  
-<note important> The binding energy value is significant only if the same basis is used for both reactants and products. </note> 
  
 +^ Basis set                  ^ Energy H [$E_h$] ^ Energy H$_2$ [$E_h$] ^ Distance H$_2$ [$Å$] ^ Binding Energy H$_2$ [$E_h$] ^
 +| mybasis (from given input) | ....             | ....              | ....            | ....                      |
 +| basis try 1                | ....             | ....              | ....            | ....                      |
 +| basis try 2                | ....             | ....              | ....            | ....                      |
 +| ....                       | ....             | ....              | ....            | ....                      |
 +| pc-0                       | ....             | ....              | ....            | ....                      |
 +| pc-1                       | ....             | ....              | ....            | ....                      |
 +| pc-2                       | ....             | ....              | ....            | ....                      |
 +| ....                       | ....             | ....              | ....            | ....                      |
  
-You can now update your table: +<note important>The binding energy is only significant if all terms were calculated with the same basis-set.</note>
- +
-^ Basis set                     ^ Energy  H (Eh)       ^ Energy H2 (Eh) ^ Binding Energy H2 (Eh) ^ +
-| mybasis (from given input)    | .... | .... | ....  | +
-| basis try 1                   | .... | .... | ....  | +
-| basis try 2                   | .... | .... | ....  | +
-| ....                          | .... | .... | ....  | +
-| pc-0                          | .... | .... | ....  | +
-| pc-1                          | .... | .... | ....  | +
-| pc-2                          | .... | .... | ....  | +
-| ....                          | .... | .... | ....  |+
  
 ===== Part III: Questions ===== ===== Part III: Questions =====
- +  - What is the effect of changing the exponents in a basis set? 
-- What is the effect of:\\ +  - What is the effect of adding p- and d-function to the basis set? Do H and H$_2$ respond differently?
-  * increasing/decreasing the value of the exponents for the given basis? +
-  * adding sets with p,d symmetry to the basis? You have the same effect in H and H2? +
- +
- +
-===== Part IV: Additional basis sets =====+
  
  
 +===== Appendix: Literature Basis Sets =====
  
 +<code>
 H  pc-0 H  pc-0
   2   2
Line 162: Line 157:
           1.25000000          1.00000000           1.25000000          1.00000000
  
- +</code>
  
  
exercises/2015_ethz_mmm/basis_sets.1429193376.txt.gz · Last modified: 2020/08/21 10:14 (external edit)