Guest Editorial: The Importance of Interdisciplinary Science: When Chemistry Needs Physics

Home > Journal Archive > Guest Editorial: The Importance of Interdisciplinary Science: When Chemistry Needs Physics

Johnson Matthey Technol. Rev., 2020, 64, (2), 101

Introduction

Johnson Matthey has over 200 years of history, creating sustainable technologies, shaped around customers’ needs. Our ambition is to research, develop and innovate solutions to make the world cleaner and healthier, today and for future generations. Much of the underpinning science behind these technologies relies on a knowledge of chemistry and its application. Like most successful organisations, Johnson Matthey reflects on the scientific capabilities that are key to developing these solutions today but also looks to the future to plan which capabilities will be required to meet future challenges and opportunities. Much of this learning comes from external insight by looking at what is happening both within the markets and scientific disciplines we are familiar with but also in parallel disciplines. Today, our core scientific capabilities can be grouped into nine key areas covering catalysis, characterisation and modelling, chemical synthesis, materials design and engineering, electrochemistry, platinum group metal and specialist metallurgy, process optimisation, product formulation, surface chemistry and coatings. Pulling these together forms a powerful toolbox to develop solutions for our customers’ needs.

When looking beyond these capabilities, one useful ‘lens’ to look through is the overlap between scientific disciplines. For Johnson Matthey this might be to look at the interface between chemistry, one of our key underpinning strengths, and other sciences. For example, the interface between chemistry and physics; chemistry and biology or with other enablers such as the digital transformation that is enabling different ways of exploring science. Following this premise, Johnson Matthey Technology Review has devoted this issue to focus on physics and a future edition will look at biology.

Figure 1 shows how Johnson Matthey’s core science capabilities today may overlap with physics and biology. Further insights can then be drawn by mapping how we apply these capabilities to provide customer solutions into our existing markets (pink text) or where they may be aligned with global drivers and world challenges. Such techniques can be a valuable tool to help discuss and identify opportunities and needs for an organisation.

Fig. 1.

Johnson Matthey’s research and development (R&D) drivers and core capabilities

Johnson Matthey’s research and development (R&D) drivers and core capabilities

Two areas increasingly dependent on capabilities bridging chemistry and physics are characterisation and modelling of materials and processes and the development of functional surfaces and coatings. These topics feature heavily in this edition of the Johnson Matthey Technology Review. For Johnson Matthey, characterisation and modelling are key capabilities to help develop new technology.

Characterisation, Modelling, Coatings and Surfaces

Characterisation provides insights into composition, structure and property-performance relationships at all length scales. The latter includes in situ and operando analysis, which is important to understanding how materials may respond in their intended application.

Modelling also encompasses all length scales and includes statistical, empirical and physical models. Modelling has been used for a long time in chemical engineering to design reactors, systems and processes. Examples include designing a new reactor for a chemical reaction, an aftertreatment system for a vehicle or a process flow sheet for recycling waste materials. More recently, advances in modelling are permitting chemists to be more predictive, to be able to design materials, reactions and their performance with far fewer experiments. For example, in this edition, the need for computational modelling methods to replace incremental experimental development to meet the need to design complex new advanced materials is explained (1).

The application of nuclear magnetic resonance (NMR) to characterising activated carbons leads to insights into kinetic exchange of solvent molecules (2). The technique makes use of the magnetic shielding properties of the carbon structure to give insights into molecular level mechanisms which can give information to the chemist about where adsorbed species are in the material’s structure. These techniques enable the industrial chemist to gain a better understanding of the materials being used which leads to faster development and better understood technology. Equally important is the fundamental understanding of new materials and their properties both at the atomic and molecular scales which in time can lead to advances in existing or new technology.

Coatings and surface properties is another area at the interface of chemistry and physics. Johnson Matthey has many examples of products which rely on the functionality of particles deposited onto a surface. Examples include precious and base metal catalysts, advanced energy materials and medical components. As the coating thickness reduces from micron to atomic, the chemist’s traditional toolbox to deposit layers of formulated slurries, pastes and inks changes towards different deposition techniques such as chemical vapour deposition (CVD) and physical vapour deposition (PVD). The ability to design and deposit functional particles of a controlled size and shape onto a surface can find application in many disciplines such as transparent or reflective coatings, semiconductor devices, energy harvesting and sensing. Typically, these applications harness a combination of electronic, optical and chemical functionality. Further examples of applications in areas such as sensing, electronics and renewable energy are explored within this edition (3, 4).

Summary

Looking forward, global drivers such as climate change, the energy transition, population growth and longevity and resource challenges will drive the need for new technologies in areas such as more sustainable products, low carbon operations, clean energy and improved health and medical care. To meet these challenges chemists will increasingly need to reach out to adjacent disciplines to develop innovative solutions. In this edition of Johnson Matthey Technology Review, we welcome you to look at some of the advances in physics and explore how they are being used to drive forward R&D.

By |2020-03-31T10:24:42+00:00March 31st, 2020|Weld Engineering Services|Comments Off on Guest Editorial: The Importance of Interdisciplinary Science: When Chemistry Needs Physics

“Solid-State NMR in Zeolite Catalysis”

Johnson Matthey Technol. Rev., 2020, 64, (2), 197

Introduction

“Solid-State NMR in Zeolite Catalysis” was written by four professors from the Wuhan Institute of Physics and Mathematics, Chinese Academy of Sciences, who have tremendous expertise in the fields of solid-state nuclear magnetic resonance (solid-state NMR) and heterogeneous catalysis. This book is Volume 103 in the series ‘Lecture Notes in Chemistry’ published by Springer. It can be organised into three sections: (a) principles of NMR theory (Chapter 1); (b) solid-state NMR studies of zeolites and related materials (Chapters 2 and 3); and (c) surface chemistry and reaction chemistry over zeolites (Chapters 4–6). Each chapter offers a sufficient introduction to move smoothly into the main content. The book guides readers to familiarise themselves with solid‐state NMR and gives sufficient practical examples, detail and illustration. However, visually supported explanations by the authors were included only after Chapter 1. Adding to this could help deepen the understanding of pre-graduate readers. A few abbreviations were not explained, which could confuse readers. Advantages and disadvantages of solid-state NMR for each purpose are addressed in comparison with other techniques in each chapter. The discussion in this book further applies to materials for heterogeneous catalysts, glasses, ceramics, superconductors, lithium-ion batteries and biological systems.

Principles of Nuclear Magnetic Resonance

In Chapter 1, the authors provide a snapshot of NMR theory with the current state of methodology development, including strategies to enhance NMR sensitivity. This chapter gives sufficient background to understand the subsequent chapters without overwhelming.

Solid-State Nuclear Magnetic Resonance Studies of Zeolites and Related Materials

Highlights and lowlights of solid-state NMR for material characterisation are addressed at the beginning of Chapter 2. The chapter goes on to explain the local framework structure characterised by different chemical environments in zeolites. NMR spectroscopy provides detailed information on synthesis mechanisms in the solid and liquid phases, which may be complementary to X-ray diffraction and electron microscopy methods that are often used to gain structural information.

The application of NMR spectroscopy in zeolite crystallisation, in situ and ex situ methods is also introduced to investigate the kinetics of crystallisation for the convenience of a general audience. As an example, the structural changes in the aluminosilicate gels of zeolite faujasite (FAU) were confirmed in addition to the changes in silicon:aluminium ratios during the ageing and crystallising processes. In addition to zeolite analysis, xenon-129 chemical shifts were used to examine the pore geometry and chemical surroundings. This also fits the discussion of the guest-host interaction study. Additional visual aids could improve understanding in this chapter. The basics of characterisation of Brønsted acid sites may also be complemented by a review by Hunger (1).

Surface Chemistry and Reaction Chemistry Over Zeolites

NMR spectroscopy has the advantage of providing detailed information on local bonding and solid-state interactions. Investigations on host-guest interactions between organic molecules, such as aromatics to surfactants and inorganic moieties of zeolites, are well summarised in Chapter 4 with sufficient examples and explanation of local structures and dynamic behaviours between adsorbed molecules and active sites. Although several informative figures were cited from references, additional visual supports to explain the host-guest interaction would further elucidate this interaction. The content of Chapter 4 reintroduced and extended this reviewer’s understanding of Li‐ion battery and Li-ion ultracapacitor materials (2).

Chapters 5 and 6 elucidate interfacial chemistry by solid-state NMR approaches. Several representative reactions are given including in situ studies with a focus on how surface and interfacial phenomena significantly influence the catalytic performances (activity and selectivity) of heterogeneous catalysts. Many essential aspects of interfacial chemistry can be emphasised on the chemical and electrochemical reactions at interfaces. The following readings may bring additional value (3, 4), and this is also applicable for an understanding of Li-ion battery and Li-ion ultracapacitor materials (5).

Conclusions

The reviewer would like to recommend this book for anyone who would like to quickly grasp the capability of solid-state NMR as well as researchers working on nanoporous materials, nanocrystals, nanomaterials, Li-ion batteries and capacitors with surface and interfacial phenomena in addition to zeolites and zeolite related materials. This is exciting reading and reminded the reviewer of the ability of solid-state NMR to assess the mechanisms of intermolecular interactions, chemical reactions and transport phenomena. Additional visual supporting materials should be added to deepen the readers’ understanding of the authors’ explanations, particularly in Chapter 4.

“Solid-State NMR in Zeolite Catalysis”

“Solid-State NMR in Zeolite Catalysis”

The Author


Shingo Watanabe received his PhD from the Pennsylvania State University, USA, in the field of fuel, surface and materials chemistry. He worked for Johnson Matthey, USA, for eight years in heterogeneous catalysis and purification in petrochemicals, oleochemicals, biorenewable chemicals and gas and liquid purification. Currently, he is working for a start-up company, Nanoramic Laboratories, USA, as a business development director for the Japan market in the areas of Li-ion batteries, ultra- and supercapacitors and thermal interface materials.

By |2020-03-24T14:25:45+00:00March 24th, 2020|Weld Engineering Services|Comments Off on “Solid-State NMR in Zeolite Catalysis”

Manufacturing and Characterisation of Robot Assisted Microplasma Multilayer Coating of Titanium Implants

Home > Journal Archive > Manufacturing and Characterisation of Robot Assisted Microplasma Multilayer Coating of Titanium Implants

Johnson Matthey Technol. Rev., 2020, 64, (2), 180

1. Introduction

Manufacturing technology for medical implants undergoes constant improvement in order to accelerate the patient’s recovery and increase the service life of the implant. Following this trend, special attention is paid to applying biocompatible coatings on the surface of implants (17). There is a huge clinical need for advanced biomaterials with enhanced functionality to improve the quality of life of patients and reduce the burden of health care for the world’s ageing population. In recent years, Ti and HA have been widely used in medical devices due to their favourable biocompatibility (18).

The research presented here offers robot assisted MPS of Ti wires and HA powder onto Ti substrates (912). Currently, MPS technology is highly specialised with few research groups working on its development. Yet it has shown promising results in delivery of high quality and well-designed coatings (13). Among existing plasma spraying processes, MPS is particularly characterised by low plasma power (up to 4 kW), small diameter of the spray deposition spot on the surface (up to 15 mm) and the possibility of forming a laminar jet. This laminar jet can be up to 150 mm in length. It heats the refractory material in a stream of argon plasma and provides low heat input into the substrate (2). The process provides deposition of Ti or HA on small size parts and components, including those with complex geometry. This is normally unachievable with any other method. The MPS generally provides a micro-rough surface and a higher degree of porosity (~20%) that in the case of biocompatible coatings facilitates bony tissues ingrowth; in most cases the bond strength of MPS coatings with substrates is acceptable (2). However, there are still a number of challenges remaining to be addressed. The most important issue is the formation of coatings with specified structure and properties.

Currently, endoprosthesis medical practice widely uses metal implants coated with HA (18, 13, 14). HA is the calcium phosphate mineral Ca10(PO4)6(OH)2 of the apatite group, which is chemically similar to the apatite of the host bone and is a source of Ca and P for the bone-HA interface (3, 4, 6, 8). HA coatings improve osseointegration and can significantly reduce the duration of implantation of the endoprosthesis. It provides a reliable connection with the bone and increases the reliability of the implant (18, 13, 14). In the case of thermal spraying of HA powder, the chemical composition of the final HA coating is dependent on the thermal decomposition occurring during spraying. The high temperatures experienced by HA powder particles in the plasma spraying process lead to the dihydroxylation and decomposition of the particles. At temperatures above 1050°C HA decomposes to tricalcium phosphate (β-TCP, Ca3(PO4)2) and tetracalcium phosphate (TTCP, Ca4(PO4)2O), and above 1120°C β-TCP is converted to α-tricalcium phosphate (α-TCP, Ca3(PO4)2) (13, 14). Thus, the resultant coating phase composition depends on the thermal history of the powder particles. The higher the plasma jet temperature and the longer the exposure of the particles to plasma, the greater the degree of phase transformation. According to the British Standards Institution standard specification BS ISO 13779-2:2000 (15), the maximum allowable content of non-HA phases in a HA coating is 5%, and the minimum allowable percentage of crystallinity is 50%. The degree of crystallinity of the HA coating affects the process of osseointegration (16). Amorphous calcium phosphate (ACP) has a higher rate of dissolution, which reduces the recovery time of the patient, but at the same time some reduction of the reliability of fixation of the endoprosthesis in the bone is also possible. Thus, increased crystallinity provides reliable fixation of the implant in the bone (2, 3).

It is assumed that to increase the biocompatibility of the implants and rapid accretion with the bone, the implant surface should be covered with a biocompatible coating with an extensive surface morphology, with recommended pore sizes in the coating from 20 μm to 200 μm, and closed porosity of at least 30% (24, 16, 17). At the same time, the coating must be firmly connected to the implant, without transversal porosity, so that the implant material does not interact directly with the human body. Therefore in this research, the proposed composition and thickness of the coating is designed to meet these requirements. In order to characterise the coated layer structures, SEM, TEM and XRD are used. Optimum modes of plasma spraying are chosen on the basis of the results of structural characterisation.

The aim of this study is to select the modes of MPS of both Ti wire and HA powder to obtain a sufficiently thick (up to 300 μm) multilayer Ti/HA coating. The dense Ti sublayer is supposed to provide good adhesion to the substrate and the porous Ti middle layer and HA top layer can accelerate the bone ingrowth. Another objective is to clarify the relationship between various plasma spray process parameters and the resultant coating structure. This is achieved through the development of process models that relate process parameters to various coating structures.

2. Materials and Methods

The HA was synthesised in the laboratories of East Kazakhstan State Technical University. The process of synthesis of HA powder with the ratio Ca:P of 1.65 by a chemical precipitation method was described in our previous paper (10). The purity of HA powder was 99.5%, which met the purity requirement (not less than 95%) set out in the ASTM International (USA) standard ASTM F1185-03(2014) (18). The purity of the synthesised HA is determined using the XRD results which are further explained below. Before spraying, the powder was dried at a temperature of 120°C for 6 h, in order to avoid clogging during powder feeding. HA powder was used as a sprayed coating material. The particles had irregular shape with smooth edges (Figure 1).

Fig. 1.

SEM images of HA particles indicating particle size

SEM images of HA particles indicating particle size

The shapes of the particles are particularly designed with smooth edges to make sure that they do not stick together. Moreover the HA particles should not cling to each other in order to provide the necessary flowability of the powder. The flowability is essential in gradual supply of the powder to the plasma jet.

According to previous studies (9, 10), the particle size of the HA powder for MPS should be in the range of 40 μm to 90 μm. Before MPS, the dried and milled HA powder was sieved through mesh diameters 40 μm and 90 μm to obtain only the powder fraction within the desired range. The flow rate of the HA powder is in the range of 120 s 50 g–1.

Samples of medical Ti alloy of Grade 5 extra low interstitials (ELI) (Table I) were used as substrates for MPS. For deposition of Ti coatings, wires of VT1‐00 commercially pure Ti with a diameter of 0.3 mm were used (Table I). The samples of medical Ti alloy of Grade 5 ELI (Table I) were cut with thicknesses of 3 mm from rods with a diameter of 50 mm and 30 mm on CTX 510 ecoline computer numerical control (CNC) machine (DMG MORI AG, Germany). Plates of size 15 mm × 15 mm × 2 mm were cut from large sheet of Grade 5 ELI alloy.

Table I

Chemical Composition of Ti-Based Materials

Materials grade Reference composition wt% of element
Al V O C N H Fe Si Ti
Grade 5 ELI alloy (19), (20) 5.50–6.75 3.50–4.50 0.13–0.20 0.08 0.05 0.015 0.25–0.40 base
VT1-00 commercially pure Ti (21) 0.3 0.12 0.05 0.03 0.003 0.15 0.08 base

MPS of the HA powders and Ti wires was carried out by MPN-004 microplasmatron (produced by E.O.Paton Institute of Electric Welding, Kiev, Ukraine) (22). The microplasmatron was mounted on an industrial robot arm (RS010L, Kawasaki Heavy Industries, Japan). It is able to move horizontally along a computed trajectory at set speed. The thickness of the coatings was varied from 80 μm to 300 μm by changing the MPS parameters. The speed of linear movement of the plasmatron along the substrate was chosen to be 50 mm s–1. The choice of speed of the plasmatron was based on preliminary estimates of the temperature of the substrate when exposed to a plasma jet (12). This was to make sure the temperatures remain well below the melting temperature. Ar served as a plasma-forming and transporting gas for MPS; additional heating of the substrate was not carried out. Using a robotic arm allows precise spraying of coatings with a uniform speed of movement of the plasmatron along the surface of the implant, as well as moving the plasmatron along a predetermined path.

Before MPS, the surfaces of the samples were degreased with acetone and subjected to ultrasonic cleaning. To ensure proper adhesion of the coatings, it was important to pre-treat the surfaces of the substrates to increase their roughness. For surface activation, gas abrasive surface treatment was carried out on a Contracor® ECO abrasive blasting machine (Comprag Group GmbH, Germany) using normal grade A14 electrocorundum. The chemical reactivity of the substrate’s surface rapidly falls due to oxidation and the adsorption of chemical gases from the atmosphere. Therefore it is important that the time interval between the gas abrasive treatment and the coating on the surface does not exceed 2 h. Before coating, samples were stored in a tightly closed container.

To evaluate the porosity of biocompatible coatings, the images obtained by the scanning electron microscope JSM-6390LV (JEOL Ltd, Japan) were processed using MicroCapture (MustCam, Hong Kong) and ATLAS.ti (ATLAS.ti Scientific Software Development GmbH) computer-aided programs. The measurements were carried out on the polished cross-section of the coatings according to ASTM E2109-01(2014) standard (23). The surface roughness of the substrates and the as-sprayed coatings was measured in accordance with ISO 4287:1997 (24) using a MarSurf PS 10 mobile roughness measuring instrument (Mahr, Germany). Four measurements were taken for each sample and the average was determined. The adhesion strength of coatings to the substrates was measured in tension using a AG-X universal testing machine (Shimadzu, Japan) in a static experiment according to ASTM C633-13(2017) standard (25).

The crystallinity and the structure-phase compositions of HA powders and plasma sprayed HA coatings were measured by X’Pert PRO diffractometer (PANalytical, The Netherlands). The interpretation of the XRD patterns was carried out using Rietveld method and powder diffraction data from the database of the International Centre for Diffraction Data (ICDD, USA, 2003), the ASTM card file and X’Pert HighScore Plus software (Malvern Panalytical, UK). The percentage of crystallinity of the HA powder was calculated using the area of crystalline peaks in the region of 15° to 45° 2θ and the area of the amorphous diffuse background in this region. Diffraction scans of the HA powder and coatings were carried out in accordance with ASTM F2024-10(2016) (26). The purity of HA powder and HA coatings was evaluated by calculating the areas of all non‐HA peaks found in the diffraction pattern. The impurity area was determined by calculating the area in the region where the highest impurity phase peaks were present. The impurity peaks that would be expected to be present in HA powders and HA coatings were those of TTCP, α-TCP and β-TCP. The Rietveld method was used to quantitatively determine the percentage of various phases of impurity in HA coatings. X’Pert HighScore Plus software was used to calculate the impurity. Three purity measurements were carried out for each of the XRD patterns.

Electron diffraction patterns of samples of HA powder were obtained by TEM on JEM-2100 (JEOL). The structural-phase composition data obtained using TEM were compared with the data obtained using XRD analysis. TEM sample preparation techniques are described in detail in a previous publication (11).

The particles of the sprayed Ti wire after collision with the substrate were studied using splat tests (2). MPS of Ti wire onto the plates of polished Ti alloy was performed in the plane perpendicular to the axis of the plasma jet. Speed of the linear movement of plasma jet was set at 50 mm s–1. As a result, single particles of the sprayed material (splats) were fixed on the substrate and deformed upon contact with the substrate surface. The splats’ visual analysis was carried out by SEM; the splats were classified according to their appearance and their spraying modes (Table II).

Table II

MPS Deposition Parameters of Ti Wires, the Porosity and the Arithmetic Ra of Sprayed Coatings

