Acid-Catalyzed Ring Opening of Furan in Aqueous Solution - Energy

Dec 20, 2017 - School of Chemical and Biomolecular Engineering, The University of Sydney, Sydney, New South Wales 2006, Australia ... furanic compound...
0 downloads 6 Views 3MB Size
Subscriber access provided by UNIV OF TASMANIA

Article

ACID-CATALYSED RING OPENING OF FURAN IN AQUEOUS SOLUTION Xiao Liang, Brian S. Haynes, and Alejandro Montoya Energy Fuels, Just Accepted Manuscript • DOI: 10.1021/acs.energyfuels.7b03239 • Publication Date (Web): 20 Dec 2017 Downloaded from http://pubs.acs.org on December 28, 2017

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

Energy & Fuels is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 29 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

ACID-CATALYSED RING OPENING OF FURAN IN AQUEOUS SOLUTION

Xiao Liang, Brian S. Haynes, and Alejandro Montoya* School of Chemical and Biomolecular Engineering The University of Sydney Sydney, NSW 2006, Australia

ABSTRACT The mechanism and energetics of furan decomposition via the opening of the five-member ring structure in dilute acid solution were investigated using density functional Car-Parrinello molecular dynamics combined with metadynamics simulations. Diffusion of an acidic proton from the aqueous medium leading to the formation of a protonated furan is found to be the rate-limiting step of the ring opening process, with protonation at the Cα position being 7 kcal mol–1 less activated than at the Cβ position. Protonation at the Cα position results in the formation of 2,5 dihydro-2-furanol or 2,3 dihydro-3-furanol via the nucleophilic attack by a solvent molecule. Subsequent protonation of the furanols at the ring oxygen initiates the opening of the furanic ring leading to the formation of 4-hydroxy-2-butenal.

*

To whom correspondence should be addressed. E-mail: [email protected]

ACS Paragon Plus Environment

Energy & Fuels 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

1. INTRODUCTION

Furans represent a family of chemicals that can be derived from cellulose and hemicellulose. Furfural and 5–hydroxymethylfurfural (5–HMF) are valuable furanic compounds obtained from renewable resources because they are essential building blocks for other chemicals with wide applications in the automotive and chemical industry.1,2 These compounds are produced by dehydration of monosaccharides such as glucose, fructose and xylose. Furans are susceptible to decomposition via opening of their five–member ring structures.3-6 The decomposition reaction pathways become predominant especially in the presence of acid catalysts.7-10 Therefore, elucidation of the reaction mechanisms and kinetics of the furanic ring opening in acidic aqueous solution is fundamental to gaining insight into the process of sustainable production of chemicals from biomass.

The stability of furan and its derivatives in acidic solutions has been studied under various acid concentrations and temperatures. The stability of the furanic ring depends significantly on the nature of its substituents. Furanics with electron-attracting substituent groups are rather stable in acidic solution,11 while those with electron-donating groups are readily converted into their corresponding 1,4-diketones.12,13 In hot dilute aqueous acids, it has been accepted that the major reaction of furanics is hydrolytic ring opening that starts from protonation of one of the ring carbon atoms.14

The selectivity of the protonation site of a furanic ring has been debated in the literature. A rate-limiting protonation at the Cβ positions was proposed on the basis of several experimental studies of acid-catalysed ring cleavage of furan, 2–methoxy–, 2–methyl– and 2,5–dimethyl furans,15-17 Although the product distribution was subsequently revised,18 the Cβ protonation is used to explain the formation of levulinic acid and polymerisation from 5– ACS Paragon Plus Environment

Page 2 of 29

Page 3 of 29 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

HMF.19,20 On the other hand, Stamhuis et al. have proposed that protonation at the Cα position is rate limiting among the ring opening reactions of furan and 2,5–dimethyl furan in aqueous HClO4 solution at 25ºC.21 Additional kinetic isotope experiments confirmed that hydrogen exchange between solvent and Cα follows a two-step A–SE2 mechanism.16,22 Two separate experimental studies of furan in mildly acidic solution23 and 2–methyl furan in dioxane–D2O solution with dilute HClO4 acid22 supported the Cα–protonation pathway. Nevertheless, kinetic studies22,23 suggests that the pre-equilibrium of Cα–protonation is faster than the subsequent cleavage of the furanic ring, conflicting with Stamhuis’ observations.

The activation energy for α–hydrogen exchange of 2–methyl furan in a water–dimethyl sulphoxide mixture with dilute HClO4 between 70 ~ 100 ºC was determined to be 21.7 ~ 24.1 kcal mol–1.16 In contrast, no hydrogen exchange between solvent and Cβ was observed.16,23 It is proposed that the Cα positions are more readily protonated in both gas phase and aqueous solution24,25, mainly due to a higher stability of the Cα–furanium ion, which has greater charge delocalization than its isomer Cβ–furanium ion.26 The ring oxygen has lower proton affinity than the Cα and Cβ positions.27 It has been shown that protonation at the ring oxygen alters only slightly the electron density distribution of the furanic ring, thereby reducing the possibility of a subsequent nucleophilic attack.28,29

While furan is detected only in low concentration during the decomposition of monosaccharides in acid solutions, its ring configuration provides a framework to study the furanic ring opening reaction without the interference of formyl and hydroxymethyl groups associated with the more prevalent furanic products. Therefore, furan has been selected as the model furanic compound in molecular modelling studies. In the gas phase, the relevant stabilities of furanium ions were computed at the G3 level of theory, showing that the most favoured protonation position is Cα atoms, which is 10.5 and 27.7 kcal mol–1 more ACS Paragon Plus Environment

Energy & Fuels 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

exothermic than that at Cβ and hetero oxygen atoms, respectively.30 The kinetics of Cβ– protonation of furan were separately studied in the gas phase at the B3LYP/6-311G(d,p) level of theory31, yielding the result that the proton transfer from an adjacent hydronium ion to the Cβ position was 5.3 kcal mol–1 exothermic and nearly barrierless (0.2 kcal mol–1).31 The energy barrier for the proton shift from α– to β– position, assisted by one water molecule, was 16.7 kcal mol–1.31 The reaction mechanism and energetics of Brønsted acid catalysed hydrolysis of dimethylfuran in aqueous solution has been recently explored using electronic structure calculations at the M062X/6-311+G(d,p) theory of level using an implicit water solvent model.32 Although the activation energy of α– protonation was computed to be 3 kcal mol–1 lower than that of β– protonation, it was proposed that the favoured reaction pathway follows a general acid catalysis mechanism by rate-limiting Cβ protonation. This mechanism was proposed because of the absence of an energy barrier in the water binding reaction of β– furanium ion and the low energy barrier (9 kcal mol–1) of the subsequent ring opening step, compared with a much higher activation energy of 27 kcal mol–1 for the ring opening of α– furanium ion.32

