DOI: 10.5281/zenodo.16881721
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: outputs_q_0p833_enfunc_2_2PN_spin_flag_ON.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$'  and '$k$', which are are give in columns  as cosphi12, cospsi1, cospsi2, costilt1, costilt2, et and 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 and 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.