A 1965 precedent to the Dewar Pyrimidone MOST system – and text book examples of the Woodward-Hoffmann pericyclic reaction selection rules

In the previous post,[1] I noted the photochemical isomerisation of a pyrimidone into what is called the bicyclic Dewar form, being part of a solar energy storage system.[2] A colleague (thanks Alan!) has recollected a very similar example dating from 1965[3] in which a related molecule known as a diazepinone 72 (scheme below) is converted by light into a Dewar form 73.

This example was first highlighted in Woodward and Hoffmann’s (WH) famous 1971 book on the topic of the conservation of orbital symmetry in which they noted that the Dewar form of a diazepinone (73 in scheme) had been observed[3] to thermally “revert to diazepinone in the dark”. The original authors[3] also specifically noted that the Dewar diazepinone was “stable to storage” after being protonated. These two properties are the exact inverse of the recent report,[2] whereby the photochemical bicyclic form of pyrimidone was found to be thermally stable, but very rapid ring opening was induced by protonation with acid. Here I explore whether these apparently contradictory reports can be reconciled.

In discussing the reaction of 72 in their book, WH suggest that the stereochemical aspects of the thermal ring opening of 73 could be explained using their rules by prior inversion of the ring nitrogen stereochemistry to that of 73-inv, followed by conrotatory/antarafacial ring opening to 72. Here, with the help of ωB97XD/Def2-TZVPP/DCM DFT calculations,[4] I discuss whether this suggestion is viable, and also propose an alternative mechanism (72-trans, Scheme above).

Firstly, I show the calculated reaction path[5],[6] along which HTS3 and 73-inv are found, being the WH suggestion for this reaction.



Figure 1. IRC Energy plot and animation for TS3

  1. At IRC ~8, (Figure 1) the potential shows what can be called a “hidden transition state”, at which point the gradient norm is close to zero. This is the point labelled HTS3, followed soon after by a “hidden intermediate” (IRC ~4) or 73-inv. The process corresponds to inversion of the nitrogen lone pair to produce a bicyclic species with a trans ring fusion. These are both “hidden” because the gradient norm (Figure 2) does not actually reach a value of 0.0 as required for “real” transition states and intermediates, but comes very close.
  2. ΔG298 at these points relative to the starting point is ~34 kcal/mol, rather higher than would be needed for a truly thermal reaction. The CN bond length has not yet started to change (Figure 3).
  3. At IRC = 0.0 the true transition state is reached (TS3), involving WH-allowed antarafacial cleavage (Figure 5) of the bicyclic C-N bond (length @TS 2.035Å). The energy is now ~65 kcal/mol above the starting point, which makes this pathway very unlikely.
  4. The thermal reaction is exothermic by -19 kcal/mol (Figure 1), significantly less than that for Dewar pyrimidone.

Figure 2. Gradient norm plot for TS3


Figure 3. C-N bond length plot for TS3


Figure 4. Dipole moment plot for TS3, just for fun!


Figure 5. Geometry of TS3, showing C-N bond with antarafacial component (top face connecting bottom face) corresponding to conrotation (both clockwise) of the two termini.

Next, I tried an alternative mechanism, involving direct ring opening via TS1 to give a 7-ring with a trans bond, 72-trans. [7],[8] (Figure 6). Back in 1971, 7-rings with trans bonds were a rarity, so WH were probably reluctant to suggest this.


Figure 6. Energy plot for TS1

  1. The activation energy (corresponding to ΔG298 27.84 kcal/mol) is looking much better, matching to a slow (hours, days) thermal reaction at room temperatures. This value is somewhat less than the value of 32.9 kcal/mol for the analogous ring opening of Dewar pyrimidone,[1] probably because the larger 5-ring ring means less transition state strain.
  2. The reaction again occurs with conrotation/antarafacial (Figure 7), C-N 2.192Å.
  3. But it is now endothermic by about +15 kcal/mol, reflecting the relatively high energy of a 7-ring product with a trans bond (Figure 6).


Figure 7. Geometry of TS1, showing C-N bond with antarafacial component (top face connecting bottom face) corresponding to conrotation (both clockwise) of the two termini.

To complete the mechanism, a route must now be found to convert 72-trans back to 72 itself.

  1. This can be done via a linear arrangement of the C-N-N atoms[9] but the barrier to doing so is prohibitive (ΔG298 39.7 kcal/mol).
  2. An alternative is direct rotation about the C=N bond via an allylic biradical transition state (TS2),[10],[11] which yields ΔG298
    26.74 kcal/mol. This value is less than that for TS1, and so is not rate determining.

Figure 8. Energy plot for TS2

When TS1 is protonated, ΔG298 becomes 26.4 kcal/mol (Figure 9, [12],[13], C-N 2.198Å) compared to the unprotonated value of ΔG298 27.8 kcal/mol. The slight decrease in barrier upon protonation does not match the observation[3] that protonated 73 is “stable to storage”. This still leaves open the question of why computations indicate that the rate of ring opening of Dewar diazepinone is relatively unchanged by protonation, whereas that of Dewar pyrimidone is greatly accelerated – the former involves protonating a hydrazine whereas the latter involves protonating an amide. Further models will need investigating to confirm whether this accounts for the essential difference in behaviour.


Figure 9. Energy plot for TS1 upon protonation.

To conclude, WH’s suggestion of a nitrogen inversion mechanism for the slow thermal pericyclic reaction of 72 followed by conrotatory C-N ring opening is instead replaced here by one invoking the electrocylic formation of a 7-ring intermediate with a trans bond and then biradical rotation of this bond.