Spraying mode Parameters/settings Porosity of sprayed Ti coatings, % Ra, μm
I, A Vpg, slpm H, mm Vw, m min–1
1 25 3.7 120 4.3 6.2 ± 0.74 12.0 ± 0.97
2 25 3.7 40 3.0 13.8 ± 1.90 45.6 ± 4.40
3 25 2.3 120 3.0 12.0 ± 1.50 44.4 ± 3.85
4 25 2.3 40 4,3 5.7 ± 0.45 12.0 ± 1.13
5 15 3.7 120 3.0 9.6 ± 0.35 31.5 ± 3.10
6 15 3.7 40 4.3 10.6 ± 1.32 34.8 ± 3.23
7 15 2.3 120 4.3 8.7 ± 2.32 30.2 ± 3.87
8 15 2.3 40 3.0 31.0 ± 3.87 >50

3. Results and Discussion

The effect of parameters of MPS such as electric arc current (I, A), plasma gas flow rate in standard litre per minute (Vpg, slpm), spraying distance (H, mm) and wire flow rate (Vw, m min–1) or powder consumption (Ppowder, g min–1) on the surface morphology, porosity and structural‐phase transformations of the coatings have been studied. The coating experiments for MPS were accomplished in a two level fractional factorial design (24–1). The experimental conditions in fractional factorial designs have been selected to provide balanced design (27). The maximum and minimum values of the parameters for feasible processing of high-quality coatings were chosen empirically.

MPS of Ti wires on Ti alloy gas-abrasive treated substrates was carried out in various modes as shown in Table II. The key characteristics of the resultant coatings such as their porosity and roughness were examined. The data in Table II and Table III represent the averaged values for three experimental runs.

Table III

The Dependency of the CTE, Phase Composition and Crystallinity of HA Coating on the Spraying Parameters

Mode I, A Vpg, slpm H, mm Ppowder, g min–1 CTE, % CTE, % estimated HAcryst. Aph β-TCP
1 45 2.0 160 1.2 54 58 92 5 3
2 45 2.0 80 0.4 64 71 98 0 2
3 45 1.0 160 0.4 89 89 93 2 5
4 45 1.0 80 1.2 69 69 96 0 4
5 35 2.0 160 0.4 29 29 93 4 3
6 35 2.0 80 1.2 48 48 94 3 3
7 35 1.0 160 1.2 40 47 88 7 5
8 35 1.0 80 0.4 56 60 98 0 2
9 40 1.5 120 0.8 60 59 90 6 4

Examination of the Ti particles splats obtained in Modes 1, 4 and 8 (Table II) showed that the samples are completely melted and have formed a disk (Figures 2(a) and 2(b)). The thicker splats are formed in Mode 8 (Table II, Figure 2(c)). The beginning of the process of solidification of the particles upon impact with the substrate is shown in Figure 2(c). The increase in the thickness of the splats in Mode 8 (Table II) is due to the minimum velocity of the particles during their interaction with the substrate.

Fig. 2.

SEM images of Ti splats sprayed in: (a) Mode 1; (b) Mode 4; and (c) Mode 8 (according to Table II)

SEM images of Ti splats sprayed in: (a) Mode 1; (b) Mode 4; and (c) Mode 8 (according to Table II)

As can be seen in Figure 2, the splats obtained in different modes have similar area but different thickness. The increase in the thickness of the splats in Mode 8 (Table II) is due to the minimum velocity of the particles during their interaction with the substrate. The minimum particle velocity in Mode 8 is caused by the large size of the sprayed particles, low gas flow rate and low spraying distance. This also leads to a decrease in the speed of the sprayed particles. Particles of Ti wire melted by a plasma jet move towards the substrate during plasma spraying. On the one hand, the size of these particles depends on the set spraying parameters. On the other hand, both the set spraying parameters and the size of the particles affect the speed and the degree of heating of the particles in the plasma jet. Before interacting with the substrate, the particles of molten metal can either heat up, being in the high-temperature zone of the plasma jet (the initial section of the plasma), or, more likely, cool down by going to the low-temperature zone (the end of the plasma jet). The size of the sprayed particles and the degree of their melting in the plasma jet can be varied by spraying parameters. A porous coating with a high surface roughness can be achieved by large particles moving with low speed, as in Mode 8. A dense coating with a relatively low surface roughness can be obtained by high speed small and completely melted particles, as in Mode 1 or in Mode 4 (Table II).

The appearance and structure of the splats varies depending on the interactions of the sprayed particles with the substrate. The interaction is generally defined by the velocity of the particles on impact and also their degree of melting. It is possible to correlate the temperature and velocity of the particles before the collision to the substrate with the resultant structure of the coating. The profile and cracks on the surface of the splat correlate with the stress state of the coating (2, 28).

The substrate average roughness (Ra) after gas abrasive treatment was Ra = 7.0 ± 0.35 μm. The analysis of the surface morphology of Ti coatings showed the possibility of obtaining dense coatings with relatively low roughness (Ra is about 12.0 μm) using Mode 1 and Mode 4 (Table II), the remaining modes provide a high surface roughness (above 30 μm). The maximum size of open pores (up to 300 μm) is observed on the surface of coatings obtained in Mode 8 (Figure 3).

Fig. 3.

SEM images of surfaces and cross-sections of Ti coatings sprayed in Mode 4 and Mode 8 (according to Table II)

SEM images of surfaces and cross-sections of Ti coatings sprayed in Mode 4 and Mode 8 (according to Table II)

The analysis of the cross-sections of the Ti coatings showed that the coatings sprayed in Mode 8 have the highest average porosity of 31.0%, while Mode 4 makes it possible to obtain dense coatings with an average porosity of 5.7% (Figure 3 and Table II). The tensile strength test established the average adhesion strength of the coating with a thickness of 100 μm sprayed in Mode 4 (Table II) to be 38.7 MPa. This meets the requirements of ISO 13179‐1:2014 (29). According to ISO 13179‐1:2014 (29), the average static tensile strength of a Ti coating should be more than 22 MPa.

As can be seen from the results presented in Table II, the arithmetic Ra of the coatings surface correlates with the porosity of the sprayed coatings. The condition of the presence of large, incompletely molten particles with a low speed in the plasma jet leads to a high surface roughness and increased porosity of the coatings. Thus, the roughness of the coating is affected by the size of the particles involved in the formation of the coating. The main purpose of plasma spraying of the Ti layers at the initial stage of the coating is to form a porous surface with high roughness in order to spray fully molten HA particles onto it. It was previously shown (30) that partially melted HA particles were not able to form dense coatings. This leads to poor adhesion of the HA coating to the substrate. A desirable HA coating shows good adhesion and high surface roughness. The surface roughness of the HA coating affects osteoblast cell attachment and thus accelerates bone growth into the implant. Whereas fibroblasts and epithelial cells prefer smoother surfaces, osteoblasts attach and proliferate better on rough surfaces (31, 32). A relatively thin HA layer (not more than 100 μm) sprayed on a porous surface with a high roughness (above 50 μm) can ensure the reliable adhesion of HA coating to the surface, while maintaining high surface roughness and porosity. XRD analysis confirmed that the phase composition of the initial HA powder was fully crystalline Ca10(PO4)6(OH)2.

The TEM images (Figures 4(a) and 4(b)) of the HA powder particle and the corresponding microelectron diffraction pattern (Figure 4(c)) are shown in Figure 4. The results of TEM analysis are in good agreement with the results of XRD analysis (Figure 4(d)). X-ray phase analysis showed that the main phase (99.5%) is HA with the hexagonal crystal system P63/m. The electron diffraction pattern corresponds to the hexagonal phase of HA with unit cell parameters a = 0.94 nm, c = 0.68 nm.

The plasma spraying of HA powder was carried out using nine different modes (Table III). The key criteria were the phase composition, the degree of crystallinity (the proportion of the amorphous phase (Aph), the proportion of the crystalline phase (HAcryst)) and the coating transfer efficiency (CTE). In powder coating, transfer efficiency is the ratio of the quantity of powder deposited on the part to the quantity of powder directed at the part. Transfer efficiency is provided as a percentage, with 100% being most desirable. An experiment was conducted to determine the impact of process parameters such as I, Vpg, H and Ppowder on the CTE.

Fig. 4.

(a, b) TEM images of the HA powder particles; (b, c) the corresponding indexed microelectron diffraction pattern; and (d) the XRD patterns of HA powder

(a, b) TEM images of the HA powder particles; (b, c) the corresponding indexed microelectron diffraction pattern; and (d) the XRD patterns of HA powder

The application of well-known methods of fractional factorial design (27) and the design of experiment method (33) for the analysis of plasma spraying of HA powder are described elsewhere (34, 35) and in our previous paper (10). A factorial experimental design to investigate the relationship between plasma spray parameters and the microstructure of HA coatings was first used by Dyshlovenko et al. (35). Three responses were examined (35). This included the fraction of HA, the fraction of decomposition phases and the amorphous content of the coatings. In this study, the next three responses were examined: the fraction of HA, the amorphous content of the coatings and CTE. The first two responses were chosen to ensure purity and crystallinity of the HA coating. CTE was selected in order to increase the efficiency of the plasma spraying process. Moreover it allows interpretation by a linear regression model which could quite easily and reliably be measured in experimental runs. The linear regression model was chosen to estimate CTE. The coefficients in the regression Equation (i) were calculated by assigning the corresponding units of measure: 2.575 A–1; –0.246 slpm–1; –0.203 mm–1; 4.06 min g–1; –0.825.

(i)

The comparison of calculated and experimental results indicates good agreement (Table III). Therefore, Equation (i) can be used for preliminary estimation of CTE when selecting MPS modes. A more complex regression model that takes into account the mutual influence of factors should be applied to determine the dependency of other values presented in Table II and Table III on the spraying parameters. This requires further investigation using the obtained experimental data.

The results of XRD analysis presented in Table III show that the phase compositions of all coatings comply with ISO 13779-2:2000 (15). However, Mode 3 provides the highest CTE. Thus, we consider Mode 3 to be the most cost-effective. This mode allows a desired HA coating thickness (about 100 μm) to be obtained in one pass of a plasma jet.

The modes of application of multilayer coatings of Ti/HA were selected on the basis of the analysis of the tests of MPS of Ti wire and HA powders with the measurement of thickness and porosity of the deposited layers. The porosity and surface roughness of Ti coatings and the high CTE value, purity and crystallinity of the HA coating indicate the optimum coating composition and MPS modes. The coating thickness and the expected adhesion to the substrate are other parameters indicating the quality of the coating. The first relatively thin (up to 80 μm) and dense layer of Ti wire coating is applied in Mode 4 (Table II), then another layer of Ti wire up to 100 μm thick is sprayed in Mode 8 (Table II) to form a porous coating with a rough surface, then the top layer of HA powder coating is applied with thickness up to 100 μm in Mode 3 (Table III). The microstructure of microplasma sprayed multilayer Ti/HÀ coatings under the above modes is shown in Figures 5(a) and 5(b). The XRD pattern for microplasma-sprayed HA coating is presented in Figure 5(c).

Fig. 5.

SEM images of the multilayer coating with the Ti lower and middle layers and the HA upper layer: (a) the surface; (b) the cross-section; and (c) the XRD patterns of microplasma-sprayed HA coating

SEM images of the multilayer coating with the Ti lower and middle layers and the HA upper layer: (a) the surface; (b) the cross-section; and (c) the XRD patterns of microplasma-sprayed HA coating

The desired porosity was achieved in the Ti lower layer (30 vol%) (Figure 5(b)). Pore sizes in both the Ti middle layer and the HA top layer are in the range 20–50 μm (Figure 5(b)). HA coating porosity is about 20% (Figure 5(b)). It should be noted that for biocompatible coatings, open porosity is essential: that is the egress of the pore on the surface of the coating, where the bone grows. Therefore, measuring the diameters of the pore craters on the surface of the coating is an appropriate way to indicate the surface morphology and profile. In our experiment, the maximum pore diameter on the surface of the HA coating was about 150 μm (Figure 5(a)).

It was established by XRD that the mode specified for HA powder provides the required structure-phase composition in the HA coating: 93% by weight of the HAcryst, 5% by weight of β-TCP phase, and 2% by weight of the Aph. The coating purity was determined using the procedure outlined in Materials and Methods section above. The highest peaks of HA and β phases for HA coatings were located at 32° 2θ and 31° 2θ respectively. The purity of coatings was found to be 95.1%. This shows that the purity meets the 95% purity requirements of ISO 13779-2:2000 (15). The measurement error was 0.06.

The loci of Aph have been found on the XRD patterns between 18° and 38° 2θ (Figure 5(c)). All the diffraction patterns in the range of 37.3° 2θ were thoroughly investigated, but even weak peaks of calcium oxide (CaO) were not found (Figure 5(c)). It confirms that no harmful CaO compound is formed through MPS coating of HA powder.

For a multilayer coating, the porosity and adhesion of the top HA layer depends on the characteristics of the lower Ti layers such as the roughness and open porosity of the middle Ti layer. To determine the dependency of the porosity and adhesion of the HA coating on the parameters of MPS, further research is needed.

This study proves that it is possible to obtain coatings from biocompatible materials with the desired level of porosity and satisfactory adhesion to the substrate using MPS. A robot assisted MPS of coatings from biocompatible materials of Ti and HA onto Ti implants has been implemented. Also the composition and modes of microplasma deposition of multilayer coatings for Ti implants have been identified. The next stage of the research includes the study of the biocompatibility of microplasma-sprayed coatings (in vitro tests) and MPS of different materials such as tantalum and zirconium.

Among the number of works that have demonstrated the advantages of thermal spraying of biocompatible coatings for use in medical applications, three recent papers (3638) have shown promising directions for further development of the research presented here. Cizek et al. (36) have reviewed the patents concerning thermal spraying for biomedical applications for the period 2005 to 2018. They have also reported recent research and development trends in this field. Among the materials recommended for bio-applications, they have mentioned Ta. It can be noted that MPS of Ta wire onto Ti alloys using the technology presented in our paper is highly feasible. Our trials with Ta have indicated great potential. This could open the potential to apply the developed technology for other materials. Fotovvati et al. (37) have compared the results of obtaining biocompatible coatings by cold and thermal spraying in favour of thermal spraying. Fousova et al. (38) have shown the benefits of using thermal plasma spray to prepare bulk Ti for bone enlargements. However, despite the advantages and relative cost effectiveness of thermal plasma spraying, its use for the manufacture of medical implants has not yet become widespread. This is mainly due to the high temperatures of the bulk resulting from the thermal spraying process. MPS avoids the issue of overheating. It allows coatings to be obtained from materials with a high melting point, such as Ti and Ta, by a microplasma jet while introducing a very small thermal impact into the substrate.

The use of robotic MPS could be considered promising for the production of patient specific implants. Three-dimensional scanning and rapid prototyping technologies facilitate the manufacture of specifically designed complex geometry implants and robot assisted plasma coating is used for coating. This is more advantageous for the production of small endoprostheses with biocompatible coatings, such as vertebral cages and dental implants (39, 40).

4. Conclusion

It has been established that the main parameters controlling the porosity of microplasma sprayed coatings are I and Vpg. MPS parameters for the formation of porous coatings of Ti wire and HA powders with rough surfaces have been determined and are reported here. The advantages of applying SEM, TEM and XRD to analyse the structure of sprayed Ti and HA coatings to substantiate the choice of plasma spraying modes of the coatings were demonstrated. It is also proven that by using the appropriate MPS process parameters, a layer of HA with a high degree of crystallinity (93%) can be obtained, controlled by changing the deposition mode. The small size of the spraying spot (up to 8 mm) provides a significant reduction in Ppowder when depositing on implants of small size compared to conventional plasma spraying.

The composition and modes of microplasma deposition of multilayer coatings for Ti implants, including a dense Ti sublayer, porous Ti middle layer and HA top layer have been established. The total thickness of such coatings is about 300 μm. The porous middle coating layer has a porosity of 30% and a pore size varying from 20 μm to 50 μm. Moreover the upper layer of HA indicates a thickness up to 100 μm with a 95% level of HA phases and 93% crystallinity. The results of this research are of significance for a wide range of researchers developing plasma spray technologies for biocompatible coatings manufacture.

The Authors


Darya Alontseva completed her PhD in Physics at East Kazakhstan State University in 2002. She completed her postdoctoral studies at Altai State Technical University, Russia, in 2013 and received the degree of Doctor of Sciences in Physics and Mathematics. In 2016 she was awarded the academic title of Full Professor of Physics. She has 19 years of research experience in developing new material and processes and management of funded scientific projects. She is a lead researcher in her research area: physics of condensed state, material science and surface engineering.


Elaheh Ghassemieh has 25 years of research experience in the areas of advanced manufacturing including additive manufacturing, development of novel materials especially composites using range of experimental and multiscale numerical methods. After completing her PhD in simulation of micromechanics of composite materials, she worked at several universities in the UK including The University of Sheffield, Queen’s University Belfast and Loughborough University. She has obtained a number of research funding grants from UK and international research councils and has also secured and managed many industrially funded projects where she has transferred novel research outcomes to relevant industrial users in the aerospace, automotive or biomedical fields.


Sergey Voinarovych completed his postgraduate study at E.O.Paton Electric Welding Institute, Ukraine, in 2008 and received a PhD degree in the specialty “Welding and related processes and technologies”. Since 2010 he has been working as a Senior Researcher at E.O.Paton Electric Welding Institute. He has 20 years of research experience in developing new processes, materials and equipment in the area of thermal spray coatings. He has designed MPS equipment and technology for forming biocompatible coatings on parts of endoprostheses. Currently he is researching new biocompatible materials and coatings.


Oleksandr Kyslytsia completed his postgraduate study at E.O.Paton Electric Welding Institute in 2010 and received PhD degree in the specialty “Welding and related processes and technologies”. Since 2011 he has been working as a Senior Researcher at E.O.Paton Electric Welding Institute. For over 20 years he has been developing new processes, materials and equipment for producing coatings by gas thermal spraying methods. He developed equipment and technology for MPS from wire materials to obtain biocompatible coatings on parts of endoprostheses. Currently he is researching new biocompatible materials and coatings.


Yuri Polovetskyi graduated from the Chernihiv Technological Institute, Ukraine, in 1999 with a degree in Welding Technology and Equipment and entered the doctoral program at the E.O.Paton Electric Welding Institute. After completing his doctorate in 2002 and to the present, Polovetskyi is a senior researcher of the department of physical and chemical research of materials. His area of scientific expertise is the structural characteristics, chemical composition and mechanical properties of welded joints and plasma coatings for various purposes.


Nadezhda Prokhorenkova received her PhD in Technical Physics in 2014. She has nine years of research experience in developing new materials and processes. She is an accomplished researcher in her research area: material science. Her area of scientific expertise is X-ray analysis of coating structures. Currently Prokhorenkova is an associate professor at the School of Engineering at D. Serikbayev East Kazakhstan State Technical University.


Albina Kadyroldina received her BS and MS degrees from D. Serikbayev East Kazakhstan State Technical University. She is currently pursuing a PhD with the School of Engineering, D. Serikbayev East Kazakhstan State Technical University. Her research interests include automation, control and mathematical modelling.

By |2020-03-24T09:21:16+00:00March 24th, 2020|Weld Engineering Services|Comments Off on Manufacturing and Characterisation of Robot Assisted Microplasma Multilayer Coating of Titanium Implants

EuropaCat 2019

Johnson Matthey Technol. Rev., 2020, 64, (2), 192

1. Introduction

The 14th European Congress on Catalysis (EuropaCat 2019), themed ‘Catalysis without Borders’, was held on the 18th–23rd August 2019 at the Eurogress conference centre in Aachen, Germany. The conference hosted over 1500 participants from academia and industry across the world, with around 400 lectures and 800 posters presented throughout the week. There were six parallel sessions; this review is a small selection of the talks attended by the reviewers.

2. Plenary Lectures

Plenary lectures of an hour were intended to showcase the significant breadth and depth of scientific knowledge acquired by established academics over the course of their careers. Summaries of key themes presented in the lectures are provided below.

Two plenary lectures were delivered on the subject of using molecular oxygen as a reagent, one by Professor Karen Goldberg of the University of Pennsylvania, USA and another by Professor Shannon Stahl from the University of Wisconsin-Madison, USA. Professor Goldberg opened by highlighting how molecular O2 represents the ideal oxidant for chemical transformations, with the ‘jackpot’ being to achieve benign and facile routes to chemicals with atmospheric O2, avoiding O2 separation from air. Professor Stahl provided examples of work studying the control of molecular O2 reactivity including maximising the energy efficiency of O2 reduction to water in fuel cells and achieving selective oxidation of organic molecules without overoxidation to carbon dioxide or other undesirable byproducts (13).

A high-impact plenary was presented by Professor Ib Chorkendorff from the Technical University of Denmark. The talk started with a few headline data to emphasise that sectors such as aviation, long haul transport and the chemical industry will always have to rely on the same chemical building blocks we use today rather than full electrification. The speaker predicted that the chemical demands of the future should be met by ‘solar’ fuels in the form of photovoltaics providing the energy for electrochemistry. Several areas of focus within the research team were summarised including a collaboration with Haldor Topsøe, Denmark, in electrified steam reforming (4), green hydrogen generation (5) and electrochemical hydrogenation of CO2 (6).

