diff --git a/doc/pub/week48/html/week48-bs.html b/doc/pub/week48/html/week48-bs.html index 1ed7e74..b7a87ee 100644 --- a/doc/pub/week48/html/week48-bs.html +++ b/doc/pub/week48/html/week48-bs.html @@ -39,10 +39,34 @@ 'sections': [('Week 48, November 25-29, 2024', 2, None, - 'week-48-november-25-29-2024')]} + 'week-48-november-25-29-2024'), + ('CCSD with twobody Hamiltonian', + 2, + None, + 'ccsd-with-twobody-hamiltonian'), + ('Two-body normal-ordered Hamiltonian', + 2, + None, + 'two-body-normal-ordered-hamiltonian')]} end of tocinfo --> + + + + + + + diff --git a/doc/pub/week48/html/week48-solarized.html b/doc/pub/week48/html/week48-solarized.html index 3108827..476f098 100644 --- a/doc/pub/week48/html/week48-solarized.html +++ b/doc/pub/week48/html/week48-solarized.html @@ -40,10 +40,34 @@ 'sections': [('Week 48, November 25-29, 2024', 2, None, - 'week-48-november-25-29-2024')]} + 'week-48-november-25-29-2024'), + ('CCSD with twobody Hamiltonian', + 2, + None, + 'ccsd-with-twobody-hamiltonian'), + ('Two-body normal-ordered Hamiltonian', + 2, + None, + 'two-body-normal-ordered-hamiltonian')]} end of tocinfo --> + + + + + + +

Week 48: Coupled cluster theory and summary of course

@@ -59,7 +83,7 @@

Week 48: Coupled cluster theory and summary of course


-

Week 48, November 25-29

+

November 25-29, 2024


@@ -72,13 +96,71 @@

