Nickel Oxide Reduction by Hydrogen: Kinetics and ... - ACS Publications

Jun 22, 2015 - Department of Electrical Engineering, University of Notre Dame, ... Department of Chemistry, Yerevan State University, Yerevan 0025, Ar...
2 downloads 0 Views 9MB Size
Subscriber access provided by UNIV OF MISSISSIPPI

Article

Nickel Oxide Reduction by Hydrogen: Kinetics and Structural Transformations Khachatur V. Manukyan, Arpi Avetisyan, Christopher Shuck, Hakob Chatilyan, Sergei Rouvimov, Suren Kharatyan, and Alexander S. Mukasyan J. Phys. Chem. C, Just Accepted Manuscript • DOI: 10.1021/acs.jpcc.5b04313 • Publication Date (Web): 22 Jun 2015 Downloaded from http://pubs.acs.org on June 23, 2015

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.

The Journal of Physical Chemistry C 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 22

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

The Journal of Physical Chemistry

Nickel Oxide Reduction by Hydrogen: Kinetics and Structural Transformations

Khachatur V. Manukyan1,2*, Arpi G. Avetisyan2, Christopher E. Shuck3, Hakob A. Chatilyan2, Sergei Rouvimov4, Suren L. Kharatyan2,5, Alexander S. Mukasyan3* 1

Department of Physics, University of Notre Dame, Notre Dame, Indiana 46556, United States 2

3

4

Laboratory of Kinetics of SHS processes, Institute of Chemical Physics NAS of Armenia, Yerevan, 0014, Armenia

Department of Chemical and Biomolecular Engineering, University of Notre Dame, Notre Dame, Indiana 46556, United States

Department of Electrical Engineering, University of Notre Dame, Notre Dame, Indiana, 46556, United States 5

Department of Chemistry, Yerevan State University, Yerevan, 0025, Armenia

1 ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

ABSTRACT We studied the reduction kinetics of bulk NiO crystals by hydrogen and the corresponding structural transformations in the temperature range 543-1593 K. A new experimental approach allows us to arrest and quench the reaction at different stages with millisecond time resolution. Two distinctive temperature intervals are found where the reaction kinetics and product microstructures are different. At relatively low temperatures, 543-773 K, the kinetic curves have a sigmoidal shape with long induction times (up to 2000 s) and result in incomplete conversion. Low temperature reduction forms a complex polycrystalline Ni/NiO porous structure with characteristic pore size on the order of 100 nm. No induction period was observed for the high temperature conditions (1173-1593 K) and full reduction of NiO to Ni is achieved within seconds. An extremely fine porous metal structure, with pore size under 10 nm, forms during high temperature reduction by a novel crystal growth mechanism. This consists of the epitaxiallike transformation of micron-sized NiO single crystals into single crystalline Ni without any crystallographic changes, including shape, size, or crystal orientation. The Avrami nucleation model accurately describes the reaction kinetics in both temperature regimes. However, the structural transformations during reduction in both nano- and atomic levels are very complex and the mechanism relies both on nucleation and the critical diffusion length for outward diffusion of water molecules.

2 ACS Paragon Plus Environment

Page 2 of 22

Page 3 of 22

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

The Journal of Physical Chemistry

INTRODUCTION Nickel is an important industrial catalyst and is typically produced by reduction of nickel oxide (NiO).1,2,3,4,5 This metal is also a promising oxygen carrier in chemical looping combustion, which is a novel approach for green power generation.6,7,8 NiO reduction and the subsequent Ni behavior in hydrogen is of practical importance in the solid oxide fuel cells9,10, where it determines the structure of the electronic conductor on the anode side. It is also a model process for investigation of oxide reduction mechanisms2,11,12,13,14,15 due to its single-stage transformation, as compared to Fe2O316,17,18,19 and CuO.20,21,22 Despite extensive investigations of NiO reduction by hydrogen, since the pioneering work23 of Benton and Emmett in 1924, many questions still remain about its kinetics and the related structural transformation mechanism. For example, the reported activation energy varies by an order of magnitude, between 10 and 150 kJ/mol.2,11,24,25,26 This large variation may be explained by the lack of consistency between NiO samples, characterization techniques, and other experimental parameters. Furthermore, different kinetic models, including reaction-order2, Avrami27,28,29, geometrical construction26 and Szekely-Evans30,31 were proposed to describe this reaction. Several in-situ techniques were used to investigate the mechanism of NiO reduction. Richardson et al. studied2 the reduction of porous NiO particles by hydrogen using a timeresolved X-ray diffraction technique in the temperature range of 450-570 K. A three-step reaction mechanism was proposed: (i) nickel clusters form during an induction period; (ii) next, the reaction rate accelerates as the clusters grow; (iii) finally, the process settles to a pseudo-firstorder reaction with respect to nickel. Jeangros et al. (573-923 K) used an environmental transmission electron microscopy (TEM) coupled with electron energy loss spectroscopy for in situ study of hydrogen reduction of NiO,11 Based on a model-based fitting procedure, they suggested that the Avrami (nucleation and nuclei growth) model describes the reaction. This assumes that hydrogen molecules predominantly dissociate on nickel atoms that are surrounded by oxygen vacancies.32 Rodríguez et al. used synchrotron time-resolved X-ray adsorption and diffraction techniques to study the reduction of single crystalline and polycrystalline NiO between 523-620 K.1 They found a strong correlation between the rate of reduction and the amount of oxygen defects in the metal oxide structure. It was suggested that the presence of oxygen vacancies increases the adsorption energy of H2 on the oxide surface, as well as 3 ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