Professor Valentin Parmon of the Boreskov Institute of Catalysis SB RAS in Novosibirsk, Russia, presented an overview of some non-traditional approaches to achieve various endothermic thermocatalytic transformations. An example was the application of nuclear radiation to supply heat for steam reforming in an integrated chemical-nuclear reactor. A steam reforming catalyst was developed for the system; a porous uranium oxide support with nickel on the surface and titanium dioxide to protect against nucleotides (7). Solar energy was also considered to provide heating for steam reforming, temperatures >1000 K can be achieved by concentrating solar radiation (8). More recent work in the Parmon group focused on the use of microwaves to supply energy to reactions such as H2 evolution from alkanes though the cost of the electrical energy required is prohibitive for large scale application (9).

A range of fundamental structure-activity relationships were described in the lecture of Professor Christophe Copéret from ETH Zürich, Switzerland. Work in the group is focused on the understanding and control of chemistry on surfaces, with the ultimate goal to generate isolated metal sites with defined chemical environment to elucidate structure-activity relationships. The approach in the team has led to highly active and selective single-site catalysts that out-performed their homogeneous counterparts, but that also provided useful information to understand industrial catalysts. Several examples were provided in the lecture including details of research into chromium polymerisation catalysts to generate a Cr(II) species directly rather than the need to reduce Cr(VI) in situ (10). An example was also provided for copper-particle catalysed methanol synthesis from methane, investigating the effect of the support; copper on alumina displayed higher activity than copper on silica (11). Figure 1 summarises work undertaken in the Copéret team and acknowledges the ‘pioneering’ work of Denis Ballard, Imperial Chemical Industries (ICI), UK.

Fig. 1.

Molecular processes investigated by the Copéret group reproduced with kind permission from Professor Copéret (12, 13)

Molecular processes investigated by the Copéret group reproduced with kind permission from Professor Copéret (12, 13)

3. Industrial Forum

The industrial forum sessions of the conference are to present work undertaken by industrial corporations or of direct industrial application. Over the course of the conference there were 21 talks in the industrial forum by both academics and speakers from commercial organisations.

3.1 Growth of Chemicals

Joe Scheper from ExxonMobil Chemical Europe Inc, Belgium, opened the industrial forum with a talk focussed on the projected growth in global chemical demand, set to outpace gross domestic product (GDP) growth. Attention was drawn to <15% of a barrel of oil being converted to chemicals pre-1990, versus 25–40% at present, and projected to be 40–80% post 2020. Specific examples were given of technology to enable this growth, including the ExxonMobil EMTAMSM process for selective para- xylene production from toluene and methanol. Para- xylene is smaller than the meta- or ortho- variants, so can be targeted for selective transport through specifically sized zeolite pores. Avoiding reaction on the outer surfaces of zeolites minimises byproduct formation and catalyst attrition resistance was highlighted as critical for this fluid bed regenerative process. Some future developments required to sustain the necessary global growth in chemicals were presented. Processes to take plastic waste and turn this into a chemical feedstock as well as algae-based routes to biofuels with 50,000 barrels (bbl) of petroleum products being projected to be produced from algae in the not too distant future.

Andrei Parvulescu from BASF SE, Germany, then went on to discuss growing demand for highly branched C7 and C8 alkanes owing to their superior fuel burning efficiency. C4 alkene dimerisation or alkylation of isobutane with propene or butenes were presented as routes (Scheme I).

Scheme I

 

Solid acid catalysts have been proposed for both of these processes but suffer from relatively poor stability and selectivity. One successfully evaluated catalyst for alkylation was a zeolite beta made via seed-assisted synthesis without organic structure directing agent, providing a high density of acid sites and easy ability to be modified by dealumination. The catalyst was shown to have superior activity and selectivity compared to conventional zeolite beta and control of product mix was demonstrated through manipulation of reaction temperature and olefin concentration (14).

Jean-Paul Lange from Shell Global Solutions International BV, The Netherlands, continued the theme of growth in chemicals, discussing how as populations become wealthier demand for petrochemicals will rise. Fossil fuels would be able to meet this increase in demand, as they are still abundant and will have lower demand in future as the automotive sector decarbonises. Shell analysis predicts that the petrochemical industry will be affected by limited carrying capacity of the earth and societal calls for a more circular economy. Therefore, focus must be placed on transforming the linear petrochemicals industry into a circular one built on recycled plastics. There was also agreement between Shell and the earlier discussions led by ExxonMobil, with Shell predicting that as transport electrifies, oil capacity will be liberated for use in chemicals and a 4% annual growth in chemicals could be expected. Some data was also shared looking at CO2 cost versus product value, with formaldehyde and monoethylene glycol coming out as the most expensive (i.e. relatively cheap chemical product made with high CO2 emissions).

3.2 European Federation of Catalysis Societies Applied Catalysis Award

The European Federation of Catalysis Societies (EFCATS) award for 2019 was awarded to Professor Glenn Sunley of BP Plc, UK, and was accompanied by a talk titled ‘Adventures in C1 Chemistry’. Glenn talked through some of his career highlights, including work done on the BP CativaTM process for acetic acid production and the elucidation of the reaction mechanism in collaboration with the University of Sheffield, UK (15). Work was also presented about a collaboration with University of California, Berkeley, USA and the California Institute of Technology, USA, investigating the conversion and mechanistic features of methanol homologation to triptane over an indium iodide catalyst (16, 17). Special mention was also given to some recent Fischer-Tropsch technology development as a result of the long-term collaboration between BP and Johnson Matthey, UK, described as ‘stunning engineering’. The first license for this technology has already been sold and will enable Fulcrum BioEnergy Inc, USA, to convert 175,000 tonnes of household rubbish into 11 million gallons of jet fuel each year (18, 19).

3.3 SunCarbon – Creating a New Value Chain from Forest to Refineries

Christian Hulteberg from SunCarbon (both a process and Swedish based company) described their newly launched process for the transformation of lignin into vehicle fuel. This process is dual-function, removing lignin from a typical pulp mill cycle allowing debottlenecking of the pulp mill evaporator and also providing the mill with an additional value stream. Lignin is extracted by membrane filtration and then homogenised in a homogeneously catalysed process. The resultant product is then purified and mixed with vacuum gas oil to make it liquid and pumpable, before being shipped to a refinery. The adjacent pulp mill is capable to take care of any and all waste products of the SunCarbon process. Once at the refinery, the material undergoes standard hydrotreatment and is transformed into vehicle fuel. The feedstock was reported to contain up to 20% O, which would place significant additional H2 demand on the refinery hydrotreater.

4. Carbon Dioxide Utilisation

An ongoing theme of EuropaCat 2019 was the focus of many research groups on utilisation of CO2. Many talks were presented on addressing various aspects of the Sabatier reaction for CO2 hydrogenation, on the dry reforming of methane (DRM) for H2 production or for the direct synthesis of methanol from CO2.

Matteo Monai from Utrecht University, The Netherlands, presented a short talk on work done tuning metal-support interactions on supported Ni catalysts for the Sabatier reaction hydrogenating CO2 to CH4 (Equation (i)):

(i)

Both a reaction mechanism with carbide and formate intermediates are observed on Ni catalysts and both mechanisms were reported to be sensitive to particle size and local structure. Using strong metal-support interactions and a reducible support, it was possible to ‘embed’ Ni particles into the support material and suppress sintering. Enhanced C–C coupling was also observed on these materials when utilising the embedded metal particles, the mechanism for which was still subject to investigation (20).

An experimental programme investigating the mechanism for carbon formation over rhodium on alumina catalysts in the DRM was shared by Gianluca Moroni from the Politecnico di Milano, Italy (Equation (ii)):

(ii)

A variety of CH4:CO2 ratios were flowed over the catalyst at temperatures from 300–700°C in an operando-Raman annular reactor with gas analysis by Micro GC Natural Gas Analyzer (Agilent, USA). An adverse dependence between CH4 concentration and activity was observed and a presented microkinetic model was said to predict experimentally observed activity with catalyst surface area as an input. Catalyst activity was also observed to reduce with increased carbon monoxide levels, and C deposition was said to begin on the Rh sites before migration and deposition on adjacent sites on the support (21).

Shohei Tada from The University of Tokyo, Japan, gave a talk on work done investigating catalysts for the conversion of CO2 to methanol. Zirconia supported Cu catalysts were discussed and work was presented showing Cu species being metallic during reaction and methanol synthesis being performed at Cu-ZrO2 interface sites (22). Tada went on to discuss the importance of finding a good support material to suppress methanol decomposition to CO, with the methanol-support interaction being said to be key. Methanol weakly adsorbs on an amorphous surface (a-ZrO2) and strongly adsorbs on a monoclinic surface (m-ZrO2), with the stronger interaction yielding more unfavourable decomposition of methanol to CO (23, 24).

5. Conclusions

EuropaCat 2019 was a large conference with wide-ranging themes delivered in multiple sessions over six days. It has not been possible to summarise the conference in its entirety. Major topics included: alternative energy inputs (plasma, electrochemistry, light); ammonia synthesis; biomass valorisation; CO2 to chemicals; dry reforming; Fischer-Tropsch catalysis; industrial forum; mechanistic insights; nitrogen oxides (NOx) reduction; organometallic catalysis; and zeolite catalysis. More information can be found on the conference website.

The Authors


Andrew Richardson graduated from the University of Newcastle Upon Tyne, UK with an MChem in Chemistry. Andrew joined Johnson Matthey in 2012 and has worked across a number of catalysis and absorbent research areas. He is currently a Principal Researcher at Johnson Matthey’s Chilton site in Billingham, UK.


Katie Smart has a PhD in Chemistry from the University of Paul Sabatier, France, where she studied at the Laboratoire de Chimie de Coordination in Toulouse. Katie has worked at Johnson Matthey since 2013 and is currently Technical Development Manager for the high-temperature shift research and development team at the Chilton site in Billingham.

By |2020-03-19T16:16:52+00:00March 19th, 2020|Weld Engineering Services|Comments Off on EuropaCat 2019

Observing Solvent Dynamics in Porous Carbons by Nuclear Magnetic Resonance

Understanding the performance of activated carbons in their applications as energy storage materials or catalyst supports requires a description of the behaviour of adsorbates, including solvents, gases, organic molecules or ions. NMR allows adsorption of species in porous carbons to be studied due to the nucleus-independent chemical shift (NICS) arising from the aromaticity of the pore walls (1, 2). The NICS allows adsorbed species to be distinguished from ex-pore species in the bulk solution external to the carbon particles and has been shown to depend on, among other parameters, the distance between the adsorbate and the pore walls. In recent years, several studies have exploited this to use NMR as a probe of the pore structure of a range of porous carbons. Borchardt et al. showed that the NICS of adsorbed organic electrolyte species varied in accordance with the pore size of titanium carbide-derived carbons (CDCs) which have very well-defined porosity (3). This was supported in work by Forse et al. who also showed that ion adsorption in CDCs is significantly reduced when the average pore size is smaller than the solvated ion size (4). In a subsequent study, density functional theory (DFT) was used to show that for model carbon slit pores, the NICS depends upon both the pore width and the size and curvature of carbon fragments making up the pore walls (5). Xing et al. used DFT calculations to derive a relationship between the magnitude of the NICS and the pore size assuming a slit pore geometry, which was found to agree well with experimental measurements of NICS for aqueous species adsorbed on poly-ether-ether-ketone (PEEK) derived carbons (PDCs) (6).

Another important phenomenon that can be studied by NMR is the dynamics of the often highly mobile adsorbate and solvent species. Diffusion coefficients of adsorbed species can be determined experimentally via NMR using pulsed field gradient (PFG) techniques (7). A number of PFG NMR studies have shown that diffusion of species confined in carbon micropores is reduced significantly compared to bulk solution. Furtado et al. observed a broad distribution of local diffusivities for a carbon with a bimodal pore size distribution which was interpreted in terms of restricted internal diffusion between pores of different sizes (8). PFG NMR measurements on ethylene carbonate and dimethyl carbonate mixtures by Alam and Osborn Popp showed that diffusion coefficients for species adsorbed in carbon micropores were reduced by up to a factor of five compared to bulk solution (9). Forse et al. used PFG NMR to study microporous carbon supercapacitor electrodes and observed significant reductions in the diffusion coefficients of adsorbed species, although acetonitrile solvent species were found to diffuse faster than the larger electrolyte ions (10).

While PFG NMR provides significant insight into the dynamics of species confined within the porous carbon network, further information regarding the dynamic exchange of species between pore environments and also between the in-pore and ex-pore environments can be obtained from exchange spectroscopy (EXSY) measurements (11). In this experiment, for which the pulse sequence is shown in Figure 1, the first pulse and subsequent delay allows the single-quantum magnetisation to precess at a characteristic frequency. The second pulse converts this coherence into a population state, which is maintained for a duration called the mixing time, tmix. For zero mixing time it is expected that there will be no diffusion and for long mixing times it is expected that the solvent molecules will reach an equilibrium between the different pore locations. After the mixing time a third pulse is applied and the spins once again precess at their characteristic frequencies and a signal is detected. After a Fourier transformation of both time dimensions is performed, a 2D spectrum is obtained. Magnetisation that has not been exchanged will have precessed at the same characteristic frequency during both times, giving rise to a peak on the diagonal of the spectrum. Magnetisation that has been exchanged between in-pore and ex-pore locations will have precessed at one frequency initially and another during detection. A peak will appear off-diagonal in the spectrum and at the direct-dimension chemical shift of the spin where the molecule finally resided. These two cases give rise to diagonal peaks and cross-peaks, respectively. Further information can be extracted by performing the experiment with varying mixing times. Integrating the cross-peak area and plotting as a function of mixing time yields a build-up curve, from which exchange rate constants can be extracted (12, 13).

Fig. 1.

Pulse sequence of the EXSY experiment. For each tmix chosen, the sequence is repeated with increasing t1 values

Pulse sequence of the EXSY experiment. For each tmix chosen, the sequence is repeated with increasing t1 values

The EXSY approach was applied by Alam and Osborn Popp who showed that the in-pore environment is inhomogeneously broadened due to species occupying a range of pore environments, between which exchange takes place on the millisecond timescale (9). Griffin et al. showed for a commercial porous carbon saturated with an organic electrolyte that in-pore–ex-pore exchange of the anionic species also takes place on the millisecond timescale, although did not fit well to a single exchange process (14). Fulik et al. subsequently showed that build up curves in EXSY data for commercial activated carbon saturated with organic electrolyte can be interpreted in terms of two processes with different rates, i.e. a slow process attributed to diffusion from the centre of the particle to the surface, and a much faster process attributed to effective exchange between ex-pore and in-pore (15). In addition, it was shown that exchange dynamics can also affect the observed NICS and lineshape in 1D spectra, whereby the onset of in-pore–ex-pore exchange upon saturation of the micrometre-sized particles leads to a reduction of the observed NICS (4). In these typical cases, the diffusion path between the ex-pore and in-pore environment is much shorter, allowing faster exchange of species between ex-pore and in-pore which leads to partial or complete exchange averaging of the in-pore and ex-pore resonances.

However in our previous work (16) we showed that the NICS did not change upon saturation for 100 μm PDC particles soaked with water. An example of this observation is shown in Figure 2(a). This was attributed to the relatively large particle size used in the experiments, meaning that in-pore–ex-pore exchange does not take place on the timescale of the NMR experiment. Furthermore, nitrogen gas sorption analyses of PDCs revealed the presence of at least three different pore widths of 0.8 nm, 1.2 nm and 2.2 nm, (Figure 2(b)) whereas only a single in-pore resonance was observed. This was rationalised in terms of fast exchange between micro- and mesopores within the pore network, leading to motional averaging of the in-pore resonance. However, the single in-pore resonance was found to be subject to inhomogeneous broadening, as evidenced by a purely diagonal in-pore diagonal peak in 2D EXSY spectra at short mixing time. This means that within a particle, there are regions of different averaged NICS, which can be due to a variation of aromaticity as well as local average pore size. The latter depends mostly on the inhomogeneity of activation throughout the carbon particles. This was experimentally minimised by activating smaller particles. The symmetric in-pore components coming from the regions of fast averaged NICS can, as shown by Merlet et al., be broadened by slow intra-particle exchange (17). The size of these regions of averaged NICS can be estimated to be around 1 μm in diameter, assuming that the in-pore diffusion coefficient is the same as ex-pore, and that a 1 kHz frequency is enough for the fast regime. From ex-pore to in-pore however, the diffusion path was much longer in our particles. It arises that there must be regions with fast in-pore–ex-pore exchange and more isolated regions.

Fig. 2.

(a) 1H NMR spectra of a typical PDC wetted with increasing amounts of water shown in terms of multiples of the total pore volume (PV); (b) corresponding N2 gas sorption pattern

(a) 1H NMR spectra of a typical PDC wetted with increasing amounts of water shown in terms of multiples of the total pore volume (PV); (b) corresponding N2 gas sorption pattern

In this work, we systematically investigate the effects of carbon particle size, porosity and solvent properties on in-pore–ex-pore and in-pore–in-pore exchange dynamics as viewed by NMR spectroscopy. We chose PDCs as our model system due to the ease of synthesis and tunability of the porosity whereby the pore volume (PV) and average pore size vary approximately linearly with burn-off (BO) and activation time (18), in addition to these materials giving rise to generally strong NICS, facilitating structural characterisation and analyses of adsorbate behaviour (1921). We first provide a review of the effects of two-site exchange on 1D NMR spectra with specific reference to a model for a carbon particle saturated with solvent species. We then discuss the experimental results in three sections: first, we compare PDC particles of ~80 μm diameter with ~15–20 μm particles obtained from the same sample, to confirm that our samples show minimal exchange broadening due to a longer path between the two environments. Second, we compare PDC samples that were steam-activated with 20% BO and 54% BO to see the extent at which diffusion coefficients in pores of different sizes influence the exchange rate constants. Finally, we compare PDC samples saturated with different solvents, namely water, hexane and cyclohexane. These solvents were chosen because of their comparable or different polarity and viscosity, two parameters which are expected to influence the diffusion of the solvent molecules in the pores. The pore filling is assessed in unsaturated samples to determine the accessibility of the pore network. Ex-pore–in-pore exchange rate constants were determined using exchange experiments and were then used to better understand how various exchange regimes perturb 1D NMR spectra.

Averaging of NMR signals is a common phenomenon observed when individual nuclei explore several magnetic environments. In the case of a solvent adsorbed in porous carbons, the diffusion of solvent molecules between regions giving rise to different NICS affects the NMR spectrum. The perturbation depends on the average time each molecule remains in each environment. Considering a nucleus able to explore two environments, several cases can be distinguished: the slow exchange regime where both peaks may be broadened but without being shifted, the intermediate regime where both peaks have merged into a single very broad peak and the fast exchange regime where the single peak narrows. To understand these observations, we can turn to the underlying principles of the Fourier transform.

In the fast exchange regime, the resulting peak is located at the average of the chemical shifts weighted by the residence time in each environment. To illustrate this, various NMR signals were simulated using cosine waves and subsequently Fourier transformed by fitting cosine functions of variable frequencies, whereby the NMR spectrum is obtained as the integral of the product between the original and the fitted functions (full details of the simulation are given in Supplementary Information (SI)). This simulation is a simplified case of fast exchange where motional broadening is not accounted for. Its purpose is to show how the averaged chemical shift depends on the relative dwell times. Figure 3(a) shows the simulated NMR signals; the blue component oscillates at a frequency of 30 Hz and the red components oscillate at 50 Hz. The ratio of the time spent oscillating at either frequency is given: an average of approximately 32 Hz is obtained when adopting 30 Hz for 0.1 ms and 50 Hz for 0.01 ms, and an average of 40 Hz is obtained when oscillating at both frequencies for the same time. For better visibility, the angle in radians travelled as a function of time is plotted in Figure 3(b). The spectra obtained after Fourier transformation, Figure 3(c), show how the apparent frequency of the exchange-averaged peak depends on the residence time in each environment. When the residence times of the nucleus in both environments are identical, we call it a symmetric exchange. This is found for example in an acid-base mixture when pH = pKa, where by definition both acid and conjugated base are equimolar. However, asymmetric exchange is more relevant to the present study due to the large population difference between the in-pore and ex-pore environments, and also between the volume of the different connected pores.

Fig. 3.

(a) Portion of simulated cosine waves oscillating at 30 Hz (blue) or 50 Hz (red) for various periods of time; (b) angle evolution as a function of time; and (c) corresponding Fourier transformed spectra

(a) Portion of simulated cosine waves oscillating at 30 Hz (blue) or 50 Hz (red) for various periods of time; (b) angle evolution as a function of time; and (c) corresponding Fourier transformed spectra

Different chemical shifts can be attributed to different regions inside and around a porous carbon particle. Figure 4 illustrates a carbon particle (black continuous line) saturated with a solvent. Molecules in the green region adopt the ex-pore chemical shift, and molecules in the blue region an in-pore chemical shift. The deeper the colour, the more likely the molecule is to change environment within a certain time. The dashed lines represent the boundaries between molecules undergoing fast and slow exchange: within the dashed lines, the in-pore–ex-pore exchange (orange arrows) is frequent enough for the chemical shift to be averaged. Note that the ex-pore dashed line is at a constant distance from the particle surface because the ex-pore diffusion coefficient is constant. The in-pore dashed line however is distorted due to inhomogeneities of in-pore diffusion coefficients and wraps around big mesopores penetrating deep into the particles. In summary we expect two non-exchanging peaks corresponding to ex-pore solvent far from the particle and in-pore solvent in the core of the particle, and a broadened peak, corresponding to solvent molecules undergoing fast exchange.

Fig. 4.

