日本新能源产业技术综合开发机构（NEDO）于2017年制定了促进氢利用扩展的技术战略路线图，并在这一战略路线图中，提出了到2030年左右累计普及80万台FCV的目标。然而，鉴于目前FCV在日本的销售情况，实现该目标仍然遥远。本文中，丰田中央研究所燃料电池第1研究室以及丰田汽车株式会社的FC开发部的研究人员对实现燃料电池性能目标的材料课题进行了讨论和探究。
1、前言
作为实现汽车二氧化碳排放量削减和燃料多样化的技术之一，以氢为燃料的燃料电池汽车（FCV）被寄予厚望。“确保人人获得负担得起、可靠和可持续的现代能源”是可持续发展目标（SDGs）之一，而氢气可以从多种能源中制取，且储存和运输也比较容易，因此FCV的普及和由此引发的氢的大规模利用可以为该目标的实现做出巨大贡献。
日本新能源产业技术综合开发机构（NEDO）制定了促进氢利用扩展的技术战略路线图（NEDO，2017）。在这一战略路线图中，提出了到2030年左右累计普及80万台FCV的目标。鉴于目前FCV在日本的年销量约为1000台的情况（米田，2020），这个目标可以说是非常高的。表1示出为实现该普及目标的FCV目标性能。续航距离的增加有助于扩大FCV的需求，铂使用量的降低有助于削减燃料电池系统的成本和降低资源风险，电堆功率密度的增加有助于向多种车型的扩展。
但是，这些目标性能与当前性能之间的技术差距很大。下面，使用表示燃料电池I-V特性的最简单的公式（公式(1)）来简单地描述这种差距的大小。
式中，U表示氧还原反应（ORR）的热力学平衡电位，b表示ORR塔菲尔斜率，Io表示ORR交换电流密度（∝阴极催化剂表面积×面积比活性），R表示欧姆阻抗（=电子电阻+质子电阻），Ilim表示极限电流密度（∝(气体运输阻力)-1）。
由该式可知，为了提高最大负荷点电压（V）和电堆功率密度（≒I×V），有效方法是增大Io，减小R并增大Ilim（如图1）。可知，通过减少反应表面积，减少贵金属的使用量不仅使Io降低，还降低了Ilim，单电池电压和输出密度的提高之间存在权衡关系。表1的目标数值中要求将贵金属使用量减少到1/3以下，同时将电堆功率密度提高到2倍，解决这一权衡问题需要材料技术的进步。将目标续航距离从650km提高到800km可能看起来比较容易，但实际上需要大幅改良材料特性。例如，考虑通过改善I-V特性只提高5%的发电效率来延长续航距离（提高5%的发电效率无法达到800km的续航距离，因此假设同时采取增加氢装载量、提高辅机效率、降低透过氢量等措施）。
简单地说，要想提高5%的发电效率，就必须将一定电流密度下的电池电压提高5%。因此，在原电池电压为0.80V的情况下，目标单电池电压应该为0.84V。当ORR塔菲尔斜率的典型值b=70mV/decade时，该电压的提高量相当于提高约4倍的Io。将贵金属使用量减少到1/3以下，同时将Io提高到4倍，这是非常高的目标。图2示出燃料电池实用化推进协议会（FCCJ）估算的目标I-V特性。从该图中也可以看出2030年的目标与现状的巨大差距。
另一方面，近年来的研发中也提出了有希望的材料种类。例如，已知通过控制铂系催化剂的形状可以大幅提高ORR活性，特别是，已有报告指出被称为铂纳米线的催化剂的面积比活性非常高，是现有催化剂的约10倍。此外，也有报告指出，通过使用中孔碳作为载体，能够抑制由催化剂层电解质（离聚物）引起的催化剂中毒，从而实现贵金属使用量的减少和输出密度的提高。而且，通过使用高透氧性的离聚物可以大幅改善I-V特性。
如果将这些具有发展前景的技术组合在一起，会产生怎样的I-V特性呢？另外，为了实现2030年的目标，什么样的材料研究将发挥重要作用呢？本文基于I-V特性的数值模拟结果，对这些课题进行了讨论。
2、数值模拟
在本文中，对此前报道的特性值最高的材料（Early-Stage Materials）不相抵地组合，利用数值模拟模型，对组合后的I-V特性进行计算，并将计算结果与目标值进行比较。发电性能虽然也可通过式(1)计算，但该式难以准确估算反应中间体覆盖率对电位和温度的依赖性以及对平面内和平面垂直方向上环境分布的影响等。因此，本文采用了关于ORR和物质运输的更详细的模型。本节将对以计算模型为基础的材料进行简单说明。另外，操作条件是固定的（温度：60℃，入口压力：空气极与氢极在2A cm-2以上时均为200kPa.abs，在0.2A cm-2以下时均为110kPa.abs，其间呈线性互补。供给气体的阻燃比为空气1.5，氢气1.25。氢极在循环系统中的循环比为5）。
2.1 氧化还原反应（ORR）模型
ORR是控制燃料电池性能的反应，其反应速度的表现方式对模拟模型的精度具有重大影响。在许多燃料电池数值模拟模型中，ORR被认为是单一的反应过程，反应速度用巴特勒-福尔默方程式来表示。该方法基于现有催化剂的发电性能进行参数拟合，从而提供适当的结果。然而，在本文中，需要对我们身边以外的材料进行研究，并且不能实施参数拟合。因此，最好使用不需要参数拟合的模型，即尽可能使用基于理论的模型。
此前，很多研究中都对氧还原反应进行了积极的理论解析，并且已经构建了几个型。无论哪种模型，氧还原反应都是通过多个基本反应的组合来表现的。使用由第一性原理模拟得到的参数解出每种基本反应的巴特勒-福尔默方程式，从而计算出整体的ORR速率。然后通过运用与各种催化剂表面上的反应中间体的自由能和反应势垒的高度相关的互相关关系（scaling规则以及Bronsted-Evans-Polanyi规则），重现室温水溶液中的反应活性在各种催化剂表面上的趋势（由实验获得）。虽然在估计部分反应势垒的高度以及其在实际催化剂中的理论应用等方面还存在一些课题，但这些模型提供了对反应活性的理论性理解，明确了提高反应活性的方针，为催化剂的开发做出了贡献。
本文中使用了将Hansen等人的理论模型（Jinnouchi et al., 2011，Hansen et al., 2014）扩展到可以在温度和湿度变化时应用的形式而成的模型（Shibata et al., 2020）。关于温度和湿度变化时的反应中间体的自由能和反应势垒的高度的变化，目前还有许多不明之处，因此采用了几种假设进行计算（假设的详细内容请参考上述文献）。虽然ORR模型是正在开发中的模型，但这是目前能够使用的最妥当的方法。通过使用该模型，可以根据在室温溶液中进行的催化剂模型评估中获得的活性来估算在实际燃料电池环境下的活性。另外，在推测单电池性能时，还需要考虑因聚离物覆盖而导致的性能下降。于是，以氧输送阻力的电位依赖性的相关报告为参考，简单地表示出依赖电位的聚离物阴离子部位的铂覆盖率。
2.2 物质运输模型
在物质输送模型中，以偏微分方程式对分子（氧、氢、氮、水）、电荷（质子、电子）和热输送进行描述，通过与2.1节的ORR模型组合，计算出燃料电池的发电性能。从模拟隔膜和材料立体形状的3D模型，到通过少数参数表现立体形状信息而简化形状的1D模型，已经开发出多种物质输送模型。需要根据目的从这些模型中选择合适的模型。
在本文讨论的范围内，没有必要规定燃料电池和材料的3D结构，因此本文使用了一种“1+1D”模型，该模型由燃料电池隔膜流路的气体流动方向的一维运输模型和面垂直方向的一维运输模型组合而成。另一方面，燃料电池的平面和面垂直方向上发电环境的不均匀性对发电性能有很大影响，因此，考虑到这种不均匀性，使用了“1+1D”模型。计算中使用的运输参数是参考当前产品的评价结果而设定的。
2.3 前提材料
前提是将形状控制催化剂中具有最高铂重量活性的铂纳米线作为活性物质，以易接近多孔碳作为催化剂载体，以高透氧聚离物作为催化剂层电解质，以EW（等价重量）低的PFSA膜作为电解质膜，以无基板的扩散层作为气体扩散层，以低接触电阻部件作为隔膜。以与各个文献中记载的现有产品相比的特性提高量为参考来设定材料特性值。原本希望在讨论各个报告值妥当性的基础上决定使用的值。例如，铂纳米线在模型评估时测得具有高活性，但是没有报告指出其在燃料电池环境下的活性提高量有限。另外，在沟槽状流路的情况下，可能会存在一种权衡关系，即使用无基板扩散层的话，会抑制气体向肋条下的供给。
今后，有必要对这些权衡关系进行定量化。虽然本文没有充分考虑上述权衡关系，但即使是这样的分析，在寻找对提高性能具有重要作用的材料技术方面，也可以获得足够的精度。另外，在假设阴极催化剂电极的贵金属使用量为0.1mg cm-2的情况下进行