exercises:2014_ethz_mmm:basis_sets
Differences
This shows you the differences between two versions of the page.
Both sides previous revisionPrevious revisionNext revision | Previous revision | ||
exercise:basis_sets [2014/03/26 15:25] – [Part III: Questions] sclelia | exercises:2014_ethz_mmm:basis_sets [2020/08/21 10:15] (current) – external edit 127.0.0.1 | ||
---|---|---|---|
Line 3: | Line 3: | ||
In this exericse you will compare different basis sets and use them for computing the binding energy of an H2 molecule. | In this exericse you will compare different basis sets and use them for computing the binding energy of an H2 molecule. | ||
- | The cp2k basis set format is the following: | + | The cp2k basis set format is described in detail [[:basis_sets|here]]. |
- | + | ||
- | < | + | |
- | + | ||
- | Nsets | + | |
- | n1 lmin lmax Nexp Ncontr | + | |
- | exp1 contr1 contr2 contr3 | + | |
- | exp2 contr4 contr5 contr6 ... | + | |
- | Nexp.... | + | |
- | n2 lmin lmax Nexp Ncontr | + | |
- | exp1 contr1 contr2 contr3 ... | + | |
- | exp2 contr4 contr5 contr6 ... | + | |
- | Nexp.... | + | |
- | ... | + | |
- | </ | + | |
- | \\ | + | |
- | Nsets = number of sets. \\ | + | |
- | n1,n2... = main quantum number (but it gets ignored by the program!) \\ | + | |
- | lmin = min //l// (angular quantum number) of the set ( s = 0 ; p = 1 ; d = 2 ...) \\ | + | |
- | lmax = max //l// (angular quantum number) of the set ( s = 0 ; p = 1 ; d = 2 ...) \\ | + | |
- | Nexp = number of exponents of the set \\ | + | |
- | Ncontr = number of contaction coefficients per angular quantum number \\ | + | |
- | + | ||
- | As an example: | + | |
- | + | ||
- | < | + | |
- | 2 2 Sets | + | |
- | 1 0 0 1 1 Set 1: lmin=0 ; lmax=0 (→ s functions!) ; 1 exponent ; 1 contraction | + | |
- | 0.35 1 exponent1 of set 1 ; contraction | + | |
- | 1 1 1 1 1 Set 2: lmin=1 ; lmax=1 (→ p functions!) ; 1 exponenet ; 1 contraction | + | |
- | 0.6 1 | + | |
- | </ | + | |
===== Part I: Different basis sets for H and H2 ===== | ===== Part I: Different basis sets for H and H2 ===== | ||
Line 74: | Line 42: | ||
| | ||
&SUBSYS | &SUBSYS | ||
+ | & | ||
+ | & | ||
+ | &END | ||
+ | &END | ||
&CELL | &CELL | ||
ABC 10.0 10.0 10.0 | ABC 10.0 10.0 10.0 | ||
Line 90: | Line 62: | ||
& | & | ||
| | ||
- | & | + | & |
| | ||
| | ||
Line 121: | Line 93: | ||
< | < | ||
H 0 0 0 | H 0 0 0 | ||
- | H 0.8 0 0 | + | H 0.75 0 0 |
</ | </ | ||
Line 129: | Line 101: | ||
Binding energy: \\ | Binding energy: \\ | ||
- | {{ : | ||
+ | \[ \sum E_\text{products} - \sum E_\text{rectants} = E(H_2) - 2 \cdot E(H) \] | ||
<note important> | <note important> |
exercises/2014_ethz_mmm/basis_sets.1395847530.txt.gz · Last modified: 2020/08/21 10:14 (external edit)