Scheme of an activated carbon particle (black continuous line, delimiting the ex-pore/in-pore boundary) saturated with solvent, blue for adsorbed and green for bulk. The orange arrows represent fast exchange of molecules between the two environments. The colour gradient represents the probability of a solvent molecule exchange during data acquisition. Molecules between the two dashed lines may exchange fast enough to appear at the exchange-averaged chemical shift, with Vinexch and Vexexch standing for the volume of such exchanging solvent initially in-pore or ex-pore respectively

Scheme of an activated carbon particle (black continuous line, delimiting the ex-pore/in-pore boundary) saturated with solvent, blue for adsorbed and green for bulk. The orange arrows represent fast exchange of molecules between the two environments. The colour gradient represents the probability of a solvent molecule exchange during data acquisition. Molecules between the two dashed lines may exchange fast enough to appear at the exchange-averaged chemical shift, with Vinexch and Vexexch standing for the volume of such exchanging solvent initially in-pore or ex-pore respectively

The position of the exchange-averaged peak depends on the time the molecules spend in each environment. However in practice, when adsorption and desorption are at equilibrium, the solvent molecules exchange by pairs because the volumes of exchanging adsorbed solvent Vinexch and exchanging free solvent Vexexch are constant. The diffusion coefficient in the pores is smaller than in the bulk, so for the in-pore molecules to exchange in the fast regime they must reside closer to the interface. Therefore, Vinexch is smaller than Vexexch. The residence time of the exchanging solvent molecules can now be correlated to Vexexch and Vinexch and the in-pore/ex-pore ratio of diffusion coefficients. Experimentally it is possible to calculate Vexexch and Vinexch from the position and volume of the broad ex-pore peak.

Particle Size Effects

To directly observe the impact of particle size on 1D NMR spectra, two activated PDC samples were reduced from approximately 100 μm particle size to approximately 15–20 μm by sieving (see Methods in SI). The samples were named xBO_y, where x is the percentage of BO and y the median particle size D50, measured by dynamic light scattering. The samples were in the first instance wetted using a microsyringe with a defined volume of deionised water less than the PV, then with a volume greater than PV, to observe the 1H NMR spectrum before and after saturation. We have previously shown that this method allows us to compare the NICS averaged over the whole pore network without, and then including, perturbations related to in-pore–ex-pore exchange. This is because the in-pore peak before sample saturation corresponds to water located in completely filled particles that are not yet surrounded by water. Figures 5(a) and 5(b) show the spectra of sample 54BO_80 and 54BO_21, respectively. The in-pore peak shifts by approximately 0.2 ppm upon saturation, regardless of the particle size. This means that in this range of particle size, diffusion of water out of the pores has a negligible impact on the NICS for most of the adsorbed water and the width of the in-pore peak is solely due to the distribution of average NICS. However, exchange averaging has a measurable impact on the ex-pore peaks. In both samples they were fitted with a narrow and a broad component, the intensity and full width at half maximum (FWHM) of which vary. The broad component is assigned to ex-pore water having experienced the pores for a period of time and is therefore also partially homogeneously broadened by exchange-averaging; the longer the residence time in the pores, the broader and the more shifted the broad component is, i.e. the bigger the ratio Vinexch/Vexexch as per Figure 4. In big particles, the broad component amounts to 30% of the ex-pore peak with a FWHM = 0.30 ppm and is located within 0.10 ppm of bulk water chemical shift. On the other hand, with small particles the broad component represents 82% of the total ex-pore water with a much bigger FWHM = 0.70 ppm and is shifted by 0.20 ppm relative to bulk water. This indicates a slight increase of Vinexch/Vexexch when the particles are reduced; in other words, a bigger proportion of ex-pore water is able to experience the in-pore environment for a longer time.

Fig. 5.

1H NMR spectra of: (a) 54BO_80 injected with 0.2*PV (red) and 2.3*PV (blue) of water; and (b) 54BO_21 injected with 0.3*PV (red) and 2.3*PV (blue) of water. The dashed lines show the maxima of the in-pore peaks

1H NMR spectra of: (a) 54BO_80 injected with 0.2*PV (red) and 2.3*PV (blue) of water; and (b) 54BO_21 injected with 0.3*PV (red) and 2.3*PV (blue) of water. The dashed lines show the maxima of the in-pore peaks

To explain the broadening of the ex-pore peak, we must adopt the point of view of water molecules that spend the majority of their time in the ex-pore environment, where we believe the packing of the particles plays an important role. A simple model consisting of spherical particles (Figure S3 in the SI) allows exchange rates to be calculated in function of the particle diameter (Figure S4). With 80 μm particles (similar to 54BO_80) we find exchange rates of 2.5 Hz, and of 36.7 Hz with 54BO_21 particles (see SI for more details on the calculations). The NICS being around 3000 Hz, the exchange regime would be slow in both cases, which is consistent with our observations, and increases by a factor of 15 when reducing the particle size from 80 μm to 21 μm.

Examples of 2D EXSY NMR spectra are shown in Figures 6(a) and 6(b) corresponding to sample 54BO_80 and 54BO_21 saturated with water. The mixing time (tmix) was 20 ms for both spectra, and it can be seen that the cross-peaks are more pronounced with small particles, giving a first indication that the exchange is faster. Figures 6(c) and 6(d) show the build-up curves of the intensity ratio of cross-peaks over diagonal peaks. Visually, we can see that the curve for big particles reaches the maximum after long tmix intervals, whereas for small particles the build-up is complete within 0.1 s.

Fig. 6.

2D exchange spectra for: (a) 54BO_80; and (b) 54BO_21 recorded with tmix = 20 m, with (c) and (d) corresponding build-up curves of the ratio of cross- and diagonal peak intensities (Icross/Idiag) as a function of mixing time

2D exchange spectra for: (a) 54BO_80; and (b) 54BO_21 recorded with tmix = 20 m, with (c) and (d) corresponding build-up curves of the ratio of cross- and diagonal peak intensities (Icross/Idiag) as a function of mixing time

In principle, the build-up of the ratio of cross- and diagonal peak intensity (Icross/Idiag) in EXSY spectra can be described by a single dependence on tanh(ktmix), where k is the exchange rate constant. However, Fulik et al., have shown that better agreement is observed if two processes with different rates are assumed, i.e. a slow process attributed to diffusion from the centre of the particle to the surface and a much faster process attributed to effective exchange between ex-pore and in-pore (15). However, over the course of our experiments, it was observed that the ex-pore resonance reduced in intensity due to evaporation from the NMR rotor leading to a global reduction in Icross/Idiag for EXSY spectra recorded at the end of the series with long mixing times. Although the exact kinetics of the solvent evaporation are complex and were not studied in detail, we found that this could be accounted for with sufficient accuracy through the incorporation of an exponential term with a characteristic decay constant Tevap (see SI). Rate constants were therefore extracted from the EXSY data using Equation (i):

(i)

where I0 is an additional constant introduced to account for t1 noise giving rise to spurious off-diagonal low-intensity signal at zero mixing time. The fits were optimised by minimising the root mean squared deviation (RMSD) between calculated and experimental points, to converged values around 10–2. The errors on the data points were calculated from the signal-to-noise ratio of each peak for a selection of 2D spectra and propagated to Icross/Idiag, and were found to be smaller than 0.02%. The best fit for 54BO_80 was obtained with k1 = 57 Hz and k2 = 4 Hz, and for 54BO_21 with k1 = 597 Hz and k2 = 50 Hz. The carbon particles with D50 = 21 μm showed around 10 times faster in-pore–ex-pore diffusion versus the D50 = 80 μm particles. This is close to the factor of 15 which we obtained from our simple calculations based on spherical particles (see Figure S4) and shows the impact of the particle size on the rate of the exchange processes, with larger particles significantly reducing the exchange kinetics between the in-pore and ex-pore environments. Therefore, particle size is an important factor to take into account when comparing in-pore–ex-pore exchange phenomena in porous carbon samples.

Regarding the in-pore resonances, the shape and position in the NMR spectrum is largely unaffected by in-pore–ex-pore exchange despite the difference in particle size. We can therefore assume that the diffusion path out of the particle from any point of the pore network, apart from the very surface, is simply too long and exchange-averaging of the in-pore environment is in the slow regime. This would mean that if the particle size is decreased further, at some point the in-pore peak should also become affected by faster diffusion of in-pore water into the ex-pore environment. To test this hypothesis, 6 μm-sized YP50 particles were wetted with deionised water. YP50 is a commercial activated carbon, and similar to our PDCs with respect to composition, pore size distribution and average pore size (see Figure S13). The 1H NMR spectra before and after saturation are shown in Figure 7. As expected due to the small particle size, an ex-pore peak with sharp and broad components is observed. Relative to 54BO_21, the broad component is shifted five times more, consistent with smaller inter-particle voids that facilitate adsorption of ex-pore water in close proximity. The in-pore peak shifts from a NICS of 7.3 ppm to 7.0 ppm upon saturation, which is similar to 54BO_21, however the intensity decreases significantly, and we notice that the peak exhibits a tail towards the ex-pore peak. This means that after saturation, a significant proportion of adsorbed water is able to quickly diffuse into the ex-pore environment and therefore does not appear at the purely in-pore chemical shift, but rather intermediate between the ex-pore and the in-pore chemical shifts. The exchange rate constants were estimated to be k1 = 1117 Hz and k2= 165 Hz. Assuming that k1 relates to the ex-pore–in-pore exchange process that we calculate as a function of the particle size (Figure S4), we find that 4 μm particles give a similar exchange rate. This is good agreement considering that 20% of YP50 particles are smaller than 2.5 μm (Figure S14), especially because for particles smaller than 10 μm the exchange rate increases sharply. This explains why YP50 presents a relatively flat valley between the two broad peaks: even a narrow particle size distribution in the few-micrometre range gives a very broad distribution of exchange rates, therefore generating exchange-averaged peaks at a continuum of shifts.

Fig. 7.

1H NMR spectra of YP50 wetted with 0.5*PV (red) and 5.3*PV (blue) of water. The dashed lines show the maxima of the in-pore peaks

1H NMR spectra of YP50 wetted with 0.5*PV (red) and 5.3*PV (blue) of water. The dashed lines show the maxima of the in-pore peaks

One point worth addressing is whether the inhomogeneous broadening of the in-pore peak is responsible for broadening the exchange-averaged ex-pore peak, which can give an idea of the reliability of the exchange rate constants obtained. The spectrum of Figure 7 was successfully reproduced in Express software (22) in the case of a homogeneous as well as a inhomogeneous broadening of the broad ex-pore peak, so the simulation does not allow it to be determined whether the exchange-averaged peak is homogeneously or inhomogeneously broadened. Refer to SI for further detail. However in 2D EXSY spectra, for any tmix, the ex-pore peak shows symmetrical off-diagonal broadening meaning the broadening is homogeneous. This suggests that one component of the in-pore peak is much more exposed to the ex-pore environment than the other components. In consequence the values of exchange rate constants may be more reliable than if the broadening were inhomogeneous. This result is consistent with a radial distribution of NICS in the particle, in agreement with our previous observations where the centre of particles of diameter greater than 100 μm is poorly affected by steam activation. The particles employed here were smaller than 100 μm but it seems that a small gradient of activation still remains. It is possible that the valley between ex-pore and in-pore is the result of in-pore–ex-pore exchange broadening, in the fast or slow regime, of the other components of the in-pore peak, which are located further and further away from the surface, and therefore have access to smaller and smaller volumes Vexexch.

In summary, from these observations it could be deduced that particles smaller than 10 μm offer a diffusion path short enough for water to diffuse out of the pores at a rate that affects the in-pore peak as well. However, given that the diffusion coefficient of adsorbed water depends on the pore size, we cannot compare 54BO_21 and YP50 because they have different average pore sizes. Using the equation provided previously (6), these can be estimated using the NICS before saturation. We obtain for 54BO_21 and YP50, 1.17 nm and 1.10 nm respectively. Therefore, it is necessary to investigate the effect of average pore size, which can be tuned by controlling the BO, on the diffusion of water and the appearance of the spectra.

The Effect of Burn-off

The effect of BO, and thus average pore size, on the dynamics of water in PDCs can be described by comparing 54BO samples (high BO) with 20BO samples (low BO). Figure 8(a) shows samples 20BO_83 and Figure 8(b) shows 20BO_15 injected with water. The NICS of samples 20BO are 8.4–9.0 ppm, giving an average pore size less than 1.0 nm based on the previous equation (6), which is smaller than sample 54BO, in agreement with the linear dependence of the average pore size on the BO (18, 19, 23). The in-pore peaks appear to contain two components as is sometimes observed for samples with low BO. Possibly, the degree of activation was not homogeneous from the surface to the centre of the particles. In sample 20BO_83, the in-pore peak is constant regardless of the injected volume (shift < 0.10 ppm). Half the ex-pore peak was fitted with a sharp component of FWHM <0.1 ppm, and half with a broad component of FWHM = 0.16 ppm with a difference in chemical shift < 0.10 ppm. This suggests in the first instance that the in-pore–ex-pore exchange process is in a slower regime than in 54BO_80, where the broad component was broader and the in-pore peaks were shifted slightly, in agreement with the study mentioned earlier (24). With small particles 20BO_15, 90% of the ex-pore peak is broad, FWHM = 0.35 ppm but within 0.10 ppm of the bulk water, suggesting again limited exchange with negligible Vinexch.

Fig. 8.

1H NMR spectra of: (a) 20BO_83 injected with 0.7*PV (red) and 3.6*PV (blue) of water; and (b) 20BO_15 injected with 0.2*PV (red) and 3.6*PV (blue) of water. The dashed lines show the maxima of one component of the in-pore peaks

1H NMR spectra of: (a) 20BO_83 injected with 0.7*PV (red) and 3.6*PV (blue) of water; and (b) 20BO_15 injected with 0.2*PV (red) and 3.6*PV (blue) of water. The dashed lines show the maxima of one component of the in-pore peaks

The exchange rate constants provide a more quantitative description of the impact of average pore size. For the big particles 20BO_83, k1= 54 Hz and k2 = 4 Hz, which are comparable to the values for 54BO_80 (57 Hz and 4 Hz). This was perhaps not expected from the 1D spectra, where 20BO_83 appeared much less affected by diffusion. However, one must keep in mind that the exchange regime depends not only on the actual rates, but also on the chemical shift difference of the two environments involved. The NICS of samples 20BO is higher than in samples 54BO, so even identical exchange rates still situate 20BO in a slower regime than 54BO.

For the small particles 20BO_15, k1 = 162 Hz and k2 = 6 Hz, which are three times as fast when compared to 20BO_83. However, the calculated factor was 31, which is a clear discrepancy. Given that sample 20BO_15 contains particles smaller than sample 54BO_21, we can attribute the discrepancy to a BO effect and not to a particle size effect. Besides the average pore sizes, other structural parameters that could perhaps vary with BO are oxygen content and tortuosity due to the smaller amount of mesopores. These observations also raise the question of homogeneity of the diffusion coefficients within the particle: it is likely that at low BOs a bigger gradient of activation within the particle is present.

Overall, these results show that a smaller average pore size hinders exchange as opposed to a smaller particle size, which promotes it. Since YP50 has a smaller average pore size than 54BO_21 but higher exchange rate constants, we can safely deduce that in YP50, the particles are truly small enough to allow for in-pore–ex-pore exchange to affect the in-pore peak. These conclusions are in agreement with simulations (25), and also an experimental study focusing on diffusion measurements in hydrophobic slit pores by neutron-scattering (23). The material contained approximately 95% carbon with 5% oxygen atoms on the surface of the pores and average pore sizes of 1.2 nm and 1.8 nm, which is similar to our PDCs. It was found that the diffusion coefficients of water in the 1.2 nm and 1.8 nm pores are 40% and 30% smaller than that of bulk water, respectively. Furthermore, on the very surface of pore, water diffuses at 0.035 × 10–5 cm2 s–1 and 0.014 × 10–5 cm2 s–1 respectively, which is two orders of magnitude slower than in the bulk. Counterintuitively, the value on the surface of the bigger pores was found to be nearly half that in the smaller pores, which was attributed to the promotion of slightly faster concerted motion in extreme confinement, although the centre of the pore followed an overall slower regime.

The Effect of Solvent Properties

Many applications of porous carbons (such as electric double-layer capacitors) commonly employ electrolytes in organic solvents as an alternative to aqueous electrolytes. To extend the knowledge presented here to such systems, it is desirable to be able to predict how the NICS may be affected. The viscosity and the polarity are the two parameters that will be considered here as tools to predict the diffusion regime adopted by any solvent.

The effect of polarity can be estimated by comparing water and cyclohexane which have the same viscosity, but different dipole moments. YP50 is chosen for this comparison because the exchange rates are high enough to significantly impact the NMR spectra. The exchange rate constants were measured to be 1117 Hz and 165 Hz for water, and 257 Hz and 23 Hz for cyclohexane, which means that cyclohexane diffuses roughly four times slower. On one hand, this could be related to increased van der Waals interactions with the hydrophobic surface of the pores, as was also found to be the case in a xerogel using the same solvents (24). On the other hand, the bulk self-diffusion coefficient of cyclohexane is smaller than for water (1.42 × 10–5 cm2 s–1vs. 2.3 × 10–5 cm2 s–1) (26), which could also contribute. At this stage it is unclear whether other parameters such as the different molecular sizes of water and cyclohexane play a role. Figure 9 shows the 1H NMR spectrum of YP50 wetted with cyclohexane. As expected with such high exchange rate constants, we see a narrow and a broad ex-pore peak, and an in-pore peak that decreases and shifts upon saturation, much like in Figure 7 with water in YP50. With water and cyclohexane, the in-pore peaks are similar; they shift by 0.30 ppm and are broadened by a factor two upon saturation. Interestingly, we note that the NICS of cyclohexane (and hexane) is smaller than water by 0.2–0.3 ppm, which means that the apolar solvents are on average located in slightly bigger pores than water. The diffusion being faster in big pores, this is in contrast with the slower exchange kinetics measured, showing that solvent parameters are prevalent over pore size effects. The broad ex-pore peak is narrower and less shifted in cyclohexane (FWHM = 0.64 ppm and shift = 0.27 ppm) than in water (FWHM = 1.69 ppm and shift = 0.98 ppm), consistent with smaller exchange rate constants.

Fig. 9.

1H NMR spectra of YP50 injected with 0.3*PV (red) and 2.9*PV (blue) of cyclohexane. The dashed lines show the maxima of the in-pore peaks

1H NMR spectra of YP50 injected with 0.3*PV (red) and 2.9*PV (blue) of cyclohexane. The dashed lines show the maxima of the in-pore peaks

The effect of viscosity can be observed by comparing cyclohexane (0.89 cP) and hexane (0.30 cP) which have the same very low polarity. Figure 10 shows the 1H NMR spectrum of YP50 wetted with hexane. Unlike water and cyclohexane, hexane shows two ex-pore peaks, which are simply due to the two visible proton environments on the molecule; 1.28 ppm are the methylene (CH2) and 0.88 ppm are the terminal methyl (CH3) protons. The in-pore peak before saturation can be fitted with two broad peaks 0.30 ppm apart, which is similar to the difference between the CH2 and CH3 chemical shifts (refer to SI for details on the fits). This suggests that all protons of adsorbed hexane molecules roughly experience the same NICS, therefore the molecule is either tumbling isotropically or aligned with the pore wall. The exchange rate constants are k1 = 784 Hz and k2 = 104 Hz for hexane in YP50, which is approximately three times the rates seen for cyclohexane. Therefore in this case, the exchange rate was roughly proportional to the viscosity, as hexane is three times less viscous than cyclohexane.

Fig. 10.

1H NMR spectra of YP50 injected with 0.3*PV (red) and 3.0*PV (blue) of hexane. The dashed lines show the maxima of the in-pore peaks

1H NMR spectra of YP50 injected with 0.3*PV (red) and 3.0*PV (blue) of hexane. The dashed lines show the maxima of the in-pore peaks

More interestingly, the broad ex-pore peak is shifted by 0.61 ppm relative to the average between the two narrow peaks, and the FWHM = 1.30 ppm. These values are between the values for cyclohexane and water. This is consistent with the shift and width being proportional to the exchange rate constants as they are even higher for water in YP50. Overall, pronounced solvent exchange affects the ex-pore and in-pore peaks and the valley between the two peaks, but the ex-pore peak seems to be the most reliable indicator to estimate at a glance perturbations due to exchange effects.

Further discussion is required about parameters that could not be assessed with this set of experiments, for example the kinetic diameter of the solvent molecules. This factor was not considered separately because it was first assumed to be reflected in the viscosity parameter, although in pores of similar size, viscosity and diameter may have independent contributions to the exchange kinetics. In addition, to identify the impact of this factor alone, two solvents with similar viscosity and polarity but different molecular sizes should be compared. The kinetic diameter of water and cyclohexane perhaps contributes to the exchange rate difference.

