# Is it possible to work out the colour of *any* molecule?

Above are some spectral emission lines for Hydrogen, Helium and Neon.

Using the Schrodinger equation, it's possible to derive the colours that Hydrogen will emit when light is fired at it. A little more complicated (although beyond my rudimentary knowledge of quantum mechanics), the same can be done for Helium, and so on.

Of course as we progress to molecules, even the relatively small $C_5H_{11}Br$, it's (I'm told) very hard to predict the colours emitted from first principles, and even then incorrect assumptions must be made (like the Born-Oppenheimer approximation), so upon progressing to larger molecules, we find that these simplifications diverge from experiment.

My question is this: with the advent of computers, are we able to (in practice, not theoretically in a few decades) accurately compute the colours given off by a molecule by numerically solving the laws of quantum mechanics for complicated systems? Much like nonlinear differential equations can now be 'solved' with computers without any simplifying assumptions (like $\sin(x) \approx x$), is it in humanity's capacity today to approximate colours from first principles, so that computing time is the only restraint? Note the difference between an incorrect model of reality (Born-Oppenheimer), and one that converges to the correct one (even if slowly).

-
Why is the Born-Oppenheimer approximation incorrect? –  Greg Apr 28 at 11:26

See, for instance, this paper by Jacquemin et al., where they calculate $\lambda_{max}$ for a wide range of dyes using TDDFT and get a MAE of 22 nm using CAM-B3LYP (a long-range hybrid functional), reducing that to 10 nm with an empirical calibration curve. The paper only really cares about accurately reproducing the location of $\lambda_{max}$, however this paper by Quartarolo and Russo computes complete spectra with fairly reasonable accuracy. Not perfect, but the general shape of the spectra are reproduced adequately. Part of the divergence between the experimental and computed spectra probably comes down to higher order broadening effects not reflected in the simple Gaussian convolution of the computed data, but I'm not a spectroscopist.