INTEGRATED PROCESS FOR γ-butyrolactone PRODUCTION

Similar documents
Design and Control of an Integrated Toluene-Aniline Production Plant a Preliminary Study

Optimization of Batch Distillation Involving Hydrolysis System

Lecture 25: Manufacture of Maleic Anhydride and DDT

Process Design Decisions and Project Economics Prof. Dr. V. S. Moholkar Department of Chemical Engineering Indian Institute of Technology, Guwahati

HW Help. How do you want to run the separation? Safety Issues? Ease of Processing

regressing the vapor-liquid equilibrium data in Mathuni et al. and Rodriguez et al., respectively. The phase equilibrium data of the other missing pai

A Definition for Plantwide Controllability. Process Flexibility

Simulation and Process Integration of Clean Acetone Plant

(1) This reaction mechanism includes several undesired side reactions that produce toluene and benzene:

Reactors. Reaction Classifications

DME(10 TPD) Process Simulation Using Aspen Plus Release Dr. Jungho Cho, Professor Department of Chemical Engineering Dong Yang University

DESIGN AND CONTROL OF BUTYL ACRYLATE REACTIVE DISTILLATION COLUMN SYSTEM. I-Lung Chien and Kai-Luen Zeng

Process Design Decisions and Project Economics Prof. Dr. V. S. Moholkar Department of Chemical Engineering Indian Institute of Technology, Guwahati

TOPIC: Conceptual Flowsheet for Production of Benzene from Toluene. Proposed Solution:

Process design decisions and project economics Dr. V. S. Moholkar Department of chemical engineering Indian Institute of Technology, Guwahati

Integrated Knowledge Based System for Process Synthesis

Exercise 1. Material balance HDA plant

Simulation of Butyl Acetate and Methanol Production by Transesterification Reaction via Conventional Distillation Process

Systems Engineering Spring Group Project #1: Process Flowsheeting Calculations for Acetic Anhydride Plant. Date: 2/25/00 Due: 3/3/00

INTRODUCTION TO CHEMICAL PROCESS SIMULATORS

Level 2: Input output structure

Ammonia Synthesis with Aspen Plus V8.0

Shortcut Distillation. Agung Ari Wibowo, S.T., M.Sc Politeknik Negeri Malang Malang - Indonesia

Novel Control Structures for Heterogeneous Reactive Distillation

CHEMICAL REACTORS - PROBLEMS OF REACTOR ASSOCIATION 47-60

Modeling, Simulation and Control of a Tubular Fixed-bed Dimethyl Ether Reactor

System Identification and Analysis to Optimize Plant Design of Aniline

Lecture (9) Reactor Sizing. Figure (1). Information needed to predict what a reactor can do.

"Thermodynamic Analysis of Processes for Hydrogen Generation by Decomposition of Water"

Exercise 1. Material balance HDA plant

Effect of the Reaction Bath Temperature in a Fixed-bed Reactor for Oxidation of o-xylene over V 2 O 5 /TiO 2 Catalysts

To increase the concentration of product formed in a PFR, what should we do?

Innovative Design of Diphenyl Carbonate Process via One Reactive Distillation with a Feed-Splitting Arrangement

CHE 404 Chemical Reaction Engineering. Chapter 8 Steady-State Nonisothermal Reactor Design

Chemical Reaction Engineering Prof. Jayant Modak Department of Chemical Engineering Indian Institute of Science, Bangalore

CHEMICAL THERMODYNAMICS. Nature of Energy. ΔE = q + w. w = PΔV

Implementation of Cleaner Production Principles in Formaldehyde Production

Aspen Plus PFR Reactors Tutorial using Styrene with Pressure Drop Considerations By Robert P. Hesketh and Concetta LaMarca Spring 2005

Analysis of processing systems involving reaction and distillation: the synthesis of ethyl acetate

Chemical Reactions and Kinetics of the Carbon Monoxide Coupling in the Presence of Hydrogen

Chemical Kinetics and Reaction Engineering

Steady State Multiplicity and Stability in a Reactive Flash

Chapter 4. Fundamentals of Material Balance

Modeling of the Unsteady State Methanol Synthesis at the Level of Catalyst Pellet

Simulation of Methanol Production Process and Determination of Optimum Conditions

Name Energy Test period Date

Lecture 8. Mole balance: calculations of microreactors, membrane reactors and unsteady state in tank reactors

