Posts Tagged ‘Hydrogen’
Sunday, March 24th, 2019
There is a predilection amongst chemists for collecting records; one common theme is the length of particular bonds, either the shortest or the longest. A particularly baffling type of bond is that between the very electronegative F atom and an acid hydrogen atom such as that in OH. Thus short C-N…HO hydrogen bonds are extremely common, as are C-O…HO.‡ But F atoms in C-F bonds are largely thought to be inert to hydrogen bonding, as indicated by the use of fluorine in many pharmaceuticals as inert isosteres.[1] Here I do an up-to-date search of the CSD crystal structure database, which is now on the verge of accumulating 1 million entries, to see if any strong C-F…HO hydrogen bonding may have been recently discovered.
The search query uses the CF…HO distance as one variable, and the C-F-H angle as the second. The first diagram shows just intermolecular interactions, up to a distance of 2.7Å which is the sum of the van der Waals radii of the two elements. The hot spot occurs at this value, and an angle of ~95°.
The intra-molecular plot shows a similar value for the most common F…H distance, with the interesting variation that the angle subtended at F is about 80°.
The outlier at the short end of the spectrum (arrow) was observed in 2014[2] with the structure shown below. It is indeed the current record holder by some margin! This length by the way is however a great deal longer than the shortest O…HO hydrogen bonds, which can be in the region of 1.2Å (with the proton sometimes symmetrically disposed between the two oxygen atoms). The value is also very similar to the record holder for the shortest C-H…H-C interaction.
It is always useful to check up on crystallographic hydrogen atom positions using a quantum calculation, so here is one at the ωB97XD/Def2-TZVPP level (Data DOI: 10.14469/hpc/5131) which replicates the values nicely.

ωB97XD/Def2-TZVPP Calculation
A QTAIM analysis of the critical points shows that the F…H BCP has a high value of ρ(r) (most hydrogen bonds only reach about 0.03 au).

NBO analysis indicates the E(2) perturbation energy for donation from an F lone pair into the H-O σ* orbital is 21.2 kcal/mol, which indicates a strong H-bond (typical C-O…HO values are 18-22 kcal/mol). The F…H bond order is 0.05.
This molecule has another interesting property, also noted in the original article;[2] the shift in wavenumber of the O-H stretching vibration. Most hydrogen bonds are characterised by the shift (mostly red and recently discovered blue shifts) that occurs in the OH group when it hydrogen bonds. These shifts are typically 100-200 cm-1 but in this molecule there is no shift, which is described as “exceptional”.
The 1H NMR shift of the OH proton is observed at δ 4.8 ppm, with the value calculated here (ωB97XD/Def2-TZVPP) being 4.75 ppm. A very large H-F coupling was observed of 68 Hz, again a very high value for a “through space” hydrogen bond.
So another record for the molecule makers to try to break!
‡Respectively 7142 and 31428 intermolecular (3859 and 10602 intra) examples using the same search parameters as above, with the shortest values being ~1.28 and ~1.2Å.

References
- S. Purser, P.R. Moore, S. Swallow, and V. Gouverneur, "Fluorine in medicinal chemistry", Chem. Soc. Rev., vol. 37, pp. 320-330, 2008. https://doi.org/10.1039/b610213c
- M.D. Struble, C. Kelly, M.A. Siegler, and T. Lectka, "Search for a Strong, Virtually “No‐Shift” Hydrogen Bond: A Cage Molecule with an Exceptional OH⋅⋅⋅F Interaction", Angewandte Chemie International Edition, vol. 53, pp. 8924-8928, 2014. https://doi.org/10.1002/anie.201403599
Tags:Chemical bond, chemical bonding, Chemical elements, Chemistry, Fluorine, Hydrogen, Hydrogen bond, Intermolecular forces, Natural sciences, perturbation energy, pharmaceuticals, Physical sciences, Refrigerants, search parameters, search query, Supramolecular chemistry
Posted in crystal_structure_mining | No Comments »
Wednesday, July 25th, 2018
Consider the four reactions. The first two are taught in introductory organic chemistry as (a) a proton transfer, often abbreviated PT, from X to B (a base) and (b) a hydride transfer from X to A (an acid). The third example is taught as a hydrogen atom transfer or HAT from X to (in this example) O. Recently an article has appeared[1] citing an example of a fourth fundamental type (d), which is given the acronym cPCET which I will expand later. Here I explore this last type a bit further, in the context that X-H bond activations are currently a very active area of research.