A broad energy mapping of possible reaction pathways initiating at the Cβ and Cα position of a furanic ring in solution using explicit water molecules is required to reveal the important dynamic effect of water molecules in the protonation of furanics, and therefore, to obtain additional information of the rate–limiting step of the ring opening process. In this work, the mechanisms and energetics of the acid–catalysed ring cleavage of furan is explored using density functional Car-Parrinello molecular dynamics (CPMD) combined with metadynamics. This computational approach has successfully been applied to some organic reaction systems and has been proved to provide accurate energetic information for aqueous phase reaction systems.33-37 The ring opening of furan will be explored in two steps. The first step involves

ACS Paragon Plus Environment

Page 4 of 29

Page 5 of 29 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

the protonation of furan and hydration of the furanium species. Protonation at the α– position leads to the formation of 2,5 dihydro-2-furanol and 2,3 dihydro-3-furanol, while protonation at the β– position leads to the formation of 2,3 dihydro-2-furanol as shown in Scheme 1. The second step comprises the acid–catalysed ring opening via C–O bond dissociation of the furanol species to acyclic products.

Scheme 1

2. COMPUTATIONAL METHODOLOGY

The molecular dynamics simulations were carried out using the CPMD programme.38 The Goedecker pseudopotentials39 and the Becke, Lee, Yang, and Parr (BLYP) functional40,41 were used to describe the potentials exerted by the core electronic shells of ions and the valence electrons, respectively. A plane wave cut-off of 70 Ry was used for the Kohn-Sham orbitals. The equations of motion were integrated with a time step of 4 au (0.1 fs) and a fictitious mass of 800 amu for electrons was used. An NVT ensemble was used at 300 K with a Nose–Hoover chain thermostat.42 A periodic unit cell with cubic dimensions of 11 Å containing 35 H2O molecules, one furanic species, one hydronium ion (H3O+) and one chloride anion was used. The corresponding density of the solution is 939 g L–3 and the concentration of H3O+ is 1.25 mol L–1. The cell with the H2O molecules was pre-equilibrated using a Monte Carlo method.43 Then, an optimised furan or furanol molecule, together with one hydronium and one chloride species were inserted into the box containing H2O molecules. The isolated furanic species were optimised using Gaussian 0944 at B3LYP/6-311++G(d,p) level; a ten–picosecond CPMD simulation was then conducted on the created supercell to equilibrate the system before performing the metadynamics simulations. The equilibrium of

ACS Paragon Plus Environment

Energy & Fuels 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 6 of 29

the systems was examined by plotting the root mean square deviation (RMSD) of each trajectory against simulation time, as presented in Figure S1 in Supporting Information. The RMSD was observed to fluctuate in a small range (± 0.15 Å) after four picoseconds in all cases.

In the first set of calculations, furan was allowed to evolve into dihydrofuranols via acidcatalysed furan hydration reactions. To map the free energy landscape of this process, we have used two collective variables (CV) involving the protonation C–H (CV1) and the reacting C–OH bonds (CV2) as presented in Figure 1. In the second set of calculations, the dihydrofuranols were allowed to evolve into the corresponding acyclic products via ring opening reactions. The free energy landscape of the acid-catalysed ring openings was evaluated using two CVs that describe the protonation of the heteroatom O1 (CV3) and the cleavage of the cyclic ether bond (CV4), as shown in Figure 1. All hydrogen atoms of the solvent system are considered in mapping the changes of CV1 and CV3.

Figure 1

In each simulation, the CV is described as the sum of the Coordination Numbers (CN) between relevant atoms. A CN between two atoms i and j is determined from equation E.1:45

,  =

1 −  ⁄  





1 −  ⁄  

(E.1)

Here, dij is the distance between the corresponding atoms; d0 is the reference distance with the defined values of 2.0 Å for the distance between oxygen and carbon atoms, and 1.4 Å for the distance between oxygen and hydrogen atoms. The high power constants p = q = 6 are

ACS Paragon Plus Environment

Page 7 of 29 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

used to distinguish between the coordinated and non-coordinated states. The CN value between two atoms indicates whether a covalent bond exists (1 = bond, 0 = no bond). The parameters of fictitious masses (m) and force constants (k), which control the dynamics of the fictitious CVs, are selected to be 100 amu and 3.0 au for CV1 and CV3, and 600 amu and 8.0 au for CV2 and CV4 in all metadynamics computations. These parameters have been shown to be able to adiabatically separate the ionic and electronic degrees of freedom and restrain the CVs close to their coupled fictitious variables (auxiliary particles) of organic reaction systems46,47, and have been previously used to study similar reaction systems.36,37

A normal Gaussian bias potential was used in our metadynamics simulations with the width of Gaussian potential set to 0.05. The heights of the Gaussian hills were determined by the shape of the underlying energy surface and set to be at most 0.6 kcal mol–1 in all simulations except for 2,5 dihydro-2-furanol ring opening reaction. The average hill heights in these simulations were in the order of 0.4 kcal mol–1. The maximum height of Gaussian hills for the simulation of 2,5 dihydro-2-furanol ring opening was set to be 0.2 kcal mol–1 due to the relatively low energy barrier of this reaction. The addition of bias potential was allowed with a minimum time separation of 100 molecular dynamics (MD) steps (~ 10 fs) and with the displacements in CVs being larger than 1.5 times the width of the Gaussian function. The addition of consecutive potentials was forced after 300 MD steps (~ 30 fs) if the displacement had not exceeded the threshold.

An average of the free energies over several independent simulations of the metadynamics biasing potential yields very accurate free energies.48 However, this strategy is impractical for a large molecular system like the ones presented here due to the computational resources required. In our study, one metadynamics simulation was computed for each reaction step. Several earlier studies have shown that the free energies computed from a single