SKF4153 PLANT DESIGN PROCESS SYNTHESIS & CREATION. Prof. Dr. Zainuddin Abdul Manan Ir. Dr. Sharifah Rafidah Wan Alwi

Technical Resource Package 1

Stoichiometric Reactor Simulation Robert P. Hesketh and Concetta LaMarca Chemical Engineering, Rowan University (Revised 4/8/09)

Pre GATE Pre-GATE 2018

Optimal Operation of Batch Reactive Distillation Process Involving Esterification Reaction System

METHYL ACETATE REACTIVE DISTILLATION PROCESS MODELING, SIMULATION AND OPTIMIZATION USING ASPEN PLUS

Ch. 6 Enthalpy Changes

Dividing wall columns for heterogeneous azeotropic distillation

Ways to Improve the Energy Efficiency of AN Solution Plants

ChE 344 Winter 2013 Mid Term Exam II Tuesday, April 9, 2013

A novel design of reactive distillation configuration for 2-methoxy-2-methylheptane process

Multiple Choice Identify the letter of the choice that best completes the statement or answers the question.

Reactor Design within Excel Enabled by Rigorous Physical Properties and an Advanced Numerical Computation Package

Study Of Molar Flow, Mole Fraction, Partial Pressure Variation In Process Of Pure Hydrogen Production In Multichannel Membrane Reactor

Chapter 12 Intermolecular Forces of Attraction

IJSRD - International Journal for Scientific Research & Development Vol. 1, Issue 8, 2013 ISSN (online):

Chemistry Chapter 16. Reaction Energy

Name Chemistry / / SOL Questions Chapter 9 For each of the following, fill in the correct answer on the BLUE side of the scantron.


Chemistry Heat Review. Heat: Temperature: Enthalpy: Calorimetry: Activation energy:

Control Study of Ethyl tert-butyl Ether Reactive Distillation

THERMAL INTEGRATION OF A DISTILLATION COLUMN THROUGH SIDE-EXCHANGERS

Fundamentals of Material Balances

Dynamics of forced and unsteady-state processes

Steam Reforming in a Packed Bed Membrane Reactor

Industrial Applications of Microreactor Technology

1. A. Define the term rate of reaction. The measure of the amount of reactants being converted into products per unit amount of time

Acetone Process Energy Recovery by Means of Energy Analysis

Integrated Catalyst System for Removing Buildup-Gas in BWR Inert Containments During a Severe Accident

A Review Paper on Design and simulation of tubular heat exchanging reactor for coupling exothermic and endothermic reactions

CHE 404 Chemical Reaction Engineering. Chapter 8 Steady-State Nonisothermal Reactor Design

Unit 7 Kinetics and Thermodynamics

CATALYTIC COMBUSTION OF HYD RO GEN/ AIR IN MICROCHANNEL REACTOR

PROCEEDINGS of the 4th International Conference on Chemical Technology 4th International Conference on Chemical Technology

10.02 PE Diagrams. 1. Given the equation and potential energy diagram representing a reaction:

Synthesis of Ethylene Glycol Evaluation of Standard Process and Possible Alternatives

ChE 201 August 26, ChE 201. Chapter 8 Balances on Nonreactive Processes Heat of solution and mixing

Zhongwei Fu, Yunyun Zhong, Yuehong Yu, Lizhen Long, Min Xiao, Dongmei Han, Shuanjin Wang *,

2. Review on Material Balances

MODELING AND SIMULATION OF AN AUTOTHERMAL REFORMER

Studies on Mo/HZSM-5 Complex catalyst for Methane Aromatization

Metals Production in The Future

How sulphur really forms on the catalyst surface

