Abstract. This contribution summarizes our current findings in the field of calculating and predicting the linear and second-order nonlinear electric susceptibility tensor components of organic crystals. The methodology used for this purpose is based on a combination of the electrostatic interaction scheme developed by Munn and his coworkers (RLFT) with high-level electronic structure calculations. We compare the results of calculations with available experimental data for several examples of molecular crystals. We show the quality of the final results is influenced by i) the chromophore geometry, ii) the method used for molecular properties calculations and iii) the partitioning scheme used. In conclusion we summarize further plans to improve the reliability and predictability of the method.