Svoboda | Graniru | BBC Russia | Golosameriki | Facebook
Perspective The following article is Open access

Future prospects for the design of 'state-of-the-art' solid oxide fuel cells

, and

Published 3 July 2020 © 2020 The Author(s). Published by IOP Publishing Ltd
, , Citation Toshiyuki Mori et al 2020 J. Phys. Energy 2 031001 DOI 10.1088/2515-7655/ab8f05

2515-7655/2/3/031001

Export citation and abstract BibTeX RIS

Original content from this work may be used under the terms of the Creative Commons Attribution 4.0 license. Any further distribution of this work must maintain attribution to the author(s) and the title of the work, journal citation and DOI.

Solid oxide fuel cells (SOFCs) are the clean and efficient power sources for generating electricity from a variety of fuels (i.e. hydrogen, natural gas, and biogas) [13]. Also, SOFCs have no corrosive components and do not require precious-metal electrocatalysts due to the high operation temperatures (800 °C–1000 °C) [3]. In the recent research, the major target focuses on design of functional interfaces in SOFCs device and system for lowering the operation temperature of SOFCs to intermediate temperatures (IT-SOFCs, 400 °C–700 °C) in order to increase the operation durability, improve the thermal compatibility, and thermal cycle capability and reduce the fabrication and materials costs by using metallic interconnectors in the SOFC stack cell. In this endeavor, it is well-known that two major strategies have been adopted to increase the performance of single cells in the intermediate temperature region. One strategic challenge is the fabrication of thin electrolyte (thickness: less than 10 μm) with high density to minimize the ohmic loss of the electrolyte. And another one is the reduction of excess overpotentials observed for the O2 reduction reaction at the cathode in the single cells [4].

However, to develop the 'state-of-the-art IT-SOFCs' with quality of exceeding the limit of conventional IT-SOFCs, an advanced strategy for the design of key interfaces in SOFC devices and systems is required.

SOFC device and system contain a number of interface regimes which play key roles in the performance of SOFC device and system. The gas-solid interfaces on both cathode and anode assemblies control the charge transfer phenomena in the fuel cell reaction. Also, the three-phase boundary (TPB) at the electrode-electrolyte interface is a key active site area where both the electron exchange and formation of conduction ions take place. As a consequence of this, the key microstructure at the interface in SOFC devices has been discussed previously in the published literature [5].

A schematic representation of key interface regimes is illustrated in figure 1 [5]. The key interfaces of SOFC devices include (a) the second phase inclusion, (b) possible formation of ordered micro-domains, (c) liquid phase accumulation induced by impurity, (d) segregation effect, (f) precipitation in the matrix, or (g) second phase at grain boundary, as demonstrated in figure 1. In each case, the control of the interface structures at cathode/electrolyte, in the electrolyte and in the cermet anode layer can influence the performance of the IT-SOFC device. Therefore, we would like to summarize key points for development of IT-SOFCs in following sections.

Figure 1. Schematic illustration of key interface regimes in SOFC device.

Standard image High-resolution image

1. Cathode and cathode/electrolyte interface

For the development of IT-SOFC devices, the first key endeavor is the design of the interface between cathode and electrolyte. Among the various kinds of previously proposed electrolyte, lanthanum strontium manganite (LSM) perovskite-based material is the most important cathode material in a SOFC device. It is because the LSM has an excellent track record of microstructural stability and long-term performance stability. However, the polarization performance and electrochemical activity of LSM based cathodes are a complex function of composition, the microstructure of cathode/electrolyte interface and activation effect. The key points for usage of LSM cathode in SOFC devices are summarized in the previously published review paper [6]. Usually, the LSM-based cathode is used for high-temperature SOFC devices, the electrocatalytic activity of LSM-based cathode can be substantially enhanced by the design of microstructure in nano-scale (e.g. by infiltration and impregnation of catalytic active and 'mixed ionic and electronic conducting or MIEC' nanoparticles) [7]. The performance of nanostructured LSM-based cathode is comparable to that of typical MIEC cathodes such as cobaltite-based perovskite cathodes (e.g. Ba0.5Sr0.5Co0.8Fe0.2O3−δ (BSCF), La0.6Sr0.4Co0.2Fe0.8O3−δ (LSCF) and PrBa0.5Sr0.5Co1.5Fe0.5O3−δ (PBSCF)). Among cobaltite-based perovskite cathodes, LSCF-based cathode is the most popular and representative MIEC cathode material for the development of IT-SOFC devices [8]. However, it was reported that the performance of IT-SOFC with LSCF cathode deteriorated by the surface segregation of Sr species and formation of Sr-rich reaction layer at the interface of the cathode/electrolyte of the single cell [9]. The previously published data suggests that the functional interfacial layer can be formed using MIEC highly ionic conducting nano-particles such as Er0.4Bi1.6O3 (ESB) decorated LSCF-based MIEC cathodes on barier-layer-free YSZ electrolyte under SOFC operation condition [10]. The polarization induced ESB-rich interface is highly active for the O2 reduction reaction and inhibits the interfacial reaction between LSCF and YSZ by supressing the Sr surface segregation.

