1. Many-body Hamiltonians, basic linear algebra and Second Quantization

1.1. Definitions and notations

Before we proceed we need some definitions. We will assume that the interacting part of the Hamiltonian can be approximated by a two-body interaction. This means that our Hamiltonian is written as the sum of some onebody part and a twobody part

(1)H^=H^0+H^I=i=1Ah^0(xi)+i<jAv^(rij),

with

(2)H0=i=1Ah^0(xi).

The onebody part uext(xi) is normally approximated by a harmonic oscillator potential or the Coulomb interaction an electron feels from the nucleus. However, other potentials are fully possible, such as one derived from the self-consistent solution of the Hartree-Fock equations to be discussed here.

Our Hamiltonian is invariant under the permutation (interchange) of two particles. Since we deal with fermions however, the total wave function is antisymmetric. Let P^ be an operator which interchanges two particles. Due to the symmetries we have ascribed to our Hamiltonian, this operator commutes with the total Hamiltonian,

[H^,P^]=0,

meaning that Ψλ(x1,x2,,xA) is an eigenfunction of P^ as well, that is

P^ijΨλ(x1,x2,,xi,,xj,,xA)=βΨλ(x1,x2,,xi,,xj,,xA),

where β is the eigenvalue of P^. We have introduced the suffix ij in order to indicate that we permute particles i and j. The Pauli principle tells us that the total wave function for a system of fermions has to be antisymmetric, resulting in the eigenvalue β=1.

In our case we assume that we can approximate the exact eigenfunction with a Slater determinant

(3)Φ(x1,x2,,xA,α,β,,σ)=1A!|ψα(x1)ψα(x2)ψα(xA)ψβ(x1)ψβ(x2)ψβ(xA)ψσ(x1)ψσ(x2)ψσ(xA)|,

where xi stand for the coordinates and spin values of a particle i and α,β,,γ are quantum numbers needed to describe remaining quantum numbers.

1.1.1. Brief reminder on some linear algebra properties

Before we proceed with a more compact representation of a Slater determinant, we would like to repeat some linear algebra properties which will be useful for our derivations of the energy as function of a Slater determinant, Hartree-Fock theory and later the nuclear shell model.

The inverse of a matrix is defined by

A1A=I

A unitary matrix A is one whose inverse is its adjoint

A1=A

A real unitary matrix is called orthogonal and its inverse is equal to its transpose. A hermitian matrix is its own self-adjoint, that is

A=A.
Relations Name matrix elements
$A = A^{T}$ symmetric $a_{ij} = a_{ji}$
$A = \left (A^{T} \right )^{-1}$ real orthogonal $\sum_k a_{ik} a_{jk} = \sum_k a_{ki} a_{kj} = \delta_{ij}$
$A = A^{ * }$ real matrix $a_{ij} = a_{ij}^{ * }$
$A = A^{\dagger}$ hermitian $a_{ij} = a_{ji}^{ * }$
$A = \left (A^{\dagger} \right )^{-1}$ unitary $\sum_k a_{ik} a_{jk}^{ * } = \sum_k a_{ki}^{ * } a_{kj} = \delta_{ij}$

Since we will deal with Fermions (identical and indistinguishable particles) we will form an ansatz for a given state in terms of so-called Slater determinants determined by a chosen basis of single-particle functions.

For a given n×n matrix A we can write its determinant

det(A)=|A|=|a11a12a1na21a22a2nan1an2ann|,

in a more compact form as

|A|=i=1n!(1)piP^ia11a22ann,

where P^i is a permutation operator which permutes the column indices 1,2,3,,n and the sum runs over all n! permutations. The quantity pi represents the number of transpositions of column indices that are needed in order to bring a given permutation back to its initial ordering, in our case given by a11a22ann here.

A simple 2×2 determinant illustrates this. We have

det(A)=|a11a12a21a22|=(1)0a11a22+(1)1a12a21,

where in the last term we have interchanged the column indices 1 and 2. The natural ordering we have chosen is a11a22.

1.1.2. Back to the derivation of the energy

The single-particle function ψα(xi) are eigenfunctions of the onebody Hamiltonian hi, that is

h^0(xi)=t^(xi)+u^ext(xi),

with eigenvalues

h^0(xi)ψα(xi)=(t^(xi)+u^ext(xi))ψα(xi)=εαψα(xi).

The energies εα are the so-called non-interacting single-particle energies, or unperturbed energies. The total energy is in this case the sum over all single-particle energies, if no two-body or more complicated many-body interactions are present.

Let us denote the ground state energy by E0. According to the variational principle we have

E0E[Φ]=ΦH^Φdτ

where Φ is a trial function which we assume to be normalized

ΦΦdτ=1,

where we have used the shorthand dτ=dx1dr2drA.

In the Hartree-Fock method the trial function is the Slater determinant of Eq. (3) which can be rewritten as

Φ(x1,x2,,xA,α,β,,ν)=1A!P()PP^ψα(x1)ψβ(x2)ψν(xA)=A!A^ΦH,

where we have introduced the antisymmetrization operator A^ defined by the summation over all possible permutations of two particles.

It is defined as

(4)A^=1A!p()pP^,

with p standing for the number of permutations. We have introduced for later use the so-called Hartree-function, defined by the simple product of all possible single-particle functions

ΦH(x1,x2,,xA,α,β,,ν)=ψα(x1)ψβ(x2)ψν(xA).

Both H^0 and H^I are invariant under all possible permutations of any two particles and hence commute with A^

(5)[H0,A^]=[HI,A^]=0.

Furthermore, A^ satisfies

(6)A^2=A^,

since every permutation of the Slater determinant reproduces it.

The expectation value of H^0

ΦH^0Φdτ=A!ΦHA^H^0A^ΦHdτ

is readily reduced to

ΦH^0Φdτ=A!ΦHH^0A^ΦHdτ,

where we have used Eqs. (5) and (6). The next step is to replace the antisymmetrization operator by its definition and to replace H^0 with the sum of one-body operators

ΦH^0Φdτ=i=1Ap()pΦHh^0P^ΦHdτ.

The integral vanishes if two or more particles are permuted in only one of the Hartree-functions ΦH because the individual single-particle wave functions are orthogonal. We obtain then

ΦH^0Φdτ=i=1AΦHh^0ΦHdτ.

Orthogonality of the single-particle functions allows us to further simplify the integral, and we arrive at the following expression for the expectation values of the sum of one-body Hamiltonians

(7)ΦH^0Φdτ=μ=1Aψμ(x)h^0ψμ(x)dxdr.

We introduce the following shorthand for the above integral

μ|h^0|μ=ψμ(x)h^0ψμ(x)dx,

and rewrite Eq. (7) as

(8)ΦH^0Φdτ=μ=1Aμ|h^0|μ.

The expectation value of the two-body part of the Hamiltonian is obtained in a similar manner. We have

