Here we derive a completely quantum mechanical expression for linear susceptibility using a density matrix formalism where the electromagnetic field is also quantized. A zeroth order susceptibility emerges naturally from the treatment. By studying the dependence of susceptibility on material asymmetry and on the quantum state of light, we predict new and intriguing phenomena such as photon number dependent refraction index.