ACS Paragon Plus Environment

Energy & Fuels 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 8 of 29

metadynamics profile are reasonably accurate if the control parameters are selected properly, and the studied CV space is well sampled.49-53 Before reconstructing the free energy landscapes with the history-dependent bias, the evolution of each CV was examined to ensure that the CV spaces of most interest are well sampled. Graphs of the CV variation as a function of metadynamics steps and the sampling of CV trajectories for each metadynamics simulation can be found in the Supporting Information (Figure S2 to Figure S11).

The statistical error from a single metadynamics sampling can be estimated following Bussi et al54:

*+ ,+ -+

. exp −  S  w δ! √2π % /+  〈ϵ 〉 =   & βDτ S π k 

(E.2)

,1

Where, S is the size of the free energy well, w is the Gaussian hill height, δs is the hill width, τg is the metadynamics time step, D is the diffusion coefficient of the CV and can be calculated from the integration of the velocity autocorrelation function of the CV. The smaller diffusion coefficient of the two CVs in each metadynamics simulation is used for the determination of error. d is the dimensionality of the CV space, k is a possible d dimensional vector of integers with non-zero norm, and β is equal to (kBT)–1, in which kB is the Boltzmann constant and T is the temperature. Using these relevant parameters determined from each metadynamics simulation, the statistical uncertainties in our computations are estimated to be in the range of 0.4 ~ 1.5 kcal mol–1. Some intermediate states in this study are identified based on small free energy differences that are close to the corresponding metadynamics errors, raising the possibility that these intermediates may not be real stable states during reactions.

ACS Paragon Plus Environment

Page 9 of 29 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

In addition to the statistical errors, it is known that the BLYP exchange-correlation functional systematically underestimates the reaction energy barriers. In our simulated molecular systems, the energy barriers are associated with proton transfer and formation/dissociation of carbon-oxygen bonds. The BLYP uncertainties in computations of the energy barriers of formation/dissociation of carbon-oxygen bonds were reported to be 2 ~ 5 kcal mol–1.55,56 However, this functional predicts the energy barriers of proton transfer comparable to highlevel ab initio computations (MP2 and MP4 methods).57-59 Other errors involved in our CPMD computations may be related to the pseudopotential for the interaction between core and valence electrons and the control parameters, such as plane wave cut-off, fictitious electronic mass and time step. These uncertainties were estimated to be 2 ~ 3 kcal mol–1 in earlier studies with the same pseudopotential and similar control parameters as used in our present work.60,61 In the present work, we are interested in comparing the relative free energy difference between α– and β– protonation of furan to determine the preferred reaction pathway of acid–catalysed ring opening of furan rather than a direct comparison with experimental activation energies.

The molecular graphics of selected snapshots were rendered with VMD.62 A hydrogen bond between a molecular pair is identified when the oxygen-oxygen (O···O) distance is less than 3.5 Å and, at the same time, the angle ∠O–H···O is greater than 120º. These geometric criteria have been used in previous studies of aqueous phase systems.63,64 The average bond lengths, and their corresponding standard deviations, of the important states along each studied reaction pathway, are determined from all the configurations falling into the region where the free energy is within ~ 1.5 kcal mol–1 of each minimum.

ACS Paragon Plus Environment

Energy & Fuels 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

3 RESULTS AND DISCUSSIONS

3.1 Acid–catalysed hydration of furan

The free energy landscape of protonation of furan at the Cα position leading to the formation of 2,5 dihydro-2-furanol, together with selected molecular configurations along the reaction pathway is presented in Figure 2. Selected bond lengths and corresponding standard deviations (in Å) are provided in each configuration. The process involves a proton transfer from the solvent to the Cα position, followed by a nucleophilic attack of H2O and regeneration of the acid proton.

Figure 2

Furan is solvated by two layers of water molecules with the acid proton in the second solvation shell in the form of an H5O2+ species. The molecular configuration R1 shows that the acid proton connects with O1 via five hydrogen bonds. However, a direct interaction between carbon atoms of furan and solvent molecules is not observed. The hydrophobicity of furan is consistent with its low solubility in water (1 wt% at 25 ºC).65 The closest distance between the hydrogen of H2O and the ring oxygen O1 is 2.62 Å indicating a weak interaction. The C2–O1 and C5–O1 bond lengths are 1.39 and 1.40 Å, slightly longer than the reported gas phase experimental value of 1.364 Å.66

The protonation of furan proceeds via migration of the acid proton from the solvent towards the Cα position, leading to the formation of α–protonated furan. The energy basin R2 in Figure 2 represents a state of α–furanium and TS1 represents the transition state configuration

ACS Paragon Plus Environment

Page 10 of 29

Page 11 of 29 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

leading to this energy basin. The free energy change in moving from R1 to R2 (R1→R2) is 18.1 kcal mol–1 endergonic, and the free energy barrier R1→TS1 is 25.1 kcal mol–1.

The bond length H1–C5 is 1.21 Å in TS1 and 1.10 Å in R2, significantly shorter than in R1. The H1–C5 bond length in R2 is in close agreement with the reported gas phase value of 1.09 Å for α-furanium.25 The C5–O1 bond increases from 1.40 Å in R1 to 1.50 Å (7.1%) in TS1 and further elongates to 1.53 Å (9.3%) in R2, while the C2–O1 decreases to 1.32 Å (5.0%) in TS1 and R2. The small differences in geometrical parameters between TS1 and R2 indicate that α-protonation of furan proceeds via a “late” transition state.

There are considerable geometrical changes of the furanic ring upon protonation; mainly due to an electron redistribution in the furanic ring, thereby activating the second Cα (C2) and the Cβ position (C4) for a nucleophilic attack of H2O. The nucleophilic attack on the second Cα position and regeneration of the acid proton produces 2,5 dihydro-2-furanol, corresponding to state PI1 in Figure 2. Inspection of the free energy landscape shows that the formation of the furanol is 22.2 kcal mol–1 activated (R1→TS2) and 4.3 kcal mol–1 endergonic (R1→PI1). The average H···OH distance in TS2 of the reacting H2O is 1.17 Å, much longer than the corresponding length in pure liquid water.67 This indicates that the nucleophilic attack occurs concurrently with the regeneration of the acid proton. The C2–O1, C5–O1 and H1–C2 bond lengths in TS2 are very close to those in R2, suggesting an “early” transition state during the hydration process.