ΦH^IΦdτ=A!ΦHA^H^IA^ΦHdτ,

which reduces to

ΦH^IΦdτ=ij=1Ap()pΦHv^(rij)P^ΦHdτ,

by following the same arguments as for the one-body Hamiltonian.

Because of the dependence on the inter-particle distance rij, permutations of any two particles no longer vanish, and we get

ΦH^IΦdτ=i<j=1AΦHv^(rij)(1Pij)ΦHdτ.

where Pij is the permutation operator that interchanges particle i and particle j. Again we use the assumption that the single-particle wave functions are orthogonal.

We obtain

(9)ΦH^IΦdτ=12μ=1Aν=1A[ψμ(xi)ψν(xj)v^(rij)ψμ(xi)ψν(xj)dxidxj
(10)ψμ(xi)ψν(xj)v^(rij)ψν(xi)ψμ(xj)dxidxj].

The first term is the so-called direct term. It is frequently also called the Hartree term, while the second is due to the Pauli principle and is called the exchange term or just the Fock term. The factor 1/2 is introduced because we now run over all pairs twice.

The last equation allows us to introduce some further definitions.
The single-particle wave functions ψμ(x), defined by the quantum numbers μ and x are defined as the overlap

ψα(x)=x|α.

We introduce the following shorthands for the above two integrals

μν|v^|μν=ψμ(xi)ψν(xj)v^(rij)ψμ(xi)ψν(xj)dxidxj,

and

μν|v^|νμ=ψμ(xi)ψν(xj)v^(rij)ψν(xi)ψμ(xj)dxidxj.

1.2. Preparing for later studies: varying the coefficients of a wave function expansion and orthogonal transformations

It is common to expand the single-particle functions in a known basis and vary the coefficients, that is, the new single-particle wave function is written as a linear expansion in terms of a fixed chosen orthogonal basis (for example the well-known harmonic oscillator functions or the hydrogen-like functions etc). We define our new single-particle basis (this is a normal approach for Hartree-Fock theory) by performing a unitary transformation on our previous basis (labelled with greek indices) as

(11)ψpnew=λCpλϕλ.

In this case we vary the coefficients Cpλ. If the basis has infinitely many solutions, we need to truncate the above sum. We assume that the basis ϕλ is orthogonal.

It is normal to choose a single-particle basis defined as the eigenfunctions of parts of the full Hamiltonian. The typical situation consists of the solutions of the one-body part of the Hamiltonian, that is we have

h^0ϕλ=ϵλϕλ.

The single-particle wave functions ϕλ(r), defined by the quantum numbers λ and r are defined as the overlap

ϕλ(r)=r|λ.

In deriving the Hartree-Fock equations, we will expand the single-particle functions in a known basis and vary the coefficients, that is, the new single-particle wave function is written as a linear expansion in terms of a fixed chosen orthogonal basis (for example the well-known harmonic oscillator functions or the hydrogen-like functions etc).

We stated that a unitary transformation keeps the orthogonality. To see this consider first a basis of vectors vi,

vi=[vi1vin]

We assume that the basis is orthogonal, that is

vjTvi=δij.

An orthogonal or unitary transformation

wi=Uvi,

preserves the dot product and orthogonality since

wjTwi=(Uvj)TUvi=vjTUTUvi=vjTvi=δij.

This means that if the coefficients Cpλ belong to a unitary or orthogonal trasformation (using the Dirac bra-ket notation)

|p=λCpλ|λ,

orthogonality is preserved, that is α|β=δαβ and p|q=δpq.

This propertry is extremely useful when we build up a basis of many-body Stater determinant based states.

Note also that although a basis |α contains an infinity of states, for practical calculations we have always to make some truncations.

Before we develop for example the Hartree-Fock equations, there is another very useful property of determinants that we will use both in connection with Hartree-Fock calculations and later shell-model calculations.

Consider the following determinant

|α1b11+α2sb12a12α1b21+α2b22a22|=α1|b11a12b21a22|+α2|b12a12b22a22|

We can generalize this to an n×n matrix and have

|a11a12k=1nckb1ka1na21a22k=1nckb2ka2nan1an2k=1nckbnkann|=k=1nck|a11a12b1ka1na21a22b2ka2nan1an2bnkann|.

This is a property we will use in our Hartree-Fock discussions.

We can generalize the previous results, now with all elements aij being given as functions of linear combinations of various coefficients c and elements bij,

|k=1nb1kck1k=1nb1kck2k=1nb1kckjk=1nb1kcknk=1nb2kck1k=1nb2kck2k=1nb2kckjk=1nb2kcknk=1nbnkck1k=1nbnkck2k=1nbnkckjk=1nbnkckn|=det(C)det(B),

where det(C) and det(B) are the determinants of n×n matrices with elements cij and bij respectively.
This is a property we will use in our Hartree-Fock discussions. Convince yourself about the correctness of the above expression by setting n=2.

With our definition of the new basis in terms of an orthogonal basis we have

ψp(x)=λCpλϕλ(x).

If the coefficients Cpλ belong to an orthogonal or unitary matrix, the new basis is also orthogonal. Our Slater determinant in the new basis ψp(x) is written as

1A!|ψp(x1)ψp(x2)ψp(xA)ψq(x1)ψq(x2)ψq(xA)ψt(x1)ψt(x2)ψt(xA)|=1A!|λCpλϕλ(x1)λCpλϕλ(x2)λCpλϕλ(xA)λCqλϕλ(x1)λCqλϕλ(x2)λCqλϕλ(xA)λCtλϕλ(x1)λCtλϕλ(x2)λCtλϕλ(xA)|,

which is nothing but det(C)det(Φ), with det(Φ) being the determinant given by the basis functions ϕλ(x).

In our discussions hereafter we will use our definitions of single-particle states above and below the Fermi (F) level given by the labels ijklF for so-called single-hole states and abcd>F for so-called particle states. For general single-particle states we employ the labels pqrs.

The energy functional is

E[Φ]=μ=1Aμ|h|μ+12μ=1Aν=1Aμν|v^|μνAS,

we found the expression for the energy functional in terms of the basis function ϕλ(r). We then varied the above energy functional with respect to the basis functions |μ. Now we are interested in defining a new basis defined in terms of a chosen basis as defined in Eq. (11). We can then rewrite the energy functional as

(12)E[ΦNew]=i=1Ai|h|i+12ij=1Aij|v^|ijAS,

where ΦNew is the new Slater determinant defined by the new basis of Eq. (11).

Using Eq. (11) we can rewrite Eq. (12) as

(13)E[Ψ]=i=1AαβCiαCiβα|h|β+12ij=1AαβγδCiαCjβCiγCjδαβ|v^|γδAS.

1.3. Second quantization

We introduce the time-independent operators aα and aα which create and annihilate, respectively, a particle in the single-particle state φα. We define the fermion creation operator aα

(14)aα|0|α,

and

(15)aα|α1αnAS|αα1αnAS