decreases the energy barrier for hydrogen dissociation. H atoms then diffuse to the NiO reaction sites, resulting in a rupturing of the Ni-O bond with simultaneous desorption of a water molecule. Studies using optical microscopy by Utigard et al33 led to the conclusion that the geometric construction, or shrinking core model represents the reduction kinetics of NiO granules in the range of 100-500 µm. This model assumes that, after the reduction reaction has begun, the solid particles consist of a non-reacted core encased in a uniform layer of the reaction product. These two regions are sharply separated from each other by a well-defined interface. As the reaction proceeds, diffusion of hydrogen through the Ni layer becomes the rate limiting step. Later Plascencia and Utigard showed34 that the reduction of large porous NiO particles can be described by the grain model proposed by Szekely and Evans.35 In this empirical model, a particle is considered to be an agglomerate of individual grains and each of them undergoes a microscopic shrinking core reduction. Depending on the process conditions and characteristics of NiO, either chemistry (low temperatures, high porosities) or internal diffusion (high temperatures, low porosities) controls the reaction kinetics. Recently, Hidyat et al. identified several distinctive porous microstructures36, during the reduction of bulk dense NiO samples. Jeangros et al. also reported11 the macropore (50-100 nm) formation during early stages of NiO nanoparticle reduction by hydrogen. However, towards the end of reduction, the pores disappeared. The severe collapse of pores has also been reported during the reduction of hollow NiO nanoparticles at 550-870 K.37,38,39 Evidently, the hydrogen reduction of nickel oxide is a complex process, which includes a variety of topological and morphological factors. It is also clear that the microstructural evolution during reduction and its relation to the kinetic models is not well understood. Herein, we report the results of reduction kinetics of NiO by hydrogen in a wide temperature range (5431593 K). A novel experimental method permits arresting and rapidly quenching (with 104 K/s cooling rate) the reaction at different stages within milliseconds. This approach, coupled with ex situ microscopy techniques allows us to measure the reaction rate and precisely track the ongoing structural transformation. The mechanism of NiO reduction by hydrogen is suggested based on direct correlations between reaction kinetics and the observed structural transformations at the micro-, nano-, and atomic levels.

4 ACS Paragon Plus Environment

Page 4 of 22

Page 5 of 22

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

The Journal of Physical Chemistry

■ EXPERIMENTAL SECTION The experiments were performed using partially oxidized nickel wires (Supporting Figure S1). To prepare these specimens an initial Ni wire (purity 99.995%, 100 µm and 8.5 cm in length) was heated by direct current to 1573 K and retained at this temperature for 150 s in air. Scanning electron microscopy (SEM) and transmission electron microscopy (TEM) images (Figure 1A, B) illustrate that the formed NiO layer is polycrystalline with grain size in the range 1-3 µm. Analysis of the select area electron diffraction pattern (Figure 1C) and high-resolution TEM images (Figure 1D) reveal that the individual NiO grains are pore-free single crystals.

Figure 1: Microstructure of NiO: Surface of wire (A), bright field TEM image (B) of a thin slice of oxide layer, diffraction pattern (C) and a high resolution TEM image (D) A reaction chamber equipped with infrared-transparent windows, high-speed solar cell-based temperature sensors, a power source, and a PC controlled unit (Figure 2S), was used in this hydrogen reduction experiment (Figure S2). The reaction setup allows rapid (up to 4.5·105 K/s) controllable heating of the wires (Figure S3) and continuous data (temperature, electrical power, resistivity of the wire) acquisition with frequency of 10 kHz40,41 (see details in Supporting Information). These NiO/Ni wires were inserted into the chamber, which was then evacuated to 10-3 kPa, purged with pure hydrogen and finally filled with H2 to the desired pressure from 1.0 to 7.0 kPa. Following each experiment, no pressure change was observed. Utilizing the internal Ni metallic core, the wires were controllably heated by direct current to the desired temperature (543-1593 K). The termination of electrical current through the thin wire allows rapid arresting of the process with a quenching rate of ~ 104 K/s (Figure S3). This feature permits both arresting of the reaction and preservation of the sample structure with millisecond time resolution.

5 ACS Paragon Plus Environment

The Journal of Physical Chemistry

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 22

A Mettler-Toledo balance (measuring accuracy of 10-5 g) was used to obtain the sample weight at different stages of the process to determine the degree of conversion (α) during reduction reaction as follows: α = −∆m/∆mo

(1)

