Mapping the Kinetic and Thermodynamic Landscape of Formaldehyde Oligomerization under Neutral Conditions

Exploring the Core Formose Cycle: Catalysis and Competition

The core autocatalytic cycle of the formose reaction may be enhanced or eroded by the presence of simple molecules at life's origin. Utilizing quantum chemistry, we calculate the thermodynamics and kinetics of reactions both within the core cycle and those that deplete the reactants and intermediates, such as the Cannizzaro reaction. We find that via disproportionation of aldehydes into carboxylic acids and alcohols, the Cannizzaro reaction furnishes simple catalysts for a variety of reactions. We also find that ammonia can catalyze both in-cycle and Cannizzaro reactions while hydrogen sulfide does not; both, however, play a role in sequestering reactants and intermediates in the web of potential reactions.

Prebiotic Synthesis of Glycolaldehyde and Glyceraldehyde from Formaldehyde: A Computational Study on the Initial Steps of the Formose Reaction

In this work, the initial steps of the mechanism of the Formose reaction (FR) is computationally studied using DFT methods. The FR has been considered to be a relevant process in the prebiotic evolution leading to several types of sugars or carbohydrates. These molecules are some of the basic building blocks of the life. The dimerization of formaldehyde was found to take place via an intramolecular deprotonation reaction, leading to the formation of an intermediate which, after an isomerization, forms a Ca-complex of the cis-enediol tautomer of glycolaldehyde. The aldol reaction of this complex with additional formaldehyde gave glyceraldehyde, the simplest aldotriose. The catalyst Ca(OH)2 plays a dual role in the reaction, acting as a base (in the intramolecular deprotonation) and as Lewis acid (activating the carbonyl group) in the aldol addition.

Matrix isolation and photorearrangement of cis- and trans-1,2-ethenediol to glycolaldehyde

1,2-Ethenediols are deemed key intermediates in prebiotic and interstellar syntheses of carbohydrates. Here we present the gas-phase synthesis of these enediols, the high-energy tautomers of glycolaldehyde, trapped in cryogenic argon matrices. Importantly, upon photolysis at λ = 180-254 nm, the enols rearrange to the simplest sugar glycolaldehyde.

Preliminary Free Energy Map of Prebiotic Compounds Formed from CO2, H2 and H2S

What kinds of CHOS compounds might be formed in a prebiotic milieu by reducing CO2 in the presence of H2 and H2S? How might the presence of sulfur influence the chemical composition of the mixture? We explore these questions by using first-principles quantum chemistry to calculate the free energies of CHOS compounds in aqueous solution, by first generating a thermodynamic map of one- and two-carbon species. We find that while thiols are thermodynamically favored, thioesters, thioacids, and thiones are less favorable than their non-sulfur counterparts. We then focus on the key role played by mercaptoacetaldehyde in sulfur analogs of the autocatalytic formose reaction, whereby the thiol group introduces asymmetry and potential thermodynamic selectivity of some compounds over others.

Thermodynamics of Potential CHO Metabolites in a Reducing Environment

How did metabolism arise and evolve? What chemical compounds might be suitable to support and sustain a proto-metabolism before the advent of more complex co-factors? We explore these questions by using first-principles quantum chemistry to calculate the free energies of CHO compounds in aqueous solution, allowing us to probe the thermodynamics of core extant cycles and their closely related chemical cousins. By framing our analysis in terms of the simplest feasible cycle and its permutations, we analyze potentially favorable thermodynamic cycles for CO₂ fixation with H₂ as a reductant. We find that paying attention to redox states illuminates which reactions are endergonic or exergonic. Our results highlight the role of acetate in proto-metabolic cycles, and its connection to other prebiotic molecules such as glyoxalate, glycolaldehyde, and glycolic acid.

Identification of Glycolaldehyde Enol (HOHC═CHOH) in Interstellar Analogue Ices

Glycolaldehyde is considered the entry point in the aqueous prebiotic formose (Butlerow) reaction although it mainly exists in its unreactive hydrated form in aqueous solution. The characterization of the more reactive nucleophilic enol form under interstellar conditions has remained elusive to date. Here we report on the identification of glycolaldehyde enol (1,2-ethenediol, HOHC═CHOH) in low temperature methanol-bearing ices at temperatures as low as 5 K. Exploiting isotope labeling and isomer-selective photoionization coupled with reflectron time-of-flight mass spectrometry, our results unravel distinct reaction pathways to 1,2-ethenediol, thus demonstrating the kinetic stability, availability for prebiotic sugar formation, and potential detectability in deep space.

An Upgrade of Apparatus and Measurement Systems for Generation of Gaseous Formaldehyde: A Review

