Eigene Flugzeit Schwarzschild-Metrik ist endlich

Beim Studium der metrischen Schwarzschild-Geodäte kann man leicht auf die folgende Differentialgleichung kommen

D R D τ = C 2 ( 1 2 G M R )
die die radiale Koordinate und die Eigenzeit außerhalb des Ereignishorizonts in Beziehung setzt R H = 2 G M (Ich benutze natürlich C = 1 ).
Stellen Sie sich vor, ein Beobachter fällt aus der Ausgangsposition R 0 = 2 G M + H (mit H > 0 um außerhalb des Ereignishorizonts zu sein) bis hinunter zum Ereignishorizont R H = 2 G M . Die Eigenflugzeit ist dank der oben gegebenen Differentialgleichung durch
Δ τ = τ 0 τ H D τ = R 0 = 2 G M + H R H = 2 G M D R C 2 ( 1 2 G M R )
Die Koordinatenflugzeit ist jedoch durch die Relation gegeben D T = C ( 1 2 G M / R ) D τ (Wo C > 0 , siehe beispielsweise http://gfm.cii.fc.ul.pt/events/lecture_series/general_relativity/gfm-general_relativity-lecture4.pdf ), resultierend
Δ T = T 0 T H D T = R 0 = 2 G M + H R H = 2 G M C 1 2 G M R D R C 2 ( 1 2 G M R )
Daraus (siehe zum Beispiel die oben angegebene Referenz) heißt es: „ Es ist leicht durch direkte Auswertung der Integrale zu erkennen, dass die Eigenflugzeit endlich ist, während die Koordinatenflugzeit unendlich ist, die Umlaufbahn Beobachter wird niemals sehen, dass der einfallende Beobachter den Ereignishorizont erreicht, außer asymptotisch ''.

Ich habe versucht, das erste Integral mit Mathematica zu berechnen, und das Ergebnis, das ich erhalte, ist absurd und gibt eine komplexe Eigenzeit an. Hier ist der Code, den ich verwendet habe:
Integrate[1/Sqrt[C^2 - (1 - 2*G*M/r)], {r, 2*G*M, 2*G*M + h}]
Und die Ausgabe ist:
ConditionalExpression[(G M (2 Sqrt[C^2 - C^4] - I Log[2] - I Log[((-I + 2 I C^2 + 2 Sqrt[C^2 - C^4]) G M)/Sqrt[1 - C^2]]))/(1 - C^2)^(3/2) + (-Sqrt[1 - C^2] (h + 2 G M) Sqrt[C^2 - h/(h + 2 G M)] + I G M Log[2] + I G M Log[(-I h + I C^2 h - I G M + 2 I C^2 G M + Sqrt[1 - C^2] h Sqrt[C^2 - h/(h + 2 G M)] + 2 Sqrt[1 - C^2] G M Sqrt[C^2 - h/(h + 2 G M)])/Sqrt[1 - C^2]])/(1 - C^2)^(3/2), ((G M)/h != 0 && Re[(GM)/h] >= 0) ||Re[(G M)/h] < -(1/2) || (G M)/h \[NotElement] Reals]

Was fehlt mir hier oben? Danke

Der beste Weg, sich Integralen zu nähern, besteht darin, einige nette Substitutionen durchzuführen, um die Integranden so weit wie möglich zu vereinfachen (wie in der zweiten Antwort vorgeschlagen), und so viele Variablen zu belassen, neigt dazu, das CAS zu vermasseln, da es Annahmen über diese Größen treffen muss. Es weiß nicht, dass M eine Masse ist und zum Beispiel positiv sein muss, wie Sie im Bedingungsausdruck sehen können. Eine Möglichkeit, dies zu unterstützen, besteht darin, Annahmen mithilfe von Refine[Ausdruck, Annahmen->{...}] festzulegen, wobei die geschweiften Klammern alle Annahmen enthalten, die zum Bereinigen des Ausdrucks erforderlich sind.
Das Analyseergebnis ist hier angegeben: physical.stackexchange.com/questions/426143/427025#427025 - Das zeigt auch die zweite Grafik T (in blau) divergiert dabei am Horizont τ (in grün) nicht.
@Triatticus Danke! Ich habe es mit Ihrer Antwort versucht, aber das angegebene Ergebnis stimmt nicht mit dem erwarteten Wert überein, obwohl es tatsächlich real und sinnvoll ist.
@safesphere Danke! Die Ergebnisse, die Sie angeben, sind korrekt, das sind die Ergebnisse, die ich mit Mathematica überprüfen wollte. Ich weiß jedoch immer noch nicht, warum der Befehl nicht funktioniert.
Gibt es hier eine Physikfrage?

Antworten (2)

Ich weiß auch nicht so recht, wie das "einfach zu sehen" ist.
Sie können das Integral jedoch nach der richtigen Zeit lösen, indem Sie zuerst das für notieren C 2 = konst. wir haben (mit R S = 2 G M )

C 2 = 1 R S R 0

was aus dem Einstecken folgt ( 4.2 ) hinein ( 4.1 ) und Betrachten eines einfallenden Teilchens, das bei ruht R = R 0 , Bedeutung R ˙ | R = R 0 = 0 .
Dies können Sie verwenden, um Ihre erste Gleichung zu schreiben als

R ˙ = R S 1 / 2 ( R 0 R R 0 R ) 1 / 2

Daraus ergibt sich das Eigenzeitintegral

Δ τ = R S 1 / 2 R 0 R S ( R R 0 R 0 R ) 1 / 2 D R

die Sie mit Mathematica oder auch durch Einführung der Parametrisierung lösen können

R ( η ) = R 0 2 ( 1 + cos η )

mit η [ 0 , π ] .

Danke! Ja, wissen Sie, wissenschaftliche Arbeiten enthalten viele Wörter wie „trivial“, „offensichtlich“, „sehr klar“, … in Momenten, in denen es nicht einmal annähernd trivial oder einfach ist, haha. Danke für Ihre Antwort, aber setzen Sie den Befehl in Mathematica mit Ihrem Ausdruck für ein Δ τ das Programm gibt kein Ergebnis, ich weiß nicht, was passiert.
Hast du es auch schon mit der Parametrierung probiert? Damit ist es eigentlich ganz einfach und kann mit Stift und Papier erledigt werden. Ich habe selbst keinen Zugriff auf Mathematica und kann Ihnen daher leider nicht wirklich mit dem Softwareteil helfen.
Ja, ich habe es auch mit der Parametrisierung gemacht und es wird kein Ergebnis berechnet. Mach dir keine Sorge! Danke für deine Antwort. Ich weiß, dass es analytisch gemacht werden kann, aber ich wollte wissen, was mit der Software los ist.

Das Ergebnis, das ich erhalte, ist absurd und gibt eine komplexe Eigenzeit an.

Ihr Ergebnis enthält Ies, von dem ich annehme, dass es sich um die Notation von Mathematica handelt ich = 1 , aber das bedeutet nicht, dass das Ergebnis komplex ist. Bei bestimmten Integrationstechniken ist es üblich, Ergebnisse zu erhalten, die komplex aussehen, sich aber tatsächlich als real erweisen, wenn man sie auswertet.

Sie machen Ihre Ergebnisse unnötig kompliziert, indem Sie nehmen G Und M als Parameter. Bei dieser Art von Integral gilt es, so viele Variablen wie möglich in einheitenlose Größen umzuwandeln, was im Allgemeinen dazu führt, dass unnötige Konstanten wie diese eliminiert werden. Hier wollen Sie auf die Variable wechseln X = R / 2 G M . Geben Sie es nicht einfach in das Computer Algebra System (CAS) ein, ohne diese Vorbereitungen zu treffen. CASs sind dumm und erzeugen komplizierte Ausgaben, geben Sie ihnen also eine Chance, indem Sie dasselbe Setup durchführen, das Sie normalerweise tun würden, wenn Sie ohne CAS arbeiten würden.

Ich verwende eher die Open-Source-CAS-Maxima als Mathematica. Es hat kein Problem, ein Ergebnis ohne explizite Ergebnisse zu erzeugen ich ist drin.

$ maxima -q --batch-string="assume(h>0 and (c^2-1)<0 and c>0); integrate(-(c^2-(1-1/x))^(-1/2),x,1+h,1);"

(%i1) assume(h > 0 and c^2-1 < 0 and c > 0)
                                     2
(%o1)                       [h > 0, c  < 1, c > 0]
(%i2) integrate(-(c^2-(1-1/r))^((-1)/2),r,1+h,1)
                                   2                     2           2
                 2       sqrt(1 - c ) sqrt(h + 1) sqrt((c  - 1) h + c )
(%o2) (sqrt(1 - c ) atan(----------------------------------------------)
                                        2           2
                                      (c  - 1) h + c  - 1
     2                         2           2     4      2
 + (c  - 1) sqrt(h + 1) sqrt((c  - 1) h + c ))/(c  - 2 c  + 1)
                                 2
             2       c sqrt(1 - c )     3
   sqrt(1 - c ) atan(--------------) + c  - c
                          2
                         c  - 1
 - ------------------------------------------
                  4      2
                 c  - 2 c  + 1

Vielleicht möchten Sie versuchen, es in Mathematik auf die gleiche Weise einzurichten. Nach der Änderung der Variablen wird es sicherlich einfacher und leichter zu interpretieren sein.

Eigentlich gibt es keine Garantie dafür, dass dieses Ergebnis wirklich herauskommt, nur weil es keine gibt ich ist drin. Es enthält eine Quadratwurzel, die aussieht, als könnte sie imaginär sein. Wir wissen, dass es reell sein sollte, weil es ein bestimmtes Integral eines reellen Integranden ist. Wenn Sie dann überprüfen möchten, ob es wirklich echt ist, versuchen Sie es einfach mit einigen Zufallszahlen C Und H . Wenn sich herausstellt, dass das Ergebnis echt ist, dann ist das kein Zufall. Wenn es komplex wird, kann eine der folgenden Aussagen zutreffen: (1) Sie haben es falsch codiert, (2) es gibt einen Fehler im CAS oder (3) es gibt ein Problem mit Zweigschnitten für die umgekehrte Tangente.

Ich habe versucht, Ihre Ersetzung in dimensionslose Einheiten zu verwenden, aber das Ergebnis ist dasselbe wie oben angegeben und hat immer noch die schreckliche Explizitheit ICH = ich = 1 . Ich verstehe nicht, was passiert, vielleicht versuche ich, Maxima oder andere Software zu verwenden und es dort zu überprüfen. Danke!