where ∆mo – mass change of Ni wire after initial oxidation and ∆m is the subsequent mass change of the Ni/NiO wire during the reduction reaction. Obtained kinetic curves were analyzed using custom Matlab™ nonlinear least squares regression analysis. A variety of models, totaling 18, were tested, including reaction-order, geometrical contraction (R1, R2, R3), one-, two- or three-dimension diffusion-limited (D1, D2, D3), Avrami-Erofeev (AE) nucleation models, and others.26,42 The fitting parameters including pre-exponential factor (ko), activation energy (Ea), and exponential term (n) were varied for each model; with the two parameter models only having ko and Ea fit. All models were analyzed in series; with the varying parameters ranging over 20 log units of perturbation over the local attempt-space. The results with the highest goodness of fit (R2) were compared across models. Each kinetic dataset was analyzed both individually and in every inclusive range. The inclusive datasets with the best goodness of fits were interpreted to be within the same kinetic regime. SEM and TEM were employed to characterize the composition and morphology of the reaction products, as well as the atomic structure of the materials. SEM analysis was conducted on a Magellan 400 (FEI, USA) with resolution of 0.6 nm. The Magellan 400 is equipped with an energy dispersive X-ray spectrometer (EDS, Bruker) with energy resolution of 123 eV. A Titan 80-300 (FEI, USA) transmission electron microscope, with resolution of 0.136 nm in scanning TEM mode and about 0.1 nm information limit in high resolution TEM mode was used. The Titan is equipped with an energy dispersive X-ray spectroscopy (EDS, Oxford Inca) system with spectral energy resolution of 130 eV. A Helios NanoLab 600 system (FEI) with dual electron/ion beam was used to produce cross-sections of reacted wires by ion milling with a gallium-ion beam. First, a ~2 µm thick layer of platinum in a 15µm × 1µm rectangular area was deposited onto the selected sample. Next, a trench with ~10 µm width, ~15 µm length, ~7.5 µm depth and a 45° base angle was milled on the particle surface under an accelerating voltage of 10 kV with a milling current of 27 nA. Then, a 40 nm thick layer was milled into the side wall of the trench under an accelerating voltage of 10 kV and a milling current of 700 pA to produce a clean cross-section without milling artifacts. 6 ACS Paragon Plus Environment

Page 7 of 22

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

The Journal of Physical Chemistry

Thereafter, the electron-beam was used to image the side wall of the trench at an angle of 52° parallel to the side wall. The TEM samples were also prepared using a Helios NanoLab 600 system by making cross-sectional slices from the top surfaces of the wires. The slice and view (S&V, FEI) software package was used on the Helios NanoLab 600 to collect a series of sequential images. A voltage of 30 kV with a milling current of 9.7 pA was used to make a series of 100 images with ~10 µm width and ~7.5 µm depth. These images were taken in series after milling 10 nm in between each image, leading to a particle volume of ~75 µm3 imaged throughout and prepared for pore characterization. ■ RESULTS Below we report the kinetics for NiO reduction by hydrogen at different temperatures (e.g. 543 K; 608 K; 1503 K and 1593 K) at constant H2 pressure (e.g. 1.33 kPa and 6.67 kPa), as well as the characteristic structures of the reaction media on micro-, nano-, and atomic scales at different stages of the process. Reaction Kinetics Kinetics curves of NiO reduction by hydrogen, which represent the degree of conversion (α) vs time (t) at different temperatures and constant H2 pressure, are presented in Figure 2. It can be seen that, at 543 K (Figure 2A) and 608 K (Figure 2B), the kinetic curves have a sigmoidal shape with low reaction rates at initial stages, followed by a significant acceleration of the process after some critical time (induction period). Supporting Table 1 shows that the induction period depends on the temperature and pressure, and can be as long as 2000 s. Another important observation is that full conversion cannot be achieved. The kinetics curves at 1503 K (Figure 2C) and 1593 K (Figure 2D) possess different behavior. No induction period is detected and the degree of conversion monotonically increases with time in a nonlinear (parabolic) fashion. Figure 1D also demonstrates that the hydrogen pressure has essentially no influence on the reaction rates. Finally, complete reduction of NiO could be achieved within 10-30 s.

7 ACS Paragon Plus Environment

The Journal of Physical Chemistry

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 22

Figure 2: Kinetics of NiO reduction at 543 K (A), 608 K (B), 1503 K (C) and 1593 K (D) at 6.67 kPa (solid point) and 1.33 kPa (open points) hydrogen pressures The results of the kinetic nonlinear least squares analysis indicate that only the nth order Avrami-Erofeev (AEn) model is relevant for the low temperature reduction (1173 K) kinetics. The AEn model has following form: αNiO = 1 − exp(−koexp(-Ea/(RT))tn)

(2)

where the Avrami exponent (n) varies from 0.5-4, with each value describing a specific type of kinetic limitation; values below 1 indicate a diffusion-limited reaction, with integers 1 and above corresponding to growth in n that number of dimensions. The D1 model takes the form:

αNiO = (koexp(-Ea/(RT))t)(1/2)