G T = G f (product) G f (reactant) = kj/mol Ka (550 O C) = exp(8840/(8.314*( )) = 3.369

SIMULATION ANALYSIS OF FULLY THERMALLY COUPLED DISTILLATION COLUMN

Process Intensification for Ethyl Lactate Production Using Reactive Distillation

PROCEEDINGS of the 5 th International Conference on Chemical Technology 5 th International Conference on Chemical Technology

COPYRIGHTED MATERIAL INTRODUCTION CHAPTER 1

A comparative study on the recovery of 1,2-dichloroethane and the removal of benzene contained in the byproducts of VCM process

FDE 211-MATERIAL AND ENERGY BALANCES: MATERIAL BALANCES ON REACTIVE SYSTEMS. Dr. Ilgın PakerYıkıcı Fall 2015

2 nd Semester Study Guide 2017

Third Quarter Cumulative Review Questions. 1. Which factor distinguishes a metallic bond from an ionic bond or a covalent bond?

Thermodynamics revisited

Transcription:

U.P.B. Sci. Bull., Series B, Vol. 76, Iss. 3, 214 ISSN 1454 2331 INTEGRATED PROCESS FOR γ-butyrolactone PRODUCTION Ahtesham JAVAID 1, Costin Sorin BILDEA 2 An integrated process for the production of γ-butyrolactone by coupling maleic anhydride hydrogenation and 1, 4-butanediol dehydrogenation in a fixed bed adiabatic reactor is studied. The reaction conditions realizing optimal hydrogen utilization and better energy efficiency are determined. Compared to stand-alone processes, the integrated process has several advantages, such as simpler reactor, easy temperature control, improved γ-butyrolactone yield, good energy efficiency and optimal hydrogen utilization. The integrated plant shows robust behavior with respect to disturbances and allows production rate change. Keywords: γ-butyrolactone, reaction coupling, process design and control. 1. Introduction γ-butyrolactone (GBL) is known as a versatile solvent and a raw material for the synthesis of pyrrolidone, N-methyl pyrrolidone, N-vinyl pyrrolidone, herbicides, and rubber additives. GBL is manufactured via two methods: hydrogenation of maleic anhydride (MA) and dehydrogenative cyclization of 1,4- butanediol (BDO). Both processes are typically performed in multi-tubular fixedbed reactors [1]. The catalytic hydrogenation of MA or its derivatives such as maleic acid, succinic anhydride or succinic maleate to GBL has frequently been reported in the literature [2]. The reaction equation of MA hydrogenation to GBL can be represented as follows: CHO+ 3 H CHO + HO H = -211 kj/mol (1) 4 2 3 2 4 6 2 2 It is obvious from equation (1) that production of 1 mole of GBL requires 1 mole of MA and 3 moles of H 2, releasing 211 kj of heat. Owing to the strong exothermic nature of the reaction, it is difficult to control the temperature of the process. In a tubular fixed-bed reactor, this might cause thermal runaway, lead to apparent hotspots and low selectivity to the desired product, GBL. A large excess of hydrogen is needed. 1 PhD student, Faculty of Applied Chemistry and Materials Science, University POLITEHNICA of Bucharest, e-mail: a_javaid@chim.upb.ro 2 Prof., Faculty of Applied Chemistry and Materials Science, University POLITEHNICA of Bucharest, e-mail: s_bildea@upb.ro

34 Ahtesham Javaid, Costin Sorin Bildea The vapor phase catalytic dehydrogenation of BDO to GBL has also been described in the literature [3]. This reaction can be expressed as follows: CH O CHO + 2H H = +61.6 kj/mol (2) 4 1 2 4 6 2 2 From equation (2) it is evident that production of 1 mole of GBL requires 1 mole of BDO, releasing 2 moles of H 2 and requiring 61.6 kj of heat. The hydrogen released from reaction cannot be used properly in a single dehydrogenation process. The new catalytic process proposed in this work combines reactions (1) and (2) into one coupled system which can significantly improve the yield of the catalytic hydrogenation of MA, apart from the better thermal balance and the effective usage of hydrogen through hydrogen transfer between the two reactants. In addition, this new integrated process produces one desired product (GBL) with a substantially increased hydrogenation yield. The new combined reaction can be expressed as follows: CHO+ 1.5 CH O 2.5 CHO + HO H = -119 kj/mol (3) 4 2 3 4 1 2 4 6 2 2 The MA to BDO ratio used in the process is 1:1.5 to balance the hydrogen consumed and produced in the system. Due to this perfect hydrogen balance it was expected to easier control of the reactor temperature due to less heat released as compared to individual MA hydrogenation. From equation (3), it can be noted that 2.5 mole of GBL requires 1 mole of MA and 1.5 moles of BDO. The system can also be coupled on the basis of balancing heat. The combined reaction can be expressed as follows: CHO+ 3 CH O 4 CHO + HO + 3 H H = -25 kj/mol (4) 4 2 3 4 1 2 4 6 2 2 2 From equation (4) it can be noticed that 4 moles of GBL requires 1 mol of MA and 3 moles of BDO. The MA to BDO ratio used in the process is 1:3 to balance the heat consumed and produced in the system. Due to this almost perfect heat balance the adiabatic operation of the chemical reactor becomes an attractive option. The combined system for hydrogenation of MA and the dehydrogenation of BDO could be carried out over the same Cu-based catalyst, and under similar reaction conditions. The reaction kinetics used for the simulations are extracted from the literature based on Cu-Zn-Al catalyst: Maleic anhydride hydrogenation [4]: k 1 kmol m s Pa (5) = 5 4.68 1 exp (-58 / T -3-1 -1 )

Integrated process for γ-butyrolactone production 35 1, 4-butanediol dehydrogenation [5]: k 1 = 17 2.57 1 exp (-262 / T -3-1 -1 ) kmol m s Pa (6) The coupled reactions can be performed at temperatures which are 3 5 C less as compared to individual reaction, because the active hydrogen species on the catalyst surface from BDO dehydrogenation promotes the MA hydrogenation reaction [6], [7]. The coupling may lead to advantages, such as good energy efficiency, the easy temperature control in an adiabatic reactor due to its moderate heat release compared with the single reaction processes which can avoid the formation of apparent hotspots or cold spots [6], [7]. The previous studies [6], [7] considered the possibility of coupling these reactions taking into account only the reactor. However, no plantwide design and control results have been reported. The goal of this paper is to fill this gap by presenting details of an integrated process for 1,4-butanediol (BDO) dehydrogenation and maleic anhydride (MA) hydrogenation, to produce one valuable product, γ-butyrolactone (GBL). The advantages of optimum hydrogen and energy utilization, reduced parametric sensitivity, dynamic stability and safe operation are emphasized. Our conclusions are supported by rigorous steady state and dynamic simulations performed in Aspen Plus and Aspen Dynamics. Conclusion and final remarks end the paper. 2. Integration by Hydrogen Balance This section presents a process which couples dehydrogenation of 1,4- butanediol and hydrogenation of maleic anhydride in a single adiabatic reactor. The benefit of coupling the exothermic and endothermic reactions (hydrogenation and dehydrogenation) is that the heat generated by the first reaction is taken up by the second one, and no additional heating or cooling arrangements are required in the system. This leads to a much simpler reactor and the parametric sensitivity of the reactor is greatly reduced. The conversion of both maleic anhydride and 1, 4- butanediol is > 99.5%, eliminating recycling of reactants. Fig. 1 presents the flowsheet of the integrated maleic anhydride hydrogenation and 1, 4-butanediol dehydrogenation process. Reactants (MA 2 kmol/h, BDO 3 kmol/h and hydrogen 5 kmol/h) are mixed and brought to reaction temperature (21 ºC). The reactions take place in a tubular reactor (5 m diameter, 3 m length) which is adiabatically operated. The reactor is designed for total conversion of reactants. Stream table of the process is given as Table 1. Fig. 2 shows mole fraction and temperature profiles along the reactor. The conversion of reactants into products is uniformly distributed along the reactor length and almost constant temperature should be remarked.

36 Ahtesham Javaid, Costin Sorin Bildea The reactor outlet is cooled to 2 ºC and vapor liquid separation is performed. As the complete conversion (> 99.5%) of reactants is achieved, there is no recycling of reactants simplifying technical difficulties. A high boiling point difference between water and GBL results in easier separation by distillation. GBL is obtained as a product at the rate of 49.6 kmol/h with high purity > 99.5%. In theory, the plant is best operated at MA: BDO = 1: 1.5 with very high ratio of hydrogen: feed = 1:1 to :1. Perfect balancing of hydrogen consumption / generation in the two reactions is impossible in practice. Therefore, a small amount of fresh and purge hydrogen streams are provided for better control purposes. 1 PURGE COMPRESR H2 H2MIXER SPLITER H2REC H2IN 13 PURG MA 1 BDO RIN ROUT WATER MIXER HEATER 2 SEPIN GBL REACTOR COOLER H2SEP FLASH Fig. 1. Integrated Process 25 23 21 Mole Fraction BDO MA GBL.5 1 1.5 2 2.5 3.5 1 1.5 2 2.5 3 Fig. 2 Temperature and Molar Profile

Integrated process for γ-butyrolactone production 37 Table 1 Stream Table (Hydrogen Balance) MA BDO GBL H2REC RIN ROUT SEPIN WATER Temperature ( o C) 25 25 2 2 21 23.3 2 2 Pressure (bar) 1.5 1.5 1.2 1.2 1.5 1.5 1.2 1.2 Vapor Fraction 1 1 1.1 Mole Flow (kmol/hr) 2 3 49.62 4548 598 5118 64.56 14.94 Mole Fraction Hydrogen.989.981.977.1 Water.9 3.228.984 MA 1.7 GBL 1.769.3 BDO 1.1.5 3. Integration by Heat balance Reactants (MA 2 kmol/h, BDO 6 kmol/h and hydrogen 3 kmol/h) are mixed and brought to reaction temperature ( ºC). The reactions take place in a tubular reactor (4 m diameter, 3 m length) which is adiabatically operated. The reactor is designed for the total conversion of reactants. Stream table of the process is given as Table 2. Fig. 3 shows mole fraction and temperature profiles along the reactor. The conversion of reactants into products is uniformly distributed along the reactor length and almost constant temperature should be remarked. 28 26 18 1 2 3 Mole Fraction.25.2 BDO 5.5 MA GBL 1 2 3 Fig. 3 Temperature and Molar Profile

38 Ahtesham Javaid, Costin Sorin Bildea Table 2 Stream Table (Heat Balance) MA BDO GBL H2REC RIN ROUT SEPIN WATER Temperature ( o C) 25 25 23 2 235 2 45 Pressure (bar) 1.5 1.5 1 1.2 1.5 1.5 1.2 1 Vapor Fraction 1 1 1 Mole Flow (kmol/hr) 2 6 79.84 3 331 339 64.56 19.22 Mole Fraction Hydrogen.991.967.962 Water.9.9 4.228.996 MA 1 GBL 1.24.769 BDO 1.18.1 The reactor outlet is cooled to 2 ºC, and vapor liquid separation is performed. The hydrogen is recycled. The amount of hydrogen generated by the reactions is withdrawn and can be stored and used for other purposes. As the complete conversion (> 99.5%) of reactants is achieved, there is no recycling of reactants simplifying technical difficulties. A high boiling point difference between water and GBL result in easier separation by distillation. GBL is obtained as a product at the rate of 79.8 kmol/h with high purity > 99.5%. In heat coupled system, the plant can be operated at much lower ratio, hydrogen : feed (BDO + MA) = 4:1 as compared to 1:1 in hydrogen coupling. Even at this low ratio, the temperature is under control. 4. Parametric Sensitivity The chemical reactor deserves a special discussion. To remove the risk of run-away in hydrogenation of maleic anhydride as an individual process a large excess of hydrogen to feed ratio (1 : 1) is provided to the system. Despite this, the reaction is very difficult to control. The operation of the integrated reactor is simulated, at various process conditions, in order to emphasize the stability of reactor with reduced parametric sensitivity. Fig. 4 and Fig. 5 show the effect of changing the reactor-inlet temperature and the hydrogen to feed ratio on the products mole fractions and temperature respectively.

Integrated process for γ-butyrolactone production 39 Mole Fraction / [MA] 28 26 Feed Tempearture: o C 21 o C o C 18 1 2 3.3.1 o C 21 o C o C 1 2 3 2 1 2 3 Fig. 4 Effect of Feed Temperature (Hydrogen Balance) Mole Fraction / [GBL] Mole Fraction / [BDO].5.3 o C 21 o C o C.1 o C 21 o C o C 1 2 3 Mole Fraction / [MA] 28 26 H 2 :Feed = 8:1 12:1 18 1 2 3 H 2 :Feed = 8:1 12:1 1 2 3 4 2 1 2 3 Fig. 5 Effect of Hydrogen to Feed Ratio (Hydrogen Balance) Fig. 6 and 7 shows the effect of changing the reactor inlet feed temperature and hydrogen to feed ratio on molar and temperature profiles for the heat-balance setup. The reduced reactor sensitivity with respect to disturbances is observed. Mole Fraction / [GBL] Mole Fraction / [BDO] H 2 :Feed = 8:1 12:1 H 2 :Feed = 8:1 1:1 12:1 1 2 3

4 Ahtesham Javaid, Costin Sorin Bildea 28 26 Feed Tempearture: 23 o C o C 21 o C o C 18 1 2 3 Mole Fraction / [GBL] 25.75.5 5 1 2 3 23 o C o C Mole Fraction / [MA] 5.15.1 o C.5 23 o C 1 2 3 Mole Fraction / [BDO] Fig. 6 Effect of Feed Temperature (Heat Balance) o C 23 o C 1 2 3 28 26 H 2 :Feed = 4:1 12:1 18 1 2 3 Mole Fraction / [GBL].25.2 5.5 H 2 :Feed = 4:1 12:1 1 2 3 Mole Fraction / [MA].75 5 H 2 :Feed = 4:1.3.15 12:1 1 2 3 6 2 12:1 1 2 3 Fig. 7 Effect of Hydrogen to Feed Ratio (Heat Balance) Mole Fraction / [BDO].2 H 2 :Feed = 4:1

Integrated process for γ-butyrolactone production 41 5. Plantwide Dynamics and Control The plantwide control system fixes the flow rate of maleic anhydride, 1,4- butanediol and hydrogen at reactor inlet. Control of heat exchangers, two-phase separator and distillation column is standard. A dynamic model was built in Aspen Dynamics and flexibility in changing the production rate was investigated. For hydrogen-balance and heat-balance designs, starting from steady state, the MA and BDO feed rates were increased or decreased, as shown in Fig. 8. The change in feed rate is followed by change in production rate of GBL while purity remains the same. Molar Flow /[kmol/hr] 6 1 5 GBL.99 4 BDO.98 3 2 MA.97 1 1% 1% 1% 1%.96.95 1 2 3 4 5 Time / [hr] Purity / [GBL] Molar Flow / [kmol/hr] 1 8 6 4 2 GBL BDO MA 5% 5% 5% 5%.99.98.97.96.95 1 2 3 4 5 Time / [hr] 1 Purity / [GBL] Fig. 8 Dynamic Simulation (Hydrogen Balance - left, Heat Balance - right) 6. Conclusions The coupling of exothermic maleic anhydride hydrogenation with the endothermic 1,4-butanediol dehydrogenation leading to γ-butyrolactone can be achieved with overall reactants conversion exceeding 99.5% and product purity above 99%. The coupling can be done according to hydrogen balance or energy balance, the later one being preferred. Compared to a conventional process, where the reactions are carried on in two separate reactors and energy is transferred from one reactor to another by means of a heat carrier, the coupled process offers several advantages: - reaction coupling (energy balance) allows reduction of the hydrogen excess by at least a factor 2.5 (from 1- to 4). - the reaction system is much simpler, namely one adiabatic reactor instead of two multi-tubular reactors with heating and cooling, respectively. There is no need of a heat carrier, and the hydrogen recycle system is much simpler. The plant is controllable, allowing for changes in production rate while keeping high product purity.

42 Ahtesham Javaid, Costin Sorin Bildea R E F E R E N C E S [1]. R. E. Kirk, F. Othmer, J. I. Kroschwitz, M. Howe-Grant, Kirk Othmer s Encylopedia of Chemical Technology, 4 th ed., John Wiley & Sons, New York, 1991, pp. 211. [2]. M. J. Bergfeld, K. Uihlein, Process for the preparation of gamma-butyrolactone and the use thereof, in U.S. Patent 675153, June. [3]. I. Tatsumi, K. Kenji, S. Sadakatsu, U. Hiroshi, Process for the preparation of gammabutyrolactone, Eur. Patene EP58448, 1994. [4]. N. Ichikawa, S. Sato, R. Takahashi, T. Sodesawa, K. Inui, Dehydrogenative cyclization of 1,4-butanediol over copper-based catalyst, in Journal of Molecular Catalysis A: Chemical, vol. 212, April 4, pp. 197-23. [5]. G. W. Zhao, Y. L. Zhu, B. T. Teng, H. Y. Zheng, H. W. Xiang, Y. W. Li, Kinetics Study on Maleic Anhydride Gas-Phase Hydrogenation to γ-butyrolactone, in Journal of Chemical Engineering of Chinese Universities, vol. 2, no. 5, October 6, pp. 74-744, [6]. Y. L. Zhu, H. W. Xiang, G. S. Wu, L. Bai, Y. W. Li, A novel route for synthesis of γ- butyrolactone through the coupling of hydrogenation and dehydrogenation, in Chem. Communication, January 2, pp. 254 255. [7]. Zhu, Y.L., Yang, J., Dong, G.Q., Zheng, H.Y., Zhang, H.H., Xiang, H.W., Li, Y.W., An environmentally benign route to γ-butyrolactone through the coupling of hydrogenation and dehydrogenation, in Applied Catalysis, vol. 57, May 5, pp. 183 19.