- #1
Luqman Saleem
- 18
- 3
- TL;DR Summary
- As an example, calculating the ground state energy and wave function using canonical formalism of a simple fermionic model.
For my own understanding, I am trying to computationally solve a simple spinless fermionic Hamiltonian in Quantum Canonical Ensemble formalism . The Hamiltonian is written in the second quantization as
$$H = \sum_{i=1}^L c_{i+1}^\dagger c_i + h.c.$$
In the canonical formalism, the density matrix is given as
$$\rho = \frac{e^{-\beta H}}{Tr [e^{-\beta H}]}$$
where ##\beta = 1/k_B T## and the expectation value of any operator ##A## is given as
$$\langle A \rangle = Tr[\rho A]$$
Things that I want to calculate:
1. ground state energy at temperature ##T##
2. ground state wave-function at ##T##
My attempt to solution:
1. I first write ##H## in matrix form by using particle number basis i.e. 1100, 1010 ##\cdots## for 4 sites and 2 particles
2. then I diagonalize ##H## and find energy eigenvalues ##E_i## and eigenvectors ##\psi_i##
3. then I write ##\rho## by using above given formula
4. finally, I try to calculate ground state energy using ##\langle H \rangle = Tr[\rho H]##
Problems I am facing:
As Boltzmann constant ##k_B## is ##1.36\times 10^{-23}##, so even at room room temperature ##\beta## is very large i.e. ~##10^{21}##. Due to which I can't calculate ground state energy at very low temperatures? Computer just give 'Infinity' values. And at very large temperatures, the energy seems to be decreasing with increasing ##T##. ... totally wrong
Questions:
1. Is my algorithm correct? How can we calculate ground state energy at ##T=0## and at small temperatures?
2. How to calculate ground state wave function at ##T##?
$$H = \sum_{i=1}^L c_{i+1}^\dagger c_i + h.c.$$
In the canonical formalism, the density matrix is given as
$$\rho = \frac{e^{-\beta H}}{Tr [e^{-\beta H}]}$$
where ##\beta = 1/k_B T## and the expectation value of any operator ##A## is given as
$$\langle A \rangle = Tr[\rho A]$$
Things that I want to calculate:
1. ground state energy at temperature ##T##
2. ground state wave-function at ##T##
My attempt to solution:
1. I first write ##H## in matrix form by using particle number basis i.e. 1100, 1010 ##\cdots## for 4 sites and 2 particles
2. then I diagonalize ##H## and find energy eigenvalues ##E_i## and eigenvectors ##\psi_i##
3. then I write ##\rho## by using above given formula
4. finally, I try to calculate ground state energy using ##\langle H \rangle = Tr[\rho H]##
Problems I am facing:
As Boltzmann constant ##k_B## is ##1.36\times 10^{-23}##, so even at room room temperature ##\beta## is very large i.e. ~##10^{21}##. Due to which I can't calculate ground state energy at very low temperatures? Computer just give 'Infinity' values. And at very large temperatures, the energy seems to be decreasing with increasing ##T##. ... totally wrong
Questions:
1. Is my algorithm correct? How can we calculate ground state energy at ##T=0## and at small temperatures?
2. How to calculate ground state wave function at ##T##?