Mass (top) and density (bottom) profiles from the three simulations: Q1D (left), ANI (middle) and SYM (right) normalized w.r.t. to the mass M_ta and radius r_ta of the turnaround region using the best-fit M₀(ϵ) at different snapshots, color coded from blue to red in increasing order of expansion factor. The dashed black curve denotes the self-similar predictions for the mass and density profiles corresponding to ϵ = 0.8, 0.7, 0.5 for Q1D, ANI, and SYM, respectively.