In addtion, a kinetics investigation of the oxygen reduction reaction on LSCF-Gd2O3 doped CeO2 (GDC) composie cathode suggested the origin of distortions in measured impedance spetraum arc observed for LSCF-GDC cathode [11]. Also, the model analysis of the reactive pathway at TPB on the basis of kinetics and thermodynamics is useful to understand the role of functional interface of LSCF-Gd2O3 doped CeO2 composite cathode [12]. The conventional modelling on the basis of ideal interface structure of MIEC cathode would be useful to optimize the complicated composition of MIEC cathodes. However, it would be difficult to develop design concept of key active interface of cathode for societal implementation of 'state-of-the-art' IT-SOFC device and system. To maximize the effect of MIEC particle addition and minimize the affection of negative factors such as Sr segregation effect for design of high quality cathode in IT-SOFC devices, the actual interface structure on the cathode [13] should be taken into account for development of new fabrication process such as decoration. Besides new fabrication process and modelling on the basis of model interface strucutre, the modelling based on microanalysis data should be used for the development of the active interface structure and useful production concept. As a consequnce of this, the authors believe that the fusion of fabrication, characterization and modelling based on microanalysis data is the only way to develop the 'state-of-the-art' IT-SOFC device and system.

2. Electrolyte interface and electrolyte microstrucutre

The second key endeavor for the development of IT-SOFC is the control and stability of the microstructure in the solid electrolyte. In the previously published work, the formation of microdomain with the ordered structure of oxygen vacancies was clearly observed for CaO stabilized ZrO2 and Sc2O3 stabilized ZrO2 electrolytes using analytical transmission electron microscopy (TEM) [5]. Since the defect interface change in the dense ZrO2 based solid electrolytes such as CaO stabilized ZrO2, Sc2O3 stabilized ZrO2 and Y2O3 stabilized ZrO2 is subtle in the microanalysis, the influence of microdomain formation on the electrolytic properties has been disregarded in the R & D field of high-temperature SOFC devices. However, M2O3 doped CeO2 (M: Gd2O3, Y2O3, Sm2O3, La2O3 and so on) electrolytes which mainly consist of fluorite structure, as well as Y2O3 stabilized ZrO2 based solid electrolyte expose a unique feature of defect structure in TEM microanalysis. Especially, it was considered that Gd2O3 doped CeO2 can be used as a solid electrolyte for IT-SOFC whose operation temperature is below 500 °C in 2000 [14]. As a consequence of this, the microdomain structure in M2O3 doped CeO2 (M: Gd2O3, Y2O3, Sm2O3, La2O3 and so on) dense sintered bodies were systematically observed using analytical TEM in the previously published works [1524].

The key points of microdomain formation processes are summarized in figure 2. In general, easily sinterable powders are prepared by using a soft chemical processes such as sol-gel method, co-precipitation method and others. In this soft chemical process, the composition heterogeneity in nanoscale remains in the prepared nano-sized powders. During the sintering process to fabricate dense sintered bodies, nano-level dopant segregation, which is formed in the grain boundaries of electrolyte sintered bodies, seriously affects the size and distribution of the microdomains in the electrolytes [1519]. While Gd2O3 doped CeO2 electrolyte includes smaller microdomains and lower observed activation energy as compared with all of other CeO2 based solid electrolytes, the observed oxide ion diffusion becomes low level due to the formation of microdomains in the bulk. As a consequence of this, the microstructure design should be performed for stabilized ZrO2 electrolyte on the basis of microanalysis results obtained on CeO2 based electrolytes which mainly consist of fluorite structure.

