Dear all, In order to get density profile of a pure system the box has to get extended in one drection (e.g Z) as we do for surface tension calculations or density profile can be also obtained from the usual box filled up with the molecules without need for empty space in Z Thanks for any comments and advise in advance :) Best E.