DOI: 10.5281/zenodo.16881720
title: Data release: Evolution of precessing binary black holes  on eccentric orbits using orbit-averaged evolution equations

description: Data behind the figures appearing in Khun Sang Phukon, Nathan K. Johnson-McDaniel, Amitesh Singh, and Anuradha Gupta, "Evolution of precessing binary black holes on eccentric orbits using orbit-averaged evolution equations", Phys. Rev. D 112, 044067, 2025,  arXiv:2504.20543.

Figure 2:



Data file:  fig2_data.zip 

Description: This zip file contains files with the timeseries of the parameters omega in  M = 1 unit, l, lamb, S1x, S1y, S1z, S2x, S2y, S2z, costilt1, costilt2, cosphi12, cospsi1, cospsi2, LNhatx, LNhaty, LNhatz, Phatx, Phaty, Phatz, et, k for 16 binaries.  Parameters of each binary are encoded in the filenames, for instance, the q_0.833_chi1_0.95_chi2_0.95_tilt1_0.99_tilt2_0.62_phi1_3.61_phi2_1.31_M_20.csv gives the binary's mass ratio q=0.833, dimensional spin magnitudes (chi1=0.95, chi2=0.95), initial spin angles (tilt1=0.99, tilt2=0.62, phi1=3.61, phi2=1.31) and total mass M=20. 


Figure 3: 



Data file: evolved_outputs_and_diff_bw_pairwise_enfunc.csv

Description of file: The columns with headers q, chi1, chi2 give binaries' mass ratio, dimensionless magnitude of spin1 and dimensionless magnitude of spin2, while columns with headers parameter_vN, with N $\in \{0, 1, 2\}$ represent final values of these parameters obtained with the Arun et al. and Loutrel and Yunes superasymptotic and hyperasymptotic enhancement functions for 0, 1, 2, respectively. The columns with header parameter_NvM_abs give absolute differences between final values of corresponding parameters obtained with enhancement functions N $\in \{0, 1,2 \}$ and M $\in \{0, 1, 2\}$.  The parameters are $\cos \phi_{12}^{}$, $\cos \psi^{}_1$, $\cos \psi^{}_2$, $\cos \theta_1^{}$, $\cos \theta_2^{}$, $e_t$, $k$, which are are give in columns  as cosphi12, cospsi1, cospsi2, costilt1, costilt2, et, k, respectively. 


Figure 4:



Data file: 'eccentricity_evolution_monotonic_increase.txt' for monotonic increase of eccentricity,  'eccentricity_evolution_monotonic_decrease.txt'  for monotonic decrease of eccentricity and 'eccentricity_evolution_nonmonotonic.txt' for non-monotonic evolution of eccentricity. 

Description: Each file contains semimajor axis values in units of M, median eccentricity and eccentricity of each binary at the corresponding radius given by columns  radius, eccentricity_median, ecc_N, where the integer N is the index of each binary.


Figure 5:



Data files:  fraction_of_binaries_a_1000M.txt,  fraction_of_binaries_a_50M.txt,  fraction_of_binaries_a_20M.txt,  fraction_of_binaries_a_10M.txt.

Description of files:  Files contain values of instantaneous eccentricities at the upper-bin edges of binaries at semimajor axes 1000M, 50M, 20M, 10M and the corresponding fraction of binaries in C morphology, L0 morphology, and Lpi morphology, given by the f_c_morph, f_lo_morph, f_lpi_morph columns, respectively.


Figure 6:



Data files: c_to_l0_transition.txt, c_to_lpi_transition.txt

Description of files: Files contain initial eccentricities of binaries (given as e_0) and corresponding orbital separation or semimajor axis at which binaries transition to either L0 or Lpi morphology from C morphology.