Ab initio dynamics of cellulose pyrolysis: Nascent decomposition pathways at 327 and 600 °c

Vishal Agarwal, Paul J. Dauenhauer, George W. Huber, Scott M. Auerbach

Research output: Contribution to journalArticle

88 Citations (Scopus)

Abstract

We modeled nascent decomposition processes in cellulose pyrolysis at 327 and 600 °C using Car-Parrinello molecular dynamics (CPMD) simulations with rare events accelerated with the metadynamics method. We used a simulation cell comprised of two unit cells of cellulose Iβ periodically repeated in three dimensions to mimic the solid cellulose. To obtain initial conditions at reasonable densities, we extracted coordinates from larger classical NPT simulations at the target temperatures. CPMD-metadynamics implemented with various sets of collective variables, such as coordination numbers of the glycosidic oxygen, yielded a variety of chemical reactions such as depolymerization, fragmentation, ring opening, and ring contraction. These reactions yielded precursors to levoglucosan (LGA)-the major product of pyrolysis-and also to minor products such as 5-hydroxy-methylfurfural (HMF) and formic acid. At 327 °C, we found that depolymerization via ring contraction of the glucopyranose ring to the glucofuranose ring occurs with the lowest free-energy barrier (20 kcal/mol). We suggest that this process is key for formation of liquid intermediate cellulose, observed experimentally above 260 °C. At 600 °C, we found that a precursor to LGA (pre-LGA) forms with a free-energy barrier of 36 kcal/mol via an intermediate/transition state stabilized by anchimeric assistance and hydrogen bonding. Conformational freedom provided by expansion of the cellulose matrix at 600 °C was found to be crucial for formation of pre-LGA. We performed several comparison calculations to gauge the accuracy of CPMD-metadynamics barriers with respect to basis set and level of theory. We found that free-energy barriers at 600 °C are in the order pre-LGA < pre-HMF < formic acid, explaining why LGA is the kinetically favored product of fast cellulose pyrolysis.

Original languageEnglish (US)
Pages (from-to)14958-14972
Number of pages15
JournalJournal of the American Chemical Society
Volume134
Issue number36
DOIs
StatePublished - Sep 12 2012

Fingerprint

Cellulose
formic acid
Pyrolysis
Decomposition
Energy barriers
Molecular Dynamics Simulation
Free energy
Molecular dynamics
Depolymerization
Railroad cars
Formic acid
Hydrogen Bonding
Gages
Chemical reactions
Hydrogen bonds
Oxygen
Temperature
Computer simulation
Liquids
1,6-anhydro-beta-glucopyranose

Cite this

Ab initio dynamics of cellulose pyrolysis : Nascent decomposition pathways at 327 and 600 °c. / Agarwal, Vishal; Dauenhauer, Paul J.; Huber, George W.; Auerbach, Scott M.

In: Journal of the American Chemical Society, Vol. 134, No. 36, 12.09.2012, p. 14958-14972.

Research output: Contribution to journalArticle

Agarwal, Vishal ; Dauenhauer, Paul J. ; Huber, George W. ; Auerbach, Scott M. / Ab initio dynamics of cellulose pyrolysis : Nascent decomposition pathways at 327 and 600 °c. In: Journal of the American Chemical Society. 2012 ; Vol. 134, No. 36. pp. 14958-14972.
@article{23b5a2406cd948f6b21a5bcbfdfd1cc5,
title = "Ab initio dynamics of cellulose pyrolysis: Nascent decomposition pathways at 327 and 600 °c",
abstract = "We modeled nascent decomposition processes in cellulose pyrolysis at 327 and 600 °C using Car-Parrinello molecular dynamics (CPMD) simulations with rare events accelerated with the metadynamics method. We used a simulation cell comprised of two unit cells of cellulose Iβ periodically repeated in three dimensions to mimic the solid cellulose. To obtain initial conditions at reasonable densities, we extracted coordinates from larger classical NPT simulations at the target temperatures. CPMD-metadynamics implemented with various sets of collective variables, such as coordination numbers of the glycosidic oxygen, yielded a variety of chemical reactions such as depolymerization, fragmentation, ring opening, and ring contraction. These reactions yielded precursors to levoglucosan (LGA)-the major product of pyrolysis-and also to minor products such as 5-hydroxy-methylfurfural (HMF) and formic acid. At 327 °C, we found that depolymerization via ring contraction of the glucopyranose ring to the glucofuranose ring occurs with the lowest free-energy barrier (20 kcal/mol). We suggest that this process is key for formation of liquid intermediate cellulose, observed experimentally above 260 °C. At 600 °C, we found that a precursor to LGA (pre-LGA) forms with a free-energy barrier of 36 kcal/mol via an intermediate/transition state stabilized by anchimeric assistance and hydrogen bonding. Conformational freedom provided by expansion of the cellulose matrix at 600 °C was found to be crucial for formation of pre-LGA. We performed several comparison calculations to gauge the accuracy of CPMD-metadynamics barriers with respect to basis set and level of theory. We found that free-energy barriers at 600 °C are in the order pre-LGA < pre-HMF < formic acid, explaining why LGA is the kinetically favored product of fast cellulose pyrolysis.",
author = "Vishal Agarwal and Dauenhauer, {Paul J.} and Huber, {George W.} and Auerbach, {Scott M.}",
year = "2012",
month = "9",
day = "12",
doi = "10.1021/ja305135u",
language = "English (US)",
volume = "134",
pages = "14958--14972",
journal = "Journal of the American Chemical Society",
issn = "0002-7863",
publisher = "American Chemical Society",
number = "36",

}