(3)

The results of model-based fitting analysis for 608 K and 1503 K temperatures are shown in Figure 2B and Figure 2C, respectively. Kinetic data at 608 K fits with R2=0.98 by using AEn model with a value of exponent n=0.58 (Figure 2B). The kinetic data for 1503 K can be also

8 ACS Paragon Plus Environment

Page 9 of 22

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

The Journal of Physical Chemistry

accurately (R2=0.99) fit by the AEn model with n = 0.95 (Figure 2C), as well as by the D1 model (R2=0.95). All other models show less fitting accuracy. Figure 3 shows the n parameter as a function of temperature over the entire studied range. In low temperatures (543-773 K), the parameter n does not vary. This would imply that both nucleation and a diffusion-related processes slow the reaction rate. At high temperatures (1173-1593 K), the AEn model fits each individual kinetic curve. However, in this case, the

Avrami

exponent

increases

with

temperature, approaching 1. As the temperature increases, nucleation becomes less important, but the diffusion limitation remains. The activation energies of nickel oxide reduction determined by model-based fitting

Figure 3: Avrami exponent vs. temperature

analysis are ~84 and ~109 kJ/mol for the low and high temperature ranges, respectively. These values are comparable to previous reported data in the literature11,26. In order to validate the kinetic modeling, the structural transformations were studied.

Microstructure transformation We selected samples heated at 608 K and 1503 K (under 6.67 kPa) with different reduction times to investigate the relationship between structural transformations and the reduction kinetics. Figure 4 summarizes the microstructure analysis of the samples reduced at early stages of the reaction for both temperatures. A porous structure is predominantly formed on the crystal vertices close to the boundaries between the NiO grains after heating the sample at 608 K for 50 s (Figure 4A). However, this structure accounts for less than 1% the total wire surface. EDS analysis of the cross-sections (Figure 4B) show approximately 35 at.% O and 65 at.% Ni in these porous regions. Kinetic data in Figure 2B suggests that the degree of conversion for this sample is essentially zero, indicating that Ni nucleation occurs on the edges of the NiO grains forming porous Ni-rich areas.

9 ACS Paragon Plus Environment

The Journal of Physical Chemistry

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 10 of 22

Figures 4C illustrates the wire surface after heating to 1503 K for only 0.007 s. The facets that are present on the initial NiO grains are not apparent after this short reduction time. The surface, however, involves two characteristic regions with different contrast. EDS analysis

Figure 4: Microstructure of surfaces (A and C) and cross-sections (B and D) of reacted samples at early stages of reduction. reveals that the darker areas contain oxygen and nickel, while the lighter ones are pure Ni. The cross-section (Figure 4D) of this surface shows a relatively uniform, ~300 nm thick layer, which consists of a primarily pure Ni phase with some NiO residue. This indicates that, at 1503 K, nucleation is essentially instantaneous, and the following growth of the reduced layer is very fast (at a rate of ~40 µm/s). SEM analysis of the samples with significant degrees of conversion is shown in Figure 5. Figure 5A shows that the grain size on the surface of a wire, after reduction at 608 K for 300 s, is similar to the initial NiO (see Figure 1A). Reduction at 608 K leads to the formation of macropores, which range from 50-100 nm and are uniformly distributed along the grain surfaces.

10 ACS Paragon Plus Environment

Page 11 of 22

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

The Journal of Physical Chemistry

The cross-section of this sample contains two distinct layers, one with a complex porous structure and the other pure NiO (Figure 5B). Utilizing sequential cross-sections (Video 1S) it is observed that the pores form a continuous network of channels throughout the reduced layer, originating from grain boundaries or the outer wire surface (Figure 4S, Video 1S). A different structure is observed after reduction at 1503 K for 8 s (Figure 5C). The topology of the sample surface is much smoother and contains no macropores compared to the low temperature reduction. A cross-section of the sample shows that a uniform, oxygen-free Ni layer forms on the surface of the wire, with a sharp interface separating this and the NiO layer (Figure 5D).

Figure 5: Microstructure of surfaces (A and C) and cross-sections (B and D) of reacted samples with significant degrees of conversion. High resolution EDS profiles of the sample, heated at 608 K for 300 s, measured across the reacted layer show approximately 25 at.% O and 75 at.% Ni (Figure 6A), indicating ~0.5 degree of conversion. This analysis is in good agreement with the degree of conversion measured by the sample weight change (Figure 2B). EDS profiles also show that, closer to the non-reduced NiO grains, the concentration of Ni gradually decreases to 50 at.% with a simultaneous increase 11 ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

of O to 50 at.%. The profiles of the reduced layer heated at 1503 K for 8 s (Figure 6B) indicate that the Ni concentration is ~100 at.%. An abrupt decrease of nickel and simultaneous increase of

Figure 6: EDS analysis of wire cross-sections heated at 608 K for 300 s (A) and 1503 K for 8 s (B) O to ~50 at.% takes place at the interface between the reduced and non-reduced layers. TEM analysis for both samples heated at 608 K and 1503 K are shown in Figure 7 and Figure 8, respectively. A bright field TEM image (Figure 7A) of the sample reduced at 608 K shows a transitional area between the reduced layer and original NiO grain. The