Figure 2. Schematic diagram of the formation process of micro-domain with the ordered structure of oxygen vacancy in doped CeO2 solid electrolytes.

Standard image High-resolution image

3. Anode and anode/electrolyte interface

A third key endeavor for the development of IT-SOFCs is microstructure control in the anode layer. In previously published works, the 'mixed ionic and electronic conducting (MIEC)' particles such as Y-doped SrTiO3 was used in the cermet anode as well as LSM-based cathode [25]. In this case, both of nano-sized CeO2 and Ru (0.4 vol%) as co-catalyst were added into the anode with MIEC particles using the inflitration method, and chemical stability against H2S in H2 fuel and performance of SOFC were improved. Also, it was reported that addition of La0.3Sr0.7TiO3 together with infiltration of CeO2 and Pd (0.5 wt%) [26], ceria doped with Gd2O3 or Sm2O3 infiltration, and CeOx-Cu infiltration methods showed improvement of anode performance and chemical stability of anode to small amout of impurity H2S in the fuel [27]. The previously reported infiltration technique using MIEC nano-particle and catalysts such as CeO2 and precious metal (less than 1 wt%) or Cu additives were useful for a lowering of excess overpotential on the anode and enhancemnt of the resistance to H2S poisioning and carbon deposition (i.e. 'coking'). However, it is not sufficient in order to minimize the precious metal or metal content and maximize the effect of improvement of anode performance. In addition, a distributed charge-transfer model analysis of Cu- and Ce- based IT-SOFCs well simulated the polarization and impedance spectra observed for Cu–Pd-mixture of CeO2 and ZrO2/Sm doped CeO2/LSCF-Gd doped CeO2 [28]. However, it would not be an easy task to develop design concept of active interface using experimental data simulation together with mass and charge conservation equations.

To design the active interface in the 'state-of-the-art' IT-SOFCs with the quality of exceeding the limit of conventional IT-SOFCs, the actual interfacial defect structure on the anode should be taken into account. Also, modelling on the basis of microanalysis data should be performed in order to propose the design paradigm for fabrication of 'state-of-the-art' IT-SOFC devices and systems. In the anode layer, the TPB areas are formed at the interface between Ni and electrolyte particle such as CeO2 doped with Sm2O3, Gd2O3 or Y2O3 stabilized ZrO2. STEM-EDS analysis examined the TPB area of an anode layer which consists of M2O3 doped CeO2 (M: Y2O3 or Gd2O3) and Ni as shown in figure 3. In the narrow interface (i.e. width: approximately 200 nm) between doped CeO2 and Ni, C-type rare earth related superstructure formation was clearly observed by STEM-EDS analysis [24, 29, 30]. To conclude the defect structure of C-type rare earth related superstructure well, the lattice statics modelling (Source code: GULP, empirical potential: Buckingham potential) was performed on the basis of analytical TEM microanalysis results (i.e. EELS analysis data and EDS analysis data) using workstation. This modelling indicates the small defect cluster which consists of dopant cation and oxygen vacancy is easily formed in the bulk of CeO2 based electrolyte [31] and at the interface of between Ni and CeO2 based electrolyte particle in the anode layer. This small defect cluster is easily combined with each other to form larger defect cluster in the narrow interface between Ni and CeO2 based particles. In addition, the oxygen vacancy position of large defect cluster agrees with that of C-type rare earth structure [24, 2931]. It means that the lattice statics modelling on the basis of microanalysis data is quite useful to conclude the defect interface structural feature of SOFC device well.

Figure 3. Schematic diagram of the formation area of C-type rare earth related superstructure in the anode of IT-SOFC single cell device.

Standard image High-resolution image

In general, the simulation of fuel cell performance in ideal modelling field (i.e. the first-principles simulation) is performed for the design of high-quality IT-SOFC device. Recently, first-principles simulation in combination with already available data predicted an innovative breakthrough on the design of high-quality anode in IT-SOFCs [32]. This first-principles simulation work suggests that the surface oxide ion diffusion on YSZ in the cermet anode is the rate-determining step for the formation of water molecules in the fuel cell anode reaction below approximately 800 °C, although the dissociation rate constant of H2 on Ni is known to be high. In addition, this first principles simulation clearly indicates that the formation reaction of water molecules on the anode will be clearly enhanced, if we design the new active site for surface oxide ion diffusion on cubic ZrO2. While the guide from the theoretical analysis is important, it is hard to create new active sites on the cubic ZrO2 surface using only data of first-principles simulations. This is because the real anode surface reaction is much more complicated as compared with the ideal anode reaction and additional factors which cannot be covered by theoretical analysis should be considered to design the new active sites on the real anode materials in IT-SOFCs.