In Eq. (14) the operator aα acts on the vacuum state |0, which does not contain any particles. Alternatively, we could define a closed-shell nucleus or atom as our new vacuum, but then we need to introduce the particle-hole formalism, see the discussion to come.

In Eq. (15) aα acts on an antisymmetric n-particle state and creates an antisymmetric (n+1)-particle state, where the one-body state φα is occupied, under the condition that αα1,α2,,αn. It follows that we can express an antisymmetric state as the product of the creation operators acting on the vacuum state.

(16)|α1αnAS=aα1aα2aαn|0

It is easy to derive the commutation and anticommutation rules for the fermionic creation operators aα. Using the antisymmetry of the states (16)

(17)|α1αiαkαnAS=|α1αkαiαnAS

we obtain

(18)aαiaαk=aαkaαi

Using the Pauli principle

(19)|α1αiαiαnAS=0

it follows that

(20)aαiaαi=0.

If we combine Eqs. (18) and (20), we obtain the well-known anti-commutation rule

(21)aαaβ+aβaα{aα,aβ}=0

The hermitian conjugate of aα is

(22)aα=(aα)

If we take the hermitian conjugate of Eq. (21), we arrive at

(23){aα,aβ}=0

What is the physical interpretation of the operator aα and what is the effect of aα on a given state |α1α2αnAS? Consider the following matrix element

(24)α1α2αn|aα|α1α2αm

where both sides are antisymmetric. We distinguish between two cases. The first (1) is when α{αi}. Using the Pauli principle of Eq. (19) it follows

(25)α1α2αn|aα=0

The second (2) case is when α{αi}. It follows that an hermitian conjugation

(26)α1α2αn|aα=αα1α2αn|

Eq. (26) holds for case (1) since the lefthand side is zero due to the Pauli principle. We write Eq. (24) as

(27)α1α2αn|aα|α1α2αm=α1α2αn|αα1α2αm

Here we must have m=n+1 if Eq. (27) has to be trivially different from zero.

For the last case, the minus and plus signs apply when the sequence α,α1,α2,,αn and α1,α2,,αn+1 are related to each other via even and odd permutations. If we assume that α{αi} we obtain

(28)α1α2αn|aα|α1α2αn+1=0

when α{αi}. If α{αi}, we obtain

(29)aα|α1α2αn+1α=0

and in particular

(30)aα|0=0

If {ααi}={αi}, performing the right permutations, the sequence α,α1,α2,,αn is identical with the sequence α1,α2,,αn+1. This results in

(31)α1α2αn|aα|αα1α2αn=1

and thus

(32)aα|αα1α2αn=|α1α2αn

The action of the operator aα from the left on a state vector is to to remove one particle in the state α. If the state vector does not contain the single-particle state α, the outcome of the operation is zero. The operator aα is normally called for a destruction or annihilation operator.

The next step is to establish the commutator algebra of aα and aβ.

The action of the anti-commutator {aα,aα} on a given n-particle state is

aαaα|α1α2αnα=0
(33)aαaα|α1α2αnα=aα|αα1α2αnα=|α1α2αnα

if the single-particle state α is not contained in the state.

If it is present we arrive at

aαaα|α1α2αkααk+1αn1=aαaα(1)k|αα1α2αn1
=(1)k|αα1α2αn1=|α1α2αkααk+1αn1
(34)aαaα|α1α2αkααk+1αn1=0

From Eqs. (33) and (34) we arrive at

(35){aα,aα}=aαaα+aαaα=1

The action of {aα,aβ}, with αβ on a given state yields three possibilities. The first case is a state vector which contains both α and β, then either α or β and finally none of them.

The first case results in

aαaβ|αβα1α2αn2=0
(36)aβaα|αβα1α2αn2=0

while the second case gives

aαaβ|βα1α2αn1α=|αα1α2αn1α
aβaα|βα1α2αn1α=aβ|αββα1α2αn1α
(37)=|αα1α2αn1α

Finally if the state vector does not contain α and β

aαaβ|α1α2αnα,β=0
(38)aβaα|α1α2αnα,β=aβ|αα1α2αnα,β=0

For all three cases we have

(39){aα,aβ}=aαaβ+aβaα=0,αβ

We can summarize our findings in Eqs. (35) and (39) as

(40){aα,aβ}=δαβ

with δαβ is the Kroenecker δ-symbol.

The properties of the creation and annihilation operators can be summarized as (for fermions)

aα|0|α,

and

aα|α1αnAS|αα1αnAS.

from which follows

|α1αnAS=aα1aα2aαn|0.

The hermitian conjugate has the folowing properties

aα=(aα).

Finally we found

aα|α1α2αn+1α=0,in particular aα|0=0,

and

aα|αα1α2αn=|α1α2αn,

and the corresponding commutator algebra

{aα,aβ}={aα,aβ}=0{aα,aβ}=δαβ.

1.4. One-body operators in second quantization

A very useful operator is the so-called number-operator. Most physics cases we will study in this text conserve the total number of particles. The number operator is therefore a useful quantity which allows us to test that our many-body formalism conserves the number of particles. In for example (d,p) or (p,d) reactions it is important to be able to describe quantum mechanical states where particles get added or removed. A creation operator aα adds one particle to the single-particle state α of a give many-body state vector, while an annihilation operator aα removes a particle from a single-particle state α.

Let us consider an operator proportional with aαaβ and α=β. It acts on an n-particle state resulting in

