Surface density map in code units (1 code unit = 200 AU) in the x − y plane for the simulations with e = 0, q = 1 (upper panel); e = 0.45, q = 0.7 (middle panel); and e = 0.9, q = 1 (bottom panel) at three different moments: t = 1,  900,  1300 P_B for the equal mass cases and t = 1,  900,  1600 P_B for the unequal mass case. In each panel, the first and the second row show results for P_tot = P_gas and P_tot = P_gas + P_rad, respectively. In the eccentric simulations, the initial cavity radius is 3a₀, while in the circular binary case, it is 2a₀. When the radiation pressure is included, the binary begins accreting at the very early stages of its evolution; otherwise, the disc experiences a transient phase lasting approximately 900 P_B, 700 P_B, and 600 P_B respectively for the cases e = 0, q = 1; e = 0.45, q = 0.7; and e = 0.9, q = 1. The spatial resolution at the edge of the cavity (r = 3a) is Δx = 0.025a.