To overcome this challenge, we tried to modify the defect structure at TPB on the Ni-cubic ZrO2 cermet anode by a combination of the lattice statics simulation and microanalysis of defect interface structure as well as combination analysis of Ni-cubic CeO2 cermet anode in the IT-SOFCs. As the authors mentioned in part 2 (electrolyte interface and electrolyte microstructure) and in the early part of the present part 3 (Anode and anode/electrolyte interface), the ordered structure of oxygen vacancy is formed as non-equiribrium defect at the interface of SOFC devices. This non-equiribrium defects is formed in the slow cooling process of fabrication of SOFC devices. In both Y2O3 stabilized ZrO2 and CeO2 doped with Sm2O3, Gd2O3 or Y2O3 systems, similar micro-domains with ordered structure of oxygen vacancy which can be assigned by C-type rare earth related structure were formed during fabriction process of SOFC devices. Also, its formation width at TPB area on the anode is quite narrow as we demonstrated in figure 3 such as 200 nm or less. In this narrow region, the amount of inactive lattice defect with C-type rare earth related defect structure would be in extremely low level such as 10 ppm (10 mg kg−1) level. To modify this trace amount of inactive lattice defect for design of active site, the conventional infiltration process is not suitable. If the doping amount level is in too high level, different types clustering will occur in the defect interface region on the anode. Then, it is hard to expect the maximization of anode performance by using conventional infiltration method.

To activate the inactive defects at TPB on Ni-cubic ZrO2 cermet, we used trace amount of PtOx whose content is approximately 10 to 100 ppm using sputtering technique [33]. The PtOx was deposited on the prepared NiOx-8YSZ, conventional cermet anode, by means of magnetron sputtering as illustrated in figure 4.

Figure 4. Schematic diagram of PtOx sputtering process for activation of useless defects in superstructure at TPB on Ni-Y2O3 stabilized cubic ZrO2 (YSZ) cermet anode in IT-SOFC.

Standard image High-resolution image

The inductively coupled plasma mass spectrometry (ICP-MS) quantitatively analysed the Pt content in whole anode layer after reduction of anode layer at 800 °C in an H2 gas flow (i.e. the conventional conditioning process of the anode layer). After the conditioning process of the anode with trace amount of PtOx (Pt content: approx. 10 to 100 ppm), we confirmed that the surface deposited Pt species were diffused into the anode layer by using analytical TEM, as we expected in figure 4.

The current density observed for the IT-SOFC single cell in which the PtOx was deposited on the anode is approximately seven times higher than that of the IT-SOFC single cell without PtOx at 700 °C, refer to figure 5 [33]. Also, the combination analysis of analytical TEM and lattice statics surface simulations (source code: GULP, empirical potential: Buckingham potential) confirmed that the defect interfacial structure located on partially oxidized Ni surface which is very close to TPB mainly consists of Pt cation (i.e. Pt2+ (Pti¨)) and lattice vacancies such as oxygen vacancy (VO¨) and interstitial oxygen (Oi''), refer to table 1 [33].

Figure 5. Current density vs. cell potential (IR-free) observed for trace amount of PtOx (Pt content: 10 ppm) sputtered anode and conventional anode. Anode: YSZ-Ni cermet anode (YSZ/Ni = 1/4), anode gas: humidified H2 gas (80 ml min−1), cathode: commercially available LSM (LSM20-I, Fuel Cell Materials. Com, USA), cathode gas: O2 gas (80 ml min−1), operation temperature: 700 °C.

Standard image High-resolution image

Table 1. Assignment of extra reflections in selected area electron diffraction pattern recorded from the surface of Ni in the cermet anode with trace amount of Pt species and calculated binding energy on the basis of microanalysis data.

Lattice spacing of extra spots / nm Assigned materials and lattice
0.152 PtO (200)
0.230 Pt (200)
0.164 unknown
0.199 Pt (311), NiO (222)
0.200 Pt (200), Pt (111)
0.204 Pt(200), Pt(111), Ni(111), NiO(200)

