Framework for ab initio calculation of nonlinear optical susceptibilities on the level of Real Time Time Dependent Density Functional Theory was developed and implemented in high performance numerical codes. Methodology ...