You signed in with another tab or window. Reload to refresh your session.You signed out in another tab or window. Reload to refresh your session.You switched accounts on another tab or window. Reload to refresh your session.Dismiss alert
# Project #4: The second-order Moller-Plesset perturbation (MP2) energy
2
2
Unlike the Hartree-Fock energy, correlation energies like the MP2 energy are usually expressed in terms of MO-basis quantities (integrals, MO energies).
3
-
The most expensive part of the calculation is the transformation of the two-electron integrals from the AO to the MO basis representation.
4
-
The purpose of this project is to understand this transformation and fundamental techniques for its efficient implementation.
3
+
The most expensive part of the calculation is the transformation of the two-electron integrals from the AO to the MO basis representation.
4
+
The purpose of this project is to understand this transformation and fundamental techniques for its efficient implementation.
5
5
The theoretical background and a concise set of instructions for this project may be found [here](./project4-instructions.pdf).
6
6
7
7
## Step #1: Read the Two-Electron Integrals
8
8
The Mulliken-ordered integrals are defined as:
9
9
10
10
<imgsrc="./figures/eri.png"height="45">
11
11
12
-
Concise instructions for this step can be found in
This approach is easy to implement (hence the word "[noddy](http://www.hackerslang.com/noddy.html)" above), but is expensive due to its N<sup>8</sup> computational order.
26
-
Nevertheless, you should start with this algorithm to get your code working, and run timings (use the UNIX "time" command)
27
-
for the test cases below to get an idea of its computational cost.
23
+
This approach is easy to implement (hence the word "[noddy](http://www.hackerslang.com/noddy.html)" above), but is expensive due to its N<sup>8</sup> computational order. Nevertheless, you should start with this algorithm to get your code working, and run timings (use the UNIX "time" command) for the test cases below to get an idea of its computational cost.
28
24
29
25
* Hint 1: Noddy transformation code
30
26
### The Smarter Algorithm
@@ -34,8 +30,7 @@ Notice that none of the *C* coefficients in the above expression have any indice
This means that each summation within brackets could be carried out separately,
37
-
starting from the innermost summation over <html>σ</html>, if we store the results at each step.
38
-
This reduces the N<sup>8</sup> algorithm above to four N<sup>5</sup> steps.
33
+
starting from the innermost summation over <html>σ</html>, if we store the results at each step. This reduces the N<sup>8</sup> algorithm above to four N<sup>5</sup> steps.
39
34
40
35
Be careful about the permutational symmetry of the integrals and intermediates at each step:
41
36
while the AO-basis integrals have eight-fold permutational symmetry, the partially transformed integrals have much less symmetry.
@@ -53,10 +48,10 @@ where *i* and *j* denote doubly-occupied orbitals and *a* and *b* unoccupied orb
53
48
* Hint 3: Occupied versus unoccupied orbitals
54
49
55
50
## Test Cases
56
-
The input structures, integrals, etc. for these examples are found in the [input directory](https://github.com/CrawfordGroup/ProgrammingProjects/tree/master/Project%2304/input).
51
+
The input structures, integrals, etc. for these examples are found in the [input directory](./input).
57
52
58
-
* STO-3G Water | [output](https://github.com/CrawfordGroup/ProgrammingProjects/blob/master/Project%2304/output/h2o/STO-3G/output.txt)
59
-
* DZ Water | [output](https://github.com/CrawfordGroup/ProgrammingProjects/blob/master/Project%2304/output/h2o/DZ/output.txt)
60
-
* DZP Water | [output](https://github.com/CrawfordGroup/ProgrammingProjects/blob/master/Project%2304/output/h2o/DZP/output.txt)
0 commit comments