- #1

- 624

- 14

- Thread starter LagrangeEuler
- Start date

- #1

- 624

- 14

- #2

DrDu

Science Advisor

- 6,092

- 803

- #3

- 624

- 14

I can't get trace in case without circular boundary condition.

- #4

DrDu

Science Advisor

- 6,092

- 803

So what? The term <S_1|T|S_N> is missing. Ok.I can't get trace in case without circular boundary condition.

You get [itex] \sum_{S_1} \sum_{S_N} \langle S_1| T^{N-1} |S_N\rangle [/itex].

What happens when you decompose |S_1/N> into eigenfunctions of T?

- #5

- 624

- 14

I'm not quite sure what you think? One eigenvector is just ##|+1\rangle##.

- #6

- 624

- 14

\langle +|\hat{T}|+ \rangle & \langle +|\hat{T}|- \rangle \\

\langle -|\hat{T}|+ \rangle & \langle -|\hat{T}|- \rangle \end{array} \right)[/tex]

For ##|S_1/_N\rangle## there are two possibilities ##|+\rangle## or ##|-\rangle##.

But I don't know how to write that in some normal way and calculate.

- #7

DrDu

Science Advisor

- 6,092

- 803

- #8

- 624

- 14

[tex] \left( \begin{array}{ccc}

e^{j+h} & e^{-j} \\

e^{-j} & e^{j-h}

\end{array} \right) [/tex]

eigenvectors are

[tex] \left( \begin{array}{ccc}

\frac{1}{2}e^{-h}(-e^{2j}+e^{2h+2j}\pm \sqrt{4e^{2h}+e^{4j}-2e^{2h+4j}+e^{4h+4j}}) \\

1

\end{array} \right) [/tex]

- #9

DrDu

Science Advisor

- 6,092

- 803

What I meant is that you can write

[itex] U^T T U=\Lambda[/itex] with [itex]\Lambda[/itex] being a diagonal matrix of the eigenvalues [itex]\lambda_1[/itex] and [itex] \lambda_2[/itex]. I assume [itex]\lambda_1>\lambda_2[/itex] and call V_0 the matrix V with j=0.

Hence [itex]Z=Tr(V^{N-1}V_0)=Tr(U^T\Lambda^{N-1}UV_0)=Tr(\Lambda^{N-1} UV_0U^T)\approx

\lambda_1^{N-1} (UV_0U^T)_{11}[/itex].

Hence Z differs from the Z for periodic boundary conditions only by the factor [itex] (UV_0U^T)_{11}/\lambda_1[/itex]. As we are only interested in Log Z, this makes a term which is 1/N times smaller than Z. This is of the same order as the error resulting from the replacement of [itex]\Lambda^{N-1}[/itex] by [itex]\lambda_1^{N-1}[/itex].

[itex] U^T T U=\Lambda[/itex] with [itex]\Lambda[/itex] being a diagonal matrix of the eigenvalues [itex]\lambda_1[/itex] and [itex] \lambda_2[/itex]. I assume [itex]\lambda_1>\lambda_2[/itex] and call V_0 the matrix V with j=0.

Hence [itex]Z=Tr(V^{N-1}V_0)=Tr(U^T\Lambda^{N-1}UV_0)=Tr(\Lambda^{N-1} UV_0U^T)\approx

\lambda_1^{N-1} (UV_0U^T)_{11}[/itex].

Hence Z differs from the Z for periodic boundary conditions only by the factor [itex] (UV_0U^T)_{11}/\lambda_1[/itex]. As we are only interested in Log Z, this makes a term which is 1/N times smaller than Z. This is of the same order as the error resulting from the replacement of [itex]\Lambda^{N-1}[/itex] by [itex]\lambda_1^{N-1}[/itex].

Last edited:

- Replies
- 1

- Views
- 1K

- Replies
- 30

- Views
- 10K

- Last Post

- Replies
- 16

- Views
- 2K

- Last Post

- Replies
- 1

- Views
- 832

- Replies
- 0

- Views
- 1K

- Replies
- 5

- Views
- 3K

- Replies
- 0

- Views
- 2K

- Replies
- 1

- Views
- 2K

- Replies
- 0

- Views
- 1K

- Replies
- 11

- Views
- 4K