To help understand these four types, I have colour-coded the electron pair constituting the X-H covalent bond in red.
- In mechanism (a), this electron pair stays with X, thus liberating a proton which is captured by the base.
- The hydride transfer (b) is so-called because in fact this electron pair travels together with the proton, hence the term hydride or H–.
- Hydrogen atom transfers as in (c) in effect transfer both a proton and one electron to another atom (oxygen in the example above), leaving behind one electron on X. The electron and the proton are said to travel together as a “true” hydrogen atom.
- The fourth mechanism (d) is fundamentally different from (c) in that whilst the electron and the proton travel in concert (at the same time), they do not travel together. In this example the proton travels to the oxygen, whilst the electron travels to the iron centre, in the process reducing its oxidation state. This mode is now called a concerted proton-coupled electron transfer, or cPCET as above.
The tool employed to distinguish between mechanisms (c) and (d) is the IBO or intrinsic bond orbital localisation scheme.[2] One practical advantage of such a scheme over better known localisation methods such as NBO (Natural bond orbitals) is that IBOs can be made to transform smoothly during the course of a reaction, as followed by say an IRC (Intrinsic reaction coordinate). NBOs may instead show discontinuous behaviour along a reaction IRC. Klein and Knizia have located transition states for examples of both (c) and (d) above and studied the IBOs along such IRCs. The two IBO reaction transformations are very different, as illustrated below (used, with permission, from the article itself). For the HAT type (X=C above), an α-spin IBO morphs from a C-H bond into a H-O bond, whilst the β-spin counterpart morphs from being located on the C-H bond into a carbon-centered radical. For the cPCET mode, the α-spin IBO morphs from C-H to a C-centered radical, but the β-spin region grows onto an iron d-orbital. It is in fact even more complex than the diagram above implies, since some reorganisation of the O-Fe region occurs and the H…:O region is still anti-bonding at the transition state.

We can see from this that mechanistic reaction analysis is starting to track the “curly arrows” we conventionally use to represent reactions in some detail, as well as informing us about the relative detailing timing of the various curly arrows used. Of course this latter aspect cannot be easily represented by conventional curly arrows. It seems timely to revisit the vast corpus of organic and organometallic “curly arrow pushing” to starting adding such information!
References
- J.E.M.N. Klein, and G. Knizia, "cPCET versus HAT: A Direct Theoretical Method for Distinguishing X–H Bond‐Activation Mechanisms", Angewandte Chemie International Edition, vol. 57, pp. 11913-11917, 2018. https://doi.org/10.1002/anie.201805511
- G. Knizia, "Intrinsic Atomic Orbitals: An Unbiased Bridge between Quantum Theory and Chemical Concepts", Journal of Chemical Theory and Computation, vol. 9, pp. 4834-4843, 2013. https://doi.org/10.1021/ct400687b
Tags:chemical reactions, Chemistry, Deprotonation, Hydride, Hydrogen, Hydrogen atom abstraction, Proton, proton travel, Proton-coupled electron transfer, Technology/Internet
Posted in Interesting chemistry | No Comments »
Wednesday, April 18th, 2018
The molecules below were discussed in the previous post as examples of highly polar but formally neutral molecules, a property induced by aromatisation of up to three rings. Since e.g. compound 3 is known only in its protonated phenolic form, here I take a look at the basicity of the oxygen in these systems to see if deprotonation of the ionic phenol form to the neutral polar form is viable.

The equilibrium being considered is shown below for compound 2:
The energetics of this equilibrium shown below, computed at the ωB97XD/Def2-TZVPPD/SCRF=water level and for which the FAIR data DOI is 10.14469/hpc/4073
For 1: X=Cl, the energy is shown below as a function of the O….H distance. Proton abstraction from HCl is exothermic by ~25 kcal/mol.
For 2: X=Cl, the exothermicity increases by only ~5 kcal/mol , despite the apparent aromatisation of a further ring. It is also worth noting that this is greater basicity than that of e.g. water, where around 4-5 water molecules acting in concert are required to deprotonate HCl.
For 1: X=OH, the proton abstraction from water is mildly endothermic by about 13 kcal/mol; indeed there is no energy minimum for carbonyl protonation and instead a relatively strong hydrogen bond to the water is formed instead.

