Articles you may be interested inAlthough it was proposed some time ago that ͑hyper͒polarizabilities might be estimated from the results of x-ray charge density refinements, early results were unconvincing. In this work we show that the one particle density obtained from the usual multipole refinement model does not contain sufficient information to determine these response properties and instead pursue the "constrained wave function" approach of fitting to x-ray structure factors. Simplified sum-over-states expressions are derived for determining the dipole polarizability from these wave functions, and these clearly show that the earlier work ignored important two-electron expectation values for the dipole polarizability, and two-and three-electron terms for , etc. Correction factors for the simplified sum-over-states polarizability tensors from the constrained wave function are obtained by calibration against coupled Hartree-Fock ab initio results to yield in-crystal effective polarizability tensors. Results obtained for benzene, urea, and 2-methyl-4-nitroaniline demonstrate that the effective molecular polarizabilities clearly include the effects of intermolecular interactions and electron correlation, especially for urea where the effects on the polarizability are known to be quite large. We also carefully consider the way in which the linear bulk susceptibility, ͑1͒ , and refractive indices are determined from the x-ray fitted polarizabilities, employing three models based on a rigorous treatment of the local field. Incorrect results are obtained for the sort of molecules that are of interest in nonlinear optical applications if the molecules are approximated by single point dipoles. In contrast, the use of Lorentz-factor tensors averaged over several sites yields excellent results, with refractive indices obtained using this model in remarkably good agreement with optical measurements extrapolated to zero frequency.