(41)aαaα|α1α2αn={0α{αi}|α1α2αnα{αi}

Summing over all possible one-particle states we arrive at

(42)(αaαaα)|α1α2αn=n|α1α2αn

The operator

(43)N^=αaαaα

is called the number operator since it counts the number of particles in a give state vector when it acts on the different single-particle states. It acts on one single-particle state at the time and falls therefore under category one-body operators. Next we look at another important one-body operator, namely H^0 and study its operator form in the occupation number representation.

We want to obtain an expression for a one-body operator which conserves the number of particles. Here we study the one-body operator for the kinetic energy plus an eventual external one-body potential. The action of this operator on a particular n-body state with its pertinent expectation value has already been studied in coordinate space. In coordinate space the operator reads

(44)H^0=ih^0(xi)

and the anti-symmetric n-particle Slater determinant is defined as

Φ(x1,x2,,xn,α1,α2,,αn)=1n!p(1)pP^ψα1(x1)ψα2(x2)ψαn(xn).

Defining

(45)h^0(xi)ψαi(xi)=αkψαk(xi)αk|h^0|αk

we can easily evaluate the action of H^0 on each product of one-particle functions in Slater determinant. From Eq. (45) we obtain the following result without permuting any particle pair

(ih^0(xi))ψα1(x1)ψα2(x2)ψαn(xn)
=α1α1|h^0|α1ψα1(x1)ψα2(x2)ψαn(xn)
+α2α2|h^0|α2ψα1(x1)ψα2(x2)ψαn(xn)
+
(46)+αnαn|h^0|αnψα1(x1)ψα2(x2)ψαn(xn)

If we interchange particles 1 and 2 we obtain

(ih^0(xi))ψα1(x2)ψα1(x2)ψαn(xn)
=α2α2|h^0|α2ψα1(x2)ψα2(x1)ψαn(xn)
+α1α1|h^0|α1ψα1(x2)ψα2(x1)ψαn(xn)
+
(47)+αnαn|h^0|αnψα1(x2)ψα1(x2)ψαn(xn)

We can continue by computing all possible permutations. We rewrite also our Slater determinant in its second quantized form and skip the dependence on the quantum numbers xi. Summing up all contributions and taking care of all phases (1)p we arrive at

H^0|α1,α2,,αn=α1α1|h^0|α1|α1α2αn
+α2α2|h^0|α2|α1α2αn
+
(48)+αnαn|h^0|αn|α1α2αn

In Eq. (48) we have expressed the action of the one-body operator of Eq. (44) on the n-body state in its second quantized form. This equation can be further manipulated if we use the properties of the creation and annihilation operator on each primed quantum number, that is

(49)|α1α2αkαn=aαkaαk|α1α2αkαn

Inserting this in the right-hand side of Eq. (48) results in

H^0|α1α2αn=α1α1|h^0|α1aα1aα1|α1α2αn
+α2α2|h^0|α2aα2aα2|α1α2αn
+
+αnαn|h^0|αnaαnaαn|α1α2αn
(50)=α,βα|h^0|βaαaβ|α1α2αn

In the number occupation representation or second quantization we get the following expression for a one-body operator which conserves the number of particles

(51)H^0=αβα|h^0|βaαaβ

Obviously, H^0 can be replaced by any other one-body operator which preserved the number of particles. The stucture of the operator is therefore not limited to say the kinetic or single-particle energy only.

The opearator H^0 takes a particle from the single-particle state β to the single-particle state α with a probability for the transition given by the expectation value α|h^0|β.

It is instructive to verify Eq. (51) by computing the expectation value of H^0 between two single-particle states

(52)α1|h^0|α2=αβα|h^0|β0|aα1aαaβaα2|0

Using the commutation relations for the creation and annihilation operators we have

(53)aα1aαaβaα2=(δαα1aαaα1)(δβα2aα2aβ),

which results in

(54)0|aα1aαaβaα2|0=δαα1δβα2

and

(55)α1|h^0|α2=αβα|h^0|βδαα1δβα2=α1|h^0|α2

1.5. Two-body operators in second quantization

Let us now derive the expression for our two-body interaction part, which also conserves the number of particles. We can proceed in exactly the same way as for the one-body operator. In the coordinate representation our two-body interaction part takes the following expression

(56)H^I=i<jV(xi,xj)

where the summation runs over distinct pairs. The term V can be an interaction model for the nucleon-nucleon interaction or the interaction between two electrons. It can also include additional two-body interaction terms.

The action of this operator on a product of two single-particle functions is defined as

(57)V(xi,xj)ψαk(xi)ψαl(xj)=αkαlψαk(xi)ψαl(xj)αkαl|v^|αkαl

We can now let H^I act on all terms in the linear combination for |α1α2αn. Without any permutations we have

(i<jV(xi,xj))ψα1(x1)ψα2(x2)ψαn(xn)
=α1α2α1α2|v^|α1α2ψα1(x1)ψα2(x2)ψαn(xn)
+
+α1αnα1αn|v^|α1αnψα1(x1)ψα2(x2)ψαn(xn)
+
+α2αnα2αn|v^|α2αnψα1(x1)ψα2(x2)ψαn(xn)
(58)+

where on the rhs we have a term for each distinct pairs.

For the other terms on the rhs we obtain similar expressions and summing over all terms we obtain

HI|α1α2αn=α1,α2α1α2|v^|α1α2|α1α2αn
+
+α1,αnα1αn|v^|α1αn|α1α2αn
+
+α2,αnα2αn|v^|α2αn|α1α2αn
(59)+

We introduce second quantization via the relation

aαkaαlaαlaαk|α1α2αkαlαn
=(1)k1(1)l2aαkaαlaαlaαk|αkαlα1α2αnαk,αl
=(1)k1(1)l2|αkαlα1α2αnαk,αl
(60)=|α1α2αkαlαn

Inserting this in (59) gives

HI|α1α2αn=α1,α2α1α2|v^|α1α2aα1aα2aα2aα1|α1α2αn
+
=α1,αnα1αn|v^|α1αnaα1aαnaαnaα1|α1α2αn
+
=α2,αnα2αn|v^|α2αnaα2aαnaαnaα2|α1α2αn
+
(61)=α,β,γ,δαβ|v^|γδaαaβaδaγ|α1α2αn

Here we let indicate that the sums running over α and β run over all single-particle states, while the summations γ and δ run over all pairs of single-particle states. We wish to remove this restriction and since

(62)αβ|v^|γδ=βα|v^|δγ

we get

(63)αβαβ|v^|γδaαaβaδaγ=αββα|v^|δγaαaβaδaγ
(64)=αββα|v^|δγaβaαaγaδ

where we have used the anti-commutation rules.

Changing the summation indices α and β in (64) we obtain

(65)αβαβ|v^|γδaαaβaδaγ=αβαβ|v^|δγaαaβaγaδ

From this it follows that the restriction on the summation over γ and δ can be removed if we multiply with a factor 12, resulting in

(66)H^I=12αβγδαβ|v^|γδaαaβaδaγ

where we sum freely over all single-particle states α, β, γ og δ.

With this expression we can now verify that the second quantization form of H^I in Eq. (66) results in the same matrix between two anti-symmetrized two-particle states as its corresponding coordinate space representation. We have

(67)α1α2|H^I|β1β2=12αβγδαβ|v^|γδ0|aα2aα1aαaβaδaγaβ1aβ2|0.

Using the commutation relations we get

aα2aα1aαaβaδaγaβ1aβ2
=aα2aα1aαaβ(aδδγβ1aβ2aδaβ1aγaβ2)
=aα2aα1aαaβ(δγβ1δδβ2δγβ1aβ2aδaδaβ1δγβ2+aδaβ1aβ2aγ)
=aα2aα1aαaβ(δγβ1δδβ2δγβ1aβ2aδ
(68)δδβ1δγβ2+δγβ2aβ1aδ+aδaβ1aβ2aγ)

The vacuum expectation value of this product of operators becomes

0|aα2aα1aαaβaδaγaβ1aβ2|0
=(δγβ1δδβ2δδβ1δγβ2)0|aα2aα1aαaβ|0
(69)=(δγβ1δδβ2δδβ1δγβ2)(δαα1δβα2δβα1δαα2)

Insertion of Eq. (69) in Eq. (67) results in

α1α2|H^I|β1β2=12[α1α2|v^|β1β2α1α2|v^|β2β1
α2α1|v^|β1β2+α2α1|v^|β2β1]
=α1α2|v^|β1β2α1α2|v^|β2β1
(70)=α1α2|v^|β1β2AS.

The two-body operator can also be expressed in terms of the anti-symmetrized matrix elements we discussed previously as

H^I=12αβγδαβ|v^|γδaαaβaδaγ
=14αβγδ[αβ|v^|γδαβ|v^|δγ]aαaβaδaγ
(71)=14αβγδαβ|v^|γδASaαaβaδaγ

The factors in front of the operator, either 14 or 12 tells whether we use antisymmetrized matrix elements or not.

We can now express the Hamiltonian operator for a many-fermion system in the occupation basis representation as

(72)H=α,βα|t^+u^ext|βaαaβ+14αβγδαβ|v^|γδaαaβaδaγ.

This is the form we will use in the rest of these lectures, assuming that we work with anti-symmetrized two-body matrix elements.

1.6. Particle-hole formalism

Second quantization is a useful and elegant formalism for constructing many-body states and quantum mechanical operators. One can express and translate many physical processes into simple pictures such as Feynman diagrams. Expecation values of many-body states are also easily calculated. However, although the equations are seemingly easy to set up, from a practical point of view, that is the solution of Schroedinger’s equation, there is no particular gain. The many-body equation is equally hard to solve, irrespective of representation. The cliche that there is no free lunch brings us down to earth again.
Note however that a transformation to a particular basis, for cases where the interaction obeys specific symmetries, can ease the solution of Schroedinger’s equation.

But there is at least one important case where second quantization comes to our rescue. It is namely easy to introduce another reference state than the pure vacuum |0, where all single-particle states are active. With many particles present it is often useful to introduce another reference state than the vacuum state|0. We will label this state |c (c for core) and as we will see it can reduce considerably the complexity and thereby the dimensionality of the many-body problem. It allows us to sum up to infinite order specific many-body correlations. The particle-hole representation is one of these handy representations.

In the original particle representation these states are products of the creation operators aαi acting on the true vacuum |0. Following Eq. (16) we have

(73)|α1α2αn1αn=aα1aα2aαn1aαn|0
(74)|α1α2αn1αnαn+1=aα1aα2aαn1aαnaαn+1|0
(75)|α1α2αn1=aα1aα2aαn1|0

If we use Eq. (73) as our new reference state, we can simplify considerably the representation of this state

(76)|c|α1α2αn1αn=aα1aα2aαn1aαn|0

The new reference states for the n+1 and n1 states can then be written as

(77)|α1α2αn1αnαn+1=(1)naαn+1|c(1)n|αn+1c
(78)|α1α2αn1=(1)n1aαn|c(1)n1|αn1c

The first state has one additional particle with respect to the new vacuum state |c and is normally referred to as a one-particle state or one particle added to the many-body reference state. The second state has one particle less than the reference vacuum state |c and is referred to as a one-hole state. When dealing with a new reference state it is often convenient to introduce new creation and annihilation operators since we have from Eq. (78)

(79)aα|c0

since α is contained in |c, while for the true vacuum we have aα|0=0 for all α.

The new reference state leads to the definition of new creation and annihilation operators which satisfy the following relations

(80)bα|c=0
{bα,bβ}={bα,bβ}=0
(81){bα,bβ}=δαβ

We assume also that the new reference state is properly normalized

(82)c|c=1

The physical interpretation of these new operators is that of so-called quasiparticle states. This means that a state defined by the addition of one extra particle to a reference state |c may not necesseraly be interpreted as one particle coupled to a core. We define now new creation operators that act on a state α creating a new quasiparticle state

(83)bα|c={aα|c=|α,α>Faα|c=|α1,αF

where F is the Fermi level representing the last occupied single-particle orbit of the new reference state |c.

The annihilation is the hermitian conjugate of the creation operator

bα=(bα),

resulting in

(84)bα={aαα>FaααFbα={aαα>FaααF

With the new creation and annihilation operator we can now construct many-body quasiparticle states, with one-particle-one-hole states, two-particle-two-hole states etc in the same fashion as we previously constructed many-particle states. We can write a general particle-hole state as

(85)|β1β2βnpγ11γ21γnh1bβ1bβ2bβnp>Fbγ1bγ2bγnhF|c

We can now rewrite our one-body and two-body operators in terms of the new creation and annihilation operators. The number operator becomes

(86)N^=αaαaα=α>Fbαbα+ncαFbαbα

where nc is the number of particle in the new vacuum state |c.
The action of N^ on a many-body state results in

(87)N|β1β2βnpγ11γ21γnh1=(np+ncnh)|β1β2βnpγ11γ21γnh1

Here n=np+ncnh is the total number of particles in the quasi-particle state of Eq. (85). Note that N^ counts the total number of particles present

(88)Nqp=αbαbα,

gives us the number of quasi-particles as can be seen by computing

(89)Nqp=|β1β2βnpγ11γ21γnh1=(np+nh)|β1β2βnpγ11γ21γnh1

where nqp=np+nh is the total number of quasi-particles.

We express the one-body operator H^0 in terms of the quasi-particle creation and annihilation operators, resulting in

H^0=αβ>Fα|h^0|βbαbβ+α>F,βF[α|h^0|βbαbβ+β|h^0|αbβbα]
(90)+αFα|h^0|ααβFβ|h^0|αbαbβ

The first term gives contribution only for particle states, while the last one contributes only for holestates. The second term can create or destroy a set of quasi-particles and the third term is the contribution from the vacuum state |c.

Before we continue with the expressions for the two-body operator, we introduce a nomenclature we will use for the rest of this text. It is inspired by the notation used in quantum chemistry. We reserve the labels i,j,k, for hole states and a,b,c, for states above F, viz. particle states. This means also that we will skip the constraint F or >F in the summation symbols. Our operator H^0 reads now

H^0=aba|h^|bbabb+ai[a|h^|ibabi+i|h^|abiba]
(91)+ii|h^|iijj|h^|ibibj

The two-particle operator in the particle-hole formalism is more complicated since we have to translate four indices αβγδ to the possible combinations of particle and hole states. When performing the commutator algebra we can regroup the operator in five different terms

(92)H^I=H^I(a)+H^I(b)+H^I(c)+H^I(d)+H^I(e)

Using anti-symmetrized matrix elements, bthe term H^I(a) is

(93)H^I(a)=14abcdab|V^|cdbabbbdbc

The next term H^I(b) reads

(94)H^I(b)=14abci(ab|V^|cibabbbibc+ai|V^|cbbabibbbc)

This term conserves the number of quasiparticles but creates or removes a three-particle-one-hole state. For H^I(c) we have

H^I(c)=14abij(ab|V^|ijbabbbjbi+ij|V^|abbabbbjbi)+
(95)12abijai|V^|bjbabjbbbi+12abiai|V^|bibabb.

The first line stands for the creation of a two-particle-two-hole state, while the second line represents the creation to two one-particle-one-hole pairs while the last term represents a contribution to the particle single-particle energy from the hole states, that is an interaction between the particle states and the hole states within the new vacuum state. The fourth term reads

H^I(d)=14aijk(ai|V^|jkbabkbjbi+ji|V^|akbkbjbiba)+
(96)14aij(ai|V^|jibabj+ji|V^|aiji|V^|iabjba).

The terms in the first line stand for the creation of a particle-hole state interacting with hole states, we will label this as a two-hole-one-particle contribution. The remaining terms are a particle-hole state interacting with the holes in the vacuum state. Finally we have

(97)H^I(e)=14ijklkl|V^|ijbibjblbk+12ijkij|V^|kjbkbi+12ijij|V^|ij

The first terms represents the interaction between two holes while the second stands for the interaction between a hole and the remaining holes in the vacuum state. It represents a contribution to single-hole energy to first order. The last term collects all contributions to the energy of the ground state of a closed-shell system arising from hole-hole correlations.

1.7. Summarizing and defining a normal-ordered Hamiltonian

ΦAS(α1,,αA;x1,xA)=1AP^(1)PP^i=1Aψαi(xi),

which is equivalent with |α1αA=aα1aαA|0. We have also

2 0 3

< < < ! ! M A T H _ B L O C K

δpq={ap,aq},

and

2 0 5

< < < ! ! M A T H _ B L O C K

|Φ0=|α1αA,α1,,αAαF

2 0 7

< < < ! ! M A T H _ B L O C K

{ap,aq}=δpq,p,q>αF

with i,j,αF,a,b,>αF,p,q,any

ai|Φ0=|Φi,aa|Φ0=|Φa

and

ai|Φ0=0aa|Φ0=0

The one-body operator is defined as

F^=pqp|f^|qapaq

while the two-body opreator is defined as

V^=14pqrspq|v^|rsASapaqasar

where we have defined the antisymmetric matrix elements

pq|v^|rsAS=pq|v^|rspq|v^|sr.

We can also define a three-body operator

V^3=136pqrstupqr|v^3|stuASapaqarauatas

with the antisymmetrized matrix element

(98)pqr|v^3|stuAS=pqr|v^3|stu+pqr|v^3|tus+pqr|v^3|ustpqr|v^3|sutpqr|v^3|tsupqr|v^3|uts.

1.8. Operators in second quantization

In the build-up of a shell-model or FCI code that is meant to tackle large dimensionalities is the action of the Hamiltonian H^ on a Slater determinant represented in second quantization as

|α1αn=aα1aα2aαn|0.

The time consuming part stems from the action of the Hamiltonian on the above determinant,

(αβα|t+u|βaαaβ+14αβγδαβ|v^|γδaαaβaδaγ)aα1aα2aαn|0.

A practically useful way to implement this action is to encode a Slater determinant as a bit pattern.

Assume that we have at our disposal n different single-particle orbits α0,α2,,αn1 and that we can distribute among these orbits Nn particles.

A Slater determinant can then be coded as an integer of n bits. As an example, if we have n=16 single-particle states α0,α1,,α15 and N=4 fermions occupying the states α3, α6, α10 and α13 we could write this Slater determinant as

ΦΛ=aα3aα6aα10aα13|0.

The unoccupied single-particle states have bit value 0 while the occupied ones are represented by bit state 1. In the binary notation we would write this 16 bits long integer as

α0α1α2α3α4α5α6α7α8α9α10α11α12α13α14α150001001000100100

which translates into the decimal number

23+26+210+213=9288.

We can thus encode a Slater determinant as a bit pattern.

With N particles that can be distributed over n single-particle states, the total number of Slater determinats (and defining thereby the dimensionality of the system) is

dim(H)=(nN).

The total number of bit patterns is 2n.

We assume again that we have at our disposal n different single-particle orbits α0,α2,,αn1 and that we can distribute among these orbits Nn particles. The ordering among these states is important as it defines the order of the creation operators. We will write the determinant

ΦΛ=aα3aα6aα10aα13|0,

in a more compact way as

Φ3,6,10,13=|0001001000100100.

The action of a creation operator is thus

aα4Φ3,6,10,13=aα4|0001001000100100=aα4aα3aα6aα10aα13|0,

which becomes

aα3aα4aα6aα10aα13|0=|0001101000100100.

Similarly

aα6Φ3,6,10,13=aα6|0001001000100100=aα6aα3aα6aα10aα13|0,

which becomes

aα4(aα6)2aα10aα13|0=0!

This gives a simple recipe:

  • If one of the bits bj is 1 and we act with a creation operator on this bit, we return a null vector

  • If bj=0, we set it to 1 and return a sign factor (1)l, where l is the number of bits set before bit j.

Consider the action of aα2 on various slater determinants:

aα2Φ00111=aα2|00111=0×|00111aα2Φ01011=aα2|01011=(1)×|01111aα2Φ01101=aα2|01101=0×|01101aα2Φ01110=aα2|01110=0×|01110aα2Φ10011=aα2|10011=(1)×|10111aα2Φ10101=aα2|10101=0×|10101aα2Φ10110=aα2|10110=0×|10110aα2Φ11001=aα2|11001=(+1)×|11101aα2Φ11010=aα2|11010=(+1)×|11110

What is the simplest way to obtain the phase when we act with one annihilation(creation) operator on the given Slater determinant representation?

We have an SD representation

ΦΛ=aα0aα3aα6aα10aα13|0,

in a more compact way as

Φ0,3,6,10,13=|1001001000100100.

The action of

aα4aα0Φ0,3,6,10,13=aα4|0001001000100100=aα4aα3aα6aα10aα13|0,

which becomes

aα3aα4aα6aα10aα13|0=|0001101000100100.

The action

aα0Φ0,3,6,10,13=|0001001000100100,

can be obtained by subtracting the logical sum (AND operation) of Φ0,3,6,10,13 and a word which represents only α0, that is

|1000000000000000,

from Φ0,3,6,10,13=|1001001000100100.

This operation gives |0001001000100100.

Similarly, we can form aα4aα0Φ0,3,6,10,13, say, by adding |0000100000000000 to aα0Φ0,3,6,10,13, first checking that their logical sum is zero in order to make sure that orbital α4 is not already occupied.

It is trickier however to get the phase (1)l. One possibility is as follows

  • Let S1 be a word that represents the 1bit to be removed and all others set to zero.

In the previous example S1=|1000000000000000

  • Define S2 as the similar word that represents the bit to be added, that is in our case

S2=|0000100000000000.

  • Compute then S=S1S2, which here becomes

S=|0111000000000000
  • Perform then the logical AND operation of S with the word containing

Φ0,3,6,10,13=|1001001000100100,

which results in |0001000000000000. Counting the number of 1bits gives the phase. Here you need however an algorithm for bitcounting. Several efficient ones available.

1.9. Exercise 1: Relation between basis functions

This exercise serves to convince you about the relation between two different single-particle bases, where one could be our new Hartree-Fock basis and the other a harmonic oscillator basis.

Consider a Slater determinant built up of single-particle orbitals ψλ, with λ=1,2,,A. The unitary transformation

ψa=λCaλϕλ,

brings us into the new basis.
The new basis has quantum numbers a=1,2,,A. Show that the new basis is orthonormal. Show that the new Slater determinant constructed from the new single-particle wave functions can be written as the determinant based on the previous basis and the determinant of the matrix C. Show that the old and the new Slater determinants are equal up to a complex constant with absolute value unity. (Hint, C is a unitary matrix).

Starting with the second quantization representation of the Slater determinant

Φ0=i=1naαi|0,

use Wick’s theorem to compute the normalization integral Φ0|Φ0.

1.10. Exercise 2: Matrix elements

Calculate the matrix elements

α1α2|F^|α1α2

and

α1α2|G^|α1α2

with

2 4 1

< < < ! ! M A T H _ B L O C K

2 4 2

< < < ! ! M A T H _ B L O C K

2 4 3

< < < ! ! M A T H _ B L O C K

G^=12αβγδαβ|g^|γδaαaβaδaγ,

and

αβ|g^|γδ=ψα(x1)ψβ(x2)g(x1,x2)ψγ(x1)ψδ(x2)dx1dx2

Compare these results with those from exercise 3c).

1.11. Exercise 3: Normal-ordered one-body operator

Show that the onebody part of the Hamiltonian

H^0=pqp|h^0|qapaq,

can be written, using standard annihilation and creation operators, in normal-ordered form as

H^0=pqp|h^0|q{apaq}+ii|h^0|i.

Explain the meaning of the various symbols. Which reference vacuum has been used?

1.12. Exercise 4: Normal-ordered two-body operator

Show that the twobody part of the Hamiltonian

H^I=14pqrspq|v^|rsapaqasar,

can be written, using standard annihilation and creation operators, in normal-ordered form as

H^I=14pqrspq|v^|rs{apaqasar}+pqipi|v^|qi{apaq}+12ijij|v^|ij.

Explain again the meaning of the various symbols.

This exercise is optional: Derive the normal-ordered form of the threebody part of the Hamiltonian.

H^3=136pqrstupqr|v^3|stuapaqarauatas,

and specify the contributions to the twobody, onebody and the scalar part.

1.13. Exercise 5: Matrix elements using the Slater-Condon rule

The aim of this exercise is to set up specific matrix elements that will turn useful when we start our discussions of the nuclear shell model. In particular you will notice, depending on the character of the operator, that many matrix elements will actually be zero.

Consider three N-particle Slater determinants |Φ0, |Φia and |Φijab, where the notation means that Slater determinant |Φia differs from |Φ0 by one single-particle state, that is a single-particle state ψi is replaced by a single-particle state ψa. It is often interpreted as a so-called one-particle-one-hole excitation. Similarly, the Slater determinant |Φijab differs by two single-particle states from |Φ0 and is normally thought of as a two-particle-two-hole excitation. We assume also that |Φ0 represents our new vacuum reference state and the labels ijk represent single-particle states below the Fermi level and abc represent states above the Fermi level, so-called particle states. We define thereafter a general onebody normal-ordered (with respect to the new vacuum state) operator as

F^N=pqp|f|β{apaq},

with

p|f|q=ψp(x)f(x)ψq(x)dx,

and a general normal-ordered two-body operator

G^N=14pqrspq|g|rsAS{apaqasar},

with for example the direct matrix element given as

pq|g|rs=ψp(x1)ψq(x2)g(x1,x2)ψr(x1)ψs(x2)dx1dx2

with g being invariant under the interchange of the coordinates of two particles. The single-particle states ψi are not necessarily eigenstates of f^. The curly brackets mean that the operators are normal-ordered with respect to the new vacuum reference state.

How would you write the above Slater determinants in a second quantization formalism, utilizing the fact that we have defined a new reference state?

Use thereafter Wick’s theorem to find the expectation values of

Φ0|F^N|Φ0,

and

Φ0G^N|Φ0.

Find thereafter

Φ0|F^N|Φia,

and

Φ0|G^N|Φia,

Finally, find

Φ0|F^N|Φijab,

and

Φ0|G^N|Φijab.

What happens with the two-body operator if we have a transition probability of the type

Φ0|G^N|Φijkabc,

where the Slater determinant to the right of the operator differs by more than two single-particle states?

1.14. Exercise 6: Program to set up Slater determinants

Write a program which sets up all possible Slater determinants given N=4 eletrons which can occupy the atomic single-particle states defined by the 1s, 2s2p and 3s3p3d shells. How many single-particle states n are there in total? Include the spin degrees of freedom as well.

1.15. Exercise 7: Using sympy to compute matrix elements

Compute the matrix element

α1α2α3|G^|α1α2α3,

using Wick’s theorem and express the two-body operator G in the occupation number (second quantization) representation.

1.16. Exercise 8: Using sympy to compute matrix elements

The last exercise can be solved using the symbolic Python package called SymPy. SymPy is a Python package for general purpose symbolic algebra. There is a physics module with several interesting submodules. Among these, the submodule called secondquant, contains several functionalities that allow us to test our algebraic manipulations using Wick’s theorem and operators for second quantization.

from sympy import *
from sympy.physics.secondquant import *

i, j = symbols('i,j', below_fermi=True)
a, b = symbols('a,b', above_fermi=True)
p, q = symbols('p,q')
print simplify(wicks(Fd(i)*F(a)*Fd(p)*F(q)*Fd(b)*F(j), keep_only_fully_contracted=True))
  Input In [1]
    print simplify(wicks(Fd(i)*F(a)*Fd(p)*F(q)*Fd(b)*F(j), keep_only_fully_contracted=True))
          ^
SyntaxError: invalid syntax

The code defines single-particle states above and below the Fermi level, in addition to the genereal symbols pq which can refer to any type of state below or above the Fermi level. Wick’s theorem is implemented between the creation and annihilation operators Fd and F, respectively. Using the simplify option, one can lump together several Kronecker-δ functions.

1.17. Exercise 9: Using sympy to compute matrix elements

We can expand the above Python code by defining one-body and two-body operators using the following SymPy code

# This code sets up a two-body Hamiltonian for fermions
from sympy import symbols, latex, WildFunction, collect, Rational
from sympy.physics.secondquant import F, Fd, wicks, AntiSymmetricTensor, substitute_dummies, NO

# setup hamiltonian
p,q,r,s = symbols('p q r s',dummy=True)
f = AntiSymmetricTensor('f',(p,),(q,))
pr = NO((Fd(p)*F(q)))
v = AntiSymmetricTensor('v',(p,q),(r,s))
pqsr = NO(Fd(p)*Fd(q)*F(s)*F(r))
Hamiltonian=f*pr + Rational(1)/Rational(4)*v*pqsr
print "Hamiltonian defined as:", latex(Hamiltonian)

Here we have used the AntiSymmetricTensor functionality, together with normal-ordering defined by the NO function. Using the latex option, this program produces the following output

fqp{apaq}14vsrqp{apaqaras}

1.18. Exercise 10: Using sympy to compute matrix elements

We can now use this code to compute the matrix elements between two two-body Slater determinants using Wick’s theorem.

from sympy import symbols, latex, WildFunction, collect, Rational, simplify
from sympy.physics.secondquant import F, Fd, wicks, AntiSymmetricTensor, substitute_dummies, NO, evaluate_deltas
# setup hamiltonian
p,q,r,s = symbols('p q r s',dummy=True)
f = AntiSymmetricTensor('f',(p,),(q,))
pr = NO((Fd(p)*F(q)))
v = AntiSymmetricTensor('v',(p,q),(r,s))
pqsr = NO(Fd(p)*Fd(q)*F(s)*F(r))
Hamiltonian=f*pr + Rational(1)/Rational(4)*v*pqsr
c,d = symbols('c, d',above_fermi=True)
a,b = symbols('a, b',above_fermi=True)

expression = wicks(F(b)*F(a)*Hamiltonian*Fd(c)*Fd(d),keep_only_fully_contracted=True, simplify_kronecker_deltas=True)
expression = evaluate_deltas(expression)
expression = simplify(expression)
print "Hamiltonian defined as:", latex(expression)

The result is as expected,

δacfdbδadfcbδbcfda+δbdfca+vcdab.

1.19. Exercise 11: Using sympy to compute matrix elements

We can continue along these lines and define a normal-ordered Hamiltonian with respect to a given reference state. In our first step we just define the Hamiltonian

from sympy import symbols, latex, WildFunction, collect, Rational, simplify
from sympy.physics.secondquant import F, Fd, wicks, AntiSymmetricTensor, substitute_dummies, NO, evaluate_deltas
# setup hamiltonian
p,q,r,s = symbols('p q r s',dummy=True)
f = AntiSymmetricTensor('f',(p,),(q,))
pr = Fd(p)*F(q)
v = AntiSymmetricTensor('v',(p,q),(r,s))
pqsr = Fd(p)*Fd(q)*F(s)*F(r)
#define the Hamiltonian
Hamiltonian = f*pr + Rational(1)/Rational(4)*v*pqsr
#define indices for states above and below the Fermi level
index_rule = {
     'below':  'kl',
     'above':  'cd',
     'general': 'pqrs'
     }
Hnormal = substitute_dummies(Hamiltonian,new_indices=True, pretty_indices=index_rule)
print "Hamiltonian defined as:", latex(Hnormal)

which results in

fpqaqap+14vqpsrasarapaq

1.20. Exercise 12: Using sympy to compute matrix elements

In our next step we define the reference energy E0 and redefine the Hamiltonian by subtracting the reference energy and collecting the coefficients for all normal-ordered products (given by the NO function).

from sympy import symbols, latex, WildFunction, collect, Rational, simplify
from sympy.physics.secondquant import F, Fd, wicks, AntiSymmetricTensor, substitute_dummies, NO, evaluate_deltas
# setup hamiltonian
p,q,r,s = symbols('p q r s',dummy=True)
f = AntiSymmetricTensor('f',(p,),(q,))
pr = Fd(p)*F(q)
v = AntiSymmetricTensor('v',(p,q),(r,s))
pqsr = Fd(p)*Fd(q)*F(s)*F(r)
#define the Hamiltonian
Hamiltonian=f*pr + Rational(1)/Rational(4)*v*pqsr
#define indices for states above and below the Fermi level
index_rule = {
     'below':  'kl',
     'above':  'cd',
     'general': 'pqrs'
     }
Hnormal = substitute_dummies(Hamiltonian,new_indices=True, pretty_indices=index_rule)
E0 = wicks(Hnormal,keep_only_fully_contracted=True)
Hnormal = Hnormal-E0
w = WildFunction('w')
Hnormal = collect(Hnormal, NO(w))
Hnormal = evaluate_deltas(Hnormal)
print latex(Hnormal)

which gives us

fii+fpqaqap14viiii14viiii+14vqpsrarasaqap,

again as expected, with the reference energy to be subtracted.

1.21. Exercise 13: Using sympy to compute matrix elements

We can now go back to exercise 7 and define the Hamiltonian and the second-quantized representation of a three-body Slater determinant.

from sympy import symbols, latex, WildFunction, collect, Rational, simplify
from sympy.physics.secondquant import F, Fd, wicks, AntiSymmetricTensor, substitute_dummies, NO, evaluate_deltas
# setup hamiltonian
p,q,r,s = symbols('p q r s',dummy=True)

v = AntiSymmetricTensor('v',(p,q),(r,s))
pqsr = NO(Fd(p)*Fd(q)*F(s)*F(r))
Hamiltonian=Rational(1)/Rational(4)*v*pqsr
a,b,c,d,e,f = symbols('a,b, c, d, e, f',above_fermi=True)

expression = wicks(F(c)*F(b)*F(a)*Hamiltonian*Fd(d)*Fd(e)*Fd(f),keep_only_fully_contracted=True, simplify_kronecker_deltas=True)
expression = evaluate_deltas(expression)
expression = simplify(expression)
print latex(expression)

resulting in nine terms (as expected),

δadvefcbδaevfdcb+δafvedcbδbdvefacδbevfdac+δbfvedac+δcdvefab+δcevfdabδcfvedab

1.22. Exercise 14: Diagrammatic representation of Hartree-Fock equations

What is the diagrammatic representation of the HF equation?

αk|uHF|αi+j=1n[αkαj|v^|αiαjαkαj|v|αjαi]=0

(Represent (uHF) by the symbol X .)

1.23. Exercise 15: Derivation of Hartree-Fock equations

Consider the ground state |Φ of a bound many-particle system of fermions. Assume that we remove one particle from the single-particle state λ and that our system ends in a new state |Φn. Define the energy needed to remove this particle as

Eλ=n|Φn|aλ|Φ|2(E0En),

where E0 and En are the ground state energies of the states |Φ and |Φn, respectively.

  • Show that

Eλ=Φ|aλ[aλ,H]|Φ,

where H is the Hamiltonian of this system.

  • If we assume that Φ is the Hartree-Fock result, find the

relation between Eλ and the single-particle energy ελ for states λF and λ>F, with

ελ=λ|t^+u^|λ,

and

λ|u^|λ=βFλβ|v^|λβ.

We have assumed an antisymmetrized matrix element here. Discuss the result.

The Hamiltonian operator is defined as

H=αβα|t^|βaαaβ+12αβγδαβ|v^|γδaαaβaδaγ.