Temperature-driven protonated glycolaldehyde chemistry

The question

Reaction mechanisms are often built from optimized reactants, intermediates and transition states. This gives precise structural and energetic information and works especially well when the important states are already known. The modelling becomes more difficult when several intermediates and competing channels may be relevant and their importance changes with temperature.

In that situation, the useful question is not only whether a pathway is possible. It is also which pathways carry most of the reaction under the conditions being studied.

What this approach adds

The study starts from finite-temperature ab initio molecular dynamics and uses enhanced sampling to explore the reactive configurational space. Structures visited by the trajectories are grouped into chemical states, and the observed transitions between those states are used to construct temperature-dependent reaction networks.

This does not replace stationary-point and barrier calculations. It adds a complementary view: the relevant states and connections are inferred from the sampled dynamics rather than requiring the complete reaction network to be specified in advance. The analysis can then distinguish state stability, network connectivity and the fraction of reactive flux carried by different pathways.

The dominant pathway changes with temperature

For protonated glycolaldehyde formation, the pathway distribution changes markedly across the temperature range. At 400 K, the direct reactant-to-product route contributes less than one percent of the total reaction flux, and intermediate-containing routes dominate. From 514 to 2343 K, the direct transition contributes roughly half of the total flux. At higher temperature, additional competing product channels become accessible and divert part of the flux away from protonated glycolaldehyde.

The overall formation rate therefore does not simply increase with temperature. In the model studied in the paper, the reported forward rate reaches a maximum near 662 K and then decreases as the balance between productive and competing pathways changes.

A visible intermediate is not necessarily the bottleneck

The network also separates frequent pathway participation from kinetic control. One proton-bound intermediate exchanges rapidly with the reactant, while the subsequent step toward protonated glycolaldehyde is much slower. A state can therefore be prominent in the network without being the step that controls the overall formation rate.

Why this matters for computational R&D

For an experiment or process-development question, changing temperature or environment may alter more than the speed of one fixed mechanism. It can change which states are populated, which pathways dominate and where competing products begin to matter.

A trajectory-derived network provides a way to compare those changes within one finite-temperature framework. It is especially useful as a complement to high-accuracy electronic-structure calculations when a system contains several plausible intermediates or competing routes and the experimentally relevant mechanism may depend on operating conditions.

Connection to an experimental observable

A complementary 2025 collaboration measured gas-phase infrared spectra of protonated glycolaldehyde, acetic acid and methyl formate and compared the measurements with several levels of spectral calculation, including AIMD-based spectra. The reaction-network study addresses which protonated structures and pathways become relevant; the spectroscopy study provides an observable that can help distinguish the corresponding isomers experimentally.