electron

diffraction

pattern (Figure 7B) of the reduced layer contains two distinct

phases;

fine

Figure 7: TEM analysis of wire heated at 608 K for 300 s: Bright field image of NiO/Ni interface (A), diffraction pattern (B) and structure of porous reacted layer (C), as well as highangle annular dark field scanning TEM image of that layer (D) 12 ACS Paragon Plus Environment

Page 12 of 22

Page 13 of 22

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

The Journal of Physical Chemistry

polycrystalline Ni (scattered rings) and single-crystalline NiO (big bright spots). Bright field TEM (Figure 7C) of the reduced layer shows that, pore-free NiO grains are transformed into a complex NiO-Ni porous structure, with both phases in well-defined regions. The fine polycrystalline Ni occurs along the pore walls within a single crystal NiO skeleton. High-angle annular dark field scanning TEM (Figure 7D) together with EDS analysis suggests that the brighter phase, closer to the surface of pores, is Ni, while the darker one is NiO TEM studies reveal that, at high temperatures, NiO grains convert into Ni without changing the grain geometry, including shape, size, or crystal orientation (Figure 8A and B). The reaction forms a network of mesoporous channels, ranging from 5-15 nm, that end at the grain boundaries (Figure 8C) or at the wire surface. Figure 8D shows the nanosized pores at the Ni/NiO interface and the corresponding channels.

Figure 8: TEM analysis of a wire heated 1503 K for 8 s: Bright field image of reduced (A), diffraction pattern (B), structure of porous single crystalline Ni layer (C), and high resolution TEM image NiO/Ni interface (D)

■ DISCUSSION 13 ACS Paragon Plus Environment

The Journal of Physical Chemistry

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 14 of 22

The kinetics of nickel oxide reduction by hydrogen were extensively studied both using supported and unsupported (bulk) samples. In this work we use unsupported polycrystalline samples with 1-3 µm single crystalline NiO (Figure 1A and B). Different kinetic models were suggested to describe the hydrogen reduction of NiO with similar characteristics. Formal fitting of the kinetic data, however, is not sufficient to characterize the reaction mechanism due to the complex nature of the process. We combined kinetic modeling with ex situ microscopy analysis to understand the direct correlations between the reaction kinetics and structural transformations. Our results show that, at 543-773 K, the Avrami nucleation model accurately describes the reaction kinetics. The kinetic data also indicates that there is a long induction period where the degree of conversion is zero (Figure 2A and B). At this stage, the Ni nuclei appeared on the outer surface of the NiO grains. This process is slow and takes hundreds to thousands of seconds to reach a measurable degree of conversion. Initial nucleation leads to the formation of small channels that propagate into the bulk of the NiO grains (Figure 4B and Figure 9). We have estimated that approximately 41% of the NiO volume will be converted into a porous structure upon complete reduction assuming that the overall crystal volume of the initial pore-free NiO remains constant. The reduction is accompanied by generation of water vapor that diffuses out through the channels. Channel growth proceeds both into the bulk and radially outward, forming the polycrystalline Ni nano-layer that coats the inner structure of the channels.

Figure 9: Schematic illustration of NiO reduction mechanism at different temperatures

14 ACS Paragon Plus Environment

Page 15 of 22

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

The Journal of Physical Chemistry

TEM observations suggest that Ni nano-layers grow extremely slowly (Figure 7). The low growth rate occurs due to the slow outward diffusion of water vapor, resulting from the pore-free Ni nano-layer (Figure 7D). Therefore, the thickness (δ) of the Ni nano-layer is a critical parameter and the growth of this layer becomes the limiting stage of the reaction (Figure 9). Due to this diffusion limitation, conversion remains low at temperatures below 773 K even for prolonged experiment duration. To avoid this limitation, the average distance between neighboring channels must not exceed 2δ. This can be accomplished by increasing the temperature, which increases the nucleation density. To summarize, at 543-773 K, the following steps occur: (i) Ni nucleates predominately on the surface vertices and grain boundaries of the NiO crystals, (ii) pore formation occurs at the nucleation sites, growing into the NiO bulk, (iii) additional Ni nucleation sites form on the inner surface of these channels, resulting in the formation of Ni nano-layers; (iv) channels begin branching, leading to an interconnected network. At 1173-1593 K, complete conversion of the NiO samples occurs. A sharp interface is evident between the unreacted NiO and completely reduced Ni. This may indicate that the geometrical construction or diffusion-limited (D1) models are suitable to describe the kinetics. However, the Avrami model best fits the kinetic data, suggesting instantaneous nucleation and one-dimensional growth. This implies that the Ni nucleation density on the outer surface of the NiO crystals is much higher than for low temperature reduction. The increased nucleation density may be due to the higher equilibrium concentration of oxygen vacancies, which act as nucleation sites. As a result, a continuous layer of porous reduced Ni covers the entire wire surface and rapidly moves toward the core (Figure 5D and Figure 9). The number of channels is much higher than the low temperature case and the distance between them is less than 10 nm, this explains the absence of an induction period. Because of the remarkably high reduction rates (Figure 2C and D) and small distances (less than 2δ) between the channels, the NiO grains proceed to full conversion (Figure 9). Similar to the low temperature reduction case, the channels are formed due to the shrinking of the cubic NiO lattice into cubic Ni. There are many similarities in the reduction process for both temperature ranges, while the main difference is the nucleation density. Due to this, the AEn model better describes both kinetics. Owing to high nuclei density, the process only proceeds in the forward direction,