On partially oxidized Ni surface at near the three-phase boundary;

ΔEb of VNi" - Pti•• - Vo•• - Oi" : 1.1 / eV

The lattice statics surface simulation on the partially oxidized Ni surface along the guide of TEM microanalysis result shows high binding energy (ΔEb) such as 1.1 eV. Note that high positive ΔEb value indicates the preferable cluster structure in the defect formation reaction. Also, the crystallographic defect cluster model indicates that buffer area for minimization of the lattice mismatch between cubic ZrO2 surface and partially oxidized Ni surface is formed at TPB area. And this buffer area would promote surface oxide ion diffusion (i.e. rate-determining step of anode reaction below 800 °C) from YSZ surface to Ni surface at TPB area. Then, it is concluded that the excess overpotential on the anode becomes low by modification of inactive non-equilibrium lattice defect at TPB area on the Ni surface. In contrast, the defect interfacial structure on cubic ZrO2 surface has not been identified due to the difficulty of observation of the defect structure with coherent interface for TEM sample. To overcome this difficulty for characterization of coherent interface defect structure, the first-principles surface simulation (source code: ABINIT, pseudopotential: PAW method) was performed using coherent interface model (i.e. Schottky defect model), refer to figure 6. This first-principles surface simulation provides us useful hint for design of active site on YSZ side. Our first-principles surface simulation suggests that di-valence (i.e. Pt2+, Fe2+) and tri-valence (i.e. Rh3+, Fe3+) cations create the useful defect cluster structure as active site at TPB area. In addition, our simulation suggests that defect cluster formation energy (ΔEf) which is calculated by defect cluster with Fe3+ cation shows high negative value such as −5.8 eV. Note that high negative ΔEf value indicates the preferable cluster structure in the defect formation reaction.

Figure 6. Defect formation energies (ΔEf) calculated by coherent interface model and cell potential as a function of current density observed for small amount of FeOx sputtered anode and conventional anode. Operation condition is same as figure 5.

Standard image High-resolution image

In addition, XPS analysis clearly detected Fe3+ species from the interface of Ni and cubic ZrO2 after conventional conditioning process (i.e. reduction of anode at 800 °C in H2 gas flow) (unpublished work).

This combination work of surface microanalysis and predictive first-principles surface simulation indicates that active site is created by Fe3+ species doping on cubic ZrO2 at TPB. Namely, the anode performance can be improved using Fe3+ species instead of Pt2+ species. Along the guide of this combination work, our multidisciplinary team successfully observed the improved anode performance using small amount of FeOx sputtered sample (i.e. Fe content: approximately 1000 ppm) as well as PtOx sputtered sample, refer to figure 6) [34].

4. Summary

The data which are described in sections 1 to 3 clearly suggest the importance of multidisciplinary research work as shown in figure 7. To see new breakthroughs for the design of 'state-of-the-art' next generation IT-SOFC devices, we must carefully characterize key defect interface structure in our devices. Also, we must perform the lattice statics surface and bulk simulation to conclude the key interface defect structure at TPB area. Then we must modify the interfaces on the basis of suggestion of combination analysis of microanalysis and surface modelling. If our model is reasonable, we will see improved properties from our IT-SOFC devices. After we confirm such the improvement of IT-SOFC device performance, more careful microanalysis such as operando-microanalysis and high precision modelling (i.e. first-principles simulation) for the design of high-quality TPB in the interfaces of IT-SOFCs are required. Through such an interactive design concept by fusing fabrication, micro-characterization and modelling knowledge, we have a much greater opportunity to purpose (target) design new electrode materials for GREEN innovations.

Figure 7. Schematic diagram of key concept on the basis of multidisciplinary research in which interactive combination of fabrication, microanalysis and modelling lead to the ultimately optimised fuel cell stack for the design of 'state-of-the-art' IT-SOFC systems.

Standard image High-resolution image

In general, the GREEN innovation is defined as technical innovations that consist of new or modified processes, practices, systems and products which benefit the environment and so contribute to environmental sustainability [35]. To establish the sustainable society through the formation of hydrogen society, R & D challenges for high quality polymer electrolyte membrane fuel cell (PEMFC) device and system or high power all solid state Li ion battery is very popular in the world. Certainly, hydrogen which is created by primary energy can be directly used for PEMFC system. In this case, however, there may be little chances to use the stored greenhouse gas (i.e. CO2 gas) which is recovered from exhausted gas well.