For 2: X=OH the endothermicity is reduced to ~9 kcal/mol.
For 1: X=CH3, deprotonation of methane is now strongly endothermic by ~40 kcal/mol.

So the molecules 1 – 2 above are clearly not superbases, which perhaps augers well for being able to deprotonate the ionic phenols into these neutral but highly polar molecules.
Tags:Antiseptics, Aromatization, Chemistry, energy, energy minimum, Hydrogen, Molecule, Neurotoxins, Science
Posted in Interesting chemistry | No Comments »
Friday, February 16th, 2018
Last year, this article[1] attracted a lot of attention as the first example of molecular helium in the form of Na2He. In fact, the helium in this species has a calculated‡ bond index of only 0.15 and it is better classified as a sodium electride with the ionisation induced by pressure and the presence of helium atoms. The helium is neither valent, nor indeed hypervalent (the meanings are in fact equivalent for this element). In a separate blog posted in 2013, I noted a cobalt carbonyl complex containing a hexacoordinate hydrogen in the form of hydride, H–. A comment appended to this blog insightfully asked about the isoelectronic complex containing He instead of H–. Here, rather belatedly, I respond to this comment!
The complex [HCo6(CO)15]– has a calculated bond index at the hydrogen of 0.988 and a calculated NMR chemical shift of 21.6 ppm (ωB97XD/Def2-TZVPPD calculation) compared to a measured value of 23.2 ppm. Despite being six-coordinate, the hydride has a bond index that does not exceed one (it is not hypervalent).
So here is the neutral helium analogue. The He bond index emerges as 0.71 at the geometry of the hydride complex. Compare this with the bond index of 0.15 calculated for Na2He and it would be fair to say that at this geometry, the helium in [HeCo6(CO)15] would have a greater claim to be a molecular compound. Back in 2010, extrapolating from a series of posts here, I had speculated[2] about other molecular species of He, including the di-cation below. This has a He bond index of 0.54, rather less than that in [HeCo6(CO)15] but much more than in Na2He. It is also vibrationally stable.
But now, [HeCo6(CO)15] goes “pear-shaped” (why do pears have such a bad press?). I started a process of optimizing the geometry of this complex (ωB97Xd/Def2-TZVPPD). Slowly, the He started to creep out of the centre of the complex and emerge from the cavity. After about 100 steps it reached the geometry shown below, at which point the Wiberg bond index has dropped to 0.62 and still going down. I think it might take a few more steps to be completely expelled, but I have stopped the geometry optimisation at this stage.

So helium appears not to be valent in [HeCo6(CO)15]. However, I have yet to try Ne, which is both larger and softer. I will post results here.
‡All data at 10.14469/hpc/3587.
References
- X. Dong, A.R. Oganov, A.F. Goncharov, E. Stavrou, S. Lobanov, G. Saleh, G. Qian, Q. Zhu, C. Gatti, V.L. Deringer, R. Dronskowski, X. Zhou, V.B. Prakapenka, Z. Konôpková, I.A. Popov, A.I. Boldyrev, and H. Wang, "A stable compound of helium and sodium at high pressure", Nature Chemistry, vol. 9, pp. 440-445, 2017. https://doi.org/10.1038/nchem.2716
- H.S. Rzepa, "The rational design of helium bonds", Nature Chemistry, vol. 2, pp. 390-393, 2010. https://doi.org/10.1038/nchem.596
Tags:chemical bonding, Chemical elements, chemical shift, Chemistry, helium, Hydride, Hydrogen, Hypervalent molecule, Matter, Metal hydrides, Reducing agents, Transition metal hydride
Posted in Hypervalency | 4 Comments »
Thursday, September 21st, 2017
A recent article reports, amongst other topics, a computationally modelled reaction involving the capture of molecular hydrogen using a substituted borane (X=N, Y=C).[1] The mechanism involves an initial equilibrium between React and Int1, followed by capture of the hydrogen by Int1 to form a 5-coordinate borane intermediate (Int2 below, as per Figure 11).‡ This was followed by assistance from a proximate basic nitrogen to complete the hydrogen capture via a TS involving H-H cleavage. The forward free energy barrier to capture was ~11 kcal/mol and ~4 kcal/mol in the reverse direction (relative to the species labelled Int1), both suitably low for reversible hydrogen capture. Here I explore a simple variation to this fascinating reaction.∞

