EoM in Schwarzschild geometry: geodesic v Hamilton formalism

Click For Summary
SUMMARY

This discussion centers on the discrepancies observed between two Mathematica scripts designed to simulate the trajectory of a timelike geodesic in Schwarzschild geometry using both the geodesic equations of motion and the Hamiltonian formalism. The first script produces a precessing orbit, while the second script results in a trajectory that falls directly into the event horizon. Key issues identified include potential errors in the Hamiltonian approach, particularly regarding the definition of momentum and the setup of initial conditions. The participants suggest verifying the equivalence of equations of motion derived from both methods to resolve the discrepancies.

PREREQUISITES
  • Understanding of Schwarzschild geometry and its implications for geodesics
  • Familiarity with Mathematica for coding and simulation
  • Knowledge of Hamiltonian mechanics and Lagrangian formalism
  • Ability to interpret differential equations and numerical solutions
NEXT STEPS
  • Review the derivation of Hamilton's equations from Lagrangian mechanics
  • Learn about the Hamilton-Jacobi method and its application in geodesic motion
  • Investigate the numerical methods for solving differential equations in Mathematica
  • Explore the properties of the Schwarzschild metric and its geodesic equations
USEFUL FOR

Researchers and students in theoretical physics, particularly those focusing on general relativity, computational physics, and numerical simulations of dynamical systems.

  • #31
vanhees71 said:
It becomes the same equation, given that for an affine parameter A=constA=constA=\text{const}.

Ah, the missing fact, thanks.
 
Physics news on Phys.org
  • #32
pervect said:
Possibly I'm confused, but something seems wrong here.

Let's take the flat space case when ##g_{\mu\nu} = \eta_{\mu\nu}## = diag(-1,1,1,1), so we have a -+++ metric signature.

We know that the relativistic Lagrangian for a free particle in the above flat space-time should be ##L = -m\,c^2 \sqrt{1-\beta^2}##, where ##\beta = v/c##. If we use geometric units so that c=1, we can write this as:

##L(t,\textbf{x},\frac{d \textbf{x}}{dt}) = -m \, \sqrt{ -g_{\mu\nu} \frac{dx^\mu}{dt} \frac{dx^\nu}{dt}}##

It's one of those lucky miracles that you get the same equations of motion using:

L = -m \sqrt{g_{\mu \nu} \frac{dx^\mu}{ds} \frac{dx^\nu}{ds}}

or

L =\frac{m}{2} g_{\mu \nu} \frac{dx^\mu}{ds} \frac{dx^\nu}{ds}

(Actually, the first gives more solutions than the second, but every solution to the first can be converted into a solution of the second by reparametrizing s)
 
  • Like
Likes   Reactions: vanhees71

Similar threads

  • · Replies 11 ·
Replies
11
Views
2K
  • · Replies 8 ·
Replies
8
Views
4K
Replies
3
Views
2K
  • · Replies 23 ·
Replies
23
Views
3K
Replies
9
Views
6K