TY - JOUR

T1 - Ab initio dynamics of cellulose pyrolysis

T2 - Nascent decomposition pathways at 327 and 600 °c

AU - Agarwal, Vishal

AU - Dauenhauer, Paul J.

AU - Huber, George W.

AU - Auerbach, Scott M.

PY - 2012/9/12

Y1 - 2012/9/12

N2 - We modeled nascent decomposition processes in cellulose pyrolysis at 327 and 600 °C using Car-Parrinello molecular dynamics (CPMD) simulations with rare events accelerated with the metadynamics method. We used a simulation cell comprised of two unit cells of cellulose Iβ periodically repeated in three dimensions to mimic the solid cellulose. To obtain initial conditions at reasonable densities, we extracted coordinates from larger classical NPT simulations at the target temperatures. CPMD-metadynamics implemented with various sets of collective variables, such as coordination numbers of the glycosidic oxygen, yielded a variety of chemical reactions such as depolymerization, fragmentation, ring opening, and ring contraction. These reactions yielded precursors to levoglucosan (LGA)-the major product of pyrolysis-and also to minor products such as 5-hydroxy-methylfurfural (HMF) and formic acid. At 327 °C, we found that depolymerization via ring contraction of the glucopyranose ring to the glucofuranose ring occurs with the lowest free-energy barrier (20 kcal/mol). We suggest that this process is key for formation of liquid intermediate cellulose, observed experimentally above 260 °C. At 600 °C, we found that a precursor to LGA (pre-LGA) forms with a free-energy barrier of 36 kcal/mol via an intermediate/transition state stabilized by anchimeric assistance and hydrogen bonding. Conformational freedom provided by expansion of the cellulose matrix at 600 °C was found to be crucial for formation of pre-LGA. We performed several comparison calculations to gauge the accuracy of CPMD-metadynamics barriers with respect to basis set and level of theory. We found that free-energy barriers at 600 °C are in the order pre-LGA < pre-HMF < formic acid, explaining why LGA is the kinetically favored product of fast cellulose pyrolysis.

AB - We modeled nascent decomposition processes in cellulose pyrolysis at 327 and 600 °C using Car-Parrinello molecular dynamics (CPMD) simulations with rare events accelerated with the metadynamics method. We used a simulation cell comprised of two unit cells of cellulose Iβ periodically repeated in three dimensions to mimic the solid cellulose. To obtain initial conditions at reasonable densities, we extracted coordinates from larger classical NPT simulations at the target temperatures. CPMD-metadynamics implemented with various sets of collective variables, such as coordination numbers of the glycosidic oxygen, yielded a variety of chemical reactions such as depolymerization, fragmentation, ring opening, and ring contraction. These reactions yielded precursors to levoglucosan (LGA)-the major product of pyrolysis-and also to minor products such as 5-hydroxy-methylfurfural (HMF) and formic acid. At 327 °C, we found that depolymerization via ring contraction of the glucopyranose ring to the glucofuranose ring occurs with the lowest free-energy barrier (20 kcal/mol). We suggest that this process is key for formation of liquid intermediate cellulose, observed experimentally above 260 °C. At 600 °C, we found that a precursor to LGA (pre-LGA) forms with a free-energy barrier of 36 kcal/mol via an intermediate/transition state stabilized by anchimeric assistance and hydrogen bonding. Conformational freedom provided by expansion of the cellulose matrix at 600 °C was found to be crucial for formation of pre-LGA. We performed several comparison calculations to gauge the accuracy of CPMD-metadynamics barriers with respect to basis set and level of theory. We found that free-energy barriers at 600 °C are in the order pre-LGA < pre-HMF < formic acid, explaining why LGA is the kinetically favored product of fast cellulose pyrolysis.

UR - http://www.scopus.com/inward/record.url?scp=84866414245&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=84866414245&partnerID=8YFLogxK

U2 - 10.1021/ja305135u

DO - 10.1021/ja305135u

M3 - Article

C2 - 22889121

AN - SCOPUS:84866414245

VL - 134

SP - 14958

EP - 14972

JO - Journal of the American Chemical Society

JF - Journal of the American Chemical Society

SN - 0002-7863

IS - 36

ER -