pyseobnr.eob.dynamics.integrate_ode_prec.compute_dynamics_quasiprecessing

pyseobnr.eob.dynamics.integrate_ode_prec.compute_dynamics_quasiprecessing(omega_ref: float, omega_start: float, H: Hamiltonian, RR: Callable, m_1: float, m_2: float, chi_1: ndarray, chi_2: ndarray, params: EOBParams, rtol: float = 1e-12, atol: float = 1e-12, step_back: float = 250, y_init=None, initial_conditions: Literal['adiabatic', 'postadiabatic'] = 'adiabatic', initial_conditions_postadiabatic_type: Literal['analytic', 'numeric'] = 'analytic')[source]

Compute the dynamics starting from omega_start, with spins defined at omega_ref.

First, PN evolution equations are integrated (including backwards in time) to get spin and orbital angular momentum. From that we construct splines either in time or orbital frequency for the PN quantities. Given the splines we now integrate aligned-spin EOB dynamics where at every step the projections of the spins onto orbital angular momentum is computed via the splines.

Parameters:

  • omega_ref (float) – Reference frequency

  • omega_start (float) – Starting frequency

  • H (Hamiltonian) – Hamiltonian to use

  • RR (Callable) – RR force to use

  • m_1 (float) – Mass of primary

  • m_2 (float) – Mass of secondary

  • chi_1 (np.ndarray) – Dimensionless spin of the primary

  • chi_2 (np.ndarray) – Dimensionless spin of the secondary

  • params (EOBParams) – Container of additional inputs

  • rtol (float, optional) – Relative tolerance for EOB integration. Defaults to 1e-12.

  • atol (float, optional) – Absolute tolerance for EOB integration. Defaults to 1e-12.

  • step_back (float, optional) – Amount of time to step back for fine interpolation. Defaults to 250.

  • y_init (np.ndarray, optional) – Initial condition vector (r,phi,pr,pphi).

  • initial_conditions (str, optional) – Type of initial conditions for the ODE evolution (‘adiabatic’ or ‘postadiabatic’).

  • initial_conditions_postadiabatic_type (str, optional) – Type of postadiabatic initial conditions for the ODE evolution (‘analytic’ or ‘numeric’).

Returns:

Aligned-spin EOB dynamics, PN time, PN dynamics, PN splines

Return type:

tuple