Liouville equation with Dirac delta as probability density

Click For Summary
SUMMARY

The discussion centers on solving the Liouville equation, ∂ρ/∂t = -{ρ,H}, with the initial condition ρ(t=0) = δ(q,p), where δ(q,p) represents a Dirac delta function in phase space. The proposed solution is of the form ρ(t) = δ(q(t),p(t)), with q(t) and p(t) being trajectories derived from Hamilton's equations. Participants emphasize the need to verify this solution by substituting it back into the Liouville equation and addressing the complexities introduced by the Dirac delta function in the context of distributions.

PREREQUISITES
  • Understanding of the Liouville equation in statistical mechanics
  • Familiarity with Hamiltonian mechanics and Hamilton's equations
  • Knowledge of Dirac delta functions and their properties in distributions
  • Basic proficiency in calculus, particularly in handling partial derivatives
NEXT STEPS
  • Study the properties of Dirac delta functions in the context of distributions
  • Learn how to apply Hamilton's equations to derive trajectories in phase space
  • Explore the method of characteristics for solving partial differential equations
  • Investigate the implications of the Liouville theorem in statistical mechanics
USEFUL FOR

Mathematicians, physicists, and students studying statistical mechanics or dynamical systems, particularly those interested in the application of the Liouville equation and Dirac delta functions in phase space analysis.

andresB
Messages
625
Reaction score
374
I would like to know the solution to Liouville equation
∂ρ/∂t=-{ρ,H}

given the initial condition

ρ(t=0)=δ(q,p)

where δ(q,p) is a dirac delta centered in some point (q,p) in phase space.

I have the feeling, but I'm not sure, that the solution is of the form

ρ(t)=δ(q(t),p(t))

where q(t) and p(t) are the trajectories from Hamilton equations.

Any help?
 
Physics news on Phys.org
Did you try inserting your suspicion into the Liouville equation to see if it does solve it?

Note that your suspicion is actually
$$
\rho(q,p,t) = \delta(q-q(t),p-p(t)),
$$
where ##q(t)## and ##p(t)## are just functions of time (that satisfy the Hamilton equations of motion), when expressed properly.
 
Orodruin said:
Did you try inserting your suspicion into the Liouville equation to see if it does solve it?

I tried, but it is confusing to me due to the nature of Dirac delta. Liouville equation read
##∂δ/∂t= -(∂δ/∂q)(∂H/∂p) + (∂δ/∂p)(∂H/∂q)## (1)

Equality should be understood in the distribution sense. Using a test function ##F(q,p)##, I get for the left hand side of (1)

##∫F(q,p) ∂δ/∂t dqdp= ∂/∂t∫F(q,p)δ(q-q(t),p-p(t) dqdp= ∂F(q(t),p(t))/∂t##

But I'm unsure how to prove the equality to the right hand side of (1).
 
You do not need to use test functions. Just apply the standard rules of derivatives that hold also for distributions. Take it step by step and show your work.
 

Similar threads

  • · Replies 4 ·
Replies
4
Views
1K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 3 ·
Replies
3
Views
1K
  • · Replies 1 ·
Replies
1
Views
3K
  • · Replies 3 ·
Replies
3
Views
2K
  • · Replies 4 ·
Replies
4
Views
2K
  • · Replies 5 ·
Replies
5
Views
2K
  • · Replies 2 ·
Replies
2
Views
2K
  • · Replies 4 ·
Replies
4
Views
2K
  • · Replies 1 ·
Replies
1
Views
2K