Figure 3 shows the free energy landscape and selected molecular configurations of protonation of furan at the Cα position leading to the formation of 2,3 dihydro-3-furanol. The mechanism of the formation of the α-furanium is similar to that in the formation of 2,5 dihydro-2-furanol discussed above, the process being 16.9 kcal mol–1 endergonic (R1'→R2')

ACS Paragon Plus Environment

Energy & Fuels 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

and 23.8 kcal mol–1 activated (R1'→TS1'). These two values are respectively 1.2 and 1.3 kcal mol–1 lower than those determined during the formation of 2,5 dihydro-2-furanol, but these differences are within the estimated statistical errors (0.4 ~ 1.5 kcal mol–1), as discussed in Section 2. Inspection of R1', TS1' and R2' in Figure 3 shows that these molecular configurations are essentially the same as those observed during the formation of 2,5 dihydro2-furanol, i.e. the furan configuration (R1') connects via five hydrogen bonds with the acid proton, with a distance between the closest H2O and O1 in the order of 2.73 Å. The closely similar free energy profiles and molecular configurations of the Cα–protonation process in the two distinct simulations indicate a good consistency and energy convergence in these calculations.

Figure 3

A nucleophilic attack of H2O on the adjacent Cβ position of the α–furanium (R2') regenerates the acid proton forming 2,3 dihydro-3-furanol (PI2). The nucleophilic attack is concerted with the regeneration of the acid proton, i.e. the average H···OH distance in TS3 is 1.19 Å. This process is 23.5 kcal mol–1 activated (R1'→TS3) and 6.1 kcal mol–1 endergonic (R1'→PI2). The free energy barrier corresponding to the formation of 2,3 dihydro-3-furanol (R1'→TS3) is only 1.3 kcal mol–1 greater than that predicted for the formation of 2,5 dihydro2-furanol (R1→TS2), indicating that there is little kinetic preference for either reaction pathway. Thermodynamically, the formation of 2,3 dihydro-3-furanol (R1'→PI2) is slightly less favoured because it is 1.8 kcal mol–1 more endergonic than the formation of 2,5 dihydro2-furanol (R1→PI1). Considering that the energy differences between these two reaction pathways are close to the uncertainties of the computational method, the reaction pathways will depend on the stabilities of the various furanols.

ACS Paragon Plus Environment

Page 12 of 29

Page 13 of 29 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

Figure 4 presents the free energy landscape and selected molecular configurations corresponding to the protonation of furan at the Cβ position, followed by a nucleophilic attack of water at an adjacent Cα position forming 2,3 dihydro-2-furanol.

Figure 4

There is no clear basin state for the formation of the β–furanium species. The distance between C2 in TS4 and the closest water molecule is in the order of 2.53 Å, significantly shorter than the corresponding lengths between water and α–furanium (3.22 Å in R2, Figure 2). The formation of the β–furanium species (R1''→TS4) is 30.5 kcal mol–1 activated, which is 7 kcal mol–1 more activated than the formation of α–furanium. The energy gain after passing TS4 is 24.7 kcal mol–1, which is larger than the energy gain from α–furanium to the corresponding furanols as seen in Figure 2 and 3, further indicating that β–furanium is less stable than α–furanium in solution. A ring cleavage mechanism via β-protonation and the possibility of an internal proton shift from α– to β–position is unfavourable due to the instability of the β–furanium configuration. Therefore, we predict that 2,5 dihydro-2-furanol and 2,3 dihydro-3-furanol are the main products of furan reactions with acid protons and water.

3.2 Acid-catalysed ring opening of furanol intermediates

The free energy landscapes and some selected geometrical parameters of the acid–catalysed ring openings of 2,5 dihydro-2-furanol and 2,3 dihydro-3-furanol as a function of proton transfer from the H2O to the furanic ring are presented in Figures 5 and 6, respectively.

ACS Paragon Plus Environment

Energy & Fuels 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 5

The interaction of H2O with 2,5 dihydro-2-furanol (PI1' in Figure 5) is significantly stronger than with furan, i.e. the closest distance is in average 0.72 Å shorter than the corresponding values in furan. PI1' is connected with the acid proton via four hydrogen bonds, equilibrated in the second solvation shell. The ring opening initiates with an increased degree of interaction between O1 and the proton. This interaction leads to the protonation of the ring oxygen O1 and the dissociation of C2–O1 bond (TS5), which is 13.4 kcal mol–1 activated. Inspection of TS5 shows that a new HP–O1 bond is almost formed (1.07 Å), indicating that the proton has been transferred from bulk H2O to the ring. The C2–O1 bond in TS5 increased from 1.48 Å in PI1' to 2.05Å, suggesting that TS5 is formed via a concerted process involving dissociation of C2–O1 and protonation of O1. The system gains 10.8 kcal mol–1 after passing TS5 and reaches the energy basin of 4-hydroxy-2-butenal. The average distance of C2–O in 4-hydroxy-2-butenal is 1.28 Å, indicating that an aldehyde group has been completely formed. Formation of 4-hydroxy-2-butenal is nearly thermo–neutral, being only 2.6 kcal mol–1 endergonic.

Figure 6 shows the free energy landscape of the ring opening of 2,3 dihydro-3-furanol (PI2' in Figure 6). The distance between the first water shell and the furanol is 1.93Å, which is only 0.07 Å shorter than in 2,5 dihydro-2-furanol. However, the acid proton is equilibrated at five hydrogen bonds from 2,3 dihydro-3-furanol, one more hydrogen bond apart than in the case of 2,5 dihydro-2-furanol. The ring opening proceeds via a single step PI2'→TS6→P2. The energy barrier is 34.2 kcal mol–1 and the process is 14.2 kcal mol–1 endergonic. The ring opening occurs via the C2–O1 bond dissociation. The bond length of the new HP–O1 in TS6 is 1.12 Å and the length C2–O1 is 2.06 Å, indicating a concerted protonation and ring

ACS Paragon Plus Environment

Page 14 of 29

Page 15 of 29 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

opening process. An H2O molecule at a point near C2 with an average distance of 2.35 Å is required to stabilise the TS6 configuration.

The direct reaction product is 3,4-dihydroxy-butenol with the proton equilibrated in the second solvation shell. During our metadynamics simulations, it was observed that 3,4dihydroxy-butenol is rapidly converted to 4-hydroxy-2-butenal after it was formed. The conversion proceeded by losing the –OH group connected to C3 and the HP connected to O1. The conversion process from 3,4-dihydroxy-butenol to 4 hydroxy-2-butenal is not shown in the free energy landscape in Figure 6 since this was not clearly distinguished based on the currently selected CV3 and CV4. This reaction was not studied further since the formation of 3,4-dihydroxy-butenol is less thermodynamically favoured than the formation of 4-hydroxy2-butenal from 2,5 dihydro-2-furanol.

Figure 6

Figures 2-6 indicate that furan decomposition in aqueous solution is driven by the formation of the α-protonated furan. The computed free energy barrier of ~24 kcal mol–1 is consistent with reported activation energies of acid–catalysed furanic ring reactions, e.g. 22.8 kcal mol–1 for 2,5–dimethyl furan and 26.0 kcal mol–1 for furan in 0.2 mol L–1 HClO4 solution15, 17 kcal mol–1 for furan in 4.61 mol L–1 HClO4 solution,23 29.9 kcal mol–1 for furfural in dilute H2SO4 solution,8 and 22.8 ~ 26.4 kcal mol–1 for 5–HMF in various dilute acid solutions.9,68-70 These different activation energies suggest an important effect of acid concentration or changes in the reaction mechanisms due to the side chains of the furanic ring.

ACS Paragon Plus Environment

Energy & Fuels 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

It should be noted that the furan and acid proton are both located inside a finite unit cell in our simulations. The free energy landscapes thus inevitably include the interactions between furan and proton species solvated in a second or third solvation shell. Previous metadynamics studies of acid-catalysed sugar reactions have shown that the proton transportation from bulk solvent to the proximity of sugar molecules is endergonic, therefore contributing to the overall reaction energy barriers.36,37,60 Therefore, additional static calculations were carried out to estimate the energy difference between an equilibrated furan in the presence of a solvated H3O+ with a system in which furan and proton are individually solvated. The gas phase equilibrium conformations were computed using B3LYP/6-311++G(d,p) level of theory and the free energies of solvation were calculated at the same level of theory using the SMD solvation model.71 The detailed description of the method for computing the free energies of infinite separation can be found in earlier publications.36,37 The energy of bringing together a solvated proton with furan is endothermic by 7 kcal mol–1, attributed to a greater proton attraction ability of bulk water compared with furan. The global free energy barrier in our simulations is therefore 7 + 23.8 (Figure 2) = ~30 kcal mol–1 referenced from an infinite separation of a solvated furan from the acid proton.

4. CONCLUSION

Molecular dynamics simulations were carried out to determine the reaction pathways and corresponding free energy landscapes for the acid-catalysed ring opening of furan in the aqueous phase. We observed that furan interacts with water only via a weak hydrogen bond linked to its ring oxygen. Protonation of furan is endergonic, indicating a greater proton attraction by bulk water molecules than furan. Protonation at the Cα position is energetically

ACS Paragon Plus Environment

Page 16 of 29

Page 17 of 29 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

preferred to that at the Cβ position. The α-protonation leads to the formation of 2,5 dihydro-2furanol and 2,3 dihydro-3-furanol. The resultant 2,5 dihydro-2-furanol is subject to the protonation on the ring oxygen triggering a ring opening to produce 4-hydroxy-2-butenal. In contrast, the ring opening of 2,3-dihydro-3-furanol is energetically less favoured due to its high free energy barrier. We predict that a primary product of decomposition of furan in acid solution is 4-hydroxy-2-butenal.

ACKNOWLEDGEMENT This work has been funded by the Australian Research Council through Discovery Grant DP1096802. The authors acknowledge the University of Sydney HPC resources and the National Computational Infrastructure Facility (NCI) supported by the Australian Government for providing resources that have contributed to the results reported in this paper.

REFERENCES

1 2

3

4

Lange, J.-P., van der Heide, E., van Buijtenen, J. & Price, R. Furfural-a Promising Platform for Lignocellulosic Biofuels. ChemSusChem 5, 150-166, doi:10.1002/cssc.201100648 (2012). Rosatella, A. A., Simeonov, S. P., Frade, R. F. M. & Afonso, C. A. M. 5Hydroxymethylfurfural (HMF) as a building block platform: Biological properties, synthesis and synthetic applications. Green Chem. 13, 754-793, doi:10.1039/c0gc00401d (2011). Kabyemela, B. M., Adschiri, T., Malaluan, R. M. & Arai, K. Glucose and fructose decomposition in subcritical and supercritical water: detailed reaction pathway, mechanisms, and kinetics. Industrial & Engineering Chemistry Research 38, 2888-2895 (1999). Chuntanapum, A., Yong, T. L.-K., Miyake, S. & Matsumura, Y. Behavior of 5-HMF in Subcritical and Supercritical Water. Industrial & Engineering Chemistry Research 47, 29562962 (2008).

ACS Paragon Plus Environment

Energy & Fuels 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

5 6 7 8

9

10 11 12 13 14 15 16

17 18 19 20

21

22

23 24

25

26

27

Aida, T. M. et al. Reactions of D-fructose in water at temperatures up to 400 DegC and pressures up to 100MPa. J. Supercrit. Fluids 42, 110-119 (2007). Aida, T. M. et al. Dehydration of D-glucose in high temperature water at pressures up to 80MPa. J. Supercrit. Fluids 40, 381-388 (2007). Girisuta, B., Janssen, L. P. B. M. & Heeres, H. J. Green chemicals a kinetic study on the conversion of glucose to levulinic acid. Chem. Eng. Res. Des. 84, 339-349 (2006). Marcotullio, G., Tavares Cardoso, M. A., De Jong, W. & Verkooijen, A. H. M. Bioenergy II: furfural destruction kinetics during sulphuric acid-catalyzed production from biomass. Int. J. Chem. React. Eng. 7, No pp given (2009). Asghari, F. S. & Yoshida, H. Kinetics of the Decomposition of Fructose Catalyzed by Hydrochloric Acid in Subcritical Water: Formation of 5-Hydroxymethylfurfural, Levulinic, and Formic Acids. Industrial & Engineering Chemistry Research 46, 7703-7710 (2007). Rose, I. C., Epstein, N. & Watkinson, A. P. Acid-Catalyzed 2-Furaldehyde (Furfural) Decomposition Kinetics. Industrial & Engineering Chemistry Research 39, 843-845 (2000). Dunlop, A. P. Furfural formation and behavior. Journal of Industrial and Engineering Chemistry (Washington, D. C.) 40, 204-209 (1948). Skrabal, A. & Skrabal, R. The velocity of hydrolysis of enol ethers. Z. physik. Chem. A181, 449-468 (1938). Gilman, H. Organic Chemistry. John Wiley and Sons Inc,. New York IV, p740 (1953). Joule, J. A. & Keith, M. Heterocyclic Chemistry, Blackwell Publishing Ltd, 2010, p347. Kankaanpera, A. & Salomaa, P. Kinetics of the hydrolytic cleavage of the furan ring. Acta Chem. Scand. 21, 575-577 (1967). Salomaa, P., Kankaanpera, A., Nikander, E., Karpainen, K. & Aaltonen, R. Electrophilic attack on the furan ring. Acid-catalyzed protodedeuteriation, protodetritiation, deuteriodeprotonation, and deuteriodetritiation at the α-position. Hydrolytic cleavage of the ring. Acta Chem. Scand. 27, 153-163, doi:10.3891/acta.chem.scand.27-0153 (1973). Kankaanpera, A. & Aaltonen, R. General acid catalysis in the hydrolysis of a furan derivative. Acta Chem. Scand. 26, 2537-2540, doi:10.3891/acta.chem.scand.26-2537 (1972). Garst, J. E. & Schmir, G. L. Hydrolysis of 2-methoxyfuran. J. Org. Chem. 39, 2920-2923, doi:10.1021/jo00933a021 (1974). Horvat, J., Klaic, B., Metelko, B. & Sunjic, V. Mechanism of levulinic acid formation. Tetrahedron Lett. 26, 2111-2114 (1985). Horvat, J., Klaic, B., Metelko, B. & Sunjic, V. Mechanism of levulinic acid formation in acid catalyzed hydrolysis of 2-(hydroxymethyl)furan and 5-(hydroxymethyl)-2furancarboxaldehyde. Croat. Chem. Acta 59, 429-438 (1986). Stamhuis, E. J., Drenth, W. & van den Berg, H. Mechanism of reactions of furans. I. A kinetic study of the acid-catalyzed hydrolysis of furan and 2,5-dimethylfuran. Recl. Trav. Chim. Pays-Bas 83, 167-176 (1964). Kankaanpera, A. & Kleemola, S. Equilibrium protonation of a carbon base prior to its hydrolysis. Acid-catalyzed cleavage of the furan ring. Acta Chem. Scand. 23, 3607-3608, doi:10.3891/acta.chem.scand.23-3607 (1969). Unverferth, K. & Schwetlick, K. Mechanism of furan hydrolysis. J. Prakt. Chem. 312, 882886 (1970). Houriet, R., Schwarz, H., Zummack, W., Andrade, J. G. & Schleyer, P. V. R. The α- vs βprotonation of pyrrole, furan, thiophene and cyclopentadiene. Gas phase proton and hydrogen affinities. The bishomocyclopropenyl cation. Nouv. J. Chim. 5, 505-509 (1981). Lorenz, U. J. et al. Protonation of heterocyclic aromatic molecules: IR signature of the protonation site of furan and pyrrole. Int. J. Mass spectrom. 267, 43-53, doi:10.1016/j.ijms.2007.02.017 (2007). Belen'kii, L. I. Positional selectivity in electrophilic substitution in π-excessive heteroaromatics. Adv. Heterocycl. Chem. 99, 143-183, doi:10.1016/s0065-2725(10)09905-8 (2010). Hiraoka, K., Takimoto, H. & Yamabe, S. Stabilities and structures in cluster ions of fivemembered heterocyclic compounds containing oxygen, nitrogen, and sulfur atoms. J. Am. Chem. Soc. 109, 7346-7352 (1987).

ACS Paragon Plus Environment

Page 18 of 29

Page 19 of 29 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

28

29

30

31

32 33 34

35

36

37

38 39 40 41

42 43 44

Gubina, T., Pankratov, A., Labunskaya, V., Voronin, S. & Kharchenko, V. Study of the mechanism of recyclization of furans into thiophenes and selenophenes in conditions of acid catalysis. 6. Experiments with labeled atoms. Quantum chemical calculations of intermediates of recyclization and hydrolysis. Chem. Heterocycl. Compd. 33, 903-909 (1997). Gubina, T., Labunskaya, V., Pankratov, A., Voronin, S. & Kharchenko, V. Study of the mechanism of recyclization of furans into thiophenes and selenophenes in conditions of acid catalysis. 5. Direction of protonation of furans. Chem. Heterocycl. Compd. 33, 898-902 (1997). van Beelen, E. S. E., Koblenz, T. A., Ingemann, S. & Hammerum, S. Experimental and Theoretical Evaluation of Proton Affinities of Furan, the Methylphenols, and the Related Anisoles. J. Phys. Chem. A 108, 2787-2793, doi:10.1021/jp0375721 (2004). Zeng, K. & Cao, Z.-X. Protonation of pyrrole and furan by H3O+ and NH4+ in the gas phase: a density functional theory study. Chin. J. Chem. 24, 293-298, doi:10.1002/cjoc.200690056 (2006). Nikbin, N., Caratzoulas, S. & Vlachos, D. G. On the Brønsted Acid‐Catalyzed Homogeneous Hydrolysis of Furans. ChemSusChem 6, 2066-2068 (2013). Park, J. M., Laio, A., Iannuzzi, M. & Parrinello, M. Dissociation Mechanism of Acetic Acid in Water. J. Am. Chem. Soc. 128, 11318-11319, doi:10.1021/ja060454h (2006). Cucinotta, C. S., Ruini, A., Catellani, A. & Stirling, A. Ab initio molecular dynamics study of the keto-enol tautomerism of acetone in solution. ChemPhysChem 7, 1229-1234, doi:10.1002/cphc.200600007 (2006). Schreiner, E., Nair, N. N. & Marx, D. Influence of Extreme Thermodynamic Conditions and Pyrite Surfaces on Peptide Synthesis in Aqueous Media. J. Am. Chem. Soc. 130, 2768-2770, doi:10.1021/ja7108085 (2008). Liang, X., Montoya, A. & Haynes, B. S. Local Site Selectivity and Conformational Structures in the Glycosidic Bond Scission of Cellobiose. J. Phys. Chem. B 115, 10682-10691, doi:10.1021/jp204199h (2011). Liang, X., Montoya, A. & Haynes, B. S. Molecular Dynamics Study of Acid-Catalyzed Hydrolysis of Dimethyl Ether in Aqueous Solution. J. Phys. Chem. B 115, 8199-8206, doi:10.1021/jp201951a (2011). Car, R. & Parrinello, M. Phys. ReV. Lett. 55, 2471 (1985). Goedecker, S., Teter, M. & Hutter, J. Separable dual-space Gaussian pseudopotentials. Phys. Rev. B Condens. Matter 54, 1703-1710 (1996). Becke, A. D. Density-functional exchange-energy approximation with correct asymptotic behavior. Phys. Rev. A Gen. Phys. 38, 3098-3100 (1988). Lee, C., Yang, W. & Parr, R. G. Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density. Phys. Rev. B Condens. Matter 37, 785-789 (1988). Nose, S. A unified formulation of the constant-temperature molecular-dynamics methods. J. Chem. Phys. 81, 511-519, doi:10.1063/1.447334 (1984). Coutinho, K. & Canuto, S. DICE: A Monte Carlo program for molecular liquid simulation. University of São Paulo, Brazil, version 2.9 (2003). Gaussian 09, R. A., M. J. Frisch, G. W. Trucks, H. B. Schlegel, G. E. Scuseria, M. A. Robb, J. R. Cheeseman, G. Scalmani, V. Barone, B. Mennucci, G. A. Petersson, H. Nakatsuji, M. Caricato, X. Li, H. P. Hratchian, A. F. Izmaylov, J. Bloino, G. Zheng, J. L. Sonnenberg, M. Hada, M. Ehara, K. Toyota, R. Fukuda, J. Hasegawa, M. Ishida, T. Nakajima, Y. Honda, O. Kitao, H. Nakai, T. Vreven, J. A. Montgomery, Jr., J. E. Peralta, F. Ogliaro, M. Bearpark, J. J. Heyd, E. Brothers, K. N. Kudin, V. N. Staroverov, R. Kobayashi, J. Normand, K. Raghavachari, A. Rendell, J. C. Burant, S. S. Iyengar, J. Tomasi, M. Cossi, N. Rega, J. M. Millam, M. Klene, J. E. Knox, J. B. Cross, V. Bakken, C. Adamo, J. Jaramillo, R. Gomperts, R. E. Stratmann, O. Yazyev, A. J. Austin, R. Cammi, C. Pomelli, J. W. Ochterski, R. L. Martin, K. Morokuma, V. G. Zakrzewski, G. A. Voth, P. Salvador, J. J. Dannenberg, S. Dapprich, A. D. Daniels, Ö. Farkas, J. B. Foresman, J. V. Ortiz, J. Cioslowski, and D. J. Fox, Gaussian, Inc., Wallingford CT, 2009.

ACS Paragon Plus Environment

Energy & Fuels 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

45

46 47

48

49

50 51 52 53

54

55 56

57

58

59

60

61

62 63 64 65

Iannuzzi, M., Laio, A. & Parrinello, M. Efficient Exploration of Reactive Potential Energy Surfaces Using Car-Parrinello Molecular Dynamics. Phys. Rev. Lett. 90, 238302/238301238302/238304, doi:10.1103/PhysRevLett.90.238302 (2003). Qian, X. Mechanisms and Energetics for Acid Catalyzed β-D-Glucose Conversion to 5Hydroxymethylfurfurl. J. Phys. Chem. A 115, 11740-11748, doi:10.1021/jp2041982 (2011). Liu, D., Nimlos, M. R., Johnson, D. K., Himmel, M. E. & Qian, X. Free Energy Landscape for Glucose Condensation Reactions. J. Phys. Chem. A. 114, 12936-12944, doi:10.1021/jp1078407 (2010). Laio, A. & Gervasio, F. L. Metadynamics: a method to stimulate rare events and reconstruct the free energy in biophysics, chemistry and material science. Rep. Prog. Phys. 71, 126601/126601-126601/126622, doi:10.1088/0034-4885/71/12/126601 (2008). Petersen, L., Ardevol, A., Rovira, C. & Reilly, P. J. Mechanism of Cellulose Hydrolysis by Inverting GH8 Endoglucanases: A QM/MM Metadynamics Study. J. Phys. Chem. B. 113, 7331-7339, doi:10.1021/jp811470d (2009). Fu, C.-W. & Lin, T.-H. Theoretical Study on the Alkaline Hydrolysis of Methyl Thioacetate in Aqueous Solution. J. Phys. Chem. A 115, 13523-13533, doi:10.1021/jp204658w (2011). Barker, I. J., Petersen, L. & Reilly, P. J. Mechanism of Xylobiose Hydrolysis by GH43 βXylosidase. J. Phys. Chem. B. 114, 15389-15393, doi:10.1021/jp107886e (2010). Alfonso-Prieto, M., Biarnes, X., Vidossich, P. & Rovira, C. The Molecular Mechanism of the Catalase Reaction. J. Am. Chem. Soc. 131, 11751-11761, doi:10.1021/ja9018572 (2009). Comas-Vives, A., Stirling, A., Lledos, A. & Ujaque, G. The Wacker Process: Inner- or OuterSphere Nucleophilic Addition? New Insights from Ab Initio Molecular Dynamics. Chemistry-A European Journal 16, 8738-8747, S8738/8731-S8738/8712, doi:10.1002/chem.200903522 (2010). Bussi, G., Laio, A. & Parrinello, M. Equilibrium Free Energies from Nonequilibrium Metadynamics. Phys. Rev. Lett. 96, 090601/090601-090601/090604, doi:10.1103/PhysRevLett.96.090601 (2006). Stirling, A. & Papai, I. H2CO3 Forms via HCO3- in Water. J. Phys. Chem. B 114, 1685416859, doi:10.1021/jp1099909 (2010). Blumberger, J., Ensing, B. & Klein, M. L. Formamide hydrolysis in alkaline aqueous solution: insight from Ab initio metadynamics calculations. Angewandte Chemie, International Edition 45, 2893-2897, doi:10.1002/anie.200600283 (2006). Geissler, P. L., Van Voorhis, T. & Dellago, C. Potential energy landscape for proton transfer in (H2O)3H+: comparison of density functional theory and wave function-based methods. Chem. Phys. Lett. 324, 149-155, doi:10.1016/s0009-2614(00)00479-6 (2000). Pan, Y. & McAllister, M. A. Characterization of low-barrier hydrogen bonds. 4. Basis set and correlation effects: an ab initio and DFT investigation. J. Mol. Struct. THEOCHEM 427, 221227, doi:10.1016/s0166-1280(97)00227-3 (1998). Xu, X. & Goddard, W. A., III. Bonding Properties of the Water Dimer: A Comparative Study of Density Functional Theories. J. Phys. Chem. A 108, 2305-2313, doi:10.1021/jp035869t (2004). Dong, H., Nimlos, M. R., Himmel, M. E., Johnson, D. K. & Qian, X. The Effects of Water on β-D-Xylose Condensation Reactions. J. Phys. Chem. A. 113, 8577-8585, doi:10.1021/jp9025442 (2009). Isayev, O., Gorb, L. & Leszczynski, J. Theoretical calculations: can Gibbs free energy for intermolecular complexes be predicted efficiently and accurately? J. Comput. Chem. 28, 1598-1609, doi:10.1002/jcc.20696 (2007). Humphrey, W., Dalke, A. & Schulten, K. VDM: visual molecular dynamics. J. Mol. Graphics 14, 33-38, plates, 27-28 (1996). Lerbret, A. et al. Molecular dynamics studies of the conformation of sorbitol. Carbohydr. Res. 344, 2229-2235, doi:10.1016/j.carres.2009.08.003 (2009). Morando, M. A. et al. NMR and molecular modeling reveal key structural features of synthetic nodulation factors. Glycobiology 21, 824-833, doi:10.1093/glycob/cwr014 (2011). McKillip, W. J., Collin, G., Höke, H. & Zeitsch, K. J. Furan and derivatives. Ullmann's Encyclopedia of Industrial Chemistry (1981).

ACS Paragon Plus Environment

Page 20 of 29

Page 21 of 29 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

66

67

68 69

70 71

Liescheski, P. B. & Rankin, D. W. H. Molecular structure of furan, determined by combined analyses of data obtained by electron diffraction, rotational spectroscopy, and liquid crystal NMR spectroscopy. J. Mol. Struct. 196, 1-19, doi:10.1016/0022-2860(89)85001-x (1989). Ichikawa, K., Kameda, Y., Yamaguchi, T., Wakita, H. & Misawa, M. Neutron-diffraction investigation of the intramolecular structure of a water molecule in the liquid phase at high temperatures. Mol. Phys. 73, 79-86, doi:10.1080/00268979100101071 (1991). Heimlich, K. R. & Martin, A. N. A kinetic study of glucose degradation in acid solution. Journal of the American Pharmaceutical Association, Scientific Edition 49, 592-597 (1960). McKIBBINS, S. W., HARRIS, J. F., SAEMAN, J. F. & NEILL, W. K. Kinetics of the acid catalyzed conversion of glucose to 5-hydroxymethyl 2-furaldehyde and levulinic acid. FOREST PRODUCTS JOURNAL (1962). Girisuta, B., Janssen, L. P. B. M. & Heeres, H. J. A kinetic study on the decomposition of 5hydroxymethylfurfural into levulinic acid. Green Chem. 8, 701-709 (2006). Marenich, A. V., Cramer, C. J. & Truhlar, D. G. Universal solvation model based on solute electron density and on a continuum model of the solvent defined by the bulk dielectric constant and atomic surface tensions. The Journal of Physical Chemistry B 113, 6378-6396 (2009).

ACS Paragon Plus Environment

Energy & Fuels 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Scheme 1. Possible reaction pathways of acid-catalysed ring opening of furan

ACS Paragon Plus Environment

Page 22 of 29

Page 23 of 29 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

 

2,5 dihydro-2-furanol

2,3 dihydro-3-furanol

2,3 dihydro-2-furanol

Figure 1. Relevant nomenclatures of atoms and definitions of CV1, CV2, CV3 and CV4 in the CPMD metadynamics computations of acid-catalysed ring opening of furan. Bonds shown in red are used to define the collective variables for the acid-catalysed hydration of furan, while those in blue colour are for the collective variables used to characterise ring opening of the dihydrofuranols.

ACS Paragon Plus Environment

Energy & Fuels 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 2. Free energy landscape of acid-catalysed hydration of furan to 2,5 dihydro-2furanol in aqueous phase at 300 K. The energy scale (in kcal mol–1) is given in the inset to the left. Insets R, TS, and PI (with relative free energy in parentheses) in the energy landscape correspond to the reactant, transition state and product basins, respectively. The average lengths for several important bonds of each state and their corresponding standard deviations are presented in Å. The redundant solvent molecules of every molecular configuration are deleted to emphasise the atomic configuration of the furanic ring.

ACS Paragon Plus Environment

Page 24 of 29

Page 25 of 29 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

Figure 3. Free energy landscape of acid-catalysed hydration of furan to 2,3 dihydro-3furanol in aqueous phase at 300 K.

ACS Paragon Plus Environment

Energy & Fuels 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 4. Free energy landscape of acid-catalysed hydration of furan to 2,3 dihydro-2furanol in aqueous phase at 300 K.

ACS Paragon Plus Environment

Page 26 of 29

Page 27 of 29 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

Figure 5. Free energy landscape of acid-catalysed ring opening of 2,5 dihydro-2furanol in aqueous phase at 300 K.

ACS Paragon Plus Environment

Energy & Fuels 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 6. Free energy landscape of acid-catalysed ring opening of 2,3 dihydro-3furanol to produce 3,4-dihydroxy-butenol in aqueous phase at 300 K.

ACS Paragon Plus Environment

Page 28 of 29

Page 29 of 29 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

The following graphic will be used for the TOC on the website

ACS Paragon Plus Environment