Woodward, R. B.; Hoffmann, Roald (1971). The Conservation of Orbital Symmetry (3rd printing, 1st ed.). Weinheim, BRD: Verlag Chemie GmbH (BRD) and Academic Press (USA). pp. 1–178. ISBN 978-1483256153. The kinetics of this process were not noted, nor was the temperature.

Author

References

  1. H. Rzepa, "A breakthrough in Molecular Solar Thermal (MOST) energy storage – Dewar Pyrimidone.", 2026. https://doi.org/10.59350/jhsbq-sfs70
  2. H.P.Q. Nguyen, A.J. Maertens, B.A. Baker, N.M. Wu, Z. Ye, Q. Zhou, Q. Qiu, N. Kaur, D.B. Berkinsky, K.E. Shulenberger, K.N. Houk, and G.G.D. Han, "Molecular solar thermal energy storage in Dewar pyrimidone beyond 1.6 megajoules per kilogram", Science, vol. 392, 2026. https://doi.org/10.1126/science.aec6413
  3. W.J. Theuer, and J.A. Moore, "Heterocyclic studies. The photoisomerization of 2,3-dihydro-5-methyl-6-phenyl-4H-1,2-diazepin-4-one and derivatives", Chemical Communications (London), pp. 468, 1965. https://doi.org/10.1039/c19650000468
  4. H. Rzepa, "WH-73", 2026. https://doi.org/10.14469/hpc/15948
  5. H. Rzepa, "[Embargoed]", 2026. https://doi.org/10.14469/hpc/15975
  6. H. Rzepa, "73 TS for trans geometry N-epimer ( G =-648.883115 => G = -648.88455 DG = 65.9 IRC mirror image", 2026. https://doi.org/10.5281/zenodo.20455914
  7. H. Rzepa, "73 TS for cis geometry suprafacial, N epimer G = -648.951789 (vs G = -648.933983 for NH isomer ) DG = 27.84 IRC", 2026. https://doi.org/10.14469/hpc/15960
  8. H. Rzepa, "73 TS for cis geometry suprafacial, N epimer G = -648.951789 (vs G = -648.933983 for NH isomer ) DG = 27.84 IRC", 2026. https://doi.org/10.5281/zenodo.20279953
  9. H. Rzepa, "73 product isomerism G = -648.932844, IRC", 2026. https://doi.org/10.14469/hpc/15963
  10. H. Rzepa, "73 cis product rotation, G=-648.953525, DG = 26.74", 2026. https://doi.org/10.14469/hpc/15967
  11. H. Rzepa, "&3 cis product rotation, G=-648.953525", 2026. https://doi.org/10.5281/zenodo.20406989
  12. H. Rzepa, "73 TS for cis geometry N-protonated G = -649.363457 DG = 30.99 => NH epimer G = -649.370744 DG = 26.41 IRC", 2026. https://doi.org/10.5281/zenodo.20474465

3 Responses to “A 1965 precedent to the Dewar Pyrimidone MOST system – and text book examples of the Woodward-Hoffmann pericyclic reaction selection rules”

  1. Jacob Manaker says:

    From citation [12], it looks like you assumed the protonation occurs at the tertiary nitrogen. That seems natural enough, but I’m hung up on the missing pH dependence during decomposition.

    The azetine is a phenyl-conjugated enamine. Is it possible that the protonation actually occurs at the (phenyl’s) benzylic position? That location is at least as electron-rich and strained as the nitrogen, and carbon is a more electropositive element.

    I don’t know a good way to determine pKas computationally. A quick calculation I did says both the nitrogen lone pair and the alkene are in the HOMO, but the algorithm I used can get orbital orderings a little wrong.

    I haven’t (yet?) tried to reproduce your reaction trajectory with the alternate protonation and see if it has a larger barrier.

  2. Henry Rzepa says:

    Thanks Jacob. re protonation, I did try the alternative alpha-N as a source of protonation, which turned out to be higher in energy, but did not try the benzylic carbon. Nitrogen is almost always more basic, which is what protonation requires and hydrazines can be more basic than amines due to the “alpha” effect. It is certainly true that more possibilities need exploring, but what we want here is to “increase” the barrier to decomposition after protonation, and at first sight that is not so obvious how to do.

    Perhaps this is a challenge for AI? After all, Tim Gowers has just announced ( https://gowers.wordpress.com/2026/05/08/a-recent-experience-with-chatgpt-5-5-pro/ ) that chatGPT 5.5 Pro has solved a “PhD level” problem requiring genuinely creative reasoning. Perhaps it can solve this one too?

  3. Henry Rzepa says:

    I noted in my previous comment that LLMs such as ChatGPT 5.5 Pro seem capable of solving math problems which an eminent mathematician deems worthy of a PhD level student (citing expts he ran which took around 15 minutes of processing time!). It strikes me that the above reaction of a diazepinone might be sufficiently well defined for AI to have a go it. The WH rules have an associated large corpus of mostly reliable and well defined data, and probably an AI inference agent could itself discover at least two routes from compound 73 to 72 based on the rules. Whether it could determine which of them is lower in energy might arise out of much more data being known about eg the stability of trans stereochemistry in small rings than was available back in 1965 or 1971, Finally, are the results for protonation consistent, or could the need to perform more experiments be flagged instead?

    Is there anyone out there more familiar with ChatGPT and the like than I am who wants to have this as a challenge?

Leave a Reply