This variation involves transposing X and Y such that Y=N+ and X=C– to form a carbon ylide such that X=C becomes much more nucleophilic than the original nitrogen nucleophile. An animation of the full IRC† (intrinsic reaction coordinate computed at ωB97XD/cc-pvtz; FAIR data doi: 10.14469/hpc/2704) is shown below.

The profile shows that the reaction is concerted between the species labelled React and Prod; no sign of Int1 and Int2!
- The region IRC -12 to -5 involves B-C bond cleavage. Because the C is so very nucleophilic, the 4-ring species labelled React becomes very stable and opening it requires a high barrier.
- Between IRC -5 and 0, the BH2 group rotates, losing its original interaction with the C to slowly create an empty acceptor orbital on the boron which can then interact with the incoming hydrogen.
- At IRC= 0 (the transition state) the hydrogen has been captured by the boron to form a 5-coordinate species, in a manoeuvre that reminds one of the orbital capture of satellites by planets on their way to the outer reaches of the solar system. If the barrier to this capture is computed from IRC= -4 (the region of Int2) it is very much lower than the original system[1], again a reflection of the higher nucleophilicity of X=C–.
- The fly past continues until IRC= +7, at which point one end of the bound hydrogen has become suitably orientated to interact with the nucleophilic carbon via lone-pair donation into the acceptor H-H σ* orbital, thus helping to break it.
- By IRC= +9, the H-H cleavage is complete.
- By IRC= +13 the reaction has reached Prod, being overall ~ -12 kcal/mol exothermic.
- The overall thermochemistry is dominated by the potent carbon nucleophile in the reactant, which in turn makes this modification entirely useless for the purposes of a hydrogen-capture system!


The evolution of the dipole moment along the IRC shows very non-linear behaviour (such plots are rarely shown in most published IRC analyses; they should be!), ending of course with the ionic zwitterion that is the imminium borohydride Prod. Indeed the entire reaction coordinate is an unusually vivid example of a non-least motion path!

This simple atom transposition has given us a very instructive exercise in reaction paths, by-passing entirely both Int1 and Int2 (making them hidden intermediates), and converting React → Prod into a concerted reaction. It would be great to probe this convoluted journey using reaction dynamics!
∞Archived as DOI: 10.14469/hpc/3096
‡ Such a species can be seen as a hidden intermediate in the mechanism of reduction of a carboxylic acid by diborane.
†None were shown in the original study.[1]
References
- L. Li, M. Lei, Y. Xie, H.F. Schaefer, B. Chen, and R. Hoffmann, "Stabilizing a different cyclooctatetraene stereoisomer", Proceedings of the National Academy of Sciences, vol. 114, pp. 9803-9808, 2017. https://doi.org/10.1073/pnas.1709586114
Tags:Ammonia borane, animation, Boranes, Chemistry, Cleaning Services, Company: React Group, free energy barrier, Hydroboration, Hydrogen, Matter
Posted in reaction mechanism | 1 Comment »
Monday, March 20th, 2017
The example a few posts back of how methane might invert its configuration by transposing two hydrogen atoms illustrated the reaction mechanism by locating a transition state and following it down in energy using an intrinsic reaction coordinate (IRC). Here I explore an alternative method based instead on computing a molecular dynamics trajectory (MD).
I have used ethane instead of methane, since it is now possible to envisage two outcomes:

An animation of the IRC starting from the located transition state is shown below (DOI: 10.14469/hpc/2331). This is based purely on the computed potential energy surface of the molecule. The IRC is computed from the forces experienced on the atoms as they are displaced from an initial set of coordinates corresponding to the located transition state and then following the direction indicated by the eigenvectors of the negative force constant required of a transition state. Importantly, there is no time component; the path is based entirely on energies and forces.
Next, a molecular dynamics simulation (ωB97XD/6-31G(d,p), DOI: 10.14469/hpc/2330). This uses the ADMP method, which requests a classical trajectory calculation using the “atom-centered density matrix propagation molecular dynamics model”. This integrates kinetic energy contributions from the molecular vibrations and so explicitly now includes a time component. In this example the evolution of the system from the transition state is charted over a period of 100 femtoseconds (1000 integrated steps). As it happens this is a relatively short period of evolution; sometimes periods of picoseconds may be required to get a realistic model, especially if one is also dealing with explicit solvent molecules (of which perhaps 500 might be required).
Such explicit inclusion of the kinetic energy from molecular vibrations in effect allows the molecule to “jump” over shallow barriers. In this case, the barrier for a [1,2] hydrogen shift from the methyl group to the adjacent carbene (watch atom 8). Simultaneously, the path taken by two hydrogens no longer corresponds to their transposition but to their elimination as a hydrogen molecule. So this quite different outcome from the IRC is very probably also a much more realistic one.
If the MD method is so much more realistic than the IRC, then why is it not always used? The simple answer is computational time! For this very small molecule and using quite a modest basis set (6-31G(d,p)), the relatively short 1000 time steps took about three times as long to compute as the IRC. The factor gets worse as the size of the molecule increases and the number of time steps for a realistic result increases. Perhaps, as technology gets better and new computer architectures emerge, MD simulations of ever increasingly complex reactions will become far more common. In ten years time, I expect most of the examples on this blog will use this method!
Tags:animation, chemical reaction, Chemistry, computational chemistry, computed potential energy surface, energy, Gaseous signaling molecules, Hydrogen, kinetic energy, kinetic energy contributions, Methane, Molecular dynamics, Physical chemistry, Quantum chemistry, Reaction coordinate, simulation, Theoretical chemistry
Posted in reaction mechanism | 2 Comments »
Wednesday, February 15th, 2017
This post arose from a comment attached to the post on Na2He and relating to peculiar and rare topological features of the electron density in molecules called non-nuclear attractors. This set me thinking about other molecules that might exhibit this and one of these is shown below.
The topology of the electron density is described by just four basic types, designed formally by the notation [3,-3], [3,-1], [3,1] and [3,3] and more colloquially by the terms nuclear attractor (NNA), line (or bond) critical point, a ring critical point and a cage critical point respectively. Mostly, the nuclear critical points coincide exactly with the actual nuclear positions, but more rarely these points are not found centered at a nucleus. It was such an NNA that was suggested as a comment on the post on Na2He. There I replied that another example of an NNA is to be found in H3+ and so its a short step to take a look at H42+ in a tetrahedral arrangement (DOI: 10.14469/hpc/2165). Since only two electrons are available for bonding, it is tempting to represent it as below, with dashed partial bonds connnecting the six edges of the tetrahedron and is associated with real normal vibrational modes; ν 416, 1490 and 1861 cm-1. A brief search on Scifinder, which appears to reference this species as hydrogen, ion (H42+), does not identify any publications associated with it (there are studies on H41+ however); if any reader here knows of any discussion please alert us!

Analysing the density however gives a different result. A NNA is located at the centre of the tetrahedron and a line (bond) critical point connects this to each of the four hydrogen nuclei. This result is similar to the obtained for H3+. It is rather odd that these non-nuclear attractors have not entered into the vocabulary used to describe the bonding in simple molecules, but this picture is certainly different from the more empirical dashed lines between the four nuclei that one is instinctively drawn to use (above).

The ELF analysis (performed using multiWFN) is more interesting. The nuclear basins associated with the hydrogens reveal each has 0.425e, but the central one (green below) has its own basin with 0.301e.