The behaviour of cyclohexane was peculiar in the case of sample 20BO_83 and 20BO_15. In the big particles, exchange was too slow to be observed, and by decreasing the particle size, the exchange became visible but the rate was higher than that of hexane and also of cyclohexane in the other PDC samples, which goes against the trends. Another interesting observation was that there were two in-pore peaks in 20BO before saturation and both peaks shifted to the right upon saturation, meaning exchange-averaging takes place within different pores but not with ex-pore, while in 54BO there was only one in-pore peak (see Figures S5 and S6). We believe these observations are all related and hint towards the ability of cyclohexane to form organised structures in slit-like micropores, which affects its diffusion coefficient. The diffusion coefficient of confined cyclohexane was shown by Fomin et al. to drop to zero in slit pores smaller than 2.1 nm, and its density to increase two-fold from 2.6 nm to 1.2 nm pores (27). Another study showed that cyclohexane forms a monolayer in 0.8 nm pores and a bilayer in 1.0 nm pores with a denser hexagonal packing structure resembling the solid phase (28). Similar behaviour has also been observed for propylene carbonate which was observed to form an ordered structure upon nanoconfinement (29). In a disordered structure the packing is expected to be less efficient, nonetheless the diffusion coefficient may still be significantly lower. The exchange rate constants determined take into account the diffusion of species in all types of pores located in Vinexch. When cyclohexane forms an immobile structure in the smallest pores, only the diffusion coefficient in the biggest pores where it is still liquid contributes to the exchange rate constant, explaining why in 20BO_15 the exchange rates are much higher than in the other PDC samples. The amount of mesopore is smaller than in 54BO and YP50, so the long-range diffusion within the particles is probably significantly slower. In that regard, the gradient of NICS will be less well averaged, and even more so with slow diffusing solvents like cyclohexane. The 3.3 ppm NICS of the small peak gives a pore size of 2.2 nm and the 7.7 ppm NICS of the main peak 1.0 nm. The small in-pore peak is therefore likely to come from few isolated mesopores while the other from micropores and mesopores in contact.

By |2020-03-19T15:51:41+00:00March 19th, 2020|Weld Engineering Services|Comments Off on Observing Solvent Dynamics in Porous Carbons by Nuclear Magnetic Resonance

Plasma Catalysis: A Review of the Interdisciplinary Challenges Faced

Home > Journal Archive > Plasma Catalysis: A Review of the Interdisciplinary Challenges Faced

Johnson Matthey Technol. Rev., 2020, 64, (2), 138

Introduction

One area where physics and chemistry come together across the disciplines is in the field of plasma catalysis. A plasma can be described as a ‘soup’ of species including molecules and atoms that are charged or excited and free electrons. Depending on the energy of the plasma it can be fully ionised and have a bulk temperature of tens of thousands of kelvin (or in the case of nuclear fusion millions of kelvin) or it can be in a low temperature non-equilibrium state where only a small portion of the gas phase is energised. At the lower end of the energy spectrum a non-thermal or cold plasma will have a bulk temperature a few tens of degrees above ambient and yet still have some exceedingly high energy species present. Even though the majority of the species in the plasma are not in equilibrium, there can exist some partial equilibria among species with similar kinetic temperatures existing at localised sites (1).

This energy can be used to initiate chemical reactions in the gas phase and on the surface of solids. The term plasma catalysis can be used to describe both the use of plasma to initiate a reaction directly and the use of plasma in combination with a catalytic material. The catalyst can be positioned after a plasma zone as post plasma catalysis (PPC) or within the plasma zone as in-plasma catalysis (IPC). The PPC configuration allows the longer-lived excited species and the products from the plasma reaction to interact with the catalyst and the IPC configuration allows a greater opportunity for the catalytic material to influence the nature of the chemistry by directly interacting with the plasma excitation as well as the products from the plasma-initiated reactions.

With the introduction of a plasma to a catalytic reactor the standard models of the chemistry and the traditional understanding of the mechanisms by which reactions take place start to become less relevant due to the non-equilibrium concentration of excited species such as free atoms, electrons and radicals (2).

Plasma Generation

Within the universe around us there are many examples of plasma (3) including the sun, stars, auroras, lightning, welding arcs and fluorescent lighting tubes. A plasma can be generated in different ways, but all require energy to be applied, either as heat in the case of a thermal ionisation or through the generation of an intense electric field via the use of electrodes, radio frequency (rf) or microwave (MW).

Electrical excitation is one of the most feasible ways for producing well controlled plasma discharges at industrial scale. This type of excitation is controlled by three main parameters: (a) the applied voltage amplitude, (b) the applied frequency and (c) the waveform shape. The combination of the aforementioned parameters defines different operating regimes (4). The main types are:

  • Direct current (DC) – earliest power supplies where constant high voltage in tens of kilovolts are applied between the electrodes for creating a sustainable discharge between fixed anode and cathode electrodes. A resistor may be used to limit the current and high voltage cables are used for power delivery to the electrodes

  • Sinusoidal – high voltage (0–40 kVp) continuous power sources in the frequency range of 50 Hz to 150 kHz are mostly used in dielectric barrier discharges (DBD). Those sources are easier to manufacture and operate in a wide frequency and power range

  • rf – continuous sinusoidal sources with hundreds of volts operating at 13.56 MHz. For optimal operation a matching network is used to restrict the reflected power

  • MW – continuous sinusoidal source operating normally at 2.45 GHz. Rectangular waveguides or coaxial cables are used for delivering the energy to the load

  • Pulsed – fast pulses with well-defined pulse rise time, duration and voltage created by switching a DC high voltage power supply. The fast transition times and control capabilities allow operation at higher power per pulse.

From the regimes described above, sinusoidal and pulsed power sources present the advantage of producing non-thermal plasma without the need for noble gasses and low pressure. Moreover, the use of dielectric barrier material protects the electrodes from erosion in chemical processes and limits the power needed for initiating and maintaining plasma discharge.

It is also possible to influence the nature and stability of the plasma through external forces such as the introduction of an applied magnetic field, a vortex flow or by reducing the pressure within the system.

In scientific and industrial applications, average power delivered in the plasma discharge is critical. Accurately measuring and monitoring the power consumption in plasma discharges is not a straightforward task. Average power consumption is calculated through the measurement of voltage and current waveforms. Although the voltage signals are monitored with high accuracy through capacitive voltage dividers, current waveforms as shown in Figure 1 and Figure 2 present significant measuring difficulties. Current spikes with very short (nanosecond) duration and high amplitude and frequency are imposed over the sinusoidal low frequency and amplitude signal. Monitoring those high dynamic range waveforms can be erroneous and special care needs to be taken in the method of acquiring those signals. In the literature (5), three district ways are described for measuring directly the current or the charge in plasma.

  • Shunt resistor method – a known resistor is inserted between the reactor and the ground electrode. The voltage drop across the resistor is logged and is converted to current using Ohms law

  • Monitor capacitor method – an integrating capacitor with known value is inserted between the plasma reactor and the ground electrode. Plotting its waveform versus the voltage generates a Lissajous curve in which the average power can be calculated

  • Rogowski coil – an inductive coil is used to measure the current through the ground electrode of the reactor.

Fig. 1.

Typical signals for sinusoidal voltage, current and time

Typical signals for sinusoidal voltage, current and time

Fig. 2.

Typical pulsed sources signal for current, voltage and time

Typical pulsed sources signal for current, voltage and time

Although the above measurements have comparable accuracy in scaled plasma reactors, in industrial environments it is preferable that in situ non-invasive techniques are used. Rogowski coils have an intrinsically safe way of operation given that they are galvanically isolated from the main circuit without compromising the accuracy.

Accurate power measurement and delivery in combination with optimised mechanical and electrical design can lead to improvements in energy efficiency. This is a critical parameter that must be considered, especially when scaling up plasma systems. It is unavoidable that some part of the initial energy is lost in the electrical transformation to high voltage and in the plasma reaction as heat. By carefully selecting the electric field and reactor characteristics those losses can be minimised. Modelling can significantly enhance the understanding of the plasma transitions in short time and space frameworks (6). The next improvement step involves the synergistic effects of plasma catalyst interaction. By introducing a catalyst in the plasma region or next to it, different works have shown considerable improvement in conversion and in total efficiency (7). All these optimisations allow energy consumption of the plasma system to be decreased and make plasma technologies feasible from the view point of economic and life cycle assessment. For example, a recent paper by Rooij et al. (8) shows that the combination of plasma with renewable energy sources is an economical method even for such an expensive process as carbon dioxide (CO2) reduction. The technoeconomic solution will of course be different for each market application based upon comparison to current scaled state of the art techniques.

Mechanisms of Excitation

Even the simplest of plasma generation arrangements takes in aspects of a wide range of physics disciplines. A schematic of a simple plasma tube is shown in Figure 3. Gas is passed through the cell, where a high voltage is applied across a central cylindrical electrode and a cylindrical outer electrode. The generation of the plasma in the cell leads to a shift in chemical make-up across the cell, hence the potential to use such cells in chemical processing. Even prior to the excitation of any plasma, the gas flow through the cell has the potential to become a complex fluid mechanics problem. Once a plasma is generated the equations of magnetohydrodynamics become applicable. This apparently innocuous point massively increases the experimental phase space that needs to be controlled and understood. A general introduction into the topic of magnetohydrodynamics can be gained from reading “An Introduction to Magnetohydrodynamics” by P. A. Davidson (9). The bulk physical continuous control parameters in a traditional reactor vessel can be broadly listed as temperatures, pressures and flow rates. In a plasma cell we have in addition voltages, currents and frequencies which must be controlled and monitored. This opens huge opportunities in chemical processing as significant parts of this phase space remain largely unmapped.

Fig. 3.

Schematic of a typical continuous flow cylindrical plasma tube

Schematic of a typical continuous flow cylindrical plasma tube

One of the oldest applications of plasma is in fluorescent lighting. This fact means that there has been a significant amount of study into the behaviour of plasma cells where there is no flow. The graph shown in Figure 4 is a schematic illustration of the characteristic (current vs. voltage) curve of a typical gas discharge in neon gas at a pressure of 1 torr, between two planar electrodes separated by 50 cm. This figure has been recreated from the information presented by Gallo (10). There are broadly speaking three classes of behaviour. A dark discharge region, a glow region and an arc region. In the dark discharge region, a voltage lower than the breakdown voltage of the gas is applied to the tube. External radiation such as gamma photons and beta particles then trigger a Townsend cascade within the tube. This behaviour is utilised in a number of simple nucleonic detection devices, notably Geiger counters (11). In this region the current is limited to remain low and the discharge is reliant on external excitation. In the glow region, the current is no longer limited, and the accelerated electrons excite further electrons and the voltage then drops. The gas emits light causing a glow, an effect utilised in fluorescent lighting. If the current is high enough the neutral gas in the tube becomes heated and arcs start to form which is the final region.

Fig. 4.

Transitions occurring in a fixed glow tube

Transitions occurring in a fixed glow tube

The graph shown in Figure 4 illustrates that the behaviour of plasmas is complex and many transitions occur in a cell where all of the standard chemical processing variables are fixed i.e. pressure is constant, temperature is fixed and the flow rate is zero. In this example there was no time dependence to the applied voltage and current.

Practical Industrial Examples of Plasma Application

There are an increasing number of references in the scientific literature giving examples of laboratory scale plasma catalysis (12, 13) or industrial processes being investigated at semi-industrial level (14). However, there are currently no known large-scale industrial applications that combine plasma with a catalytic material. The following two examples show the scaled production of chemicals using plasma excitation.

Ozone Production

Plasma is used to generate ozone (O3) industrially for applications including cleaning, disinfection, deodorisation and sanitisation. Commercially this is done using ultraviolet (UV) light at a wavelength of 185 nm, electrolytically or via a DBD plasma depicted in Figure 5. Siemens, Germany, were the first to use plasma for an industrial application in 1857 to produce O3 (15). This plasma is generated as a large number of statistically distributed micro-discharges between electrodes where the potential is insufficiently high to create an arc. Diatomic oxygen is broken down through interaction with the electrons produced within these micro-discharges that have sufficient energy to split the O2 double bond. These newly separated O atoms then combine with other diatomic O2 molecules to produce O3. The micro-discharges are individually only present for a few nanoseconds each and the number of micro-discharges generated is dependent upon the gap between the electrodes, the humidity and pressure of the air, the properties of the dielectric barrier and the characteristics of the electrical supply (16). O3 production can be as high as 100 kg h–1 from a horizontal honeycomb reactor and is closely related to the specific voltage and frequency applied, with typical voltages being 7–30 kV and frequencies between 50–1000s Hz.

Fig. 5.

O3 is formed by passing oxygen through an electrical discharge that is diffused over an area using a dielectric to create a corona discharge

O3 is formed by passing oxygen through an electrical discharge that is diffused over an area using a dielectric to create a corona discharge

Acetylene Production

Thermal plasma has been used to produce acetylene (C2H2) since the 1940s in the Huels process. The original Huels plant used the low-boiling components of the motor fuel industry as raw material; however, a wide range of hydrocarbons including natural gas were shown to be suitable as process feed stocks. The equilibrium formation of C2H2 is characterised by the requirement of very high temperatures, around 3000°C. Therefore, the reaction gas should be rapidly cooled by liquid water spray injection downstream of the plasma reaction zone to avoid formation of solid carbon that is a thermodynamically preferred product between about 1000°C and 2500°C. This fast quenching prevents decomposition of the C2H2 formed in the plasma. Because the formation of C2H2 from methane (CH4) is strongly endothermic, relatively large amounts of energy are required. In industry the best energy performance was shown by the DuPont process (a modification of the Huels reaction) with the specific energy consumption 8.8 kWh kg–1 of C2H2 (17).

More recently different methods of C2H2 formation in plasma were compared (18) and it was found that pulsed spark discharges gave the highest C2H2 yield (54%) with 69% of CH4 conversion in a pure CH4 system. It was suggested that the main disadvantage of the plasma method of C2H2 production was the fact that excited species reacted with formed C2H2 and decomposed it to undesirable byproducts. As a result, the energy costs for CH4 conversion and C2H2 formation increased with CH4 conversion percentage and were found to be best in pulsed spark discharges (highest CH4 conversions 18–69%). A Korean group reported the energy cost in its systems as 9 kWh kg–1 C2H2. They calculated their theoretical minimum energy requirement as 4.03 kWh kg–1 C2H2 based on heat of reaction (19). In this publication, it was also found that hydrogen in the plasma counterintuitively increased the selectivity to C2H2 in the process, from a respectable 70% to well above 90%. The paper claims that it is possible to decrease the specific energy consumption to 6 kWh kg–1 C2H2 in the low-temperature arc with the application of argon as the recycling reactant. It is close to the theoretical limit of 4.03 kWh kg–1 of C2H2, which is 50% of the energy required for the DuPont process. These numbers show that there is space to improve energy efficiency of the industrial production of C2H2 by plasma methods.

Plasma Catalysis

It is possible to introduce a catalyst after a plasma zone and transform some of the still excited species present in the gas phase into products at the surface of the catalyst, known as PPC (20).

Plasma can also generate activated species at a surface as well as in the gas phase. There are many terms used to describe the use of a catalyst and a plasma together, including plasma catalysis, plasma enhanced catalysis, plasma assisted catalysis, plasma driven catalysis (21); this IPC uses the excitation of both the catalyst surface and the gas phase reactants to effect chemistry through reactions between:

  • Excited gas phase species interacting with a non-excited catalyst surface

  • Excited surface species interacting with non-excited gas phase reactants

  • Excited gas phase interactions with excited surface species.

These plasma and catalyst interactions can be thought of as occurring between one or more excited states and can be representative of both the surface and the gas being excited. An electrically induced surface potential of a material is equally defined as a catalyst as a material that satisfies more traditional thermal catalysis ideals of adsorption and reaction, such as the use of glass beads for CH4 conversion (22).

In addition to the above interactions, further thought should be given to the interdependency of the plasma formed on the properties of the catalytic material present, and vice versa where the plasma will have an impact on the properties of the catalyst surface, including even impacting upon the physical morphology of the material. These synergistic interactions have been proven to offer a different mechanism for chemical reaction when the excitation comes from electrical plasma rather than thermal means (23). A study of plasma activated catalytic (palladium/aluminium oxide (Pd/ Al2O3)) CH4 oxidation, conducted in a synchrotron beamline, concluded that the Pd nanoparticles are heated within the plasma but the temperature of the nanoparticles remains lower than that required to initiate the thermal CH4 oxidation reaction. Thus, an alternative reaction mechanism with a lower activation barrier must be taking place (24).

There have been many studies looking at the application of plasma catalysis for chemical synthesis. In particular, small molecules that are difficult to activate using more traditional thermal methods, such as CO2 and CH4 lend themselves towards activation using plasma techniques (13). A special issue of the journal Catalysts titled ‘Plasma Catalysis’ was recently published including papers covering the application of plasma catalysis for CO2 splitting (25), ammonia (NH3) synthesis (26) and CH4 reforming (27).

A key challenge for plasma catalysis is to design a reactor that is suitable to house a plasma and a catalyst that has low backpressure, but retain good catalyst and gas-plasma interaction similar to the ceramic monolith widely used in automotive emission treatment. Uytdenhouwen et al. (1) identify power, pressure and gap size in a reactor as key process parameters for utilisation in design of plasma reactor and then go on to discuss their effect on CO2 disassociation in a DBD microreactor.

Mizuno (28) describes multiple approaches to tackle this problem in his review: including (a) micro-discharge plasma with the process gas flowing through a catalyst coated metal plate with very narrow (micron) gaps to improve the catalyst-gas interaction, (b) a metal mesh DC powered electrode in front of a ceramic monolith and a packed bed alternating current (AC) electrode at the rear of the ceramic monolith in order to introduce surface streamers along the monolith channels and (c) a sliding three-electrode DBD system combining a negative AC electrode, a DC electrode and a ground to create more homogeneous and widely dispersed surface plasma (29).

Plasma catalysis has been studied for automotive emission control in both the PPC and IPC configurations. Some of the first studies were conducted using a packed bed for nitrogen oxides (NOx) control (30, 31) and proposed a two-stage process whereby the nitric oxide (NO) was oxidised to nitrogen dioxide (NO2) which was then subsequently selectively reduced over the catalyst by the hydrocarbons present. This system has also been proposed as a pre-particulate filter plasma reactor to attain additional particulate matter oxidation benefits from the increased NO2 generated by the plasma (32). A successful demonstration of a combined plasma and catalytic system has also taken place for CH4 removal from dual fuel engines at low temperatures (33).

Plasma Surface Treatment

If the plasma can generate excited species at a surface, then it follows that it should be possible to change the surface by plasma treatment of a material. One area where plasma surface treatment has garnered significant interest is in the treatment of plastics and polymers. Plasma treatment of plastics and polymers can have a significant effect on the chemical and physical properties of the materials, these changes occur rapidly, often within seconds (34), through the following proposed mechanisms (35):

  • Etching and stripping surface material: plasma which reacts with the surface to clean it of contaminants, for example an oxygen plasma, or in more extreme cases such as with the inclusion of tetrafluoromethane (FreonTM (Chemours, USA)) to etch the surface itself. This happens through breaking of the polymeric covalent bonds due to the bombardment of the surface with high energy particles

  • The partially destroyed polymers can react with other similarly separated units and form crosslinks, thus extending polymeric chains and mechanically stabilising the surface

  • The surface can become functionalised by including a functional molecule in the plasma gas.

Examples of surface functionalisation include the incorporation of hydroxyl groups (OH) from humidity present in the plasma or N fixation using an NH3 plasma. The introduction of polar groups such as hydroxyls allows for a significant improvement in the wettability of polymers and therefore has significant advantages for the printing and adhesives markets. This is a complex area and large bodies of work have been produced documenting the effects of different types of plasma on different polymeric materials and summarised elsewhere (34, 36).

Increasingly plasma has been used for vapour deposition (VD) processes applying a uniform thin film coating to a material (37). These films are usually within the nanometre thickness range and used for modification of optical, chemical, electronic, physical and decorative properties of the materials. The methods of plasma application for physical VD include sputtering, ion plating and cathodic arc deposition. Sputter deposition involves deposition onto a substrate of a molecule previously vaporised from a target. The target is vaporised through the mechanism of momentum transfer from gaseous ions accelerated from a plasma. The plasma ion plating process uses the material vaporised from a target (by whichever method is suitable) and bombards the depositing film with molecules produced from a reactive gas plasma as the film is deposited in order to change the properties of the depositing film. Arc VD occurs when an electrode is vaporised through the application of high current across a biased cell with the vaporised molecules being accelerated towards, and deposited on, the polarised substrate.

These treatment and deposition techniques require an understanding of physics and electrical engineering to generate, measure and optimise a plasma in order to effect the chemical change on, or within, a surface.

Plasma for Catalyst Preparation

It has been reported that using plasma as a preparative technique can improve catalyst dispersion, increase metal-support interactions and change metal particle morphology, which in turn can lead to improved catalytic activity and stability (38). An example of improved activity from plasma preparation is the supported nickel catalysts used for steam reforming; using a DBD plasma reactor to decompose the precursors such as nickel nitrate, for catalyst preparation it is possible to increase the proportion of the (111) Ni facets which show enhanced performance and coke resistance.

Fig. 6.

A catalytic DBD plasma reactor used for the demonstration of CH4 removal from an engine exhaust

A catalytic DBD plasma reactor used for the demonstration of CH4 removal from an engine exhaust

Another example of plasma used for catalyst preparation is the use of DBD plasma instead of the standard thermal calcination for silicon dioxide (SiO2) supported cobalt materials for Fischer Tropsch synthesis. Li et al. found that plasma prepared materials had enhanced activity and a greater yield of heavy hydrocarbons when compared to the thermally calcined materials. This performance was attributed to the measured increase in Co dispersion, smaller Co(II,III) oxide (Co3O4) cluster size and more even Co distribution. A byproduct of the plasma preparation is the claim that this route can be a ‘greener’ method of preparing materials: using a low temperature electron reduction instead of using H2 as a reductant removes the need for dealing with H2 in the process (39).