As illustrated in figure 8(a), H2 and CO2 gas can be converted to mixture of H2O and CH3OH (H2O/CH3OH molar ratio = 1) on the conventional catalyst. Then, this mixture gas can be used as fuel for SOFC device and system for generation of electricity in high efficiency. In this case, PEMFC can be used as backup generator to cover the peak electricity demand in summer season and in day time of industrial work. In this case, the deposition of carbon (i.e. 'coking') would be one of problems for design of 'state-of-the-art' IT-SOFC device with good balance of high performance and good durability. It is well-know that the 'coking' is depressed along the guide of previously reported the relationship between steam/carbon molar ratio in the fuel and operation temperature of SOFCs, refer to figure 8(b) [36].

Figure 8. Schematic diagram of IT-SOFC for GREEN innovation (a) and steam to carbon molar ratio vs. operation temperature of fuel cells (b).

Standard image High-resolution image

When we use the mixture of H2O and CH3OH (H2O/CH3OH molar ratio = 1) as fuel for IT-SOFC device, this fuel composition would be in boundary region between risky region of 'coking' (below the line in figure 8(b)) and the safety region (over the line in figure 8(b)). This point can be shifted to the safety region side in IT-SOFC operation by improvement of kinetics at TPB on the anode. Then, we will have great opportunity to see GREEN innovation by development of high-quality system with 'state-of-the-art' IT-SOFC in the near future.

Of course, we must make good balance between resources, production cost, high quality performance and long operation durability for economical, ecological and environmental implementation of the 'state-of-the-art' next generation IT-SOFC devices during such an endavour. It is well-known that e.g. the performance of IT-SOFC devices were degraded by many factors such as carbon deposition on the anode [37], chromium deposition/poisoning on the cathode [38] and micro-domain growth in solid electrolyte [5, 14] in the IT-SOFCs. Those important information are now available in big database. We know that the combination of artificial intelligence (AI) and big database will accelerate the conventional style of R&D innovative work (i.e. incremental technical innovation). However, there is serious weakness in this combination approach. It is 'null set' of big database. Our suggested multidisciplinary research work will input new data (i.e. something new) in the 'null set' of big database and overcome the weakness of AI/big database combination work for the design of the 'state-of-the-art' next generation IT-SOFCs.

With the uptake of this approach and the availability of AI, the authors predict a GREEN innovation wave in the solid oxide full cells development to happen up to 2030 by a combination of AI/big database and our suggested multidisciplinary reference based work in R & D, refer to figure 9.

Figure 9. Schematic diagram of future prospects of research and development work for 'state-of-the-art' IT-SOFC device and system by fusion of fabrication, microanalysis and modelling.

Standard image High-resolution image

As explained in this prospect, the conventional modelling on the basis of kinetics and thermodynamics without microanalysis data can simulate polarization experiments and impedance experiments well. It would be useful for optimization of complicated composition of composite cathode and anode. However, it would be scientifically insufficient to complete the societal implementation of the 'state-of-the-art' IT-SOFC device and system, referring to the part 1 (cathode and cathode/electrolyte interface) and part 3 (anode and anode/electrolyte interface).

After we input the data of interactive works between fabrication and modelling on the basis of microanalysis results of the defect interface into the big database, the deep learning by combination of AI and the big database which includes both data of conventional modelling and our proposed modelling based on the real defect interface will evolve the fabrication techniques and modelling techniques for the social implementation of 'state-of-the-art' IT-SOFC device and system within one decade. That is one typical example of radical innovation as schematically suggested in figure 9.

We foresee that the proposed optimised material development concept will be helpful for acceleration of R&D of the next genearation of IT-SOFCs and open the door for a next GREEN innovation wave in the fuel cell field.

Acknowledgments

The authors express our sincere gratitude to the members of our international collaboration team (i.e. PhD students, postdoctoral researchers, and research assistant staffs in each research team) for the development of the international collaboration between Japan and Australia.

Funding

This work was partially supported by JSPS Joint Research Project (Open Partnership) under bilateral program between Japan and Australia (FY 2019- FY2020, DG 1270) and Australian Research Council (DP180100731). Also, the present work was additionally supported by the Global Research Center for Environmental and Energy based on the Nanomaterials Science (GREEN), National Institute for Materials Science (NIMS), Japan.

Please wait… references are loading.