The NICS value associated with the non-nuclear attractor is -27 ppm, which is indicative of strong spherical aromaticity.
All of which goes to show that even the simplest of molecular species may still have properties that are unexpected or certainly not well-known!
Tags:Attractor, brief search, Chemistry, Electron, Electron density, Hydrogen, Molecule, Nature, Physics, Quantum chemistry
Posted in Interesting chemistry | 11 Comments »
Saturday, February 11th, 2017
On February 6th I was alerted to this intriguing article[1] by a phone call, made 55 minutes before the article embargo was due to be released. Gizmodo wanted to know if I could provide an (almost)† instant‡ quote. After a few days, this report of a stable compound of helium and sodium still seems impressive to me and I now impart a few more thoughts here.
The discovery originates from 17 authors based in 17 different institutions, an impressive illustration of global science and cooperation. I illustrate with this diagram, to be found not in the main article body but in its supporting information and for which the caption reads:

Computed charge density (eÅ-3) of Na2He at 300 GPa, plotted in the [110] plane of the conventional cell. The color bar gives the scale.
The nuclei carry of course the greatest charge density, but the density labelled “2e” is not nuclear-centered. This is typical of species known as electrides, where positive cations are associated with just electrons acting as the counter-anion and about which there was an extensive debate earlier on this blog. There is much discussion in the article[1] about the essential role of the He atoms in bringing about the formation of such an electride, an effect that is summarised in a second diagram also found in the supporting information:

I found myself thinking that it would be great to have the first diagram represented as a movie, evolving as the pressure is increased from say ambient to 300 GPa, and presumably showing the “2e” feature (which means diamagnetic electrons) forming as the pressure increases. Would their evolution be abrupt (a step change) or gradual as the pressure increases and the interatomic distances all decrease? As I understand it, this chemical phenomenon is due not so much to the usual coulombic attraction between positive nuclei and negative charge density from the electronic wavefunction leading to e.g. covalent bonds, but to electron repulsions induced by decreasing nuclear separations resulting in electride-like ionisation and hence electron localisation into the “interstitial cavities” of the lattice. Without pressure, you would just have sodium and helium atoms!
The urge to obtain this intriguing electronic wavefunction for myself now appeared (wavefunctions are rarely if ever included in supporting information). To do this you must have atom coordinates available, But such data was not to be found in the supporting information. It was eventually tracked down (by a crystallographer; thanks Andrew!) to the caption in Figure 2.

However, you probably do need to be a crystallographer to convert this data into a set of coordinates. This was done and is here deposited as a CIF file for you to play with if you wish (DOI:10.14469/hpc/2154)[2]. I have reduced the packing of the unit cell obtained from this CIF file (198 atoms) to just 60 and you can enjoy them by clicking on the diagram below. I should point out that if one uses a program that can recognise the periodic lattice such as Crystal (used in the article discussed here), there is no need to make such reductions, but in this instance I wanted to use a program such as Gaussian in discrete (non-periodic) mode, for which the calculation (B3LYP/Def2-SVPD) has DOI: 10.14469/hpc/2156[3] and where you can also find a wavefunction file to play with if you wish.

Click for 3D model
An ELF analysis for this non-periodic wavefunction looks as below. The ELF basins labelled “2e” located in the centre of the cube show an integrated electron population of ~1.9e and correspond to the localised electron pairs noted in the article above.

Click for 3D
The basins on the boundaries of this non-periodic unit show reduced integrations (red arrows below, 0.08 – 1.7e) and are artefacts of the non-periodic approximation introduced.