15 ACS Paragon Plus Environment

The Journal of Physical Chemistry

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 16 of 22

without formation of polycrystalline Ni, so that the NiO grains epitaxialy transform into Ni without losing single crystallinity or changing the orientation. Several studies11,43,44 reported on the formation of both epitaxial and non-epitaxial Ni domains (islands) at the early stages of NiO reduction. However, the number of epitaxial Ni nuclei decreases as the reaction proceeds. It was shown that the loss of epitaxy results from the rotation of the Ni domains as they grow due to the appearance of cracks around them.43 It was also shown that the non-epitaxial Ni regions take longer to nucleate than the epitaxial ones, but grow faster when formed.45 These effects always induce randomly oriented polycrystalline Ni structures, similar to our samples reduced below 773 K. At high temperatures (e.g. 1503 K) the initial epitaxial Ni nuclei do not have enough time to rotate due to instantaneous nucleation. The high initial nucleation density forms a thin epitaxial-like Ni layer on the outer surface of the NiO single crystals; this leads to an increased NiO reduction rate in the lateral direction without any crystallographic changes. ■ CONCLUSIONS Hydrogen reduction of 1-3 µm single-crystalline NiO grains was studied from 543-1593 K by scheduled rapid heating and quenching of Ni/NiO wires, coupled with ex situ electron microscopy of the reacted samples. Both kinetic data and structural observations are consistent with the Avrami nucleation model. Depending on the reaction temperature, two distinctive porous product structures were identified. At 543-773 K, non-complete reduction results in a NiO-Ni product with a complex porous networked structure due to the low nucleation rates and slow outward diffusion of water molecules. Above 1173 K, the increased nuclei density causes instantaneous nucleation of Ni on the whole NiO crystal surface, followed by rapid growth of Ni phase. A new crystal growth mechanism was observed for metal during the reduction of an oxide, where the atomic lattice of single-crystalline NiO transforms into Ni without any change in crystallography. In the 773-1173 K range, the product represents a mixed microstructure of both low and higher temperature. The kinetics of the process shows some distinctive features. These investigations currently are in progress and will be reported in a follow up article. ■ ASSOCIATED CONTENT Supporting Information 16 ACS Paragon Plus Environment

Page 17 of 22

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

The Journal of Physical Chemistry

Detailed temperature measurements and calibrations, heating and cooling temperature profiles, additional SEM images of Ni and NiO/Ni wires, and duration of induction periods. Schematic representation of experimental setup. Web enhanced video of “Slice and View” sequential SEM analysis. This material is available free of charge via the Internet at http://pubs.acs.org. ■ AUTHOR INFORMATION Corresponding Authors * [email protected] (A.M.) * [email protected] (K.M)