Another advantage of the low bulk temperature plasma treatment as a preparative technique in comparison to a standard thermal treatment is the ability to remove precursors without inducing the detrimental changes that are associated with the temperatures normally required to oxidise the precursor molecules. An example of where this is useful is in the preparation of zeolite materials through plasma template removal. Liu et al. have conducted the removal of zeolite templates at around 125°C using a DBD plasma technique with O2. They identified the two major reactions taking place as the dissociation of template molecules by active species such as electrons or excited O2 and oxidation by excited O2 or O3 molecules (40). When directly comparing template removal from ZSM-5 zeolite by thermal and plasma techniques Liu et al. found that the rate of removal was approximately eight times higher using the DBD plasma method (41).

Plasma for Catalyst Modification

Using a DBD system AZO Materials, UK, reported differences in the temperature and in the intensity of the peaks resulting from temperature programmed reduction of magnesium oxide (MgO) supported Ni catalysts compared to non-plasma treated materials (42). The differences were attributed to Ni particle morphology and dispersion.

Zhu et al. also treated supported Ni catalysts with DBD plasma and found an increase in the catalytic activity and stability for the partial oxidation of CH4 (43). The scanning electron microscopy (SEM) images support a case for enhanced dispersion and increased interaction between the 10% Ni and the Al2O3 support. As well as a measurable increase in catalytic activity (3–5%) they also report a reduction in the formation of C around the Ni. This is consistent with a change in the Ni particle morphology towards having more (111) facets, as observed by others during the plasma preparation of catalysts (38).

Plasma for Material Regeneration

The literature related to the plasma application for material regeneration is limited. The current state and perspectives of plasma applications for catalyst regeneration was discussed in a recent review (44). Plasma regeneration was successfully applied for the reduction of oxidised catalysts and removal of poisons and C deposits. The largest advantage of plasma is that it allows catalyst regeneration to be performed at temperatures lower than those of typical thermal regeneration. The supply of gaseous reactive species and alteration of the surface structure to a more energetic state were identified as prerequisites of successful low-temperature regeneration and it was also shown that plasma can supply heat in a more cost-effective way than conventional thermal treatment. The energetic species produced in non-thermal plasma can initiate diverse reactions and open up or enhance reaction pathways other than those expected for equilibrium chemistry. As thermal regeneration can result in catalyst sintering and thus a reduced number of active sites, plasma is a viable alternative to thermal treatment.

The advantage of plasma regeneration was shown in the recent work that has been undertaken at the University of Central Lancashire, UK, studying deactivated coked zeolite regenerated with the application of different techniques including thermal, MW plasma and DBD discharge plasma (45). This work showed that plasma not only removes C from the deactivated catalyst but increases the activity of the catalyst significantly. Toluene disproportionation was used as the probe reaction in this study. Unlike thermally regenerated catalysts the material regenerated by plasma shows improved catalyst performance and the activity of the regenerated catalyst is even higher than that of the virgin material. Characterisation methods including pyridine and collidine infrared studies, NH3 temperature programmed desorption and solid state nuclear magnetic resonance were used to explain the changes in catalytic activity. Results showed MW plasma regeneration extended the catalytic life of zeolite due to the destruction of Brönsted acid sites caused by dealumination, without loss of crystal structure. In the toluene disproportionation reaction, this reduces the amount of cracking which occurs, subsequently leading to less coke deposition and therefore an extended catalytic life.

Cross Disciplinary Challenges

A key challenge for the field of plasma catalysis is to assemble teams with the relevant complementary skills in electronics, physics, engineering and chemistry together to gain an understanding of the system in order to produce the desired technological progress. A subset of this challenge is the communication between the different disciplines; involving not just the different language used, but also the models derived to express the concepts and understanding of processes which are often not ideally accessible to other branches of science and engineering.

As well as the broader topic of reducing barriers for communication, there are key challenges remaining within each of the disciplines required, including but not limited to:

  • Electronic engineering: nanosecond pulsed plasma sources tend to provide the most energetic waveforms for plasma systems. For that reason, it is critical that solid state nanosecond (5–100 ns) high voltage switches with high frequency capabilities are to be commercially available. Other important components in order to improve electrical and chemical efficiency are high dielectric strength materials with different dielectric constants (6) and reduction of electromagnetic interference (EMI) presented in such fast-rising waves

  • Physics: determining the appropriate signals to measure to obtain information that can be analysed in new ways. Analytical techniques borrowed from nuclear physics may be appropriate, for example pulse height spectrometry and more familiar optical spectrometry techniques can and should be used to characterise the species that are generated in real time. The outstanding challenge is to map out the outcomes and determine what they mean in the very large potential experimental space that presents itself in these systems

  • Engineering: scaling up a plasma reactor that requires both an interelectrode gap and sufficient gas-plasma-surface interaction to take advantage of the plasma and catalyst synergies, while simultaneously having a low pressure drop

  • Chemistry: knowledge around the catalogue of catalysts for thermally activated heterogenous reactions is not fully valid for plasma activated catalytic processes and a relevant body of knowledge does not currently exist. An additional challenge is to address the limited penetration depth of plasma into structures (46) and to develop methods for treating coated components including monolithic low pressure drop structures

  • Data analysis and modelling: large amounts of data generated from multiple variables increases the complexity of the experimental design and interpretation. This therefore makes the dataset an obvious candidate for supervised machine learning with input expertise from each technical discipline.

Concluding Remarks

It is necessary to understand physics to be able to correctly use electronic engineering to generate the requisite plasma to react with the designed catalyst in order to affect the desired chemistry. Whether this chemistry change is in the gas phase, at the surface of a catalyst or within the surface, there are clearly a large number of challenges to be faced when considering the large experimental space brought by the additional variables from non-chemistry fields. The general challenges for plasma catalysis can be expressed in broad terms of, for instance, increased efficiency, targeted catalysts or improved scaling of reactors; however within each of these targets underlying incremental improvements are required from additional scientific and engineering disciplines to achieve this.

Therefore, a key challenge for the field of plasma catalysis is to assemble teams with the relevant knowledge and skills in their own area of expertise that can work together and communicate ideas to initiate and progress the technology. Within Johnson Matthey there are teams with expertise across all the relevant disciples who can overcome these barriers and work successfully together to realise the potential of plasma catalysis on a commercial scale.

The Authors


Peter Hinde obtained his PhD in photocatalysis from The University of Bradford, UK, and has been working for Johnson Matthey at their central research and development facility supporting multiple global business units since 2004. He currently leads Johnson Matthey’s plasma catalysis research group.


Vladimir Demidyuk received his MSc and PhD degrees in Physical Chemistry from Moscow State University, Russia. For the past 15 years he has been working in academia and in industry in the UK. He joined Johnson Matthey, Sonning Common, UK, in 2015. His research interests include the application of plasma technologies for industrial purposes.


Alkis Gkelios gained his Masters in Electrical Engineering from the University of Patras, Greece in 2014. He has been working on high voltage electronics and plasma generation systems. For the last three years, he has been working for Johnson Matthey. His interests are focused on efficient plasma production and control for industrial process.


Carl Tipton gained his PhD in Nonlinear Physics from the University of Manchester, UK, in 2003. Since then, he worked as a Development Physicist at Tracerco, UK. He is currently a Measurement Engineer at Johnson Matthey, Chilton, UK, where he works to optimise Johnson Matthey industrial processes.

By |2020-03-19T13:35:16+00:00March 19th, 2020|Weld Engineering Services|Comments Off on Plasma Catalysis: A Review of the Interdisciplinary Challenges Faced

Insights into Automotive Particulate Filters using Magnetic Resonance Imaging

Home > Journal Archive > Insights into Automotive Particulate Filters using Magnetic Resonance Imaging

Johnson Matthey Technol. Rev., 2020, 64, (2), 165

Introduction

Particulate filters, comprising so-called wall-flow filter substrates, are of increasing importance in reducing pollutant emissions from vehicles to the levels required by legislation. Early legislation addressed the emissions of carbon monoxide, hydrocarbons and nitrogen oxides and the removal of these pollutants was achieved using FTM catalysts. Later, when PM emissions from diesel vehicles came under scrutiny, wall-flow filter devices were added to the emissions control system: often these were uncatalysed extruded cordierite or silicon carbide filter monoliths, though sometimes a catalyst was incorporated on the filter to widen their operating window. Wall-flow filters differ in their operation from FTMs, since adjacent channels are alternately blocked meaning exhaust gas must pass through the porous monolith wall to flow from inlet to outlet; in this way PM is deposited on the inlet channel walls. Previous PM legislation was based on particulate mass emissions, however, more recently the legislation, such as in Euro 6 (1), has turned to address particulate number emissions. Due to this change in emphasis both diesel and gasoline vehicles now require filter systems to achieve compliance. In gasoline applications, a GPF is typically a wall-flow filter, similar to that used for diesel, with a catalyst coating applied: by combining the catalyst and filter devices, multifunctional emissions control systems have also been developed resulting in reduced packaging volume.

In the manufacture of catalysed filter devices, first a catalyst coating is applied to the monolith substrate in the form of a slurry. The subsequent coated monolith is then dried and finally calcined to fix the catalyst coating. The drying process is both energy intensive and known to influence the final metal distribution within the catalyst throughout the catalysed filter and therefore also influences filter performance on a vehicle. Indeed, it has previously been shown (2, 3) that non-ideal drying in FTMs can result in the macroscopic redistribution of the catalyst. In the work of Vergunst et al. (3) and Wahlberg et al. (4) it was observed that the metal phase, when not bound to the catalyst support, will migrate to the surface where evaporation occurs. This results in varying degrees of inhomogeneous catalyst distribution depending on the method of drying used. Finally, enrichment of the catalyst in the washcoat layer of ceramic monoliths during drying has been observed using MRI (5). Traditionally, drying is studied using gravimetric methods, humidity and temperature measurements from which the drying kinetics can be determined (6). Although these techniques are well established in both industry and research, they are somewhat limited in that they are only able to provide macroscopic measurements and the process itself must be treated as a ‘black box’. Spatially resolved information has most often been obtained by sectioning the sample and then weighing the individual components. To gain a greater understanding of the intrinsic water migration characteristics during drying, in the first part of this work we use MRI to image the time-resolved water distribution during the drying. The data for the filter are compared with the same data acquired during drying of the related FTM.

The second part of this work uses MRI methods to measure the gas velocity in the channels of a clean and particulate loaded filter. During the operation of an automobile, PM-laden exhaust gas passes through the particulate filter and the particulate or soot is deposited inside. However, the micro- and macroscopic distribution of this soot deposition impacts the subsequent filtration behaviour, pressure drop, regeneration behaviour and ultimately the useful lifetime of the filter. Hence, a complete understanding of the filtration process is needed in order to optimise the function of particulate filters. The regeneration behaviour of particulate filters is coupled with the gas fluid dynamics. Both heat transfer and the mass transfer of oxidative species (oxygen for active regeneration and nitrogen dioxide for passive regeneration) are affected by the gas flow fields and will impact the efficacy and safety of the regeneration process. Large thermal gradients within the filter can be formed if the soot distribution is non-uniform that may damage the filter. While modelling and macroscopic measurements of these effects have been performed by many authors, there is a lack of experimental work focusing on the relationship between the filter structure, the gas transport and the perturbations of these by the loaded soot. As with drying, the range of techniques available to non-invasively measure flow in opaque filter systems is very limited, and most studies have used models (710). Magnetic resonance (MR) can provide spatial information on the effect of the PM on the filter operation.

MR techniques have gained prominence in chemical engineering research as they provide a non-invasive method of studying the chemistry and dynamics of a range of opaque systems (11). They are particularly useful for the study of porous media such as catalysts (1216), construction materials (1719) and pharmaceuticals (20, 21). Such techniques have also been used to provide information on the drying mechanism in a range of other applications such as detergent powders (22, 23), dehydration and preparation of foodstuffs (2426), evaporation from contaminated surfaces (2729) and fired-clay brick at elevated temperature (30). Applications of MR to drying and sorption have been covered in the review by Koptyug (31). While MR has previously been used to study drying and active component distribution in FTMs by Koptyug and coworkers (5, 15), no studies have looked at the drying process in particulate filters to date. This study investigates how the structure of the filter substrate influences the water migration behaviour during drying. MR techniques are used to characterise water migration behaviour in a wall-flow filter and FTM substrate under identical drying conditions. While MRI of liquids and their transport is common, imaging studies of nuclear magnetic resonance (NMR) active gas flows are relatively few. This is mainly due to challenges associated with the low signal-to-noise ratio (SNR) presented to the experimentalist. While hyperpolarised gases (for example, xenon-129) offer a significant boost to the SNR, they are costly and unsuitable for studying porous materials. The first application of imaging thermally polarised gases was demonstrated by Koptyug et al. (3234) who acquired two-dimensional (2D) velocity images of hydrocarbon gases at atmospheric pressure flowing through a cylindrical pipe and alumina monoliths of different channel geometries. In the monolith studies (32), from the spatially resolved profiles of the axial component of the velocity vector on the individual channel scale, information pertaining to shear rates and entry lengths were obtained which enabled useful insights into the mass transfer between the bulk gas flow and the porous channel walls in the monolith to be made. Codd and Altobelli (35) have also shown the application of using thermally polarised gases as a probe of the structure of porous materials. Ramskill et al. (36) performed an early study using sulfur hexafluoride gas to image flow profiles inside a clean emissions control filter. They implemented compressed sensing (CS) techniques, which enable spectra and images to be reconstructed with sufficient accuracy from relatively few data points and allow a reduction in data acquisition times (3741). Here we show a continuation of this study with more commercially relevant samples and with particulate loading.

This study highlights two ways in which MRI can be used to gain insight into automotive particulate filters. First a comparison of water migration during drying within a wall-flow particulate filter is reported and compared with the analogous process occurring within a FTM. The mechanism of water migration during drying will influence catalyst distribution in the manufactured filter and valuable insights are gained by being able to image this process in more than one dimension. Second, MRI is used to image gas fluid dynamics within a GPF both in the clean state and following two stages of soot loading. From these flow profiles, permeabilities as a function of space and time are predicted. Using MRI to gain a greater understanding of automotive particulate filters could eventually lead to improved catalyst properties, filtration efficiency and more efficient and improved manufacturing.

Principles of Magnetic Resonance Imaging

This section provides the reader with a brief introduction to the principles of MRI but the interested reader is directed to the texts by Callaghan (42) and Haacke (43) and review articles by Mantle and Sederman (44) and Caprihan and Fukushima (45) for further detail.

When NMR active nuclei (such as 1H or 19F) are placed in an external magnetic field (B0), the nuclei will precess at a characteristic frequency known as the Larmor frequency (ω0) as given by Equation (i):

(i)

where γ is the gyromagnetic ratio which is characteristic of the nuclei under observation. Considering the system in the rotating frame at reference frequency, ω0, spatial dependence of the precession frequency is achieved by applying a spatially-dependent magnetic field gradient (G) and can be expressed as follows (Equation (ii)):

(ii)

where ω(r) is the precession frequency at position vector r. The precessing nuclei in a volume element induce a voltage in the receiver coil and the complex NMR signal, S(t), is detected in the time domain. The time-dependent signal is then Fourier transformed into the frequency domain where it is represented by the spin-density function, ρ(r). The Fourier conjugate relationship between the time and frequency domain of the NMR signal is shown below in Equations (iii) and (iv):

(iii)

(iv)

From this, a spin-density map, i.e. an image, is obtained by taking the modulus of the complex function, ρ(r). In the case of the drying experiments the spin-density map provides a spatially-resolved measurement of the water content within the monolith sample.

Magnetic field gradients can be used to make the NMR signal sensitive to nuclei displacements in addition to position. This is achieved by first encoding the nuclei positions by applying a given magnetic field gradient (G) for time (δ), then decoding by applying the same magnetic field gradient in the opposing direction after an evolution time (Δ). Any static nuclei will be unchanged but any moving nuclei will create a phase shift in the NMR signal, ΔΦ proportional to their displacement Δr (Equation (v)):

(v)

Hence the gas velocity can be measured by consideration of the signal phase at each spatial location.

Materials and Methods

Comparison of Drying in a Particulate Filter and a Flow Through Monolith

A cordierite wall-flow particulate filter and a cordierite FTM were used. These are typical of substrates used commercially and the relevant properties of the two samples are listed in Table I. Porosity and mean pore size measurements were made using an AutoPore IV system mercury porosimeter (Micromeritics Instrument Corporation, USA). Drying of pure water from the respective substrates has been investigated under identical conditions i.e. with air at 20 l min–1 ± 2 l min–1 and temperature of 19.5°C ± 0.5°C.

Table I

Properties of the Diesel Particulate Filter and FTM Substratesa

Wall-flow FTM
Material cordierite
Length (L), mm 75
Core diameter, mm 26
Channel hydraulic diameter, mm 1
Substrate porosity (ɛ), % 48±4 24±2
Mean pore size, μm 13.8±7.8 2.8±1.1
Water content (mc), g 6.1±0.1 2.8±0.1

To saturate the samples, the substrates were immersed in deionised water for two minutes and then shaken to remove any water blocking the channels. A schematic of the experimental set up used is shown in Figure 1. The filter substrate was held within a PERSPEX® cell (Perspex International Ltd, UK) above an air distributor plate used to produce a uniform flow of air over the cross-sectional area of the filter. The relative humidity (RH) and temperature of the air flow at the inlet and outlet of the drying cell were recorded at 10 s intervals over the course of the process using a Humidiprobe (Pico Technology, UK). Temperature and RH measurements were recorded with an accuracy of ±0.5°C and ±2%, respectively. These measurements allow the total uptake of moisture by the air to be determined and thus the drying rate can be calculated through conservation of the total water mass (36, 46).

Fig. 1.

Schematic of the experimental setup for the drying experiments. A = compressed air line, B = pressure regulator, C = rotameter, D = air distributor plate, E = wall-flow filter or FTM substrate sample, F = imaging region, G = MRI spectrometer and H = magnetic field gradients

Schematic of the experimental setup for the drying experiments. A = compressed air line, B = pressure regulator, C = rotameter, D = air distributor plate, E = wall-flow filter or FTM substrate sample, F = imaging region, G = MRI spectrometer and H = magnetic field gradients

Effect of Soot Loading on Gas Fluid Dynamics in a Gasoline Particulate Filter

A cordierite GPF sample was prepared in the laboratory for this study. The cordierite substrate (55% porosity) was coated with a Pd/Rh alumina three-way catalyst typical of commercial catalysts used for GPF applications. The properties of the sample are shown in Table II, with the porosity and pore size of the catalyst coated filter listed for the front, middle and rear of the sample.

Table II

Properties of the GPF Samplea

Bare filter Catalyst coated filter
Material cordierite
Length (L), mm 145
Core diameter, mm 25 mm for soot loading, 6 mm for MRI
Channel hydraulic diameter, mm 1
Substrate porosity (ɛ), % 58 28.1, 25.9, 29.2
Mean pore size, μm 21 19, 16, 18

Samples were soot loaded using a 2 l, four-cylinder gasoline direct injection (GDI) turbocharged engine, and then subsequently removed and transferred to a different sample holder for the MRI flow experiments. Such engines are typical of current passenger automobiles. The engine was run at 2100 rpm, producing a torque of 60 Nm. Three 25 mm diameter cores were bored from the centre of the GPF, each acting as a sample for the soot loading. The filter sample to be loaded was held downstream inside the exhaust manifold. The manifold was surrounded by a furnace, allowing the filter to be held at different temperatures. A temperature of 300°C was chosen as representative of real-world gasoline exhausts; the temperature was measured using a thermocouple inserted into the manifold upstream of the filter. Two pressure transducers were placed either side of the filter sample, allowing measurement of the pressure drop during the loading process. Measurements were made at 180 ms intervals at both transducers. The pressure readings were subtracted and averaged over 2 min intervals to give the transient pressure drop. Three particulate loading protocols were used:

  • Protocol I: no soot loading

  • Protocol II: normal running of the engine for 50 min

  • Protocol III: Protocol II followed by 10 min of accelerated soot loading achieved by delaying the fuel injection by a crank-shaft angle of 50 degrees.

The backpressure profiles recorded for Protocols II and III are shown in Figure 2.

Fig. 2.

Pressure drop measurements for the GPF sample subject to particulate loading Protocols II and III

Pressure drop measurements for the GPF sample subject to particulate loading Protocols II and III

Magnetic Resonance Characterisation

