From top to bottom: evolution of inverse semimajor axis and eccentricity in the pure N-body regime for physical models A, B, C, and D (left) and numerical models B1-4 (right) based on physical model B. It is worth noting that for numerical models B3 and B4 we turned on PN terms from the start of simulations. For numerical models, residuals between model B and models B1-4 are shown as colored dashed lines. Colored arrows mark the time when each system becomes bound t_b (pale color) and then hard t_h (intense color). Colors were chosen according to the main color scheme for each model.