Formaldehyde (FA) is ubiquitous in the atmospheric environment. It is generally the dominant atmospheric carbonyl compound. Due to its well-known carcinogenicity, FA is a compound that arises the attention in the scientific community. In studies concerning the toxicological effects of FA on humans, animals, and the environment, testing and calibration of air sampling systems and analytical instruments are pivotal. Therefore, the preparation of controllable standard gaseous atmospheres containing FA at levels known with precision and accuracy is essential. This review summarizes the procedures for generating the FA atmosphere, given that operative solutions have been evolving recently. Furthermore, an overview on the available system to collect and store gaseous standard is reported. The progressively implemented FA generation techniques, together with commercially-available instruments, are herein described, classified, and compared.

Early Steps of Glycolonitrile Oligomerization: A Free-Energy Map

Building on our previous free-energy map (J. Phys. Chem. A 2018, 122, 6769-6779) examining the reactions of monomeric glycolonitrile, we explore the formation of its dimers and trimers in aqueous solution under neutral conditions. While 5-membered rings are kinetically favored, open-chain oligomers with ester or amide linkages are thermodynamically favored. Accessing the 5-membered rings also provides a potential route to glyoxal that bypasses preforming glycolamide, the thermodynamic sink for monomers. However, finding a kinetically accessible route to glycine starting from glycolonitrile in the absence of added ammonia at room temperature proved challenging; the best case involved an intramolecular nucleophilic substitution reaction in a dimer containing neighboring imine and amide groups. Our free-energy map also examines routes to experimentally proposed moieties, explaining why some are observed in low yield or not at all.

Systematic microsolvation approach with a cluster-continuum scheme and conformational sampling

Solvation is a notoriously difficult and nagging problem for the rigorous theoretical description of chemistry in the liquid phase. Successes and failures of various approaches ranging from implicit solvation modeling through dielectric continuum embedding and microsolvated quantum chemical modeling to explicit molecular dynamics highlight this situation. Here, we focus on quantum chemical microsolvation and discuss an explicit conformational sampling ansatz to make this approach systematic. For this purpose, we introduce an algorithm for rolling and automated microsolvation of solutes. Our protocol takes conformational sampling and rearrangements in the solvent shell into account. Its reliability is assessed by monitoring the evolution of the spread and average of the observables of interest.

Role of Acid in the Co-oligomerization of Formaldehyde and Pyrrole