For the drying experiments, the MR experiments were performed using a 2 Tesla (85 MHz for 1H) horizontal bore magnet controlled by an AV spectrometer (Bruker Corporation, USA). An 85 mm radio frequency (rf) coil tuned to a frequency of 85.1 MHz was used for excitation and signal detection and spatial resolution was achieved with magnetic field gradients with a maximum strength of 10.7 gauss cm–1. Three MR techniques were used as follows.

  • NMR spectroscopy was used to provide a quantitative measurement of the bulk water content through calibration with gravimetric measurements. Due to the short deadtime between rf excitation and detection, negligible relaxation weighting is associated with the spectra and hence they are directly proportional to the spin density of water. A recycle time of 3.5 s and eight scans for signal averaging were used, resulting in a total acquisition time of 0.5 min per spectrum

  • 2D images were acquired over the course of drying using the rapid acquisition with relaxation enhancement (RARE) pulse sequence (47). Images were acquired in the yz plane with a slice width of 10 mm, an in-plane field-of-view (FOV) of 80 mm × 30 mm and a data matrix symbol size of 32 × 32, giving an in-plane pixel resolution of 2.5 mm px–1 × 0.94 mm px–1 in the read (z) and phase (y) directions respectively. A RARE factor of four with eight scans were used, allowing acquisition of a full image in 3.5 min

  • One-dimensional (1D) profiles in the axial (z) direction were acquired using a spin-echo profiling sequence that integrates the spin density along the x and y directions (42). A FOV of 80 mm in the z direction and a matrix size of 128 points was used, giving a spatial resolution of 0.625 mm px–1. Eight scans were used for signal averaging, giving a total acquisition time of 0.5 min. The echo time between excitation and acquisition of the NMR signal was 10 ms, giving a maximum error of 3% for the relaxation times present in the system.

The MRI method used to measure gas velocity in the filter samples is described fully in Ramskill et al. (36). In the present study SF6 has been chosen as the NMR active gas to be used for velocity imaging due to its favourable MR properties in comparison with other potential candidates such as the hydrocarbon gases (48). Eleven images were acquired along the length of the GPF samples, each with a slice width of 6 mm. An SF6 gas pressure of 5.0 barg ± 0.1 barg and mass flow rate of 16 g min−1 was used for each sample. Axial velocity profiles were acquired for the GPF samples after all three soot loading protocols. The mean volume flow for each sample agreed with the value calculated from the mass flow rate to within 8.5%. The through-wall velocities were calculated for each based on the gas mass balance. Velocity profiles inside the inlet channels were extracted from the MR velocity images through the mid-point of the channels parallel to the filter wall.

Results and Discussion

Comparison of Drying in a Particulate Filter and a Flow Through Monolith

2D images in the yz plane of the wall-flow filter and FTM have been acquired over the course of drying as shown in Figure 3 and Figure 4 respectively. The initial water distribution along the two substrates has been determined by integrating the signal intensity in the first image of each sequence (Figure 3(a) and Figure 4(a)) in the y-direction. These data are shown in Figure 5. A uniform wetting of the channels of both substrates is seen with an area of relatively high signal intensity corresponding to the higher moisture content contained in the plugs of the wall-flow filter substrate. Figure 6 shows the moisture content at the three radial positions (marked as A, B and C on Figures 3 and 4) as a function of drying time. Apart from the longer drying time associated with the wall-flow filter substrate, the behaviour of the water saturation at the three positions across the monolith as a function of time show very similar behaviour. This result confirms that the 1D profiles in the axial (z) direction are sufficient to be able to study the drying mechanism in the wall-flow filter and FTM and thereby allow the process to be studied at a higher temporal resolution than would be permitted using 2D MRI.

Fig. 3.

2D images over the course of the drying of a wall-flow filter at: (a) 3.5 min; (b) 10.5 min; (c) 17.5 min; (d) 24.5 min; (e) 30.5 min and (f) 36.5 min. The signal intensity has been normalised relative to the maximum signal intensity in the filter substrate. Air flow is from the bottom with a volumetric flow rate of 20 l min–1. These images were acquired with a FOV of 80 mm × 30 mm in the zy plane corresponding to a spatial resolution of 2.5 mm × 0.94 mm

2D images over the course of the drying of a wall-flow filter at: (a) 3.5 min; (b) 10.5 min; (c) 17.5 min; (d) 24.5 min; (e) 30.5 min and (f) 36.5 min. The signal intensity has been normalised relative to the maximum signal intensity in the filter substrate. Air flow is from the bottom with a volumetric flow rate of 20 l min–1. These images were acquired with a FOV of 80 mm × 30 mm in the zy plane corresponding to a spatial resolution of 2.5 mm × 0.94 mm

Fig. 4.

2D qualitative images over the course of the drying of a FTM at: (a) 3.5 min; (b) 7 min; (c) 10.5 min; (d) 14 min: (e) 17 min and (f) 20.5 min. The signal intensity has been normalised relative to the maximum signal intensity in the filter material. Air flow is from the bottom with a volumetric flow rate of 20 l min–1. These images were acquired with a FOV of 80 mm × 30 mm in the zy plane corresponding to a spatial resolution of 2.5 mm × 0.94 mm

2D qualitative images over the course of the drying of a FTM at: (a) 3.5 min; (b) 7 min; (c) 10.5 min; (d) 14 min: (e) 17 min and (f) 20.5 min. The signal intensity has been normalised relative to the maximum signal intensity in the filter material. Air flow is from the bottom with a volumetric flow rate of 20 l min–1. These images were acquired with a FOV of 80 mm × 30 mm in the zy plane corresponding to a spatial resolution of 2.5 mm × 0.94 mm

Fig. 5.

1D profiles of signal intensity in: (a) the wall-flow filter obtained though numerical integration of the 2D images in the radial (y) direction; (b) FTM obtained though numerical integration of the 2D images in the radial (y) direction

1D profiles of signal intensity in: (a) the wall-flow filter obtained though numerical integration of the 2D images in the radial (y) direction; (b) FTM obtained though numerical integration of the 2D images in the radial (y) direction

Fig. 6.

Relative saturation taken at three radial (y) positions (A, B and C) from the 2D images (Figure 3(a) and Figure 4(a)) plotted over the course of drying for: (a) the wall-flow filter; and (b) FTM. Drying appears to be uniform in the radial direction

Relative saturation taken at three radial (y) positions (A, B and C) from the 2D images (Figure 3(a) and Figure 4(a)) plotted over the course of drying for: (a) the wall-flow filter; and (b) FTM. Drying appears to be uniform in the radial direction

Figures 7 and 8 show the quantitative drying curves and rate of drying for the filter and FTM, respectively. The water content data shown in Figure 7(a) and Figure 8(a) are quantitative and obtained directly from the traditional measurements (RH and temperature measurement) and integration of the signal from NMR spectroscopy; the error bars were calculated based on the instrument sensitivities and the standard deviation of repeat measurements respectively. The rate of drying data shown in Figure 7(b) and Figure 8(b) are obtained from the time derivative of the data in Figure 7(a) and Figure 8(a). It is seen that for both the wall-flow filter and FTM the trends in drying behaviour appear similar; in particular, a slow falling rate period followed by a faster falling rate as drying proceeds. The only significant difference is during the induction period of both samples; the NMR data shows an increasing rate whereas the traditional measurements show a decreasing rate. This is attributed to the temperature dependence of the NMR signal. As the sensitivity of the NMR signal depends on the population difference of nuclear spin energy levels, as described by the Boltzmann distribution, a reduction in temperature causes an increase in the observed NMR signal. During the induction period, a temperature drop of up to 6 K is observed due to the heat of evaporation, resulting in an increase of up to 2% in the observed NMR signal and hence a reduction in the calculated drying rate. While this is negligible for the modest temperature changes during most of the drying process, the temperature drop during the induction period is sharp and the increase in NMR signal decreases the measured drying rate by up to 20%. Thus, from a quantification standpoint we have successfully been able to validate the MRI using simple humidity measurements; however, the MRI is able to give spatial information, as will be shown.

Fig. 7.

(a) Drying; and (b) rate of drying curves for the wall-flow filter as determined from the RH/temperature and MR measurements. An air flow rate of 20 l min–1 was used

(a) Drying; and (b) rate of drying curves for the wall-flow filter as determined from the RH/temperature and MR measurements. An air flow rate of 20 l min–1 was used

Fig. 8.

(a) Drying; and (b) rate of drying curves for the FTM as determined from the RH/temperature and MR measurements. An air flow rate of 20 l min–1 was used

(a) Drying; and (b) rate of drying curves for the FTM as determined from the RH/temperature and MR measurements. An air flow rate of 20 l min–1 was used

For the wall-flow filter and FTM respectively, Figures 9 and 10 show the time evolution of the 1D profiles during the drying process, with 15 mm slices extracted from the data to show the average water saturation time evolution at 20 mm, 40 mm and 60 mm from the front of the substrates. In the case of the filter substrate (Figure 9), it can be seen that drying proceeds uniformly in the axial (z) direction up to a critical point at which a developing drying front is present until the filter is dry. In the initial stages, the rate of drying at the three axial positions is the same until ~20 min, after which the drying front develops and the front of the filter will dry more quickly than the middle and back sections. Between 20 min and 35 min, the middle and back sections of the filter continue to dry at the same rate until the drying front reaches the middle section and begins to dry more quickly than the back section. Finally, the remaining water is removed as the drying front moves through to the back of the filter and is completely dry after 50 min. In contrast, as is seen in Figure 10, for the FTM the drying front propagates through the substrate from the very beginning of the drying process.

Fig. 9.

(a) Time series of 1D axial (z) profiles over the course of the drying of a wall-flow filter; (b) average saturation over three 15 mm slices centred at 20 mm, 40 mm and 60 mm along the length of the sample. Air flow is from the bottom with a flow rate of 20 l min–1. FOV in the z direction was 80 mm and data matrix with 128 points thus corresponding to a spatial resolution of 0.625 mm px–1. The development of the drying front at ~20 min can be clearly seen

(a) Time series of 1D axial (z) profiles over the course of the drying of a wall-flow filter; (b) average saturation over three 15 mm slices centred at 20 mm, 40 mm and 60 mm along the length of the sample. Air flow is from the bottom with a flow rate of 20 l min–1. FOV in the z direction was 80 mm and data matrix with 128 points thus corresponding to a spatial resolution of 0.625 mm px–1. The development of the drying front at ~20 min can be clearly seen

Fig. 10.

(a) Time series of 1D axial (z ) profiles over the course of the drying of a FTM; (b) average saturation over three 15 mm slices at 20 mm, 40 mm and 60 mm along the length of the sample. Air flow is from the bottom with a flow rate of 20 l min–1. The FOV in the z direction was 80 mm and data matrix with 128 points thus corresponding to a spatial resolution of 0.625 mm px–1. It is evident that a drying front is present from the start of drying

(a) Time series of 1D axial (z ) profiles over the course of the drying of a FTM; (b) average saturation over three 15 mm slices at 20 mm, 40 mm and 60 mm along the length of the sample. Air flow is from the bottom with a flow rate of 20 l min–1. The FOV in the z direction was 80 mm and data matrix with 128 points thus corresponding to a spatial resolution of 0.625 mm px–1. It is evident that a drying front is present from the start of drying

The most striking result from these studies is that whilst the spatially-unresolved data (Figures 7 and 8) might suggest similar drying characteristics between the wall-flow filter and FTM substrates, the spatially-resolved measurements obtained from 1D and 2D MR measurements (Figures 6, 9 and 10) reveal significant differences. These data therefore provide important information on the microscopic contributions to the drying process which must be reflected in any computational model. Furthermore, the differences seen imply that care must be taken when drying particulate filters, since water migration and transport could result in a final catalyst location and properties that are different from those expected.

In order to reconcile the predictions of the bulk and spatially resolved measurements, the structure and operation of each sample must be considered. FTMs operate through the axial flow of gas through the channels of the structure, with no flow through the substrate itself. Thus, drying occurs through evaporation at the interface of the channels and the substrate, with water migrating from the saturated pores to the surface through capillary flow. Hence the rate of drying is determined by the humidity gradient in the channels. Air entering the front of the FTM has the lowest humidity, and so the rate of evaporation is greatest at the front and decreases along the length of the filter. As water is removed from the front, the humidity gradient shifts down the channels and the drying front propagates until it reaches the end of the filter as seen in Figure 10(a), causing the slowly decreasing drying rate seen between 5 min and 18 min in Figure 8(b). Once most of the moisture has been removed from each axial position, the residual is controlled by diffusive mass transport in the porous medium and is slower, resulting in the quickly decreasing drying rate after 18 min in Figure 8(b). This is consistent with the findings of Koptyug et al. (3) for a segmented monolith, suggesting that axial capillary flow is not significant due to the larger pore sizes present in this study.

In particulate filters, alternate channels are blocked meaning that gas is forced to flow through the porous substrate in order to exit the filter, adding extra complexity to the drying process. The air first removes water in the largest pores in the substrate, so the air is humidified in the walls and not the channels, so there is no humidity gradient along the length of the sample. These pores form the paths of least resistance between channels, allowing air to flow through the filter walls and become fully saturated. For this air flow rate, the transverse velocity through the walls is expected to be reasonably uniform (36, 49), and so the drying rate in this regime is highly uniform across the filter. From Figure 9, this takes around 20 min, at which point the air becomes less than fully saturated and a humidity gradient can form, starting the drying front that propagates through the wall-flow filter between 20 min and 50 min in Figure 9(a). Figure 7(b) shows no change at 20 min, indicating that despite the change in drying mechanism for the wall-flow filter, the rate-limiting process has not changed and is a similar rate to the FTM shown in Figure 8(b). The rest of the moisture is removed via a similar mechanism to that in FTMs; the transition from evaporation limited to diffusion limited drying can be seen in Figure 9(b) for each axial position.

This additional mechanism in the drying of wall-flow substrates may affect the distribution of catalyst in the monolith following the drying process. The propagation of a drying front in the filter material will create similar issues as have been observed in FTMs by Vergunst et al. (3). However, the effect of the uniform drying regime is uncertain due to the limited number of drying studies on particulate filters. Spatially uniform drying methods such as microwave drying (50) produce more homogeneous distributions of the active phase in the monolith, but as with static and mobile air drying, the mobility of the catalyst and the solvent leads to non-uniform distributions. Limiting or ceasing this mobility, either through fast drying protocols, deposition-precipitation (51) or freeze-drying (3), can improve the homogeneity of catalyst distribution at the expense of economic viability.

Effect of Soot Loading on Gas Fluid Dynamics in a Gasoline Particulate Filter

The velocity profiles measured for the catalyst coated wall-flow filter are shown in Figure 11 for (a) axial and (b) through-wall velocities. Under no soot loading (Protocol I) the axial velocity profile resembles those seen previously (36), with the characteristic U-shaped form of the through-wall velocity profile. After Protocol II of soot loading, the axial profile has only changed slightly, with a more linear change in channel velocities and the cross-over of the inlet and outlet velocities occurring further forward in the filter. A more uniform through-wall velocity is seen, although there is still a parabola-like section in the filter centre. After Protocol III, the axial velocity profile shows a linear decrease and increase in the inlet and outlet channel velocity respectively. The corresponding through-wall velocity profile is much more uniform. Velocity profiles were then extracted from the MRI measurements for the central inlet channel of the sample, and show the gas velocity radially across the cross-section of the single channel to show the evolution of the radial flow profile at different axial positions with increasing soot load (Figure 12): axial positions of z/L = 0.14, 0.33, 0.52, 0.70 and 0.87, referred to as P1, P2, P3, P4 and P5 respectively. The profiles at P1 to P4 are close together for the soot-free sample (Protocol I) but become more spaced out as the soot loading increases. Some profiles show step-like features towards the channel edge, for example P5 in Protocol I, P3 in Protocol II and P1 and P3 in Protocol III. At the highest soot loading (Protocol III), the shape of the flow profiles at P4 and P5 are narrower than the expected paraboloid.

Fig. 11.

MRI measurements (markers) of: (a) the inlet and outlet channel velocities; (b) through-wall velocities for the GPF sample with Protocol I, Protocol II and Protocol III

MRI measurements (markers) of: (a) the inlet and outlet channel velocities; (b) through-wall velocities for the GPF sample with Protocol I, Protocol II and Protocol III

Fig. 12.

Axial velocity flow profiles for the GPF sample after loading: (a) Protocol I; (b) Protocol II; and (c) Protocol III. Profiles are shown at axial positions of z/L = 0.14, 0.33, 0.52, 0.70 and 0.87. Lines are only shown as a guide

Axial velocity flow profiles for the GPF sample after loading: (a) Protocol I; (b) Protocol II; and (c) Protocol III. Profiles are shown at axial positions of z/L = 0.14, 0.33, 0.52, 0.70 and 0.87. Lines are only shown as a guide

Using a numerical 1D model (7) and the MRI data, it is possible to extract information regarding the axial properties of the filter wall in the form of a permeability from Darcy’s law. The permeability calculated is an average value for the combined effect of the filter wall and the particulate. However, because the through-wall velocity is known from the MRI measurements, it is possible to calculate a spatial permeability in the axial dimension (z) of the filter. The permeability profiles are shown in Figure 13. Some variability in the permeability profile for the clean filter is observed, but the permeability is much more uniform after loading Protocol II and especially after Protocol III. The results here show that as the GPF operates, regions of higher through-wall velocity will have a greater loading of soot due to the correspondingly high mass flow of PM; this leads to the ‘self-correction’ effect predicted by Bensaid et al. (52). The axial velocity profiles inside the inlet channels (Figure 12) show changes at the filter rear with increasing soot load. The profiles become ‘narrower’ with lower velocity towards the wall. This is similar to the profiles observed by York et al. (53) also using MRI for high soot loadings in a diesel particulate filter and may be consistent with the development of a soot cake layer in these regions. These regions also correspond to the largest reduction in wall permeability (Figure 13). Two limitations of the MRI method are that it cannot quantify the soot loading and it cannot differentiate between different diameters of soot particle. However, other techniques such as gravimetric analysis or microscopy may allow these to be related to the MRI results in the future.

Fig. 13.

Simulated permeability profiles for the GPF sample for loading: (a) Protocol I; (b) Protocol II; and (c) Protocol III. Lines are only shown as a guide

Simulated permeability profiles for the GPF sample for loading: (a) Protocol I; (b) Protocol II; and (c) Protocol III. Lines are only shown as a guide

Conclusions

A range of MR and traditional techniques, namely RH and temperature measurement, have been applied to study wall-flow filter substrates used for PM vehicle emissions control. Drying of the filter material has been compared with an FTM using 2D RARE imaging to investigate the effect of the structure on the drying kinetics. Since little deviation in the radial drying profile was seen the problem was reduced to 1D in the axial (z) direction. This allowed increased temporal resolution that revealed different drying mechanisms are associated with the wall-flow and FTM substrates, attributable to differences in the physical structures of the two autocatalyst substrates. MRI and velocimetry has also been used to investigate the effect of PM on filter channel fluid dynamics. In the gasoline system studied, any inhomogeneities in the filter wall permeability are ‘self-corrected’ by the particulate loading over the initial hour of filter operation. MRI is a method that can indirectly visualise soot location in the filter, while also providing important information on the flow characteristics and substrate properties. The application of MRI to study filters used in vehicle emissions control has provided new insight into their manufacturing and operation. The greater understanding of the drying process could ultimately result in an improved and more efficient drying process, while great understanding of their operation can lead to improved final product performance, for example higher filtration efficiency.

Acknowledgements

The authors would like to thank the Engineering and Physical Sciences Research Council (EPSRC) for Cooperative Awards in Science and Technology (CASE) to N. P. Ramskill and J. D. Cooper. We would also like to thank Johnson Matthey, UK, for permission to publish.

The Authors


Jonathan Cooper has a Master’s degree in chemistry and recently completed his PhD at the Department of Chemical Engineering, University of Cambridge. His research focused on applying MRI methods to study the gas hydrodynamics in particulate filter systems. He now teaches chemistry at a school in London.


Nicholas Ramskill holds an MEng in Chemical Engineering from the University of Leeds, UK, and a PhD in Chemical Engineering from the University of Cambridge, UK. Nicholas completed his PhD at the Magnetic Resonance Research Centre in collaboration with industrial partner Johnson Matthey where his research focused on MRI studies of diesel particulate filters (DPFs). Subsequently, Nicholas was a Postdoctoral Research Associate in Cambridge where he worked in collaboration with Royal Dutch Shell, The Netherlands, on the development and application of MRI techniques to characterise enhanced oil recovery processes in the laboratory at representative reservoir conditions.


Andy Sederman is a Reader in Magnetic Resonance in Engineering at the University of Cambridge where he also gained his PhD in 1998. His research focus is on developing magnetic resonance techniques for application to engineering and materials. He has worked extensively in the area of velocity and transport measurement and methods to increase the imaging speed to be able to investigate transient systems, both by fast data acquisition and by utilising innovative reconstruction methods allied to data under-sampling. Areas of application for these methods have focused on single and multiphase flows, fluid flow in porous media and reaction and hydrodynamics in multiphase reactors.


Lynn Gladden is Shell Professor of Chemical Engineering at the University of Cambridge in the Department of Chemical Engineering and Biotechnology. She is recognised for her work on advancing magnetic resonance imaging techniques, originally developed for use in the medical environment, and using them in engineering research to gain a greater understanding of the physical and chemical phenomena that determine the performance of chemical processes and their resulting products.


Professor Athanasios Tsolakis has academic and industrial expertise in the field of low carbon energy carriers, environmental catalysts, combustion and pollutant control technologies. He works at the forefront of basic and translational research to improve fuel efficiency and reduce the environmental impact of the transportation and power generation sectors. Prior to his academic appointment at the University of Birmingham in 2005 he worked as a research scientist at Johnson Matthey in the design and characterisation of environmental catalysts for modern aftertreatment systems.


