PN inspiral quantities

[nilsvu]

Hi Mike! I'm considering to use your new PN implementation to get good initial orbital parameters for eccentricity reduction in SpECTRE. Currently I use the SpEC scripts for this. Specifically, see here for the function in SpECTRE that calls the SpEC scripts and returns initial orbital parameters. My question is: does the Julia code already expose equivalents to the omegaAndAdot and nOrbitsAndTotalTime functions in SpEC's ZeroEccParamsFromPN.py script, or would those be easy to add?

[nilsvu]

Edit: I got the number of orbits and time to merger from sxs.julia.PNWaveform. So the only remaining question is if the PN approximation for Omega(r, q, chi1, chi2) and adot(r, q, chi1, chi2) is available as well somehow, ideally in Python? Here, r is the coordinate separation of the quasicircular binary.

[moble]

The output from sxs.julia.PNWaveform also has a field containing the PN expansion parameter v as a function of time, which is defined as $v := (M \Omega)^{1/3}$, so I think Omega = w.v**3 / (w.M1 + w.M2) would be what you want.

As for adot, I don't know that I've actually ever seen a PN approximation directly for that. It's certainly not involved in the standard (quasicircular) PN computation. I'm on vacation at the moment, so I'm afraid I won't be able to help much for the next couple weeks.

[nilsvu]

Ok thank you! Have a good vacation! We can talk more when you're back. For then:

PNWaveform takes Omega_i as input and integrates the orbital_evolution. That's great and we can use it for the number of orbits and the time to merger. But I'm also looking for just the approximation of Omega(r, q, chi1, chi2) and adot(r, q, chi1, chi2) without having to do the time integration. That's at least how SpEC determines these initial orbital parameters given an initial separation of the binary. Edit: looks like SpEC uses Eq. (4.2) in https://arxiv.org/abs/1212.5520v1 for Omega and Eq. (4.12) in http://arxiv.org/abs/gr-qc/9506022 for adot, according to the comments in the code. Or Eq. (369-370) in https://arxiv.org/pdf/1310.1528v5 plus spin effects.