The ionization into an electride is brought about by the close proximity of the atoms as induced by high pressure. Releasing the pressure would allow the ionized electrons to re-attach themselves to the valence shell of the sodium atoms, thus destroying the unique properties of the system. It is certainly true that this system challenges our normal concepts of what a molecule is. The presence of He is essential and yet its electrons are hardly involved in the re-organised wavefunction. I cannot wait for more examples to be discovered!
†To meet the 55 minute deadline, I was given about 15 minutes thinking time!
‡Instant responses on social media now seem a sine qua non of the political world, so why not the scientific one?!
References
- X. Dong, A.R. Oganov, A.F. Goncharov, E. Stavrou, S. Lobanov, G. Saleh, G. Qian, Q. Zhu, C. Gatti, V.L. Deringer, R. Dronskowski, X. Zhou, V.B. Prakapenka, Z. Konôpková, I.A. Popov, A.I. Boldyrev, and H. Wang, "A stable compound of helium and sodium at high pressure", Nature Chemistry, vol. 9, pp. 440-445, 2017. https://doi.org/10.1038/nchem.2716
- H. Rzepa, "Na2He: a stable compound of helium and sodium at high pressure.", 2017. https://doi.org/10.14469/hpc/2154
- H. Rzepa, "He20Na40", 2017. https://doi.org/10.14469/hpc/2156
Tags:10.1038, Atom, Chemical elements, chemical phenomenon, Chemistry, Company: P. Acucar-CBD, Electride, Electron, Food Retail & Distribution - NEC, helium, Hydrogen, Matter, Oxygen, Physics, social media
Posted in Bond slam, crystal_structure_mining, Interesting chemistry | 11 Comments »
Saturday, December 31st, 2016
My holiday reading has been Derek Lowe’s excellent Chemistry Book setting out 250 milestones in chemistry, organised by year. An entry for 1920 entitled hydrogen bonding seemed worth exploring in more detail here.
As with many historical concepts, it can often take a few years to coalesce into something we would readily recognise today, and hydrogen bonds are no exception. Wikipedia is another source of the history and it cites a 1912 article as the origin of the term in relation to the solvation of amines[1] but also notes that the better known setting of water occurs later in 1920.[2] Here I try to capture the essence of the concept with a few diagrams taken from these two articles.
Firstly “The state of amines in aqueous solution“[1] which is mostly concerned with the measurement of ionization constants of primary, secondary and tertiary amines. It boils down to the below:


and the connection to ionization is laid out as:


Since in 1912, Lewis’ electron pair theory of the covalent bond had not yet emerged, the authors use the terms “strong union” and “weak union”, and of course it is the “weak union” that we now know of as the hydrogen bond. Some other comments about this seminal diagram:
- The article contains the very explicit and modern term stereochemical, which is used in a manner that suggests it was already common.‡ But there is only a hint at most that the nitrogen atoms might be tetrahedral, or that the “weak union” between (what we now think of as the lone pair on) the nitrogen and the hydrogen of the water is directional.
- The second weak union between the tetramethyl ammonium (which we now describe as a cation) and the hydroxide (now described as an anion; both terms are however implied by the description strong electrolyte) is probably not what we would now call a hydrogen bond, more an intimate ion-pair.
The second article in 1920 on water itself[2] is post-Lewis, but perhaps applied in a manner which we would not entirely agree with nowadays. Thus dinitrogen, N≡N is shown as below with just a single connecting bond.

Then we get the interaction between ammonia and water,† analogous to the example shown above:

and for water itself:♠

which in each case shows the central hydrogen having what we now call a valence shell of four electrons,♣ and hence more equivalent to the “strong unions” above. This shows that in 1920 chemists were rapidly adopting Lewis’ representations, but not always entirely successfully.
On balance, I think the 1912 article sets out the modern concept of a hydrogen bond representing a weak union to a hydrogen rather better than the Latimer and Rodebush attempt (at least diagrammatically).
‡Stereochemical notation is discussed in this post, and it dates from the 1930s.
†The modern take is explored here, in which the equilibrium set up between a “weak union” between ammonia and water (the weak electrolyte) and an isomeric “strong union” which represents ionization into an ammonium hydroxide ion-pair (the strong electrolyte) is favoured for the former by ΔG ~6 kcal/mol.
♠The equilibrium between a “weak union” of two water molecules and the fully ionized strong union of hydronium hydroxide favours the former by ΔG ~23 kcal/mol.
♣ This 1920 representation does imply symmetry for the hydrogen, being ~equally disposed between the two oxygens. We now know that such symmetric hydrogen bonding is not unusual (see this post for how to fine-tune a hydrogen bond into this situation) but rather than requiring four electrons as implied in the diagram above, it is now better described as a three-centre-two-electron bond instead.
References
- T.S. Moore, and T.F. Winmill, "CLXXVII.—The state of amines in aqueous solution", J. Chem. Soc., Trans., vol. 101, pp. 1635-1676, 1912. https://doi.org/10.1039/ct9120101635
- W.M. Latimer, and W.H. Rodebush, "POLARITY AND IONIZATION FROM THE STANDPOINT OF THE LEWIS THEORY OF VALENCE.", Journal of the American Chemical Society, vol. 42, pp. 1419-1433, 1920. https://doi.org/10.1021/ja01452a015
Tags:10.1021, aqueous solution, Chemical bond, chemical bonding, Chemistry, Derek Lowe, Hydrogen, Hydrogen bond, Intermolecular forces, Lowe's, Nature, Supramolecular chemistry
Posted in Historical | 2 Comments »
Saturday, December 31st, 2016
My holiday reading has been Derek Lowe’s excellent Chemistry Book setting out 250 milestones in chemistry, organised by year. An entry for 1920 entitled hydrogen bonding seemed worth exploring in more detail here.
As with many historical concepts, it can often take a few years to coalesce into something we would readily recognise today, and hydrogen bonds are no exception. Wikipedia is another source of the history and it cites a 1912 article as the origin of the term in relation to the solvation of amines[1] but also notes that the better known setting of water occurs later in 1920.[2] Here I try to capture the essence of the concept with a few diagrams taken from these two articles.
Firstly “The state of amines in aqueous solution“[1] which is mostly concerned with the measurement of ionization constants of primary, secondary and tertiary amines. It boils down to the below:


