Skip to content

Feature/admquantities#275

Merged
KAClough merged 5 commits intomainfrom
feature/admquantities
Sep 15, 2025
Merged

Feature/admquantities#275
KAClough merged 5 commits intomainfrom
feature/admquantities

Conversation

@llibert94
Copy link
Member

@llibert94 llibert94 commented Aug 26, 2025

This is a pull request to recover Tiago's one in #121. It adds the feature for computing the ADM mass and spin of a spacetime in the KerrBH example. It uses eq. (7.15) in [arxiv:gr-qc/0703035] for the ADM mass, which differs slightly from eq. (3.128) in Baumgarte & Shapiro's book used in #121, and eq. (3.191) in Baumgarte & Shapiro's book for the ADM momentum.
When extracted in the KerrBH using the CCZ4 formalism and over a box of L=256M for a BH of mass=2M and spin a=0.8M
image
image

By doing the extrapolation to infinity using FitToInfinity.py in https://github.com/GRTLCollaboration/Postprocessing_tools/tree/feat/grchombo_and_gw_tools/PythonTools/DataAnalysis/GWTools, one gets the correct results, namely J=1.603717148220815458M^2 and mass=1.992012925023880276M

@llibert94 llibert94 requested a review from KAClough August 26, 2025 06:39
@KAClough
Copy link
Member

This looks great! Can you add some evidence of a test to show that you get the right values. Ideally it would be a proper test of convergence or something, but I will accept e.g. a plot of the output over time for M = 2 a/M = 0.8.

@KAClough KAClough added the enhancement Modification of existing feature/general improvement label Aug 29, 2025
0 1 2 3 #Theta and Gamma
0 1 2 3 1 2 3 #lapse shift and B
vars_parity_diagnostic = 0 1 2 3 #Ham and Mom
0 0 #Madm and Jadm
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Can you align these with the values above?

0 1 2 3 #Theta and Gamma
0 1 2 3 1 2 3 #lapse shift and B
vars_parity_diagnostic = 0 1 2 3 #Ham and Mom
0 0 #Madm and Jadm
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Align with values above.

@llibert94 llibert94 requested a review from KAClough September 15, 2025 09:31
Copy link
Member

@KAClough KAClough left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

This is great, thanks @llibert94 for adding it to the public code! There are some features that could be added, like making the spin direction general, but this is good for now and we can upgrade it later if needed.

@KAClough KAClough merged commit 9c299bd into main Sep 15, 2025
47 checks passed
@KAClough KAClough deleted the feature/admquantities branch September 15, 2025 10:53
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

enhancement Modification of existing feature/general improvement

Projects

None yet

Development

Successfully merging this pull request may close these issues.

2 participants