The interaction of a molecule with optical fields is customarily interpreted by means of induced time-dependent electric polarizabilities, magnetizabilities and mixed electric-magnetic polarizabilities. In general, these properties can be rationalized by integrals of density functions formulated in terms of induced charge and current densities. In this perspective, we focus on what has been done so far at the theoretical level, and on what can be expected to be unveiled from the topological study of suitable density functions, endowed with the fundamental requirement of origin invariance. Densities characterized by such a property can be integrated all over the configuration space to obtain electric dipole polarizability and optical rotatory power. Corresponding maps visualize domains mainly involved in the molecular response. The diagonal components of origin-independent density tensor functions that, on integration, yield corresponding electric dipole polarizability tensor of benzene, naphthalene, phenanthrene and ovalene, have been computed, confirming the ubiquitous presence of counter-polarization regions in the proximity of the atomic nuclei. They are associated with toroidal electron currents, induced by time derivative of the electric field of impinging radiation. Electron (de)localization in these systems is readily observed and estimated. The optical rotation density of the carbonyl chromophore is studied in detail. Its essential feature is the separation in quadrants of alternating sign of density about the CO bond. The presence of an extrachromophoric perturbation determines asymmetry in the extension of the quadrant distribution, thus causing optical rotation.