Top: Simulation scheme. The e-beam impinges normally on the accretion disk surface. Bottom: Normalized density profile (see Eq. (6)) represented for the inner border of a sample disk with an accretion rate of 10⁻⁸ M_⊙ yr⁻¹. The structure is symmetric around the disk midplane, and only the upper part is represented. The density profile from Eq. (6) is plotted with a blue line, and the limits of the five layers implemented in the simulation are marked with dashed lines.