Week 48, November 25-29, 2024

  • How to write your own coupled-cluster theory code, pairing model example
  • Coupled cluster theory, singles and doubles excitations, diagrammatic expansion
  • -
  • Firday: +
  • Friday:
    1. Coupled cluster theory for singles and doubles excitations using a diagrammatic derivation
    2. Summary of course and discussion of final oral exam
  • Lecture material: Lecture notes and Shavitt and Bartlett chapters 9 and 10. See also slides at https://github.com/ManyBodyPhysics/FYS4480/blob/master/doc/pub/week48/pdf/cc.pdf
  • +









    +

    CCSD with twobody Hamiltonian

    + +

    Truncating the cluster operator \( \hat{T} \) at the \( n=2 \) level, defines +CCSD approximation to the Coupled Cluster wavefunction. The +coupled cluster wavefunction is now given by +

    + +$$ +\begin{equation*} + \vert \Psi_{CC}\rangle = e^{\hat{T}_1 + \hat{T}_2} \vert \Phi_0\rangle, +\end{equation*} +$$ + +

    where

    +$$ +\begin{align*} + \hat{T}_1 &= + \sum_{ia} + t_{i}^{a} a_{a}^\dagger a_i \\ + \hat{T}_2 &= \frac{1}{4} + \sum_{ijab} + t_{ij}^{ab} a_{a}^\dagger a_{b}^\dagger a_{j} a_{i}. +\end{align*} +$$ + + +









    +

    Two-body normal-ordered Hamiltonian

    + +$$ +\begin{align*} + \hat{H} &= \sum_{pq} f_q^p\left\{ a_p^\dagger a_q \right\} + + \frac{1}{4} \sum_{pqrs} \langle pq \vert \hat{v} \vert rs \rangle \left\{ a_p^\dagger a_q^\dagger a_s a_r \right\} \\ + & \quad + \mathrm{E}_0 \\ + &= \hat{F}_N + \hat{V}_N + \mathrm{E}_0 + = \hat{H}_N + \mathrm{E}_0, +\end{align*} +$$ + +

    where

    +$$ +\begin{align*} + f_q^p &= \langle p\vert \hat{t} \vert q\rangle + \sum_i \langle pi\vert \hat{v} \vert qi\rangle \\ + \langle pq \vert \vert rs \rangle &= \langle pq \vert \hat{v} \vert rs \rangle \\ + \mathrm{E}_0 &= \sum_i \langle i\vert \hat{t} \vert i\rangle + \frac{1}{2} \sum_{ij} \langle ij\vert \hat{v} \vert ij\rangle. +\end{align*} +$$ + + + + + + + + + +
    © 1999-2024, Morten Hjorth-Jensen. Released under CC Attribution-NonCommercial 4.0 license diff --git a/doc/pub/week48/html/week48.html b/doc/pub/week48/html/week48.html index ef7ff9a..b8c4d25 100644 --- a/doc/pub/week48/html/week48.html +++ b/doc/pub/week48/html/week48.html @@ -117,10 +117,34 @@ 'sections': [('Week 48, November 25-29, 2024', 2, None, - 'week-48-november-25-29-2024')]} + 'week-48-november-25-29-2024'), + ('CCSD with twobody Hamiltonian', + 2, + None, + 'ccsd-with-twobody-hamiltonian'), + ('Two-body normal-ordered Hamiltonian', + 2, + None, + 'two-body-normal-ordered-hamiltonian')]} end of tocinfo --> + + + + + + +

    Week 48: Coupled cluster theory and summary of course

    @@ -136,7 +160,7 @@

    Week 48: Coupled cluster theory and summary of course


    -

    Week 48, November 25-29

    +

    November 25-29, 2024


    @@ -149,13 +173,71 @@

    Week 48, November 25-29, 2024

  • How to write your own coupled-cluster theory code, pairing model example
  • Coupled cluster theory, singles and doubles excitations, diagrammatic expansion
  • -
  • Firday: +
  • Friday:
    1. Coupled cluster theory for singles and doubles excitations using a diagrammatic derivation
    2. Summary of course and discussion of final oral exam
  • Lecture material: Lecture notes and Shavitt and Bartlett chapters 9 and 10. See also slides at https://github.com/ManyBodyPhysics/FYS4480/blob/master/doc/pub/week48/pdf/cc.pdf
  • +









    +

    CCSD with twobody Hamiltonian

    + +

    Truncating the cluster operator \( \hat{T} \) at the \( n=2 \) level, defines +CCSD approximation to the Coupled Cluster wavefunction. The +coupled cluster wavefunction is now given by +

    + +$$ +\begin{equation*} + \vert \Psi_{CC}\rangle = e^{\hat{T}_1 + \hat{T}_2} \vert \Phi_0\rangle, +\end{equation*} +$$ + +

    where

    +$$ +\begin{align*} + \hat{T}_1 &= + \sum_{ia} + t_{i}^{a} a_{a}^\dagger a_i \\ + \hat{T}_2 &= \frac{1}{4} + \sum_{ijab} + t_{ij}^{ab} a_{a}^\dagger a_{b}^\dagger a_{j} a_{i}. +\end{align*} +$$ + + +









    +

    Two-body normal-ordered Hamiltonian

    + +$$ +\begin{align*} + \hat{H} &= \sum_{pq} f_q^p\left\{ a_p^\dagger a_q \right\} + + \frac{1}{4} \sum_{pqrs} \langle pq \vert \hat{v} \vert rs \rangle \left\{ a_p^\dagger a_q^\dagger a_s a_r \right\} \\ + & \quad + \mathrm{E}_0 \\ + &= \hat{F}_N + \hat{V}_N + \mathrm{E}_0 + = \hat{H}_N + \mathrm{E}_0, +\end{align*} +$$ + +

    where

    +$$ +\begin{align*} + f_q^p &= \langle p\vert \hat{t} \vert q\rangle + \sum_i \langle pi\vert \hat{v} \vert qi\rangle \\ + \langle pq \vert \vert rs \rangle &= \langle pq \vert \hat{v} \vert rs \rangle \\ + \mathrm{E}_0 &= \sum_i \langle i\vert \hat{t} \vert i\rangle + \frac{1}{2} \sum_{ij} \langle ij\vert \hat{v} \vert ij\rangle. +\end{align*} +$$ + + + + + + + + + +
    © 1999-2024, Morten Hjorth-Jensen. Released under CC Attribution-NonCommercial 4.0 license diff --git a/doc/pub/week48/ipynb/ipynb-week48-src.tar.gz b/doc/pub/week48/ipynb/ipynb-week48-src.tar.gz index 598e231..5285d30 100644 Binary files a/doc/pub/week48/ipynb/ipynb-week48-src.tar.gz and b/doc/pub/week48/ipynb/ipynb-week48-src.tar.gz differ diff --git a/doc/pub/week48/ipynb/week48.ipynb b/doc/pub/week48/ipynb/week48.ipynb index b7f2fec..989ec22 100644 --- a/doc/pub/week48/ipynb/week48.ipynb +++ b/doc/pub/week48/ipynb/week48.ipynb @@ -2,7 +2,7 @@ "cells": [ { "cell_type": "markdown", - "id": "fbc1ab02", + "id": "5b72ba9b", "metadata": { "editable": true }, @@ -14,7 +14,7 @@ }, { "cell_type": "markdown", - "id": "e687d25a", + "id": "2f04faba", "metadata": { "editable": true }, @@ -22,12 +22,12 @@ "# Week 48: Coupled cluster theory and summary of course \n", "**Morten Hjorth-Jensen**, Department of Physics and Center for Computing in Science Education, University of Oslo, Norway\n", "\n", - "Date: **Week 48, November 25-29**" + "Date: **November 25-29, 2024**" ] }, { "cell_type": "markdown", - "id": "1268050a", + "id": "d3a293b3", "metadata": { "editable": true }, @@ -41,7 +41,7 @@ "\n", "c. Coupled cluster theory, singles and doubles excitations, diagrammatic expansion\n", "\n", - "4. Firday:\n", + "4. Friday:\n", "\n", "a. Coupled cluster theory for singles and doubles excitations using a diagrammatic derivation\n", "\n", @@ -49,6 +49,131 @@ "\n", "3. Lecture material: Lecture notes and Shavitt and Bartlett chapters 9 and 10. See also slides at " ] + }, + { + "cell_type": "markdown", + "id": "b5c08389", + "metadata": { + "editable": true + }, + "source": [ + "## CCSD with twobody Hamiltonian\n", + "\n", + "Truncating the cluster operator $\\hat{T}$ at the $n=2$ level, defines\n", + "CCSD approximation to the Coupled Cluster wavefunction. The\n", + "coupled cluster wavefunction is now given by" + ] + }, + { + "cell_type": "markdown", + "id": "36a9582e", + "metadata": { + "editable": true + }, + "source": [ + "$$\n", + "\\vert \\Psi_{CC}\\rangle = e^{\\hat{T}_1 + \\hat{T}_2} \\vert \\Phi_0\\rangle,\n", + "$$" + ] + }, + { + "cell_type": "markdown", + "id": "1ada48ee", + "metadata": { + "editable": true + }, + "source": [ + "where" + ] + }, + { + "cell_type": "markdown", + "id": "a0bf1b88", + "metadata": { + "editable": true + }, + "source": [ + "$$\n", + "\\begin{align*}\n", + " \\hat{T}_1 &= \n", + " \\sum_{ia}\n", + " t_{i}^{a} a_{a}^\\dagger a_i \\\\\n", + " \\hat{T}_2 &= \\frac{1}{4} \n", + " \\sum_{ijab}\n", + " t_{ij}^{ab} a_{a}^\\dagger a_{b}^\\dagger a_{j} a_{i}.\n", + "\\end{align*}\n", + "$$" + ] + }, + { + "cell_type": "markdown", + "id": "0d82aa94", + "metadata": { + "editable": true + }, + "source": [ + "## Two-body normal-ordered Hamiltonian" + ] + }, + { + "cell_type": "markdown", + "id": "8dd39643", + "metadata": { + "editable": true + }, + "source": [ + "$$\n", + "\\begin{align*}\n", + " \\hat{H} &= \\sum_{pq} f_q^p\\left\\{ a_p^\\dagger a_q \\right\\} + \n", + " \\frac{1}{4} \\sum_{pqrs} \\langle pq \\vert \\hat{v} \\vert rs \\rangle \\left\\{ a_p^\\dagger a_q^\\dagger a_s a_r \\right\\} \\\\\n", + " & \\quad + \\mathrm{E}_0 \\\\\n", + " &= \\hat{F}_N + \\hat{V}_N + \\mathrm{E}_0 \n", + " = \\hat{H}_N + \\mathrm{E}_0,\n", + "\\end{align*}\n", + "$$" + ] + }, + { + "cell_type": "markdown", + "id": "56de3d02", + "metadata": { + "editable": true + }, + "source": [ + "where" + ] + }, + { + "cell_type": "markdown", + "id": "643501dc", + "metadata": { + "editable": true + }, + "source": [ + "$$\n", + "\\begin{align*}\n", + " f_q^p &= \\langle p\\vert \\hat{t} \\vert q\\rangle + \\sum_i \\langle pi\\vert \\hat{v} \\vert qi\\rangle \\\\\n", + " \\langle pq \\vert \\vert rs \\rangle &= \\langle pq \\vert \\hat{v} \\vert rs \\rangle \\\\\n", + " \\mathrm{E}_0 &= \\sum_i \\langle i\\vert \\hat{t} \\vert i\\rangle + \\frac{1}{2} \\sum_{ij} \\langle ij\\vert \\hat{v} \\vert ij\\rangle.\n", + "\\end{align*}\n", + "$$" + ] + }, + { + "cell_type": "markdown", + "id": "21ab8f94", + "metadata": { + "editable": true + }, + "source": [ + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "" + ] } ], "metadata": {}, diff --git a/doc/pub/week48/pdf/week48.pdf b/doc/pub/week48/pdf/week48.pdf index 12c77b3..579579f 100644 Binary files a/doc/pub/week48/pdf/week48.pdf and b/doc/pub/week48/pdf/week48.pdf differ diff --git a/doc/src/week48/week48.do.txt b/doc/src/week48/week48.do.txt index 2b71dd5..e57ad23 100644 --- a/doc/src/week48/week48.do.txt +++ b/doc/src/week48/week48.do.txt @@ -1,6 +1,6 @@ TITLE: Week 48: Coupled cluster theory and summary of course AUTHOR: Morten Hjorth-Jensen {copyright, 1999-present|CC BY-NC} at Department of Physics and Center for Computing in Science Education, University of Oslo, Norway -DATE: Week 48, November 25-29 +DATE: November 25-29, 2024 !split ===== Week 48, November 25-29, 2024 ===== @@ -8,7 +8,64 @@ o Thursday: o Short repetition from last week o How to write your own coupled-cluster theory code, pairing model example o Coupled cluster theory, singles and doubles excitations, diagrammatic expansion -o Firday: +o Friday: o Coupled cluster theory for singles and doubles excitations using a diagrammatic derivation o Summary of course and discussion of final oral exam o Lecture material: Lecture notes and Shavitt and Bartlett chapters 9 and 10. See also slides at URL:"https://github.com/ManyBodyPhysics/FYS4480/blob/master/doc/pub/week48/pdf/cc.pdf" + + +!split +===== CCSD with twobody Hamiltonian ===== + +Truncating the cluster operator $\hat{T}$ at the $n=2$ level, defines +CCSD approximation to the Coupled Cluster wavefunction. The +coupled cluster wavefunction is now given by + +!bt +\begin{equation*} + \vert \Psi_{CC}\rangle = e^{\hat{T}_1 + \hat{T}_2} \vert \Phi_0\rangle, +\end{equation*} +!et +where +!bt +\begin{align*} + \hat{T}_1 &= + \sum_{ia} + t_{i}^{a} a_{a}^\dagger a_i \\ + \hat{T}_2 &= \frac{1}{4} + \sum_{ijab} + t_{ij}^{ab} a_{a}^\dagger a_{b}^\dagger a_{j} a_{i}. +\end{align*} +!et + +!split +===== Two-body normal-ordered Hamiltonian ===== + +!bt +\begin{align*} + \hat{H} &= \sum_{pq} f_q^p\left\{ a_p^\dagger a_q \right\} + + \frac{1}{4} \sum_{pqrs} \langle pq \vert \hat{v} \vert rs \rangle \left\{ a_p^\dagger a_q^\dagger a_s a_r \right\} \\ + & \quad + \mathrm{E}_0 \\ + &= \hat{F}_N + \hat{V}_N + \mathrm{E}_0 + = \hat{H}_N + \mathrm{E}_0, +\end{align*} +!et +where +!bt +\begin{align*} + f_q^p &= \langle p\vert \hat{t} \vert q\rangle + \sum_i \langle pi\vert \hat{v} \vert qi\rangle \\ + \langle pq \vert \vert rs \rangle &= \langle pq \vert \hat{v} \vert rs \rangle \\ + \mathrm{E}_0 &= \sum_i \langle i\vert \hat{t} \vert i\rangle + \frac{1}{2} \sum_{ij} \langle ij\vert \hat{v} \vert ij\rangle. +\end{align*} +!et + +#%\input{src/ccsd_summary02} +#%\input{src/ccsd_summary03} + +#\include{src/ccsd_H2diagrams} +#\include{src/ccsd_barh_expansion} +#\include{src/ccsd_factoring} +#\include{src/ccsd_algorithm} + + +