E. Hugh Stitt is a Scientific Consultant at Johnson Matthey, Chilton, UK. He is a Visiting Professor at the University of Birmingham; Fellow of the Institution of Chemical Engineers and a Fellow of the Royal Academy of Engineering. He has 30 years of industrial research experience across a variety of themes related to catalytic reaction engineering and catalyst manufacture with over 100 refereed publications.


Andrew York is a Research Manager responsible for the Emissions Control Reaction Engineering and Modelling group. He joined Johnson Matthey, Sonning Common, in 2000, and has had a variety of roles, including in gasoline and diesel catalyst research, and leading many collaborations with academia on various engineering and catalysis related projects.

By |2020-03-16T15:39:51+00:00March 16th, 2020|Weld Engineering Services|Comments Off on Insights into Automotive Particulate Filters using Magnetic Resonance Imaging

Nanosurfaces 2019

Johnson Matthey Technol. Rev., 2020, 64, (2), 148

Introduction

Nanosurfaces 2019 took place on 15th October 2019 at the Institute of Materials, Minerals & Mining (IOM3), London, UK. This event included commercialised techniques featuring plasma, surface coatings and graphene. Applications ranged from anticorrosion and waterproofing of electronics to battery materials and biomedical applications. Scalability, reproducibility, economics and sustainability were considered by many of the speakers.

Plasma Applications

The Keynote ‘Water Protection Technologies for Smartphones and Beyond’ was presented by Angeliki Elina Siokou, P2i Ltd, UK. P2i originated as a spinoff from Durham University, UK. The original product was an oil repellent coating for textiles created in 2004. They then developed a hydrophobic coating for electronics. Initial applications were in hearing aids then, around 2011, the company moved into water protection for smartphones. They are now coating half a billion smartphones per year. Customers include Samsung, Apple, Bose, Huawei and ASUS. Around 18% of smartphones are thought to be damaged by water ingress every year. P2i’s technology uses a radio frequency (rf) plasma attachment and polymerisation of surface molecules to provide a very thin hydrophobic conformal coating. This is suitable for mass production with good quality control. If a direct current (DC) plasma is used, the molecular structure is disrupted whereas pulsed DC retains coherent molecular structure and the best liquid repellence. Photoelectron spectroscopy analysis was used to analyse bond structure. The original product coating was splash proof but not completely waterproof, so a new coating was developed which is thicker (~100 nm) and has increased crosslinking. Thickness alone does not determine the level of water repellence (standards from IPX2 to IPX8), the process parameters affect porosity and film defects (1, 2). The coating was further improved by ensuring all product pieces are placed close to a plasma electrode. The coating was claimed to be durable and is tested with methods including the adhesive tape test, thermal, chemical and abrasion. It is economic for high-value items. The maximum temperature for the coating is 200ºC, determined by the maximum temperature for the items to be coated.

‘Enhancement in the Reactive Deposition Process Using Remote Plasma Spluttering’ by James Dutson, Plasma Quest Ltd, UK, presented the use of a remote plasma to generate ions whose path when accelerated is curved by magnets and focussed onto the sputter target. Key advantages of this technology are: high target usage, control over energy of sputter ions, low temperature at the substrate. They have had good success with aluminium doped zinc oxide as an indium tin oxide (ITO) replacement. Substrate heating is not required to control the phase or crystallinity of deposited material as this can be done via sputter ions energy control. Deposition rates are typically 1–6 μm h–1, and 1 μm h–1 for the ZnO:Al materials. Plasma Quest is working with CSEM, Neuchâtel, Switzerland, on photovoltaics and transistors. The process can be carried out at ambient temperature so is suitable for coating plastics or sensitive substrates. Reactive deposition is possible, i.e. to incorporate an oxide or nitride from a metal target. The main driver is replacement of ITO for transparent conductive coatings and there is a Horizon 2020 European Union (EU) project: INREP, Towards Indium Free TCOs. The company is currently scaling up to large surface areas ≥50 cm2. The process does not require a clean room.

‘Molecular Plasma Technology’ was presented by Britta Kleinsorge, Molecular Plasma Group, Luxembourg. This company offers atmospheric plasma (50 kHz) processes for surface modification. It is lower temperature compared to corona like arc plasma processes. Organic precursors could be functionalised and adhere to a surface at low temperature. Some interesting work on biological non-specific binding and antibody immobilisation was carried out with KU Leuven, Belgium, in which the process time was reduced from 24–72 h to 10 s. Superhydrophobic barriers have been developed for anticorrosion applications, waterproofing and adhesion. They presented a case study of a process they developed for a German automotive original equipment manufacturer that could be automated with reduced chemical usage.

‘A Radically New Antimicrobial Nanosurface Formed by Plasma Processing’ was presented by Alistair Kean, NikaWorks Ltd, Watlington, UK. There is a need for more environmentally friendly coatings, for example crisp packets which currently use 30–40 nm Al on plastic film. Tungsten carbide (WC) physical vapour deposition (PVD) coatings for surgical scissors need to be durable for decontamination cycles. The BeBionic prosthetic hand includes a titanium alloy coating on one surface. Three dimensional (3D) nanomaterials (metamaterials) are inspired by nature and nanosurfaces can have surface areas of the order of ~1000 m2 g–1. SOLAMON was a 7th Framework Programme for ‘large’ NP 20 nm, 30 nm or 40 nm scale. Gencoa Ltd is a PVD company in Liverpool, UK, which is developing antimicrobial coatings with NikaWorks.

‘Virtual Cathode Deposition’ was presented by Dmitry Yarmolich, Plasma App Ltd, UK. Plasma App is based at Harwell Oxford Science and Innovation Campus, UK. The speaker described a thin film battery project with the University of Cambridge, UK, which can deposit carbon as a combination of graphite and graphene immediately followed by the lithium cobalt oxide (LCO) cathode. It is a platform technology that can deposit a thick or thin film of almost anything from nanometre scale to 50 μm. Superior adhesion and solvent-free scalability were claimed. The speaker claimed that the technique can use two source virtual cathode deposition (VCD) to design materials in minutes, not hours.

Nanoparticles and Graphene

‘Manufacture and Applications of Nanoclusters’ was presented by Richard Palmer, Swansea University, UK. He set up Grove Nanomaterials to commercialise matrix assembly cluster source (MACS) technology (3) for the solvent-free synthesis of nanoparticles (NP) for catalysis, sensors and electrodes. Atomic nanoclusters can be created to mimic enzymes with formation of grooves and pockets. Size selected cluster beam deposition (CBD) of gold or molybdenum sulfide, which can be nickel or cobalt doped to replace platinum in water splitting, was carried out. The present plan is to scale up the research system to be able to produce 1 g h–1 and the technology has the potential to be scaled to mass production (tonne scale).

‘Silica Nanoparticles for Super-Hydrophobic Ice-Repellant Coatings’ was presented by Simon Haas, Promethean Particles Ltd, UK. Promethean Particles Ltd is a spinout from the University of Nottingham, UK. Using a continuous hydrothermal process, copper, silver, ZnO, barium titanate and silica NP can be produced as dispersion in liquid with no dry powders for safer handling. The liquids can be used in further processes and agglomeration is prevented. The dispersions are highly stable and the process is scalable. A production plant has been built which can produce 1000 tonnes per year based on dry weight equivalent. Conductive inks and printed electronics are currently being developed with partners. The speaker presented an Innovate UK funded programme, ICEMART, looking at SiO2 NP to prevent icing of plane wings. The Welding Institute, Great Abington, UK, was part of this project and provided a process to functionalise the SiO2 particles to render them hydrophobic. This technology is currently being commercialised via Sharc® Matter, an Opus Materials Technology Company, UK.

‘Graphene Enhanced Products’ was presented by Thanuja Galhena, Versarien plc, UK. Versarien acquired Cambridge Graphene Ltd, UK, in 2017 and supplies proprietary grades of two-dimensional materials. There are two main products: NaneneTM (graphene) and HexoteneTM (boron nitride). Both can be manufactured at a 3 tonnes per annum scale. The company has certification from the Graphene Council, New Bern, USA and the National Physical Laboratory (NPL), Teddington, UK, European Registration, Evaluation, Authorisation and Restriction of Chemicals (REACH) registration in the EU and is applying for similar certifications in China. Potential applications are composites, fire retardants, barrier coatings, energy storage, filtration and heat dissipation. A fairly unusual product application was a railway arch supporting sensors, claiming that lighter weight reduces transport, installation time and personnel costs. Graphene based inks known as GRAPHINKS® have been developed (4). Another example is a graphene coated diaphragm for some earphones – apparently bass is clearer. FLEXIBAT is an Innovate UK funded project. The company is working on flexible batteries with Zinergy UK Ltd, Cambridge, UK, making graphene coated electrodes for corrosion resistance. There are also supercapacitors using graphene plus metal oxide NP. The speaker noted that the cost of graphene has reduced considerably, making it cost effective for many applications.

‘Enhanced Surface-Analysis Capability’ was presented by Mike Petty, Loughborough Surface Analysis Ltd, UK. Petty focussed on the new secondary ion mass spectrometry (SIMS) system (Cameca, France) which can do quite sophisticated analysis mainly on semiconductor materials. Data cubes are formed so one can go back and analyse specific coordinates within a scanned cube (Figure 1). This technique is sensitive to isotope identification.

Fig. 1.

SIMS images of 3D silicon/tantalum test object taken using the CAMECA SIMS instruments. Image courtesy of Mike Petty, Loughborough Surface Analysis Ltd, UK

SIMS images of 3D silicon/tantalum test object taken using the CAMECA SIMS instruments. Image courtesy of Mike Petty, Loughborough Surface Analysis Ltd, UK

‘Photocatalytic Antimicrobial Surfaces’ by Jeremy Ramsden, University of Buckingham, UK. Ramsden discussed the antimicrobial application of photocatalytic powders coated via spray from titanium dioxide powders in a water solution. He described the problem of hospital acquired infections, especially resistant infections such as methicillin-resistant Staphylococcus aureus (MRSA), and presented studies of deaths and disability-adjusted life year (DALY) statistics in Europe, USA and elsewhere dating from 1930s to the present. There are many vectors of microbial transfer in and around hospital environments, these are surprisingly little studied and there is much debate about the relative importance of air, walls, floors, shoes or wheels and ceilings in transmission of pathogens. In particular there are few or no controlled studies on hospital cleaning efficacy. Hand hygiene compliance seems to have peaked at ~40%. He studies non-sacrificial catalytic coatings, i.e. TiO2 which can kill microbes even under ambient lighting with continuous effectiveness. Estimated 10 colony forming units (CFU) m–2 s–1 of bacteria arriving and 6000 s–1 oxidising equivalents forming. He carried out a trial in which a coating was applied after a hospital deep clean on the high-touch surfaces, for example bed rails, table and tubing. The surfaces were tested using Agar pads (selective for MRSA and general). The coating used a commercial sol of TiO2 NP in a sol-gel process to form a hard coating within 20–30 mins. The coating is invisible (1 μm) and can be applied to mirrors and glass. Sampling was carried out for three months, the coating was durable in this timeframe. Evidence suggested the coating resulted in lower microbial growth (Figure 2) (5). Microbial resistance is unlikely as the microbial cell has no mechanism for defence against peroxide. The half-life of bacteria was ~6 h.

Fig. 2.

Stained Staphyloccocus epidermidis within a biofilm (confocal laser scanning microscopy). Green fluorescence: live bacteria; red fluorescence: membrane-compromised bacteria. (a) 1.5 h exposure to UVA only; (b) 1.5 h exposure to photocatalytic treatment; (c) 3 h exposure to UVA only; (d) 3 h exposure to photocatalytic treatment; (e) 3 h exposure to TiO2 in the dark; (f) 3 h no treatment control (no TiO2, no UVA exposure). Reprinted from (5) with permission from Elsevier

Stained Staphyloccocus epidermidis within a biofilm (confocal laser scanning microscopy). Green fluorescence: live bacteria; red fluorescence: membrane-compromised bacteria. (a) 1.5 h exposure to UVA only; (b) 1.5 h exposure to photocatalytic treatment; (c) 3 h exposure to UVA only; (d) 3 h exposure to photocatalytic treatment; (e) 3 h exposure to TiO2 in the dark; (f) 3 h no treatment control (no TiO2, no UVA exposure). Reprinted from (5) with permission from Elsevier

Conclusions

Martin Kemp, Chairman of the IOM3 Nanomaterials Conference organising committee, remarked how many of the presentations featured plasma techniques and was excited by the many examples of commercialised techniques. There will be another meeting in early 2020.

The Authors


Alistair Kean is a New Business Development Consultant at Johnson Matthey, Sonning Common, UK. In 2013 he set up NikaWorks Ltd with the aim of commercialising nanotechnology. He is a visiting professor at Manchester Metropolitan University, UK and in 2019 he accepted an appointment at University of the Highlands & Islands (UHI), UK, as professor of medical nanotechnology.


Sara Coles is Editor of Johnson Matthey Technology Review. She obtained her Master’s degree in Chemistry from Lancaster University, UK, and has previous experience in pharmaceuticals research. She is interested in all areas of science and technology relevant to Johnson Matthey’s vision of enabling cleaner air, improved health and the efficient use of natural resources.

By |2020-03-16T14:51:55+00:00March 16th, 2020|Weld Engineering Services|Comments Off on Nanosurfaces 2019

Driving home clean transportation with two new International Standards

Responsible for nearly a quarter of global carbon emissions, transport is the achilles heel of nearly every city. Yet, with energy demand on the rise, so much needs to be done to tackle the 2030 Agenda targets. New standards have just been published to help achieve smart and efficient transportation in new and existing urban areas.

Car-sharing apps, electric cars and buses and new bike lanes have all helped to stem the rise in energy demand for transport, but they are not enough. As populations grow, demand goes with it, meaning that road vehicles continue to account for the vast majority of transport-related CO2 emissions.

Smart transportation that uses less energy and has a reduced impact on the environment means everyone can breathe easier: citizens, the environment and the world. What’s more, it makes for a more attractive city. Two new International Standards have just been published to help.

ISO 37161, Smart community infrastructures – Guidance on smart transportation for energy saving in transportation services, provides guidance to transport organizations, local and national governments and any other industry player on how to reduce energy used in transportation for passengers, delivery, freight and postal services. It features energy-saving options that can be adopted, along with the maintenance and monitoring of such measures.

The standard is complemented by ISO 37162, Smart community infrastructures – Smart transportation for newly developing areas, that helps town planners of new developments and areas implement transport services that not only meet the needs of the population but address sustainability issues.

Electric bus on charging with using pure energy from windmills.

Dr Yoshiaki Ichikawa, Chair of the ISO subcommittee that developed the standards, said reducing energy consumption and implementing changes that favour the environment must take a holistic approach.

“Transportation networks can be very complex, involving not only the vehicles used but the entire operational system,” he said.

“These standards provide guidance on all aspects of transportation, including the criteria and parameters to be considered in the selection of energy-saving options, for long-term results.”

ISO 37161 and ISO 37162 are two of many ISO standards aimed at smart transportation, with others including ISO 37157 (compact cities), ISO 37158 (battery-powered buses for passenger services) and ISO 37159 (rapid transit in and between large city zones).

The standards were developed by ISO technical committee ISO/TC 268, Sustainable cities and communities, subcommittee SC 1, Smart community infrastructures, the secretariat of which is held by JISC, ISO’s member for Japan. They can be purchased from your national ISO member or through the ISO Store.

By |2020-03-12T07:46:36+00:00March 12th, 2020|Weld Engineering Services|Comments Off on Driving home clean transportation with two new International Standards

What’s in a Brand? Quite a bit, actually

An independent and objective system for evaluating brand performance is essential. Two ISO standards help to get the job done.

It is widely agreed that a brand can be a key strategic asset for a company, but exactly how do you value it? Brands are intangibles, unlike tangibles such as factories, offices or store fronts, and this presents a challenge for conventional accounting standards.

With the recent release of ISO 20671, Brand evaluation – Principles and fundamentals, Dr Bobby J. Calder, Chair of ISO technical committee ISO/TC 289 on brand evaluation, explains how it covers the technical requirements and methods involved in measuring the strength of a brand. This standard extends the scope of ISO 10668, Brand valuation – Requirements for monetary brand valuation, that focuses primarily on the financial valuation of brands. Ultimately, ISO 20671 aims to help resolve the differences of opinion around conventional accounting treatments of a brand.

Here is where the potential struggle arises, in Dr Calder’s opinion. He believes the crux of the matter is that “finance and marketing don’t speak the same language. Marketing focuses on justifying brand expenditures and finance focuses on controlling them”. Both need to work together to treat brands, not as an expense, but as a key financial asset.

So how can these two – often at-odds – functions located in the same company move toward a common understanding of brands? Dr Bobby J. Calder shares his insights.

ISOfocus: What does “brand” mean and why is brand evaluation important?
Bobby J. Calder.

Dr Bobby J. Calder, Chair of ISO technical committee ISO/TC 289 on brand evaluation.

Dr Bobby J. Calder: Brand entails distinctive images that include names, terms, logos and/or signs to help a company distinguish itself from others in the market. The brand can be commercial or not-for-profit, but the function of all brands is to create a recognizable entity in the market that, in the mind of consumers, adds value to the product. Thus, for consumers, brands are the perceptions, benefits and experiences that they associate with a good or a service. From the point of view of the company, the main purpose of a brand is to increase cash flow through price premiums, lower costs, increased volume, or greater repeat purchase loyalty.

As an intangible asset, brands, unlike machinery, buildings or products, have no physical substance. Whereas brands can have tangible value as trademarks or customer lists, the primary value of the brand is intangible. The expected economic benefit to the company flows from association in the consumer’s mind. The challenge of how to effectively recognize and value a brand as an intangible asset is what led to the creation of ISO 20671.

In general, intangible assets such as brands are becoming more and more economically important. Though not formally recognized in GDP reports, economists find that the value of intangibles is higher than tangibles in most developed economies. And intangibles increasingly distinguish successful companies from less successful ones. This macro situation underscores the need to transcend accounting and finance traditions and treat brands as part of a company’s value creation process.

So there’s a long-standing need to bridge the divide between marketing and finance. Activities such as social media, mobile apps, sponsorships and the like certainly help build brands, but we need to understand that brands exist in the mind of the consumer. The greater the “strength” of the brand in affecting consumer buying decisions, the greater the value of the brand to the company as a financial asset. Companies need to periodically evaluate the strength of the brand and its contribution to economic returns. Based on this evaluation, better decisions can be made about investing in brand-building activities.

Time Square, New York.
Why did ISO create a technical committee to address brand evaluation?

Currently, there are many marketing metrics for assessing a brand (awareness, willingness to recommend, etc.), but no shared or accepted framework for linking the value of a brand to consumers with the value it creates for the company. ISO 10668, Brand valuation – Requirements for monetary brand valuation, tackled this from a financial valuation standpoint, but it was acknowledged that this was limited in scope. It emphasized three valuation methods. The market approach values the brand against the price of a comparable brand. The income approach uses the present value of future cash flows that a company would receive when using the brand. The hybrid approach involves royalty relief through basing the royalties a company would be required to pay if it had to license the brand from another entity.

Our technical committee developed ISO 20671, Brand evaluation – Principles and fundamentals, to provide a broader framework for evaluating brands. ISO 20671 is designed to be a useful resource for companies to rationalize their treatment of brands and to be able to more accurately report on brand value to internal decision makers and external investors. Moreover, there continues to be the need for companies to reduce economic risk, so being able to identify the value of a brand is another useful tool in this regard.

ISO 20671 is the world’s first International Standard for brand evaluation. Can you explain the main principles of this standard?
Two women working in a modern office.

Well, if ISO can do it for the quality of products such as USB cables, then why not brands? Brand evaluation standards are a logical next step for ISO. These days, marketers have many different ways of analysing and communicating about brands, and companies vary greatly in their methods. Sometimes, we marketers can be our own worst enemies by inventing our own jargon, which makes it very confusing for others within the same company, and even more so for outsiders. ISO 20671 provides clear definitions for specific terms to help eliminate this confusion and an overall framework that everyone can refer to.

The main principle is for organizations to build on this framework to begin to link the marketer’s view of brands to the internal investment and governance process, and to explore ways of reporting this information to external investors. ISO 20671 is therefore a useful starting point for companies and organizations wanting to increase their brand value. It provides a universal view of non-financial and financial measures with the intention for more specific standards to be developed in the future in conjunction with companies undertaking this process.

What future aspirations do you have for these standards?

As a concept, ISO 20671 can be applied to all companies wishing to evaluate the value of their brand more effectively. It does not need to be specific to all products, industries and services just now; it is much more important to agree on the general definitions and framework first and work on more specific guidelines later.

For me, longer-term, I believe this work could tie into many issues facing the future of business. For instance, there is an ongoing debate as to whether companies should focus exclusively on maximizing shareholder and owner value or on creating value for all stakeholders. There are good arguments on both sides. I think it is possible to approach this issue from the perspective that brands could potentially allow a company to do both. Companies can build brands that incorporate sustainability and other constructive purposes into the very idea of the brand in the consumer’s mind. Such brands would be stronger and yield higher returns to shareholders as well as returns in the form of societal capital to stakeholders.

Top view of young designers discussing a colour palette.
By |2020-03-10T09:03:47+00:00March 10th, 2020|Weld Engineering Services|Comments Off on What’s in a Brand? Quite a bit, actually
Go to Top