One of my chemical heroes is William Perkin, who in 1856 famously (and accidentally) made the dye mauveine as an 18 year old whilst a student of August von Hofmann, the founder of the Royal College of Chemistry (at what is now Imperial College London). Perkin went on to found the British synthetic dyestuffs and perfumeries industries. The photo below shows Charles Rees, who was for many years the Hofmann professor of organic chemistry at the very same institute as Perkin and Hofmann himself, wearing his mauveine tie. A colleague, who is about to give a talk on mauveine, asked if I knew why it was, well so very mauve. It is a tad bright for today’s tastes!
The first thing to note about mauveine is that it is not a single compound; actual samples can contain up to 13 different forms! These all vary in the number of methyl groups present which range from none up to four, in various positions. These compounds all have absorption maxima λmax in the range 540-550nm, the colour of purple. The structure of one of these, known as mauveine A, is shown below.
You can see from this that something is missing. The so-called chromophore is a cation, and an anion needs to be provided to balance the charge. We will now attempt to predict the color of purple using purely the power of quantum mechanics (for many years, accurate prediction of colour was a holy grail amongst dye chemists for obvious reasons). The anion can be chloride, and the colour is often measured in methanol as solvent. So the first task is to calculate this ion-pair. This used to be easier said than done (and in the past, the anion was often simply neglected). But using the ωB97XD density functional procedure (to get the van der Waals interactions modelled correctly) and a 6-311++G(d,p) basis set, coupled with a smoothed-cavity continuum solvation procedure, and two molecules of water (standing in for methanol, which is a bit bigger) as explicit solvent molecules, we get the structure apparent when you click on the diagram above (DOI: 10042/to-7320). Application of time-dependent density function theory (TD-DFT) gives a measure of the UV-optical spectrum (below, loaded as a scaleable SVG image. If you are using a modern browser, it should display. If not, try the latest FireFox, Chrome, Safari etc).
This has several noteworthy aspects.
LUMO | HOMO |
---|---|
So this post ends with a bit of a mystery. The fanciest most modern computational theory gets the colour of mauveine wrong by ~100nm. Why?
In the mid to late 1990s as the Web developed, it was becoming more obvious…
I have written a few times about the so-called "anomeric effect", which relates to stereoelectronic…
The recent release of the DataCite Data Citation corpus, which has the stated aim of…
Following on from my template exploration of the Wilkinson hydrogenation catalyst, I now repeat this…
In the late 1980s, as I recollected here the equipment needed for real time molecular…
On 24th January 1984, the Macintosh computer was released, as all the media are informing…
View Comments
Is there experimental evidence for ion-pairing? Mauveine seems to retain its color despite the anion employed (chloride, acetate, sulfate, etc.) and it is difficult for me to believe that all anions would coordinate in the same manner. Because there is no chloride character in either the HOMO or LUMO, I think we can assume its effect will be electrostatic. The more localized nature of the LUMO means that it will be raised higher in energy relative to the HOMO (compared to a model w/o chloride), and thus the transition would be higher in energy. This effect could be worth a few hundred meV, so I removed the chloride / two methanols and compute, without reoptimization, the first excited state at wB7XD/6-31G(d,p). The transition energy of 530nm seems quite close to experiment but one would have to test this with a better basis set and checking for significant structural reorganization in the absence of the chloride / methanols.
My apologies, one should be more careful blogging before morning coffee!
(i) I did not mean to undermine your work - I missed that you had tested this hypothesis yourself. That being said, with both the 6-31G(d,p) and 6-311G(d,p) basis sets I find a large shift from 450nm to ~740nm (see point ii for different # than quoted earlier) for the cation w/o anion and explicit methanols. Could you elaborate on what calculations led you to conclude that the explicit Cl / MeOHs didn't shift lambda max much?
(ii) In testing new input creation software for a colleague, I seem to have stumbled across a bug that causes it to ignore user-defined charge/multiplicity. This may be quite fortuitous. The data I listed in my previous post is for the neutral radical, and the value of 530-540nm seems to be independent of basis now that a few other calculations have finished. Is enough known about the electrochemistry of Mauveine to say whether the neutral radical is plausible?
Richard,
I had indeed earlier done a calculation on the pure cation rather than the ion-pair, and this came up with an (approximately) HOMO-LUMO excitation energy corresponding to ~437 nm (see DOI: 10042/to-7340) which is very similar to the ~440nm obtained for the full system. From which I had concluded that the electrostatic effect was surprisingly small, and that it hardly changed the HOMO-LUMO gap. This for the geometry of the cation fully optimized without the chloride. Your estimate based on removing the chloride and two solvent molecules, but leaving the rest of the geometry unchanged does not, as you say, take into account structural reorganisation, although in fact it does appear to be fairly small (
As a follow up to my earlier response, a quick check of CAS suggests nothing under the search term mauveine cyclic voltametry and only one irrelevant hit under mauveine reduction potential. The way forward might be to identify a system which is well characterised in this regard, and see if the methods I used in this post at least reproduce that system well.
I don't know if it will work but a relatively simple pH-dependent colour changing system like Congo Red might be a nice test.
I think virtually any intense nitrogen-containing chromophore might do the trick (Congo red is not a quaternary nitrogen cation, so it is a little different), but it would be worth selecting one that has as much structural similarly to mauevine as possible, and one where the 1-electron reduction potential is also known (with or without solvation).
Another update to my earlier notes. I have repeated the TD-DFT calculation with another functional (CAM-B3LYP), getting the same results (λmax ~440nm). Not too sensitive to counter-ion, solvent, basis set or functional then.
About the only thing I can think of is that the system aggregates/stacks in some way.
Here is a colour table which makes it easier to relate λmax to the actual observed colour.
I am still trying to find a model which will explain why mauveine is mauve (purple). In other words, why λmax is 540nm. Perhaps if one deprotonates mauveine to form what is called the free base?
The calculated spectrum is shown below, for which λmax is 481nm. Still not right!
Hello, bit of an odd question: do you happen to know an approximate value of a piece of fabric which would have been dyed with the original mauveine (or the value of Prof Charles Rees bow tie). I would love to know if such an item would be valued as a collectible or bought at an antiques auction.
Thank you very much for your time
The short answer is no. I do not believe any such item has ever come up for sale!
Posted on behalf of Stefan Grimme
======================
Dear Henry,
A colleague pointed me to your page http://www.ch.imperial.ac.uk/rzepa/blog/?p=3507 asking about this 100 nm failure and possible reasons (because I have some experience in theoretical UV/CD). Here is my opinion on that:
1. comparisons theory-exp. should be done on a linear eV or cm-1 scale but not in nm (which is non-linear). The difference from 450 nm to 550 nm corresponds to an error of only 0.45 eV.
2. Comparisons of vertical excitation energies with exp. band maxima can easily introduce errors of 0.2 eV.
3. Your solvation model is crude and may add another uncertainty of about 0.1 eV.
4. If you add all this you arrive at a more realistic 0.2-0.3 eV TDDFT error which we also observe for smaller chromophore benchmarks compared to e.g. high-level WFT results.
5. if you convert 0.2 eV into 5 kcal/mol you will find that this is not an uncommon DFT error also for e.g. a ground state reaction energy of about 50 kcal/mol. So in summary its not special or surprising at all (its our eye which is very sensitive).
Cheers!
Stefan
--
____________________________________________________________
Prof. Dr. Stefan Grimme (grimme@thch.uni-bonn.de)
Chair of Theoretical Chemistry
Mulliken Center for Theoretical Chemistry
University of Bonn
Beringstr. 4, D-53115 Bonn, Germany