10.1002/cssc.201700394
ChemSusChem
FULL PAPER
transition state and calculate reaction barriers. The transition states have
been further verified by vibrational calculations. An HSE (Heyd-Scuseria-
Ernzerhof) hybrid functional[32] was used to calculate the total energy of
the initial, transition, and final states that had already been optimized by
DFT-D3 calculations to reduce the underestimation of the reaction
barriers caused by the charge delocalization error of the local and semi-
local exchange-correlation functionals.[22]
upgraded to longer chain molecules via different chemistries
such as aldol condensation due to the ketone functionality.[1a, 25]
Furthermore, in this reaction, once the acylium ions are
generated, the aromatic compounds with oxygen-containing
functionalities, such as hydroxyl –OH and alkoxy – OCH3, are all
active to pick them up with the activity of this C-C coupling step
decreasing in the order of m-cresol > phenol > anisole. Alkyl
aromatic compounds such as toluene are not as effective for
picking up the acylium ions, which leads to a low acylation yield
and high loss of acylium ion to coke. The presence of oxygen-
contained functional groups in the lignin derived aromatic
compounds was not considered favorable due to the
requirement for intensive hydrotreating step. However, as shown
in this study, these functional groups have the advantage of
activating the aromatic ring and make them easier in the
upgrading process, particularly for C-C bond formation reactions.
It is thus very important to perform C-C coupling reaction via
acylation with the available functionalities in bio-oil before the
hydrotreating process.
Acknowledgements
This work was supported by the U.S. Department of Energy,
DOE/EPSCOR (Grant DESC0004600). The computational
research used the supercomputer resources of the National
Energy Research Scientific Computing Centre (NERSC) and the
OU Supercomputing Centre for Education & Research (OSCER)
at the University of Oklahoma.
In the last part of the study, the Fries rearrangement of
phenyl acetate to HAPs was investigated to differentiate
between the intermolecular and intramolecular reaction pathway.
The reaction was found to happen via both mechanisms, in
which the intermolecular one is more dominant.
Keywords: acylation • acetic acid • aromatic ester • alkyl ester •
aromatic substrate
[1]
aT. N. Pham, D. C. Shi, D. E. Resasco, Appl Catal B-Environ
2014, 145, 10-23; bM. A. Gonzalez-Borja, D. E. Resasco, Aiche J
2015, 61, 598-609.
[2]
[3]
J. A. Herron, T. Vann, N. Duong, D. E. Resasco, S. Crossley, L. L.
Lobban, C. T. Maravelias, Energy Technology 2017, 5, 130-150.
aT. N. Pham, D. C. Shi, T. Sooknoi, D. E. Resasco, J Catal 2012,
295, 169-178; bT. N. Pham, D. C. Shi, D. E. Resasco, J Catal
2014, 314, 149-158.
P. A. Zapata, J. Faria, M. P. Ruiz, D. E. Resasco, Top Catal 2012,
55, 38-52.
Experimental Section
Three H-form zeolites were used as catalysts. Two of them (Beta
CP814C, and ZSM-5 CBV-5524G) were received in the NH4-form and
then converted to the H-form by calcining under air at 600oC for 5 hours,
after a linear heating ramp of 2oC/min. The third sample was an HY
zeolite (CBV 760) received in the H-form from the manufacturer (Zeolyst
International).
[4]
[5]
[6]
P. A. Zapata, J. Faria, M. P. Ruiz, R. E. Jentoft, D. E. Resasco, J
Am Chem Soc 2012, 134, 8570-8578.
aG. Sartori, R. Maggi, Chem Rev 2006, 106, 1077-1104; bG.
Sartori, R. Maggi, Chem Rev 2011, 111, Pr181-Pr214.
A. Gumidyala, T. Sooknoi, S. Crossley, J Catal 2016, 340, 76-84.
A. Gumidyala, B. Wang, S. Crossley, Science Advances 2016, 2.
D. S. Park, K. E. Joseph, M. Koehle, C. Krumm, L. M. Ren, J. N.
Damen, M. H. Shete, H. S. Lee, X. B. Zuo, B. Lee, W. Fan, D. G.
Vlachos, R. F. Lobo, M. Tsapatsis, P. J. Dauenhauer, Acs Central
Science 2016, 2, 820-824.
aD. Rohan, C. Canaff, E. Fromentin, M. Guisnet, J Catal 1998,
177, 296-305; bD. Rohan, C. Canaff, P. Magnoux, M. Guisnet,
Journal of Molecular Catalysis a-Chemical 1998, 129, 69-78; cE.
G. Derouane, G. Crehan, C. J. Dillon, D. Bethell, H. He, S. B.
[7]
[8]
[9]
The reactions were carried out in the liquid phase using a 160-mL
batch stainless steel autoclave reactor (Parr Corporation). The reactor is
equipped with a stirring impeller, temperature controller, pressure gauge,
and sampling port. The reactions were tested over a wide range of
temperatures, under N2. The initial pressure of N2 was kept constant at
500 psi, the total volume of reactants was 80 ml and the stirring speed
500 rpm. In a typical reaction run, a measured amount of zeolite was
mixed with the reactant mixture in the reactor vessel. After the reactor
was sealed and pressurized with N2, the temperature was increased to
the desired value. After reaction, the reactor was rapidly cooled down to
room temperature and the liquid product was filtered and analyzed on a
GC-MS (Shimadzu) for identification and GC-FID (Agilent HP 5890) for
quantification.
[10]
Derouane-Abd Hamid,
J Catal 2000, 194, 410-423; dE. G.
Derouane, C. J. Dillon, D. Bethell, S. B. Derouane-Abd Hamid, J
Catal 1999, 187, 209-218.
[11]
aY. D. Ma, Q. L. Wang, W. Jiang, B. J. Zuo, Applied Catalysis a-
General 1997, 165, 199-206; bO. Kresnawahjuesa, R. J. Gorte, D.
White, Journal of Molecular Catalysis A: Chemical 2004, 208, 175-
185; cK. Gaare, D. Akporiaye, Journal of Molecular Catalysis A:
Chemical 1996, 109, 177-187.
[12]
[13]
[14]
[15]
aA. P. Singh, A. K. Pandey, Journal of Molecular Catalysis a-
Chemical 1997, 123, 141-147; bC. P. Bezouhanova, Applied
Catalysis a-General 2002, 229, 127-133.
V. D. Kumari, G. Saroja, A. Ratnamala, M. Noorjahan, M.
Subrahmanyam, Reaction Kinetics and Catalysis Letters 2003, 79,
43-51.
aA. Vogt, H. W. Kouwenhoven, R. Prins, Applied Catalysis a-
General 1995, 123, 37-49; bY. Pouilloux, N. S. Gnep, P. Magnoux,
G. Perot, J Mol Catal 1987, 40, 231-233.
aV. D. Chaube, P. Moreau, A. Finiels, A. V. Ramaswamy, A. P.
Singh, Catal Lett 2002, 79, 89-94; bV. D. Chaube, P. Moreau, A.
Finiels, A. V. Ramaswamy, A. P. Singh, Journal of Molecular
Catalysis a-Chemical 2001, 174, 255-264.
V. Durga Kumari, G. Saroja, A. Ratnamala, M. Noorjahan, M.
Subrahmanyam, Reaction Kinetics and Catalysis Letters 2003, 79,
43-51.
Q. L. Wang, Y. D. Ma, X. D. Ji, H. Yan, Q. Qiu, Journal of the
Chemical Society-Chemical Communications 1995, 2307-2308.
G. M. Loudon, Organic Chemistry, Roberts and Company, 2009.
I. Neves, F. Jayat, P. Magnoux, G. Perot, F. R. Ribeiro, M.
Gubelmann, M. Guisnet, J Mol Catal 1994, 93, 169-179.
DFT calculations
Density functional theory calculations were performed using the VASP
package[26] to investigate the reaction mechanism. The PBE (Perdew-
Burke-Ernzerhof) exchange-correlation potential[27] was used, and the
electron-core interactions were treated in the projector augmented wave
(PAW) method.[28] The van der Waals interactions have been taken into
account through the so-called DFT-D3 semi-empirical method via a pair-
wise force field.[29] The calculations have been performed using an H-
beta unit cell that includes 64 Si and 128 O atoms. One Si atom was
replaced by one Al atom, so the Si/Al ratio was 63/1. The structure of the
unit cell was taken from an experimental work (a = 17.829 Å, b = 17.829
Å, c = 14.671 Å)[30] and fixed during the calculation. All the atoms were
relaxed until the atomic forces were smaller than 0.02 eV Å-1, using a
single Γ point of the Brillouin zone with a kinetic cut off energy of 400 eV.
The nudged elastic band (NEB)[31] method has been used to find the
[16]
[17]
[18]
[19]
This article is protected by copyright. All rights reserved.