Computational quantum chemistry has made fantastic strides in the last 30 years. Often deep insight into all sorts of questions regarding reactions and structures of molecules has become possible. But sometimes the simplest of questions can prove incredibly difficult to answer. One such is how accurately can the boiling point of water be predicted from first principles? Or its melting point? Another classic case is why mercury is a liquid at room temperatures? The answer to that question (along with another, why is gold the colour it is?) is often anecdotally attributed to Einstein. More accurately, to his special theory of relativity.[cite]10.1002/andp.19053231314[/cite] But finally in 2013 a computational proof of this was demonstrated for mercury.[cite]10.1002/anie.201302742[/cite] The proof was built up in three stages.
- Potential energy surfaces for the Hg2 dimer showed that inclusion of a relativistic Hamiltonian contracts the Hg-Hg distance by 0.2Å. This can be traced back to the value of the 1S0(6s2)→3P0(6s16p11/2) electronic excitation in the atom being 4.67eV, compared to a non-relativistic value of 3.40eV. This in turn is the result of the strong relativistic 6s shell contraction and hence stabilisation. But it has previously been shown that bulk mercury cannot be described by such a simple two-body interaction.
- Next many-body clusters of various sizes were built. This is a complex task, since for each size, various types of packing might be possible. The largest was a “two-layer Mackay icosahedron”, with 55 Hg atoms. These clusters however did not show any monotonic convergence to a clear melting point.
- Finally, using Monte Carlo (MC) simulations within a quantum diatomics-in-molecules (DIM) model and with periodic boundary conditions to simulate the bulk metal, it was possible to show that without a relativistic Hamiltonian, the predicted melting point was predicted as 355K (82°C) but when the relativistic effects were switched on, this decreased to 250K (-23°C). This lowering of 105K is dominated by scalar relativistic effects through many-body contributions.
- The experimental melting point is 234K. The density is well predicted as well (the non-relativistic model predicts mercury to be denser than it actually is).
What I did not get from this article is why mercury is such a very special case (i.e. why neither gold, m.p. 1337K nor thallium, m.p. 577K, are liquids at room temperature). No doubt someone will explain. In the past, gold and mercury were said to be the only two visual manifestations of Einstein’s special theory in every-day objects. I say the past, because mercury is now rarely seen in any every-day objects (digital thermometers have taken over, and the mercury barometer has long since gone). If anyone knows of other examples, do let us know.
Note added in 2019: Similar calculations have now been performed on Copernicium[cite]10.1002/anie.201906966[/cite] predicting it too is a liquid (MP 283 ± 11K) and has noble gas character (group 18) due to both relativistic and dispersion effects. Without those, it would simply be a common group 12 element.
Tags: bulk metal, pence, potential energy surfaces
Nice post.
Spin-orbit coupling is a purely relativistic effect.
It gives rise to magnetic anisotropy and hysteresis in ferromagnets.
They are “everyday” phenomena.
A colleague has kindly brought to our attention the observation that whilst mercury is unusual amongst the elemental metals, there exists a class of non-relativistic molecular metal, of which a member Li(NH3)4 is even odder (DOI: 10.1021/ja109397k) in melting at 89K! These are described as expanded metals, since “the distance between the molecular bearers of the eventually itinerant electrons becomes large”.
There can be (much) more to the term metal than meets the eye.
In the last week, the first example of a relativistic effect on the properties of a molecule (as opposed to an element) has been published, DOI: 10.1126/science.1255720. The molecule in question is Sg(CO)6. I quote here a pertinent passage from the article itself:
because of the relativistic contraction of s and p1/2 orbitals, the d-orbitals are better shielded from the atomic nucleus and expand. The overlap with the LUMO of the CO ligand becomes stronger. This may lead to a stronger contribution of π-back-donation relative to σ-donation in comparison with the corresponding ratio in the complexes of the lighter members of group 6.
Only by including full relativistic effects can quantum mechanics predict the hexacarbonyl to be both stable and volatile, which was indeed verified by experiments. Presumably without relativistic effects, the hexacarbonyl would have been predicted to be unstable, easily dissociating to a lower-coordination carbonyl and CO gas.
Someone correct me who has a better understanding, but gold is not liquid because its 6s valence shell is not filled and is therefore more susceptible to valence bonding (the single electron has a higher energy to allow overlap with other atoms while at the same time the atom is accepting of another electron to fill the electron pair). Thallium has an unfilled p-orbital, which has a very different shape than the s orbitals of the valence electrons for gold and mercury. It's shape reduces the exposure of the electron to the effects of the heavy nucleus and frees it up valence bonding with neighbors. In thallium's case, this bond does not fulfilll any specific electon pair, so it is not as strong as the valence bonds found in gold. This at least superficially explains why gold has a higher melting point than thallium.