Building on previous work (J. Phys. Chem. A 2017, 121, 8154-8166) under neutral conditions, we examined the co-oligomerization of CH₂O and pyrrole to form porphryinogen under acidic conditions using density functional theory (B3LYP//6-311G**). Thermodynamically, we found that azafulvene intermediates were significantly stabilized under highly acidic conditions. Kinetically, energy barriers were lowered for C-C bond formation, discriminating in favor of reactions that lead to porphyrinogen. However, it was challenging to satisfactorily combine our thermodynamic and kinetic profiles into a unified free-energy profile because of difficulties in optimizing transition states of cationic species involving proton hops. Instead, we used neutral carboxylic acids as a proxy to study how energy barriers changed. By combining data from both neutral and acidic conditions, we estimate a free-energy profile for the initial steps of oligomerization under milder acidic conditions more relevant to prebiotic chemistry.

Exploring Free Energy Profiles of Uracil and Cytosine Reactions with Formaldehyde

Simple polymers can be potentially formed by the co-oligomerization of pyrimidine nucleobases, uracil and cytosine, with the small molecule formaldehyde. Using density functional calculations, we have constructed a free energy map outlining the thermodynamics and kinetics for (1) the addition of formaldehyde to uracil and cytosine to form hydroxymethylated uracil (HMU) and hydroxymethylated cytosine (HMC), (2) the deamination of cytosine and HMC to uracil and HMU, respectively, and (3) the initial oligomerization of 5-HMU. For the initial formation of monomeric HMU, addition of formaldehyde to the C5 and C6 positions is thermodynamically favored over N1 and N3, but faces higher kinetic barriers, and explains why 5-HMU is the main product observed experimentally. Oligomerization of 5-HMU is thermodynamically favorable although decreasingly so at the tetramer stage. In addition, decreasing concentrations of initial monomer shifts the equilibrium disfavoring oligomer formation.

Multiscale molecular thermodynamics of graphene-oxide liquid-phase exfoliation

Liquid-phase exfoliation is one of the most feasible methods for mass-production of two-dimensional (2D) nanomaterials such as graphene, graphene-oxide (GO), etc. Assessing requirements for successful exfoliation necessitates molecular-level thermodynamic analysis that can provide quantitative measures such as free energy changes. Here we explain this methodology and apply it to the production of GO that is used as a precursor for graphene synthesis and as an ultrathin substrate for many applications. Three different routes to GO exfoliation are studied, namely parallel and perpendicular to the GO surface as well as exfoliation via edge bending, using multi-scale combination of density functional, force field, and continuum approaches. Detailed analysis of free energy variations reveals relative feasibility of different exfoliation mechanisms and their dependence on system size and surface coverage. The methodology is general and can be applied to liquid-phase exfoliation of other 2D nanomaterials.

1,3-Dioxolane-4-ol Hemiacetal Stores Formaldehyde and Glycolaldehyde in the Gas-Phase

We report the spontaneous gas-phase formation of 1,3-dioxolane-4-ol, a mixed hemiacetal resulting from the addition of glycolaldehyde to formaldehyde. It was spectroscopically characterized by matching matrix IR spectra with coupled cluster computations. The formation of the hemiacetal must be surface-catalyzed owing to the very high computed reaction barrier of 39.8 kcal mol^-1. The reaction barrier is lowered by almost 20 kcal mol^-1 when a single water molecule acts as a proton shuttle in a favorable six-membered transition state. We characterized the hemiacetal in solution via NMR spectroscopy and followed its decomposition into its constituents within a few hours; it also dissociates upon contact with water. Sugars form in the presence of Ca(OH)₂, in line with formose-type reactivity. 1,3-Dioxolane-4-ol may be considered a storage form for formaldehyde and glycolaldehyde that is rather stable in the gas-phase.

Reactions of Glycolonitrile with Ammonia and Water: A Free Energy Map

Glycolonitrile, the product of combining CH₂O and HCN, is an intermediate in the Strecker reaction leading to the synthesis of the amino acid glycine. However, besides glycine, a plethora of other compounds are also generated when CH₂O and HCN react in the presence of ammonia and water. As a starting point to analyze the possible components of this complex mixture, we have employed density functional theory to construct a free energy map of all two-carbon (C2) species that may be present when glycolonitrile participates in addition or elimination reactions with ammonia and water. By identifying thermodynamic sinks and kinetic barriers, we find that the myriad C2 species can be grouped into three broad regions across the free energy landscape. This allows us to trace possible routes to glycine and other molecules of interest in the reaction mixture. The present map also extends our previous work on one-carbon (C1) species. We had previously found one issue with our computational protocol in the C1 map; however, our present C2 map provides a larger data set that supports using an empirical correction to our original protocol for imidic acid to amide transformations, without increasing the computational cost, while retaining the original protocol for other classes of reactions.

Uncertainty quantification for quantum chemical models of complex reaction networks

For the quantitative understanding of complex chemical reaction mechanisms, it is, in general, necessary to accurately determine the corresponding free energy surface and to solve the resulting continuous-time reaction rate equations for a continuous state space. For a general (complex) reaction network, it is computationally hard to fulfill these two requirements. However, it is possible to approximately address these challenges in a physically consistent way. On the one hand, it may be sufficient to consider approximate free energies if a reliable uncertainty measure can be provided. On the other hand, a highly resolved time evolution may not be necessary to still determine quantitative fluxes in a reaction network if one is interested in specific time scales. In this paper, we present discrete-time kinetic simulations in discrete state space taking free energy uncertainties into account. The method builds upon thermo-chemical data obtained from electronic structure calculations in a condensed-phase model. Our kinetic approach supports the analysis of general reaction networks spanning multiple time scales, which is here demonstrated for the example of the formose reaction. An important application of our approach is the detection of regions in a reaction network which require further investigation, given the uncertainties introduced by both approximate electronic structure methods and kinetic models. Such cases can then be studied in greater detail with more sophisticated first-principles calculations and kinetic simulations.

Primary Formation Path of Formaldehyde in Hydrothermal Vents

Formaldehyde is abundant in the universe and one of the fundamental molecules for life. Hydrothermal vents produce a substantial amount of hydrogen molecules by serpentinization and promote reductive reactions of single carbon compounds. The abundance of formaldehyde is expected to be low due to the high Gibbs free energy in hydrothermal vents. We consider two competing formation pathways of formaldehyde: (1) the reduction of CO by H₂ and (2) the reduction of HCOOH by H₂ to form a methanediol, followed by the dehydration of the methanediol. We performed a number of quantum chemical simulations to examine the formation of formaldehyde in the gas phase as well as in aqueous solution. The energy barrier is significantly reduced by the catalytic effect of water molecules in aqueous solution and becomes lowest when a water cluster consisted of 5 water molecules catalyzes the reduction. The energy barrier to form a methanediol by the reduction of HCOOH is lower by 17.5 kcal/mol than that to form a formaldehyde by the reduction of CO. Considering the low energy barrier to dehydrate methanediol, the primary pathway to form formaldehyde in hydrothermal vents is concluded to be the reduction of HCOOH by H₂, followed by the dehydration of methanediol.