General strategy for computing nonlinear optical properties of large neutral and cationic organic chromophores in solution.