Quantitative information of solvation and transfer free energies is often needed for the understanding of many physicochemical processes, e.g the molecular recognition phenomena, the transport and diffusion processes through biological membranes and the tertiary structure of proteins. Recently, the molecular free energy surface density concept (MolFESD) has been introduced [1,2]. This model is based on the assumptions that the overall hydrophobicity can be obtained as a superposition of fragment contributions and that the corresponding free energy can be calculated as an integral of the MolFESD over the molecular solvent accessible surface. The scalar quantity 3D-FED, the three-dimensional free energy density offers a physical basis to the establishment of a new predictive model with limited empirical character. Although this volume density is accessible from more accurate methods e.g. Monte-Carlo methods, the Grid program  has proved to be suitable for a rapid evaluation of the 3D-FED. It uses empirical force field to determine the interaction energy of a particular probe molecule (e.g water) for all points of a regular three dimensional grid in which the target molecule is enclosed. In the present study, parametrization of the MolFESD is done in two steps using the same strategy described in the previous work of Jäger et al. [1,2]; some model parameters were fitted to computed data (provided by Grid program) while other parameters were fitted to experimental logP values. Unlike most fragment based predictive models for logP, the MolFESD concept calculates logP values which depend on the area and shape of molecular surfaces and in turn, on the corresponding molecular conformations. A parametrization of the MolFESD is done for three solvents (n-octanol/water, alkane/water and cyclohexane/water) followed by its application to the prediction of brain-blood barrier partition coefficients (logBB).