Partition function from the density of states

Click For Summary

Homework Help Overview

The discussion revolves around calculating the canonical partition function from a given density of states, which includes delta functions and a step function. The problem involves the use of the incomplete gamma function and the Laplace transform approach to evaluate an integral that arises from the third term of the density of states.

Discussion Character

  • Exploratory, Mathematical reasoning, Problem interpretation

Approaches and Questions Raised

  • Participants discuss the evaluation of an integral involving a step function and exponential decay. There is uncertainty about the application of the incomplete gamma function and the manipulation of the integral limits. Some participants suggest re-evaluating the integral by changing variables and considering the properties of the gamma function.

Discussion Status

Several participants are engaged in exploring the integral's evaluation, with some providing guidance on potential methods. There is an ongoing examination of the steps involved, and while some participants express doubt about their calculations, others offer alternative approaches to simplify the problem.

Contextual Notes

The discussion is framed within the constraints of a homework assignment, emphasizing the need for clarity in mathematical reasoning and the exploration of different methods without arriving at a definitive solution.

snatchingthepi
Messages
145
Reaction score
37
Homework Statement
See post
Relevant Equations
See post
I'm given the following density of states

$$ \Omega(E) = \delta(E) + N\delta(E-\Delta) + \theta(E-\Delta)\left(\frac{1}{\Delta}\right)\left(\frac{E}{N\Delta}\right)^N $$

where $ \Delta $ is a positive constant. From here I have to "calculate the canonical partition function as a function of $$ x=\beta\Delta $$ using the incomplex gamma function

$$ \Gamma(n,x) = \int_x^\infty dt e^{-t} t^{n-1} $$

I know this can be solved for a partition function by taking a Laplace transform of the density of states. I can do the first two term very easily, but for the third term

$$ z_{can} = \int_0^\infty \theta(E-\Delta)\left(\frac{1}{\Delta}\right)\left(\frac{E}{N\Delta}\right)^N exp[-\beta E] dE $$

I'm not sure how to go forward from here. I've never seen an integral like this. I am thinking the step function changes the integral lower bound, but I'm kinda strung out so near the end of term, and am not seeing where to go now. Can anyone please help out?
 
Physics news on Phys.org
Please tell us exact step where you are stuck. The integral is doable using incomplete gamma functions
 
I am unsure *how* to do this integrla with the incomplete gamma function. My thought hit a dead-end at

$$ z_{can} = \int_0^\infty \theta(E-\Delta)\left(\frac{1}{\Delta}\right)\left(\frac{E}{N\Delta}\right)^N exp[-\beta E] dE $$

$$ z_{can} = \int_\Delta^\infty \left(\frac{1}{\Delta}\right)\left(\frac{E}{N\Delta}\right)^N exp[-\beta E] dE $$

let $$ x=\beta \Delta $$

and for the incomplete gamma function let

$$ t = \frac{-xE}{\Delta}, dt = \frac{-x}{\Delta} dE $$

so

$$ z_{can} = \left(\frac{1}{N \Delta}\right)^N \int_x^\infty dt t^N exp[-t] $$

$$ z_{can} = \left(\frac{1}{N \Delta}\right)^N \Gamma(N+1, x) $$

I'm not convinced of my math in these last few steps.
 
Last edited:
snatchingthepi said:
Homework Statement:: See post
Relevant Equations:: See post

I'm given the following density of states

$$ \Omega(E) = \delta(E) + N\delta(E-\Delta) + \theta(E-\Delta)\left(\frac{1}{\Delta}\right)\left(\frac{E}{N\Delta}\right)^N $$

where $ \Delta $ is a positive constant. From here I have to "calculate the canonical partition function as a function of $$ x=\beta\Delta $$ using the incomplex gamma function

$$ \Gamma(n,x) = \int_x^\infty dt e^{-t} t^{n-1} $$

I know this can be solved for a partition function by taking a Laplace transform of the density of states. I can do the first two term very easily, but for the third term

$$ z_{can} = \int_0^\infty \theta(E-\Delta)\left(\frac{1}{\Delta}\right)\left(\frac{E}{N\Delta}\right)^N exp[-\beta E] dE $$

I'm not sure how to go forward from here. I've never seen an integral like this. I am thinking the step function changes the integral lower bound, but I'm kinda strung out so near the end of term, and am not seeing where to go now. Can anyone please help out?
You could also replace the factor ##E^N## by ##(-1)^N## times the N-th derivative of the exponential with respect to ##\beta##. The integral will be trivial and then you can apply the N-th derivative wrt ##\beta## on the result to get the final answer.
 
  • Like
Likes   Reactions: snatchingthepi

Similar threads

  • · Replies 1 ·
Replies
1
Views
2K
Replies
0
Views
2K
  • · Replies 5 ·
Replies
5
Views
2K
  • · Replies 0 ·
Replies
0
Views
2K
Replies
2
Views
2K
Replies
1
Views
2K
  • · Replies 4 ·
Replies
4
Views
2K
  • · Replies 9 ·
Replies
9
Views
6K
Replies
19
Views
3K
Replies
2
Views
2K