We have assessed six range-separated hybrids (CAM-B3LYP, LC-PBE, ωB97, ωB97X, ωB97X-D and M11) in evaluating the dynamic quadratic polarizability along the direction of the molecular dipole moment, β μ (−2ν; ν; ν; 0), by comparing with MP2(Full)-calculated data. There are 104 solution-phase organic molecules included in the set. Both geometry optimizations and β μ (−2ν; ν; ν; 0) computations have been carried out by the same functional. Calculated results showed that simultaneous modeling geometrical and electronic effects improved significantly the nonlinear optical (NLO) property of interest. More generally, our numerical assessment revealed that the CAM-B3LYP functional could be the reasonable choice for NLO property computations.