The colour of purple

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!

Charles Rees, wearing a bow tie dyed with (Perkin original) mauveine and holding a journal named after Perkin.

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.

Mauveine A. Click to load 3D

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.

  1. The visible (right hand side) part of the spectrum is very monochromatic, with λmax ~440nm. In other words, mauveine has a pure and intense colour.
  2. This λmax is hardly affected by the presence of the counterion.
  3. The electronic transition responsible for this band is a simple HOMO (highest-occupied-molecular-orbital) to LUMO (lowest-unoccupied-molecular-orbital) excitation of an electron.
  4. These orbitals are shown below.
    LUMO HOMO

    Mauveine A. LUMO. Click for 3D

    Mauveine A. HOMO. Click for 3D

  5. Note how the excitation involves the central region of the molecule, and one of the pendant aryl groups, but not the other. One might presume that tuning the colour would only work if changes are made to the first of these aryl groups.
  6. There is a real mystery about the calculated value of λmax, which differs from the observed value by about 100nm (the wrong colour, making mauveine orange rather than purple). Normally, this sort of time dependent density functional theory has errors no greater than 15-20nm. The calculated value of λmax is not sensitive to the basis set, or the presence or not of the counter ion and solvent. Clearly, a discrepancy of this magnitude must have some other explanation. Watch this space!

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?

Henry Rzepa

Henry Rzepa is Emeritus Professor of Computational Chemistry at Imperial College London.

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

  • 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

Recent Posts

Internet Archeology: reviving a 2001 article published in the Internet Journal of Chemistry.

In the mid to late 1990s as the Web developed, it was becoming more obvious…

1 month ago

Detecting anomeric effects in tetrahedral carbon bearing four oxygen substituents.

I have written a few times about the so-called "anomeric effect", which relates to stereoelectronic…

1 month ago

Data Citation – a snapshot of the chemical landscape.

The recent release of the DataCite Data Citation corpus, which has the stated aim of…

2 months ago

Mechanistic templates computed for the Grubbs alkene-metathesis reaction.

Following on from my template exploration of the Wilkinson hydrogenation catalyst, I now repeat this…

2 months ago

3D Molecular model visualisation: 3 Million atoms +

In the late 1980s, as I recollected here the equipment needed for real time molecular…

3 months ago

The Macintosh computer at 40.

On 24th January 1984, the Macintosh computer was released, as all the media are informing…

3 months ago