Notes The authors declare no competing financial interest. ■ ACKNOWLEDGMENTS This work is partially supported by a joint Armenian-Russian research project (AR 13RF-057 // 13-03-90604) of State Committee of Science of Armenia and Russian Foundation for Basic Research. Also, we like to acknowledge the Notre Dame Integrated Imaging Facility (NDIFF) for SEM and TEM results of this research.

■ REFERENCES (1)

Rodriguez, J. A; Hanson, J. C.; Frenkel, A. I.; Kim, J. Y.; Pérez, M. Experimental and Theoretical Studies on the Reaction of H2 with NiO: Role of O Vacancies and Mechanism for Oxide Reduction. J. Am. Chem. Soc. 2002, 124, 346–354.

(2)

Richardson, J. T.; Scates, R.; Twigg, M. V. X-Ray Diffraction Study of Nickel Oxide Reduction by Hydrogen. Appl. Catal. A. 2003, 246 (1), 137–150.

(3)

Munnik, P.; Velthoen, M. E. Z.; de Jongh, P. E.; de Jong, K. P.; Gommes, C. J. Nanoparticle Growth in Supported Nickel Catalysts during Methanation Reaction-Larger Is Better. Angew. Chem. 2014, 53, 1–6.

17 ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

(4)

Ji, N.; Zhang, T.; Zheng, M.; Wang, A.; Wang, H.; Wang, X.; Chen, J. G. Direct Catalytic Conversion of Cellulose into Ethylene Glycol Using Nickel-Promoted Tungsten Carbide Catalysts. Angew. Chem. 2008, 120, 8638–8641.

(5)

Du, X.; Zhang, D.; Shi, L.; Gao, R.; Zhang, J. Morphology Dependence of Catalytic Properties of Ni/CeO2 Nanostructures for Carbon Dioxide Reforming of Methane. J. Phys. Chem. C 2012, 116, 10009–10016.

(6)

Adánez, J.; De Diego, L. F.; García-Labiano, F.; Gayán, P.; Abad, A.; Palacios, J. M. Selection of Oxygen Carriers for Chemical-Looping Combustion. Energy and Fuels 2004, 18, 371–377.

(7)

Erri, P.; Varma, A. Solution Combustion Synthesized Oxygen Carriers for Chemical Looping Combustion. Chem. Eng. Sci. 2007, 62, 5682–5687.

(8)

Adanez, J.; Abad, A.; Garcia-Labiano, F.; Gayan, P.; de Diego, L. F. Progress in Chemical-Looping Combustion and Reforming Technologies. Prog. Energy Combust. Sci. 2012, 38, 215–282.

(9)

Ni, M.; Leung, M.; Leung, D. Technological Development of Hydrogen Production by Solid Oxide Electrolyzer Cell (SOEC). Int. J. Hydrogen Energy 2008, 33, 2337–2354.

(10)

Fukui, T.; Murata, K.; Ohara, S.; Abe, H.; Naito, M.; Nogi, K. Morphology Control of Ni– YSZ Cermet Anode for Lower Temperature Operation of SOFCs. J. Power Sources 2004, 125, 17–21.

(11)

Jeangros, Q.; Hansen, T. W.; Wagner, J. B.; Damsgaard, C. D.; Dunin-Borkowski, R. E.; Hébert, C.; Van herle, J.; Hessler-Wyser, A. Reduction of Nickel Oxide Particles by Hydrogen Studied in an Environmental TEM. J. Mater. Sci. 2013, 48, 2893–2907.

(12)

Hidayat, T.; Rhamdhani, M. A.; Jak, E.; Hayes, P. C. The Kinetics of Reduction of Dense Synthetic Nickel Oxide in H2-N2 and H2-H2O Atmospheres. Metall. Mater. Trans. B 2009, 40, 1–16.

(13)

Hidayat, T.; Rhamdhani, M. A.; Jak, E.; Hayes, P. C. On the Relationships between the Kinetics and Mechanisms of Gaseous Hydrogen Reduction of Solid Nickel Oxide. Metall. Mater. Trans. B 2009, 40, 474–489.

(14)

Richardson, J.; Scates, R.; Twigg, M. X-Ray Diffraction Study of the Hydrogen Reduction of NiO/α-Al2O3 Steam Reforming Catalysts. Appl. Catal. A 2004, 267, 35–46.

(15)

Doppiu, S.; Langlais, V.; Sort, J.; Suriñach, S.; Baró, M. D.; Zhang, Y.; Hadjipanayis, G.; Nogués, J. Controlled Reduction of NiO Using Reactive Ball Milling under Hydrogen Atmosphere Leading to Ni-NiO Nanocomposites. Chem. Mater. 2004, 16, 5664–5669.

18 ACS Paragon Plus Environment

Page 18 of 22

Page 19 of 22

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

The Journal of Physical Chemistry

(16)

Wimmers, O. J.; Arnoldy, P.; Moulijn, A. Determination of the Reduction Mechanism by Temperature-Programmed Reduction: Application to Small Fe2O3 Particles. J. Phys. Chem. 1986, 90, 1331–1337.

(17)

Lin, H.-Y.; Yu-We, C.; Li, C. The Mechanism of Reduction of Iron Oxide by Hydrogen. Thermochim. Acta 2003, 400, 61–67.

(18)

Huang, W.; Ranke, W.; Schloegl, R. Reduction of an α-Fe2O3(0001) Film Using Atomic Hydrogen. J. Phys. Chem. C 2007, 111, 2198–2204.

(19)

Zieliński, J.; Zglinicka, I.; Znak, L.; Kaszkur, Z. Reduction of Fe2O3 with Hydrogen. Appl. Catal. A. 2010, 381, 191–196.

(20)

Kim, J. Y.; Rodriguez, J. A.; Hanson, J. C.; Frenkel, A. I.; Lee, P. L. Reduction of CuO and Cu2O with H2: H Embedding and Kinetic Effects in the Formation of Suboxides. J. Am. Chem. Soc. 2003, 125, 10684–10692.

(21)

Svintsitskiy, D. A.; Kardash, T. Y.; Stonkus, O. A.; Slavinskaya, E. M.; Stadnichenko, A. I.; Koscheev, S. V.; Chupakhin, A. P.; Boronin, A. I. In Situ XRD, XPS, TEM, and TPR Study of Highly Active in CO Oxidation CuO Nanopowders. J. Phys. Chem. C 2013, 117, 14588–14599.

(22)

Yamukyan, M. H.; Manukyan, K. V.; Kharatyan, S. L. Copper Oxide Reduction by Hydrogen under the Self-Propagation Reaction Mode. J. Alloys Compd. 2009, 473, 546– 549.

(23)

Benton, A.; Emmett, P. The Reduction of Nickelous and Ferric Oxides by Hydrogen. J. Am. Chem. Soc. 1924, 46, 2728–2737.

(24)

Erri, P.; Varma, A. Diffusional Effects in Nickel Oxide Reduction Kinetics. Ind. Eng. Chem. Res. 2009, 48, 4–6.

(25)

Janković, B.; Adnađević, B.; Mentus, S. The Kinetic Study of Temperature-Programmed Reduction of Nickel Oxide in Hydrogen Atmosphere. Chem. Eng. Sci. 2008, 63, 567–575.

(26)

Zhou, Z.; Han, L.; Bollas, G. M. Kinetics of NiO Reduction by H2 and Ni Oxidation at Conditions Relevant to Chemical-Looping Combustion and Reforming. Int. J. Hydrogen Energy 2014, 39, 8535–8556.

(27)

Avrami, M. Kinetics of Phase Change. I General Theory. J. Chem. Phys. 1939, 7, 1103– 1112.

(28)

Avrami, M. Kinetics of Phase Change. II Transformation-Time Relations for Random Distribution of Nuclei. J. Chem. Phys. 1940, 8 (2), 212–224.

19 ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

(29)

Avrami, M. Granulation, Phase Change, and Microstructure Kinetics of Phase Change. III. J. Chem. Phys. 1941, 9, 177–184.

(30)

Sohn, H. Y.; Szekely, J. A. Structural Model for Gas—solid Reactions with a Moving boundary—IV. Langmuir—Hinshelwood Kinetics. Chemical Engineering Science, 1973, 28, 1169–1177.

(31)

Szekely, J.; Evans, J. W. A. Structural Model for Gas—solid Reactions with a Moving Boundary. Chem. Eng. Sci. 1970, 26, 1091–1107.

(32)

Furstenau, R.; McDougall, G.; Langell, M. Initial Stages of Hydrogen Reduction of NiO (100). Surf. Sci. 1985, 150, 55–79.

(33)

Utigard, T. A.; Wu, M.; Plascencia, G.; Marin, T. Reduction Kinetics of Goro Nickel Oxide Using Hydrogen. Chem. Eng. Sci. 2005, 60, 2061–2068.

(34)

Plascencia, G.; Utigard, T. The Reduction of Tokyo and Sinter 75 Nickel Oxides with Hydrogen. Chem. Eng. Sci. 2009, 64, 3879–3888.

(35)

Szekely, J.; Lin, C. I.; Sohn, H. Y. A. Structural Model for Gas-Solid Reactions with a Moving Boundary -V An Experimental Study of the Reduction of Porous Nickel-Oxide Pellets with Hydrogen. Chem. Eng. Sci. 1973, 28, 1975–1989.

(36)

Hidayat, T.; Rhamdhani, M. A.; Jak, E.; Hayes, P. C. Investigation of Nickel Product Structures Developed during the Gaseous Reduction of Solid Nickel Oxide. Metall. Mater. Trans. B 2009, 40, 462–473.

(37)

Medford, J. A.; Johnston-Peck, A. C.; Tracy, J. B. Nanostructural Transformations during the Reduction of Hollow and Porous Nickel Oxide Nanoparticles. Nanoscale 2013, 5, 155–159.

(38)

Chenna, S.; Crozier, P. A. In Situ Environmental Transmission Electron Microscopy to Determine Transformation Pathways in Supported Ni Nanoparticles. Micron 2012, 43, 1188–1194.

(39)

Nakamura, R.; Tokozakura, D.; Lee, J.-G.; Mori, H.; Nakajima, H. Shrinking of Hollow Cu2O and NiO Nanoparticles at High Temperatures. Acta Mater. 2008, 56, 5276–5284.

(40)

Kharatyan, S. L.; Chatilyan, H. A.; Arakelyan, L. H. Kinetics of Tungsten Carbidization under Non-Isothermal Conditions. Mater. Res. Bull. 2008, 43, 897–906.

(41)

Kharatyan, S. L.; Chatilyan, H. A.; Aghayan, M. A.; Rodriguez, M. A. Non-Isothermal Phenomena in Mo/Si Diffusion Couple: Reaction Kinetics and Structure Formation. Int. J. Self-Propagating High-Temperature Synth. 2013, 22, 18–26.

20 ACS Paragon Plus Environment

Page 20 of 22

Page 21 of 22

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

The Journal of Physical Chemistry

(42)

Khawam, A.; Flanagan, D. R. Solid-State Kinetic Models: Basics and Mathematical Fundamentals. J. Phys. Chem. B 2006, 110, 17315–17328.

(43)

Ostyn, K. M.; Carter, C. B. On the Reduction of Nickel Oxide. Surface science, 1982, 121, 360–374.

(44)

Buckett, M.; Marks, L. Electron Irradiation Damage in NiO. Surf. Sci. 1990, 232, 353– 366.

21 ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

Table of Contents (TOC) image

22 ACS Paragon Plus Environment

Page 22 of 22