and the connection to ionization is laid out as:


Since in 1912, Lewis’ electron pair theory of the covalent bond had not yet emerged, the authors use the terms “strong union” and “weak union”, and of course it is the “weak union” that we now know of as the hydrogen bond. Some other comments about this seminal diagram:
- The article contains the very explicit and modern term stereochemical, which is used in a manner that suggests it was already common.‡ But there is only a hint at most that the nitrogen atoms might be tetrahedral, or that the “weak union” between (what we now think of as the lone pair on) the nitrogen and the hydrogen of the water is directional.
- The second weak union between the tetramethyl ammonium (which we now describe as a cation) and the hydroxide (now described as an anion; both terms are however implied by the description strong electrolyte) is probably not what we would now call a hydrogen bond, more an intimate ion-pair.
The second article in 1920 on water itself[2] is post-Lewis, but perhaps applied in a manner which we would not entirely agree with nowadays. Thus dinitrogen, N≡N is shown as below with just a single connecting bond.

Then we get the interaction between ammonia and water,† analogous to the example shown above:

and for water itself:♠

which in each case shows the central hydrogen having what we now call a valence shell of four electrons,♣ and hence more equivalent to the “strong unions” above. This shows that in 1920 chemists were rapidly adopting Lewis’ representations, but not always entirely successfully.
On balance, I think the 1912 article sets out the modern concept of a hydrogen bond representing a weak union to a hydrogen rather better than the Latimer and Rodebush attempt (at least diagrammatically).
‡Stereochemical notation is discussed in this post, and it dates from the 1930s.
†The modern take is explored here, in which the equilibrium set up between a “weak union” between ammonia and water (the weak electrolyte) and an isomeric “strong union” which represents ionization into an ammonium hydroxide ion-pair (the strong electrolyte) is favoured for the former by ΔG ~6 kcal/mol.
♠The equilibrium between a “weak union” of two water molecules and the fully ionized strong union of hydronium hydroxide favours the former by ΔG ~23 kcal/mol.
♣ This 1920 representation does imply symmetry for the hydrogen, being ~equally disposed between the two oxygens. We now know that such symmetric hydrogen bonding is not unusual (see this post for how to fine-tune a hydrogen bond into this situation) but rather than requiring four electrons as implied in the diagram above, it is now better described as a three-centre-two-electron bond instead.
References
- T.S. Moore, and T.F. Winmill, "CLXXVII.—The state of amines in aqueous solution", J. Chem. Soc., Trans., vol. 101, pp. 1635-1676, 1912. https://doi.org/10.1039/ct9120101635
- W.M. Latimer, and W.H. Rodebush, "POLARITY AND IONIZATION FROM THE STANDPOINT OF THE LEWIS THEORY OF VALENCE.", Journal of the American Chemical Society, vol. 42, pp. 1419-1433, 1920. https://doi.org/10.1021/ja01452a015
Tags:10.1021, aqueous solution, Chemical bond, chemical bonding, Chemistry, Derek Lowe, Hydrogen, Hydrogen bond, Intermolecular forces, Lowe's, Nature, Supramolecular chemistry
Posted in Historical | 2 Comments »