Nature-inspired electrocatalysts and devices for energy conversion †

The main obstacles toward further commercialization of electrochemical devices are the development of highly eﬃcient, cost-eﬀective and robust electrocatalysts, and the suitable integration of those catalysts within devices that optimally translate catalytic performance at the nanoscale to practically relevant length and time scales. Over the last decades, advancements in manufacturing technology, computational tools, and synthesis techniques have led to a range of sophisticated electrocatalysts, mostly based on expensive platinum group metals. To further improve their design, and to reduce overall cost, inspiration can be derived from nature on multiple levels, considering nature’s eﬃcient, hierarchical structures that are intrinsically scaling, as well as biological catalysts that catalyze the same reactions as in electrochemical devices. In this review, we introduce the concept of nature-inspired chemical engineering (NICE), contrasting it to the narrow sense in which biomimetics is often applied, namely copying isolated features of biological organisms irrespective of the diﬀerent context. In contrast, NICE provides a systematic design methodology to solve engineering problems, based on the fundamental understanding of mechanisms that underpin desired properties, but also adapting them to the context of engineering applications. The scope of the NICE approach is demonstrated via this comparative state-of-the-art review, providing examples of bio-inspired electrocatalysts for key energy conversion reactions and nature-inspired electrochemical devices. biofuel cells to


Panagiotis Trogadas is Assistant
Research Professor in the Department of Chemical Engineering at UCL since 2017, working with the Centre for Nature-Inspired Engineering. Upon completion of his diploma in Chemical Engineering at National Technical University of Athens and his doctorate in Chemical Engineering at IIT, Chicago, he has undertaken postdoctoral posts in Chemical Engineering at Georgia Institute of Technology and TU Berlin, where he has worked on projects related to hydrogen powered vehicles. His research interests lie in the exploration of nature-inspired electrochemical components, devices and systems. He has won several research grants and delivered invited lectures and keynotes.

Marc-Olivier Coppens
Marc-Olivier Coppens is Ramsay Memorial Professor and Head of Department of Chemical Engineering at UCL, since 2012, after professorships at Rensselaer and TUDelft. He founded and directs the UCL Centre for Nature Inspired Engineering, which was granted EPSRC ''Frontier Engineering'' (2013) and ''Progression'' (2019) Awards. He is most recognised for pioneering nature-inspired chemical engineering (NICE): learning from fundamental mechanisms underpinning desirable traits in nature to develop innovative solutions to engineering challenges. He is Fellow of AIChE, IChemE, Corresponding Member of the Saxon Academy of Sciences (Germany), Qiushi Professor at Zhejiang University, and has delivered 450 named lectures, plenaries and keynotes.

Introduction
The genesis of electrochemistry is related to the discovery of the movement of the muscles in a frog's leg by charge, 1 entitled ''animal electricity''. The ''animal electric fluid'' was an early example of a conducting electrolyte showcasing that the elemental composition of metallic electrodes in contact with the electrolyte is crucial in determining the electrochemical response. 2 These concepts of metallic electrode and electrolyte laid the foundations for the field of electrochemistry.
The growth of this field was slow, as it took at least 200 years to explore the relationship between electricity and chemistry. [3][4][5][6][7] Eventually, in the 20th Century, relations between interfacial properties and rates of electrochemical reactions were established. 8 The term electrocatalysis was introduced in the 1930's, 9 and has been used since to describe the dependence between the nature of electrode materials and electrochemical charge transfer reactions, as the kinetics of the latter vary significantly from one electrode to another. It is a subcategory of heterogeneous catalysis, aiming to increase the rate of electrochemical reactions occurring at the surface of electrodes. It involves the study of electrode kinetics and the determination of the current at an applied electrode potential. [10][11][12] Electrocatalysis is closely related to electrochemical energy conversion, material design and synthesis, and lies at the heart of several technologies related to environmental protection and sustainable development. 13,14 During the past decades, research on the design and synthesis of electrocatalysts for energy conversion has blossomed. Materials employed as electrocatalysts have a dual role: lower the energy barrier for electrochemical reactions, and simultaneously promote the electron charge transfer on their surface. 14-21 A plethora of highly efficient noble or non-noble electrocatalysts for electrochemical devices has been developed, including platinum group metals (PGMs) [22][23][24][25][26][27][28] or alloys,  core-shell PGM alloys, 28,[50][51][52][53][54][55][56] shape controlled PGM nanocrystals, 10,35,[55][56][57][58][59][60][61][62][63][64][65][66][67][68][69][70] PGM nanoframes, [71][72][73][74][75][76][77][78][79][80] and non-precious metal electrocatalysts (NPMCs). [81][82][83][84][85][86][87][88][89][90][91][92][93][94][95][96][97][98][99][100] The activity and stability issues of PGMs and non-PGMs as well as the high cost of PGM electrocatalysts directed research to the exploration of more sophisticated designs to meet the targets of automotive and other industrial applications. These catalysts are usually supported on carbonaceous materials; as this review is not focused on support materials, readers are referred to other published reviews on this topic. 15,22,[101][102][103][104] Inspiration for the design of electrocatalysts and devices can be derived from nature, which is full of hierarchically scaling structures. This critical review demonstrates the scope of a systematic, nature-inspired chemical engineering approach (NICE) across all scales, and emphasizes the difference between bio-imitation and bio-inspiration in design. The bio-inspired approach does not aim at replicating the exact details of an enzyme's active site or the veins of a leaf by copying seemingly useful features; instead, it thoroughly investigates the structural and functional characteristics of biological organisms to generate innovative and better performing electrocatalysts and devices for the same applications. Examples of bio-inspired electrocatalysts for key energy conversion reactions, as well as nature-inspired electrochemical devices are provided in the following sections.

Nature-inspired, integral design of electrocatalysts and devices
With advances in computational tools and material synthesis techniques, the creation of electrocatalysts has become a progressively more sophisticated craft, achieving optimal sizes of nanoparticles, high surface area, and accurate control of crystal facets and the specific position of individual catalytic atoms within the crystal lattice to increase the number of active sites. 105 However, all electrocatalyst designs discussed hitherto focus on alterations at the nano-scale only. To further improve on the properties and use of these electrocatalysts in devices, a new methodology for hierarchical catalyst design and implementation, which considers all scales, is pivotal. Not only at the nanoscale, but also across larger scales, inspiration can be drawn from nature, which provides us with numerous examples of intrinsically scaling, hierarchical structures and includes bio-catalysts for the same electrochemical reactions as those employed in electrochemical devices. Thorough examination of the structure and dynamics of natural systems can reveal key mechanisms that underpin desired properties, such as efficiency, scalability and robustness. For example, the branching structure of trees, the vascular network and mammalian lungs is not arbitrary, but consists of similarly repeated divisions across many length scales; the specific geometric properties of these natural transport networks result in facile scaling, as well as in extraordinary overall thermodynamic efficiency, through the minimization of entropy production. 101,[106][107][108][109][110][111] At meso-to macroscales, hierarchical transport networks in nature obey fractal scaling relationships, some of which relate to Murray's law, [112][113][114][115] or generalizations thereof, discussed in Section 3.1. Careful observation and characterization of geometrical and dynamic patterns in nature could, therefore, help us to discover commonalities. This analysis should be carried out beyond superficial appearances, by seeking underlying physico-chemical, mechanistic principles, while being conscious of how particular features might depend on the constraints of the natural environment -constraints that could differ in several ways from the context of an engineering application. This way to study and learn from nature is the basis for the nature-inspired chemical engineering (NICE) approach. 106 One example illustrating a basis for nature-inspired design of chemical reactors is a typical tree (Fig. 1). The quest for process intensification is an important one in chemical reactor engineering, as multiphase catalytic reactors in particular are plagued by issues related to scale-up, efficiency and robustness. Addressing such problems requires integration of selective, stable catalysis and transport phenomena across scales, which trees resolve through their design. At the macro-scale, tree roots provide support, and extract water and nutrients from the soil, which are transported via its fractal network throughout its volume. The crown of the tree consists of ramified branches and twigs with increasingly smaller diameters following fractal, self-similar scaling. Twigs bear leaves employing a venation pattern at the meso-scale for efficient chemical transport, while, at the nano-scale, leaves contain molecular complexes to capture sunlight and convert CO 2 into sugars and oxygen via photosynthesis, providing the essential nutrients for the growth of the tree. All scales, and the design at and between scales, matter.
The NICE approach we advocate and widely employ at the Centre for Nature Inspired Engineering (CNIE) at UCL in areas ranging from cancer immunotherapy to water purification, bioseparations, and the built environment, allows, in the context of this review's subject, for the design of bio-inspired electrocatalysts and electrochemical devices based on two fundamental mechanisms that are particularly effective in nature, namely, hierarchical transport networks to bridge length scales, and force balancing at individual scales. It should be noted that both principles are relevant to a much broader range of applications in heterogeneous catalysis and reaction engineering.
In terms of hierarchical transport networks, we capitalize on the unique structural characteristics of the lung to design and engineer flow fields for polymer electrolyte fuel cells (PEFCs). The NICE approach to design PEFCs will be discussed in more detail in Section 5.3, but Fig. 2 summarizes the methodology. The unique features of the lung (''Nature'') derive from its ability to bridge length scales and scale-up, irrespective of size, while providing uniform distribution of oxygen into the blood cells and minimizing thermodynamic losses across its structure; this is achieved by transitioning from a fractal tree of increasingly narrowing bronchi to a more uniform channel architecture in the acini, at a bronchiole size where the Péclet number is approximately 1, corresponding to a transition between convection-driven Fig. 2 Systematic, step-by-step employment of the NICE approach for the design and engineering of lung-inspired flow fields for PEFCs. Fig. 1 A tree as an example from nature to inspire the design of chemical reactors, from its scalable, fractal architecture of transport pathways at the macroscale, via the uniform distribution of channels at the mesoscale of the porous catalyst to the structure of the active sites at the nanoscale. Image of the tree (macro-scale) and active sites during photosynthesis (nano-scale) was reprinted with permission from ref. 116. Copyright 2010 American Chemical Society. and diffusion-dominated gas transport (''Nature-inspired concept''). A mathematical model could then be developed and used in computer simulations, based on these characteristics of the lung, to calculate the required number of fractal generations in a fuel cell flow field to achieve uniform distribution of reactants (when Péclet B1) on the catalyst layer (''Nature-inspired design''), after which scalable, robust lung-inspired flow fields could be created by various manufacturing techniques, including a form of metal 3D printing, called selective laser sintering (''Nature-inspired design''). Lung-inspired flow field based PEFCs made out of stainless steel or printed circuit boards (PCBs) exhibited improved performance over commercial serpentine flow field based PEFCs (''Experimental realization'').
With regards to the catalyst itself, which this flow plate feeds, the local environment of the active site is crucial to its activity, selectivity, and stability. The NICE approach explores the nanoconfinement effects induced by local curvature, as well as the impact of chemical and geometrical heterogeneity on transport in nanoporous materials and utilizes it in the design of the electrocatalytic layers. A hierarchically porous structure with a network of broad and narrow pores is required to reduce transport limitations and increase overall rates; nature-inspired catalyst design will be discussed in Sections 4 and 5.
This review discusses a range of examples pertaining to these two fundamental mechanisms. Other ubiquitous fundamental mechanisms employed by nature to underpin attractive properties, such as resilience and adaptability, fall outside of this review and will not be discussed here, however we see these as holding great promise in future electrochemical systems engineering. These mechanisms include dynamic self-organization and self-healing properties of natural systems, and the structure of natural networks (e.g. in ecosystems and metabolic networks), control mechanisms and modularity -in the CNIE, these two additional mechanisms are already applied to address a range of other engineering challenges, and the validation of the NICE approach will help to set the stage for applications in electrochemical devices and their operation.
In summary, nature's universal, fundamental mechanisms are an excellent guide to redesign processes and catalytic materials, as they epitomize efficient, robust, and scalable hierarchical structures. However, equally important is the appreciation that almost any optimization problem, in nature and technology, is a constrained multi-objective optimization challenge; the differences in constraints need to be considered when adopting natural mechanisms to solve engineering problems. Here, we will discuss systematically when and how a thorough scientific investigation into natural processes and materials in relation to their function and structure can provide nature-inspired solutions for the design of electrocatalysts and devices.
3. Ways to connect nature to electrocatalyst design: inspiration vs. imitation Thus far, the majority of research on the development of new electrocatalysts that turns to nature for inspiration is, in fact, based on much more direct bio-imitation, which mimics isolated features of biological or non-biological structures, leading to low activity and durability issues. 101,110,111 We use the term bio-imitation to differentiate this direct approach to the systematic, stepwise, nature-inspired engineering methodology discussed in Section 2. Even though the actual mechanistic features or physical processes that govern the biological system are neglected, there are successful examples of bio-imitating electrocatalysts. 101,110,111 A series of electrocatalysts for hydrogen evolution and oxidation reactions were synthesized, following two strategies: (i) utilization of the active metal center of NiFe or FeFe hydrogenase as catalyst for electrochemical reactions (bio-imitation); 117-120 and (ii) immobilization of enzymes on carbon support (biointegration). 121 Based on the first bio-imitation approach, metallic complexes (such as FeFe or NiFe) using the same metal sites of the active center of hydrogenases were synthesized, and demonstrated activity toward the hydrogen evolution reaction. [117][118][119][120]122 One of the most effective catalysts produced by this approach was a mono-Ni bis-diphosphine complex 123 based on the characteristics of FeFe hydrogenase. Even though it could catalyze the hydrogen oxidation and evolution reaction in organic solvents 123,124 or acidic aqueous solutions (pH = 0-6), 125 its high sensitivity to oxygen hinders its use in fuel cell systems. 126 To circumvent this issue, a bio-integration approach was used with the adsorption of hydrogenase on gold and graphite electrodes. [127][128][129][130] Even though this method contributed to the understanding of the catalytic activity of the enzymes, its major drawback was the weak stability of most hydrogenases to oxygen. 131 A decrease in activity was observed over time as a result of the slow modification of the active site by residual traces of O 2 and protein reorientation at the interface. 129 Additionally, the large size of the metalloenzyme (B6 nm) dictated the development of high surface area electrodes (such as carbon nanotubes) with high loading to achieve current densities comparable to those of Pt electrocatalysts. 128 While the above examples introduce interesting catalysts and provide useful insights, it is evident that the difference in context between nature and technology, such as in size, medium, and position of active metal center, is not accounted for. Doing so could help to make further progress, because a genuinely natureinspired approach is based on the mechanistic understanding of fundamental mechanisms underlying desired traits, without necessarily copying or directly utilizing the source of inspiration; these mechanisms are incorporated into the design and synthesis of artificial systems encompassing the traits of the natural model. 101,[109][110][111] In the following sections, we provide examples of bio-inspired designs of electrocatalysts for key electrochemical reactions for energy conversion at the nanoscale (Fig. 3) and of electrochemical devices at larger scales, including our NICE approach in the design of lung-inspired flow fields for PEFCs.

Nature-inspired electrocatalysts for the oxygen reduction reaction
The overall oxygen reduction reaction (ORR) on an electrocatalyst is a multi-electron, multi-step reaction involving several reaction intermediates. The distinction between the different mechanisms of the ORR relies on the number of proton-coupled electron transfer steps preceding the breaking step of the O-O bond. 132 A four-electron (4e À ) route is desirable to directly produce H 2 O (in acidic medium) and OH À (in alkaline medium) as the final products; depending on the electrocatalytic properties, a twoelectron (2e À ) sub-route can be included producing H 2 O 2 (in acidic medium) and HO 2 À (in alkaline medium) as the intermediate species (Table 1). 27,59,[132][133][134][135][136][137][138][139][140] The fundamental understanding of the ORR mechanism is still unclear, due to the complexity of its kinetics; 132 the rate determining step, sequence of electron and proton additions, 141 as well as adsorption of intermediates are still under debate. The rate determining step is considered to be the first electron transfer from the topmost electrode surface to the adsorbed molecular oxygen (O 2ads ), 135,[142][143][144][145] while the adsorption of intermediate species, such as oxygen (O*) and hydroxyl (OH*) radicals, is suggested to be the crucial step of the ORR mechanism on metallic catalysts. 136,146,147 For metals that bind oxygen too strongly, the reaction rate is limited by the removal of the adsorbed O* and OH* species. For metal surfaces that bind oxygen too weakly, the reaction rate is limited by the dissociation of O 2 , or, more likely, the transfer of electrons and protons to adsorbed O 2 . 135,136 Successful bio-imitation approaches for the design of electrocatalysts for the oxygen reduction reaction have been employed, mimicking the active metal center of enzymes coordinated by carboxyl, amine, thiol, and imidazole groups of the amino acid side chains. 148 Metallophthalocyanines, [149][150][151][152][153][154][155][156][157][158] and metalloporphyrins 149,150,[159][160][161][162][163][164][165][166][167][168][169] have been extensively investigated as non-precious metal electrocatalysts for ORR, even though they lack long-term stability in the harsh environment of a fuel cell. 159,170,171 Iron phthalocyanines catalyze the direct reduction of oxygen to water via a four-electron pathway, promoting scission of the O-O bond. 149,150 On the contrary, the ORR activity of metalloporphyrins depends on whether their meso-positions are occupied by aryl, pyridyl, or alkyl substituents. 168,169 For example, cobalt metalloporphyrins adsorbed on graphite electrodes without meso substituents demonstrated two-electron oxygen reduction, whereas the adsorbed cobalt porphyrins with meso substituents directly oxidized oxygen to water via a four-electron pathway. 169 This catalytic behavior was attributed to the spontaneous formation of van der Waals dimers facilitating the reduction of oxygen to water via four electrons. 169 The interaction between the cobalt centers of the porphyrin dimer and the atoms of oxygen leads to a transition state where the bridging oxygen molecule is activated, accepting four electrons from the graphite electrode, which causes the scission of O-O bond. 169 Based on these attractive activity characteristics, metal organic networks employing well-defined metal surfaces decorated with porphyrins or phthalocyanines have been synthesized via supramolecular chemistry. 148,155,166,167,172,173 Iron phthalocyanines adsorbed on the surface of Au(111) enhanced its ORR activity in alkaline media (B0.030 V Tafel slope), allowing the direct reduction of oxygen to water via a four-electron route with cleavage of the O-O bond, in contrast to the observed oxygen reduction via a two-electron route in the case of bare Au (111) surface. 155 Additionally, chemisorbed cobalt porphyrin on the surface of Au(111) increased its ORR activity in acidic media, 166,167 demonstrating an increase of the reductive current. 166,167 Despite these promising preliminary results, these bio-imitating catalysts are in early development and challenges concerning long-term stability and activity, as well as the tuning of adsorbate-substrate interactions have to be resolved. 174 Recently, mononuclear iron porphyrins with distal basic residues (i.e., pyridine, aliphatic and dibenzyl aromatic amines) imitating the proton transfer pathways and hydrogen bonding groups in cytochrome c oxidase 149,175-179 demonstrated high selectivity (greater than 90%) and reactivity (rate constant greater than 10 7 M À1 s À1 ) for ORR at pH = 7 when adsorbed on an edge-plane graphitic electrode. 180 The measured rate constant of oxygen reduction was approximately two orders of magnitude higher than the heme/Cu complexes, [181][182][183][184][185][186] suggesting that the distal basic residues stabilize the Fe III -OOH intermediate species and enable heterolytic cleavage via protonation of their distal oxygen, promoting a four-electron ORR route. 180 Another promising example of a bio-imitating, porphyrinbased ORR electrocatalyst mimicking the active site of oxygenactivating heme-containing enzymes, such as cytochrome c oxidase, is an axial-imidazole coordinated iron porphyrin, covalently grafted on multi-wall carbon nanotubes (MWCNTs). 187 This electrocatalyst exhibited high activity and stability in acidic and alkaline media; its half-wave potential (E 1/2 ) in both media was B38 mV (E 1/2 B 0.88 V vs. NHE, 0.1 M HClO 4 ) and B47 mV (E 1/2 B 0.922 V vs. NHE, 0.1 M KOH) higher than the commercial Pt/C (E 1/2 B 0.842 V and B0.875 V vs. NHE in acidic and alkaline media, respectively). Koutecky-Levich plots revealed a fourelectron oxygen reduction route, with the production of a minimal amount of hydrogen peroxide (less than 1% H 2 O 2 ). The stability of this bio-imitating electrocatalyst was evaluated via US DoE's accelerated test protocol, demonstrating a B14 mV loss in E 1/2 , which was less than half of the decrease exhibited for commercial Pt/C (B37 mV loss) after 10 000 continuous cycles between 0.6 V and 1 V (vs. NHE) in an oxygen saturated acidic solution (0.1 M HClO 4 ). 187 The high activity and stability of this electrocatalyst was attributed to its low half-wave potential value for ORR and, hence, the associated low amount of H 2 O 2 produced, insufficient to deactivate and degrade this electrocatalyst. 187 Furthermore, the ORR activity of metalloporphyrins is affected by their structure and the distance between their active metal centers. 188,189 Hence, crystalline porous materials such as metal organic frameworks (MOFs) or covalent organic frameworks can be used as scaffolds to accurately control the structure of metal-metalloporphyrin frameworks, 190 resulting in efficient ORR electrocatalysts with high surface area and porosity. However, their major disadvantage is their low electronic conductivity and long-term stability, 191 and, thus, further research is needed to effectively utilize these electrocatalysts in fuel cells.
Metallocorroles have also attracted attention as electrocatalysts for the ORR due to their high activity and tunable structure. [191][192][193][194][195][196][197][198][199] Metallocoroles are tetrapyrrolic macrocyclic compounds with one less carbon atom than porphyrins. 191 Their tri-anionic charge stabilizes the metal center and affects the chemistry of the chelated transition metal ions, resulting in low valent corrole-transition metal complexes, which are much more reactive than their porphyrin analogs. 191 Several transition metal complexes with brominated corrole (b-pyrrole-brominated 5,10,15-tris-pentafluorophenyl-corrole, M(tpfcBr 8 ), where M = Mn, Fe, Co, Ni, and Cu) adsorbed on a high surface area activated carbon (BP2000) have been investigated as electrocatalysts for ORR in alkaline media. 194 Iron-and cobalt-based brominated corroles exhibited the highest activity following a four-electron pathway for the direct reduction of oxygen to water, similar to the ORR activity of commercial Pt/C (20 wt% Pt) tested under the same conditions. 194 In acidic media, the ORR activity increases in the order of Co 4 Fe 4 Ni 4 Mn 4 Cu, with brominated cobalt-corrole demonstrating the best activity with an onset potential B80 mV lower than commercial Pt/C (20 wt% Pt). 199 In addition to the metal center of the brominated corrole, its ORR activity is suggested to depend on the support, since adsorption of metallocorroles on activated carbon or carbon nanotubes decreases the reaction overpotential. 191,[199][200][201] Further research is needed to fully comprehend this phenomenon.
Despite these promising ORR activity results, the incorporation of these bio-imitating electrocatalysts into PEFCs remains a great challenge. Their activity and stability have been evaluated only at RRDE level, and not in a membrane electrode assembly under harsh fuel cell conditions. Thus far, these electrocatalysts demonstrate poor stability and improved catalyst designs would have to be introduced to meet the requirements for practical PEFC applications. 202 As this review is not focused on bio-imitating electrocatalysts, readers are referred to other published reviews for more information on this topic. 148,191,202,203 Since the reduction and evolution reactions occurring in a fuel cell comprise three phases (gas/liquid/solid), electrocatalysts with a hierarchical porous structure are desirable to promote mass transfer and improve utilization of reactant. To improve the design of electrocatalysts, the hierarchical pore networks utilized by biological organisms (such as trees, other plants and their leaves, etc.) were used as a source for inspiration. These organisms grow obeying power laws that are generalizations of Murray's law, utilizing precise diameter ratios to connect pores from macro-to meso-and micro-scale, which results in the minimization of transport resistance and fluent transfer throughout the network. 113,204 According to Murray's law, the cube of the diameter of the parent vessel (d p ) is equal to the sum of the cubes of the diameters of the daughter vessels (d i , i = 1. . .n, where n is the number of macro-, meso-, or nano-pores in each particle) at each level of bifurcation: [109][110][111] Based on eqn (1) and several assumptions presented in detail in the ESI † (Section S1), a series of new eqn (S2)-(S5) (ESI †) correlating the macro-, meso-, and micro-pores of a ''Murrayinspired'' material were suggested. 205 A layer-by-layer, evaporation-driven self-assembly process utilizing a microporous material (ZnO, in this example) as the primary building block under ambient conditions is employed for the synthesis of these so-called ''Murray materials''. 205 The synthesis procedure of these materials is a tedious, multi-step process, since various reaction conditions are utilized to produce a combined micro-, meso-, and macro-porous ZnO (Fig. 4). Microporous ZnO is prepared in inert atmosphere using a mixture of zinc acetylacetonate hydrate and oleylamine. Mesoporous ZnO nanoparticles (B30 nm size) are synthesized using the same method with additional calcination at B290 1C.
To prepare the micro-, meso-, and macro-porous material, micro-and meso-porous ZnO particles are dispersed separately in a volatile solvent (e.g. hexane, ethanol). The first microporous layer is created by drop-casting of microporous ZnO solution on a silicon wafer and evaporation of the volatile solvent. Then, mesoporous ZnO solution is drop-casted on the already formed ZnO microporous layers.
The concentration of ZnO in solution directly affects the evaporation rate of the solvent. At low ZnO concentration (B0.03 mg mL À1 ), the nanoparticles are self-assembled into isolated islands; small macropores gradually appear as the ZnO concentration is increased (B0.06 and 0.12 mg mL À1 ) due to the enhanced nucleation and restricted expansion of the holes. Further increase in ZnO concentration (B0.5 mg mL À1 ) leads to a porous network with ZnO nanoparticles forming packed structures. 205 These initial experimental results are very promising; however, the general applicability of Murray's law to the synthesis of a wide range of materials for various applications has to be thoroughly investigated. The synthesis method advocated for hierarchical Murray ZnO may not be feasible or appropriate for different materials; hence, alternative synthesis methods should be explored, such as templating, crystallization, the Kirkendall effect, etc. 101 Recently, Murray-like cobalt supported on nitrogen doped carbon (Co-N-C) electrocatalysts for oxygen reduction and evolution, and hydrogen evolution reactions in alkaline medium were synthesized using Prussian blue analogue (PBA) as a precursor. 206 Pyrolysis was used instead of the layer-by-layer drop casting method (Fig. 5). The 3D cubic porous network of PBA imitated the hierarchical structure obtained by Murray's law, while PBA acted as the carbon, nitrogen, and cobalt source, and polyvinyl pyrrolidone (PVP) as the surfactant. After thermal treatment (pyrolysis at 900 1C in inert atmosphere), a spherical porous assembly of Co nanoparticles covered by nitrogen doped graphitic carbon was produced. Acid leaching (2 M HCl solution) was employed to remove the unstable species and optimize the porosity of the material, increasing its electrochemically active surface area (ECSA) as demonstrated by double layer capacitance measurements (B8 and B7 mF cm À2 for acid and non-acid treated Co-N-C, respectively). 206 Physicochemical characterization (X-ray diffraction (XRD), Raman and X-ray photoelectron spectroscopy (XPS), scanning (SEM) and transmission (TEM) electron microscopy) of the material post acid treatment demonstrated that its structure remained intact, while the composition and valence states of Co, N, C, Zn, and O at the surface of the electrocatalyst were similar to the ones before the acid modification step.
As a result, the acid treated Murray-type Co-N-C materials exhibited similar ORR activity (diffusion-limited current density of 5.7 mA cm À2 ) and electron transfer (B4e À ) to commercial Pt/C in alkaline solution (Fig. 6a). In the case of OER and HER (Fig. 6d), acid treated Murray-type Co-N-C displayed low onset overpotentials of B150 and B350 mV, corresponding to the hydrogen evolution reaction and oxygen evolution reaction in alkaline solution (1 M KOH), respectively, highlighting its potential to be used in overall water splitting with a total splitting voltage of 1.73 V. 206 Finally, these acid treated Murray-type electrocatalysts exhibited immunity to methanol (2 M) compared to commercial Pt/C (Fig. 6b), as well as higher stability than nonacid treated Murray-type and commercial Pt/C, with B99% retention at 0.5 V (vs. RHE) after a B6 h continuous chronoamperometric measurement, compared to B93% and B84% retention for non-acid treated Murray-type Co-N-C and commercial Pt/C, respectively (Fig. 6c). 206 Hence, the ORR activity of these electrocatalysts is affected by N species, whereas the OER and HER activity is determined by the Co species, providing more sites for the adsorption of H* and OH* radicals. 59,207,208 Despite the low water splitting performance of Murray-type Co-N-C electrocatalyst compared to literature values, 209 these initial results demonstrate the potential of Murray-type electrocatalysts with optimal porous structure to be utilized in triphase electrochemical reactions.  Moreover, inspiration for the design of new electrocatalysts for ORR can be derived from the active metal center of cytochrome c oxidase and laccase, which comprises iron and copper ion complexes. 210,211 However, directly mimicking the Cu 2+ or Fe 3+ complexes is not a fruitful strategy, as it leads to low activity, due to the absence of mediators for the transfer of electrons and the steric variation of the coordination structures after these complexes are attached onto the electrode. 212 As an example of effective tuning of the electron density of such an active site (Cu 2+ based), copper nanocomposites (CPG) were synthesized via the pyrolysis of a mixture of graphene oxide (GO) and Cu(phen) 2 (Cu 2+ 1,10-phenanthroline). 212 The electron density of Cu 2+ sites was tuned by the electron donation effect from Cu 0 of copper nanoparticles and nitrogen ligand incorporated into graphene. The electron transfer in CPG was increased via the electronic connection of Cu 2+ and nitrogen, as well as Cu 2+ and Cu 0 in graphene, resulting in enhanced ORR activity. 212 Pyrolysis temperature also affected the activity of CPG, with 900 1C being the optimum in this study. 212 The exhibited ORR onset potentials of B0.85 and B0.98 V (vs. RHE) and Tafel slopes of 71 and 49 mV dec À1 in 0.5 M H 2 SO 4 and 1 M KOH, respectively, were similar to commercial Pt/C, demonstrating the high activity of this catalyst. 212

Nature-inspired electrocatalysts for the oxygen evolution reaction
The oxygen evolution reaction (OER) is the reverse of the ORR: water is oxidized to oxygen via a 4e À pathway ( Table 2). [213][214][215][216][217] The mechanism of the OER is sensitive to the structure of the electrode surface.
The symbol ''*'' represents a surface with one oxygen vacancy site in the topmost layer, while the symbols ''OH 2 *'', ''OH*'', ''O*'', and ''OOH*'' represent the surface with the corresponding chemisorbed species residing in the oxygen vacancy site. 59 The complete OER mechanism in acidic medium consists of four oxidation steps, in each of which a proton is released into the electrolyte. Water is first adsorbed onto the surface of the oxygen vacancy site, forming ''OH 2 *'' species. ''OH 2 *'' then undergoes two subsequent oxidation reactions to form ''O*'', which reacts with another water molecule to form an ''OOH*'' intermediate. In the last step, oxygen is released from ''OOH*''. 59 The oxidation of water to oxygen in nature is catalyzed by the oxygen-evolving complex (OEC) in photosystem II, located within the thylakoid membranes of plants, algae, and cyanobacteria. [218][219][220] Its active site responsible for water oxidation is the oxygen-evolving complex, which consists of four manganese ions and a calcium ion (Mn 4 CaO 5 ) surrounded by a protein. [219][220][221][222] This complex has a sophisticated structure, in which a Mn 3 CaO 4 heterocubane is tethered to a fourth manganese ion via a m-oxo bridge and a corner oxo-attachment. [223][224][225] Even though the oxygen-evolving complex is still not fully understood at mechanistic level, metal-based clusters with cubane cores have already been synthesized as efficient catalysts for oxygen evolution reaction, imitating the Mn 4 CaO 5 active site of photosystem II. 221,224,[226][227][228] To improve the design of these bio-imitating electrocatalysts, fundamental understanding of the OER mechanism by oxygenic multi-metallic clusters is required, along with a significant improvement of their stability under the harsh oxidizing conditions of OER. 229,231,232 The approach that nature utilizes to circumvent the oxidative degradation of photosystem II by continuously replacing and repairing the damaged protein subunit 220,229 in its structure cannot be applied in these artificial electrocatalysts.
The state-of-the-art electrocatalysts for the OER are ruthenium (RuO 2 ) and iridium oxide (IrO 2 ), but their high cost prohibit their application at large scale. 230 Hence, extensive research is conducted on the synthesis and design of efficient and cost-effective electrocatalysts based on transition metals and layered double hydroxide materials. 213,[231][232][233][234][235][236][237][238][239][240][241][242][243][244] The design of new electrocatalysts for OER is focused on three scale levels: (i) at the atomic scale, the alteration of oxidation state, 245 coordination, 246,247 electron density, 245 and composition of metal composites 245,248 leads to an improvement of their electronic structure and an increase in OER activity; (ii) at the nano-scale, different material combinations (such as metal oxides), 249,250 chalcogenide, 249 boride, 251 nitride, 252 phosphide 250 supported on nanostructures (e.g. nanowires, 253,254 nanosheets, 251,254 nanotubes 255,256 ) have increased activity toward OER due to their high surface area and number of active sites; and (iii) at the meso-scale, the creation of a porous architecture of the support can enhance mass transport of the electrolyte, improving the activity and conductivity of the electrocatalyst. 257 The structure of a leaf has served as a source of inspiration for the design of a 2D/1D cobalt oxide (CoO x ) electrocatalyst (Fig. 7). 258 Each leaf has a 2D morphology with optimized surface area, favoring light absorption and surface reactions, while the presence of 1D, hollow tubular structures under the leaf facilitates the transport of nutrients and water to each leaf. 258 The leaf-inspired CoO x electrocatalyst consists of CoO x nanosheets (100-200 nm average diameter and 20-50 mm length) in a nanotube (200 nm average diameter, 20 mm length), which was synthesized via in situ etching of a Cu 2 O nanowire  258 Hence, the structure of the CoO x electrocatalyst is improved across different scales; at the atomic scale, the presence of Co 2+ in octahedral symmetry increases electron transfer, while the high surface area (B371 m 2 g À1 ) of CoO x nanosheets and 3D porous framework of CoO x nanotubes provide a high number of active sites and increased ion transport, respectively. 258 The leaf-inspired CoO x electrocatalyst demonstrated an onset potential of B1.46 V (vs. RHE) similar to commercial IrO 2 , a current density of B51 mA cm À2 at 1.65 V (vs. RHE) and a Tafel slope of 75 mV dec À1 . However, its stability is a concern, since a B35% loss in its current density (1.5 mA cm À2 initial current density at 1.5 V vs. RHE) was observed after 2 h of operation in an electrolysis cell. 258

Nature-inspired electrocatalysts for hydrogen oxidation and evolution reactions
The overall reactions for hydrogen oxidation (HOR) and evolution (HER) reactions involve either protons in acidic media or hydroxide ions in alkaline media (Table 3). 59,259,260 The first step of the reaction mechanism always involves the adsorption of a H intermediate on the electrode surface via a proton and electron transfer from the acid electrolyte and the surface of the electrode, respectively. In acidic media, H + is the proton source for the initial Volmer step, whereas, in alkaline media, H 2 O constitutes the proton source producing OH À after electron transfer. Recently, it was suggested that the Volmer step of the HOR/HER is the rate determining step for noble metals. 260 There are two different reaction routes for the final step, determined by the Tafel slope values from polarization curves. 259 The Heyrovsky reaction, in which the adsorbed hydrogen (H*) combines with an electron transferred from the electrode surface and a proton from the electrolyte to form hydrogen; or the Tafel reaction, in which two adsorbed hydrogen combine to form a hydrogen molecule (Table 3).
Pt based electrocatalysts are the state-of-the-art choice for hydrogen oxidation and evolution reactions in acidic media, as earth-abundant electrocatalysts cannot survive under these highly acidic conditions and exhibit similar activity to their Pt based counterparts. 261,262 Even though the HER/HOR activity of Pt based electrocatalysts is decreased in alkaline media, their ability to catalyze these reactions at moderate overpotentials makes them the most efficient choice. 261 The utilization of alkaline media for HER and HOR reactions enables the use of non-precious metal catalysts, such as heteroatom-doped carbon, transition metal oxides, sulfides, phosphides and their alloys, [263][264][265][266][267][268] and molecular complexes based on porphyrin and corrole. 150,264,269,270 Less attention has focused on non-precious metal electrocatalysts for the HOR. [271][272][273][274][275] However, these non-Pt electrocatalysts suffer from poor kinetics and stability; modification of their electronic structure, crystal facets, 265 composition variation, 273,276 and surface area is required to optimize their hydrogen adsorption free energy, a key descriptor of their HOR/HER activity, 14,[277][278][279] and improve their catalytic performance and durability.
Inspiration for the design of new electrocatalysts for the HER and HOR can be derived from nature, and from hydrogenases, specifically. Hydrogenases are metalloenzymes widely used in nature for hydrogen evolution and oxidation in living systems, since they can achieve low overpotentials and high turnover frequencies employing Fe and/or Ni in their active site. 126,280 Hydrogenases comprise three coordination spheres (Fig. 8). 281,282 The inner coordination sphere consists of the metal center and the atoms bonded directly to it. The active site of FeFe hydrogenases has two iron atoms with CN À and CO ligands attached to each of these two atoms (Fig. 9), while the active site of NiFe hydrogenases has nickel ions deposited in cysteinate ligands, which are in contact with the dicyanocarbonyl iron center (Fig. 10). [283][284][285][286][287][288][289] The most important ligand in this inner coordination sphere is azadithiolate, which plays an important role in hydrogen cleavage and proton transfer between the metal center and surface of the enzyme. [283][284][285][286][287][288][289] The second coordination sphere contains functional groups incorporated into the ligand structure interacting directly with substrates bound to the metal during a reaction, but do not interact with the metal center. In hydrogenases, pendant bases are positioned in this coordination sphere to facilitate proton transfer between the metal center and the acid or base in the solution. 281 The outer coordination sphere encompasses the remainder of the ligand structure and the solvent surrounding the catalyst. 281 The structure of hydrogenases inspired the design of costefficient molecular catalysts oxidizing or producing hydrogen. The most notable are Dubois catalysts, synthetic nickelbis-diphosphine complexes using pendant amines as a Lewis base. 124,[290][291][292] The pendant amine is positioned near the iron center functioning as a proton relay promoting the creation or scission of the H-H bond. 124,[290][291][292] As a result, these nickel complexes exhibit high turnover frequencies for H 2 production of B33 000 and 100 000 s À1 in dry acetonitrile and water, respectively. 290 However, these catalysts suffer from severe stability issues when immobilized onto an electrode; 293-295 a B25% loss in catalytic current within 6 h was observed during chronoamperometric measurements in acidic media (0.1 M HClO 4 ), while their activity toward hydrogen oxidation was lost after 10 min upon the addition of B2% O 2 in the hydrogen gas feed. 294 Incorporating amino acids or peptides directly into the P 2 N 2 ring of hydrogen oxidation catalysts [Ni(P Cy 2 N AminoAcid 2 ) 2 ] n+ (where Cy is cyclohexyl; amino acid is either glycine with n = 3, or arginine with n = 7), results in efficient Ni based complexes for  hydrogen oxidation and production. 296 An arginine-based Ni complex exhibited a turnover frequency of B210 s À1 at acidic pH (0), which decreased at higher pH values, reaching B40 s À1 at pH = 7. 296 This arginine-based Ni complex was immobilized onto single-wall carbon nanotubes (SWCNTs) covalently modified with naphthoic acid groups and used as the anode (0.3 mg cm À2 loading) of an H 2 /O 2 PEFC. Pt/C was used in the cathode (1 mg cm À2 loading) and Nafion s as the polymer electrolyte in the MEA. PEFC with the arginine-based Ni complex exhibited an OCV of B0.9 V and B14 mW cm À2 power density at B0.5 V (Fig. 11), in comparison to an OCV of B1 V and 94 mW cm À2 power density for a PEFC employing Pt/C as the anode and cathode catalyst. These initial in situ PEFC performance results demonstrated the ability of an arginine-based Ni complex to operate at highly acidic conditions, delivering only five times less current or power density at very low loading compared to Pt/C. 297 These preliminary PEFC performance results are superior to the ones reported earlier utilizing arginine-based Ni bis-diphosphine complexes as the anode catalyst. 127,298,299 In the case of a PEFC with cobalt (Co) supported on nitrogen doped carbon as the cathode and arginine-based Ni bis-diphosphine complexes supported on multi-wall carbon nanotubes, an B0.74 V OCV and 23 mW cm À2 max power density was reported, whereas an OCV of B0.85 V and 70 mW cm À2 max power density was demonstrated when Pt/C was used as the cathode (operation at 60 1C, B0.45-0.65 mg cm À2 loading at each side). 300 The utilization of Ni bis-diphosphine complexes as hydrogen oxidation and evolution catalysts is not novel; it was proposed a decade earlier by incorporating these complexes into multi-wall carbon nanotubes via grafting. 127 However, the preliminary MEA testing results for hydrogen production and evolution reactions were disappointing, since a two orders of magnitude smaller current density was reported for Ni bis-diphosphine complexes based MEAs compared to Pt/C MEAs. This discrepancy was attributed to the different catalyst loadings used; Ni bisdiphosphine complexes based MEAs had a B0.06 mg cm À2 loading compared to B0.5 mg cm À2 loading for Pt/C MEAs. 127 Recently, a revisited design was reported, where ligands were covalently immobilized via amide coupling on the surface of MWCNTs and, then, the nickel metal centers were introduced with [Ni(CH 3 CN) 6 ] 2+ or [Ni(H 2 O) 6 ] 2+ . 298 The activity of these immobilized Ni bis-diphosphine complexes on MWCNTs were compared to commercial Pt/C (46% wt) MEAs with 0.05 mg cm À2 loading in a rotating disk electrode (RDE) setup using 0.5 M H 2 SO 4 solution (Fig. 12). The total amount of Pt present in the electrode was B2.5 Â 10 À7 mol cm À2 , whereas the total amount of Ni present in the electrode was 10-fold lower (B2.5 Â 10 À8 mol cm À2 ). At room temperature (B25 1C), the activity of the bio-inspired electrocatalyst was approximately one order of magnitude less than the activity of commercial Pt/C; at 85 1C operation, though, the bio-inspired catalyst was only B30% less active than Pt/C for hydrogen oxidation, while it outperformed the commercial catalyst by B20% for hydrogen production. 298 A chemically inert polymer P(GMA-BA-PEGMA) (poly(glycidyl methacrylate-co-butyl acrylate)-co-poly(ethylene glycol)methacrylate) was used to immobilize Dubois' Ni bis-diphosphine  complex with cyanide and glycine [Ni(P Cy 2 N Gly 2 ) 2 ] 2+ on the surface of the electrode and prevent the electrical contact between the metal complexes within the polymer and the electrode. 280 The thick polymer film creates two distinct regions, where different reactions occur. In the region close to the surface of the electrode, hydrogen is oxidized via the Ni-complex, while, in the outer region, which is electrically disconnected from the electrode, preventing re-oxidation of the catalyst, a protonated Ni 0 complex reduces oxygen to water by using the electrons produced from the hydrogen oxidation reaction in the inner layer (Fig. 13).
The thickness of the polymer film is crucial for the effective protection of the metal center from oxygen. Even though thinner films produced high current densities for the oxidation of hydrogen, they degraded immediately upon the addition of oxygen in the gas feed (Fig. 14a). For thicker polymer films, the current density for hydrogen oxidation was not affected by the presence of oxygen in the gas feed for the first 7 h; after 24 h of continuous exposure to oxygen, though, the current density decreased to B75% of its initial value, due to the oxidation of the Ni bis-diphosphine complex and oxidation of phosphines (Fig. 14b). 280,301 3.4 Nature-inspired electrocatalysts for the carbon dioxide reduction reaction Electrochemical reduction of carbon dioxide (CO 2 ) has attracted great research interest, as it could contribute to the reduction of the greenhouse effect and generate useful chemicals. It involves multi-proton and multi-electron routes to convert CO 2 to an activated form of carbon (with the exception of carbon monoxide), which have sluggish kinetics and require high reducing potentials. 302,303 Theoretically, the formation of different carbon products depends on the applied potential, 304 even though directing the reduction of CO 2 toward a specific reaction pathway is difficult to achieve. 305 The initial transfer of a single electron to CO 2 to create CO 2,ads is considered the rate determining step, due to the high energy input required to enable this step. 305 The reactivity of CO 2,ads on the surface of metal catalyst dictates the reaction pathway that will be chosen, leading to different end products. For example, in the presence of gold or silver catalyst, CO 2,ads is converted to CO, whereas copper catalyzes the reaction of CO 2 to formic acid (HCOOH), methane (CH 4 ), and ethylene (C 2 H 4 ) (Fig. 15). 305 Several electrocatalysts for the reduction of CO 2 have been developed, 306 including metal nanoparticles and oxides, chalcogenides and carbon based materials, and their activity was optimized via modification of their defect density on their surface, 302,307-309 particle size, 310 morphology, 311,312 and electrode thickness. 313 Another important issue contributing to the sluggish kinetics of CO 2 reduction is the insufficient concentration of CO 2 on the surface of the catalyst layer. The concentration of carbon dioxide near the cathode increases with increasing metal cation size, 315 but this effect is restricted by the solubility of relevant salts. The application of large potentials can also improve the adsorption of CO 2 at the expense of increased hydrogen evolution. 316 On the contrary, nanostructured metal electrodes (Au, Pd nanoneedles) produce a high local electric field at low overpotentials, which concentrates electrolyte cations and enhances the concentration of CO 2 near the cathode catalyst layer. 317 High current densities of B22 mA cm À2 at À0.35 V vs. Ag/AgCl (gold nanoneedles) and B10 mA cm À2 at À0.2 V vs. Ag/AgCl (palladium nanoneedles) for CO and formate production were reported, respectively. 317 Hence, bio-inspired approaches have been employed to improve the activity of electrocatalysts toward reduction of CO 2 , via inspiration from biological processes or the structure of enzymes, which catalyze the same reactions. These will now be reviewed.
3.4.1 Conversion of carbon dioxide to carbon monoxide. Electrochemical reduction of carbon dioxide to carbon monoxide  (CO) in the presence of oxygen is a challenging task, since oxygen reduction is thermodynamically favored over CO 2 reduction. 35 To overcome this obstacle, the stream of CO 2 is purified before entering the electrolyzer, a task which requires the addition of an additional gas separation system, increasing the overall cost of the reactor setup. 318,319 On the contrary, in nature, reduction of carbon dioxide occurs during photosynthesis, in which enzymes are located near the active sites to significantly increase the local concentration of CO 2 123,320 and, hereby, minimize the O 2 to CO 2 ratio, which enhances the conversion rate of CO 2 to glucose, despite the presence of oxygen. 321 An oxygen tolerant cathode was recently synthesized, inspired by the process leading to an increase in CO 2 concentration during photosynthesis (Fig. 16a). 322 This bio-inspired cathode consists of a carbon fiber gas diffusion electrode with a catalyst layer and a polymer with intrinsic microporosity (PIM) 323 coated on its opposite sides. 322 PIM imitates the role of enzymes; it filters oxygen and is highly permeable to CO 2 . Its size-selective pores reject oxygen molecules (B0.35 nm kinetic diameter) 324,325 and facilitate the transport of slightly smaller CO 2 molecules (B0.33 nm kinetic diameter) 324,326 effectively decreasing the O 2 to CO 2 ratio of the feed stream reaching the catalyst layer, and, hence, favoring the reduction of CO 2 to CO (Fig. 16b). 322 The incorporation of this bio-inspired electrode into the electrolyzer resulted in a B76% faradaic efficiency to CO and B27 mA cm À2 current density at À1.1 V vs. RHE (feed gas stream with 5% O 2 ). As the oxygen concentration in the gas stream increased, the faradaic efficiency of CO decreased, reaching B20% and B35 mA cm À2 current density at 35% O 2 in the gas stream, indicating that oxygen reduction dominated the electrode reaction at high O 2 concentrations in the gas stream (Fig. 16c). 322 During continuous operation over 18 h, stable operation was observed for gas streams containing 5% and 20% O 2 (Fig. 16d)   Another electrode design employed to increase the concentration of CO 2 at the surface of the catalyst layer was inspired from the structure of alveoli in a mammalian lung. 327 Alveoli consist of thin epithelial cellular membranes with low water permeability and high gas diffusivity. 107,108,328 During pulmonary circulation, inlet air rapidly penetrates through bronchii to alveoli and reaches the blood cells, where haemoglobin protein binds oxygen and releases CO 2 . 327 The structure of the alveoli-inspired electrode (B20-80 nm thickness) comprises a layer of gold nanoparticles (B0.15 mg cm À2 loading) acting as the catalyst sputtered on a polyethylene membrane, whose hydrophobicity and network of interconnected fibers (B40-500 nm pore size) allow unobstructed diffusion of CO 2 toward the catalyst layer. To fully replicate the structure of a closed alveolus with macroscopic tubes (bronchioli) for gas transport to and from the lung, the flat Au coated polyethylene membrane was rolled and sealed to form a bi-layer pouch-type structure. Electrochemical reduction of CO 2 was conducted in an H-type cell using a Selemiont anion exchange membrane as separator and CO 2 saturated potassium bicarbonate (0.5 M KHCO 3 ) as the electrolyte. A faradaic efficiency of B92% for CO production and high current density of B25 mA cm À2 at À0.6 V vs. RHE were achieved. 327 An additional nature-inspired strategy for the design of electrocatalysts for the reduction of CO 2 is derived from the structure of dehydrogenases. CO-(CODH) and formate-dehydrogenases (FDH) catalyze the reversible CO 2 reduction via their single-or multi-metallic active sites composed of earth abundant metals, such as nickel-iron or molybdenum-sulfur-copper complexes in CODHs or single metals tungsten or molybdenum in FDHs. [329][330][331][332] In addition to their active metal centers, proton relays on their outer coordination sphere play an important role in the catalytic activity of these enzymes. [329][330][331][332]  NiFe CODHs utilizing [NiFe 4 S 4 ] as the active site are highly active toward CO 2 to CO conversion with B12 s À1 turnover frequency at low overpotentials (less than 100 mV). [333][334][335][336] Ni acts as a Lewis base transporting an electron to the unoccupied molecular orbital of CO 2 , increasing the negative partial charges of oxygen atoms bound on the Fe center acting as a Lewis acid. [333][334][335][336] Moreover, a CODH-inspired catalyst containing cofacial Fe tetraphenylporphyrin dimer (o-Fe 2 DTPP) demonstrated a high faradaic efficiency (B95%) and turnover frequency (B4300 s À1 ) at high overpotential (B0.7 V) for the conversion of CO 2 to CO in a DMF-water solution (10 wt%). 337 The introduction of electronically different substituents to porphyrin positions can alter the overpotential and turnover frequency. Attachment of electron withdrawing groups, such as perfluorophenyl substituents to the meso-positions of the porphyrin rings of CODH-inspired catalyst decreased the overpotential to B0.3 V; however, the withdrawal of electrons also reduced the electron density of the active metal center, resulting in a decrease in turnover frequency. 337 On the contrary, incorporation of electron donating groups, such as mesityl on the porphyrin groups resulted in high turnover frequencies but high overpotentials. 337 3.4.2 Conversion of carbon dioxide to formic acid, formamide, and formate. In the case of CO 2 reduction to formic acid, a cobaltbased, CO 2 -reducing catalyst, CpCo P R 2 N R based on the structural characteristics of CODH/FDH. 338 It contained P R 2 N R 2 0 À Á diphosphine (1,5-diaza-3,7-diphosphacyclooctane) ligands, with two pendant amine molecules acting as proton relays during the reduction of CO 2 . Four different diphosphine ligands with phenyl or cyclohexyl substituents on phosphorus (P Cy or P Ph ) and benzyl or phenyl substituents on nitrogen (N Bn or N Ph ) were prepared and were active toward the electrochemical reduction of CO 2 to formic acid, exhibiting 90 AE 10% faradaic efficiencies. 338 The cobalt catalyst with (P Cy 2 N Bn 2 ) ligand demonstrated the highest activity with a turnover frequency (TOF) value higher than 1000 s À1 . However, these enzyme-inspired molecular catalysts were not stable, as a complete loss of their activity was observed after 1 h of electrolysis. 338 Furthermore, a distinct structural characteristic of iron hydrogenases was used as the source of inspiration for the design of molecular catalysts for the hydrogenation of CO 2 to formamide and formate. The ortho-hydroxypyridine present in iron-based hydrogenases is crucial toward hydrogen splitting as it promotes bond scission. 339 Hence, the hydrogenase-inspired catalyst consisted of a manganese complex with a nitrogen ligand containing ortho-hydroxypyridin, namely 6,6 0 -dihydroxy-2,2 0 -bipyridine. 340 Manganese was chosen as the metal center due to the presence of Mn hydride intermediates derived from the interaction between Mn and hydroxy-functionalized ligand, leading to increased activity toward CO 2 reduction to HCOOH. 341,342 Turnover numbers of B6250 (in the presence of DBU (1,8diazabicyclo [5.4]undec-7-ene)) and B588 (in the presence of secondary amine diethylamine) were achieved for the hydrogenation of CO 2 to formate and formamide, respectively. 340 These initial results show the potential of these enzymeinspired electrodes as CO 2 reduction catalysts, even though further research is needed to optimize their structure and improve their durability. Thorough reviews about molecular catalysts for reduction of CO 2 are available elsewhere. [343][344][345]

Connecting nano-to macro-scale: hierarchical transport networks in electrocatalysis
In an electrochemical reaction involving a porous electrocatalyst, reactants must move from the gas or fluid stream into the electrocatalyst layer and reach the active sites on its surface. A high specific, internal surface area is desirable to achieve a high concentration of active sites per unit catalyst mass, suggesting the employment of nanoporous electrocatalysts. However, the diffusion of reactants through electrocatalytic layers is greatly influenced by the geometry and organization of their pores. In non-hierarchical structures, there is high mass transfer resistance within the nanopores, resulting in decreased overall rates. 101,111 Hence, it is crucial to minimize the effects of transport limitations to increase yield; this is not a trivial problem, due to the various phases required for transport of molecular reactants and products, ions and electrons.
A nature-inspired approach to reduce diffusion limitations is the utilization of hierarchically structured porous materials with an optimized network of narrow and broad pores. Nature can provide the mechanistic basis to discover the optimal network structure of the pores, as hierarchical transport networks are widely employed in various natural systems in plants and animals, such as their circulatory and respiratory systems (Fig. 17). In such systems, convective flow dominates transport at large length scales, while diffusion takes over at small length scales. This transition corresponds to a Péclet number of unity and is accompanied by a quite sudden change in the channel architecture: the fractal distribution of channel diameters at large scales (pressure driven) transitions into a uniform distribution of narrow channels (concentration or chemical potential driven). This optimal structural characteristic is universal and witnessed in lungs, plants, and other biological transport networks. [107][108][109][110][346][347][348][349] A numerical model was developed to demonstrate the significance of the microstructure of the cathode catalyst layer (CCL) in the performance of a PEFC (Fig. 18), because transport limitations in the CCL are responsible for a sharp drop in efficiency when operating at high current density. 351 This CCL was assumed to consist of several spherical Pt/C agglomerates, each surrounded by Nafion ionomer. The composition and size of these agglomerates (B50-5000 nm) depend on the formulation of the catalyst ink, which is sprayed onto the polymer membrane to form the MEA.
This optimized CCL design leads to high Pt utilization, even at very low Pt loadings (B0.01 mg cm À2 ), as sufficient oxygen concentration is present along the entire radius of the agglomerate (Fig. 19a). On the contrary, in a conventional design, Pt resides at the center of the agglomerate and is exposed to very low oxygen concentrations, underutilizing the catalyst (Fig. 19b). When this optimized CCL design is incorporated into a lung-inspired PEFC (described in detail in Section 5.3), the DoE target for platinum utilization of B8 kW g Pt À1 is surpassed 352 at N = 4 generations of the fractal flow field and reaches a maximum of B36 kW g Pt À1 at N = 6 generations (Fig. 19c). Hence, these modeling results demonstrate the significance of a rational design based on hierarchically structured electrocatalysts, which can significantly reduce the cost of the MEA. However, the mitigation of degradation of Pt must be considered into the design, as, at such low catalyst loadings, the electrochemically active surface area can rapidly decrease leading to significant PEFC losses.
In summary, it is evident that there is a need for a new methodology for the design of highly active and cost-efficient electrocatalysts for key, energy conversion related, electrochemical reactions, which can be successfully fulfilled by the implementation of a nature-inspired approach, using the NICE methodology. Even though challenges regarding long-term stability, cost, scale-up, and fundamental understanding of kinetic and reaction barriers in key elementary steps of electrochemical reactions still need to be addressed, the nature-inspired approach is a powerful avenue towards the rapid growth of effective electrocatalyst designs.

Bio-inspired electrochemical devices
In the following section, we provide examples of the nature-inspired approach for the design of nature-inspired electrochemical devices (Fig. 20). In the case of bio-fuel cells, the term ''bio-inspired'' is limited to the different living enzymes or microbes employed as electrodes, and is, in fact, bio-integration, as biological components are explicitly used. In contrast, in the case of other types of fuel cells and batteries, the nature-inspired approach does not require biological components; in addition, the approach aims to redesign electrochemical devices across different scales, encompassing not only the electrodes but the entire structure of the battery or the flow fields of, e.g., polymer electrolyte fuel cells (PEFCs).

Biofuel cells
Biofuel cells (BFCs) are fuel cells using a biogenous fuel or catalyst, or their combination, to convert chemical energy stored in biodegradable substances directly into electricity. 354,355 Biofuel cells are categorized into enzymatic (EFCs) and microbial fuel cells (MFCs), depending on the nature of the biocatalyst employed, and offer a clean energy alternative to fossil fuels due to the utilization of cost efficient, environmentally friendly, and renewable fuels (such as sugars and components in wastewater) to produce electrical energy. However, the main obstacle toward the commercialization of this technology is its short lifetime and low power density, 355 limiting the range of applications to powering microelectronic systems, such as sensors and actuators, [356][357][358] and for environmental remediation. 359,360 The different biocatalysts incorporated in these two types of biofuel cells leads to large differences between them in volumetric size, power output and targeted applications. 355,361 EFCs and  MFCs use redox enzymes and living microorganisms (e.g. bacteria, yeast cells), respectively, as biocatalysts for power generation. As a result, EFCs can produce high power density and are well-suited to miniaturization, sensors, and wearable or implantable devices fueled by endogenous substances (such as glucose in the blood stream). [362][363][364][365][366] On the contrary, MFCs are mainly used as bioreactors for purification of wastewater or as long-term power generators for small devices in remote locations. Despite the promising potential of biofuel cells, the commercialization of these systems is limited by their low power density and stability. 355 To circumvent these issues, research has been focused on the development of nanostructured materials as electrodes with tailored pore size (B40-300 nm), high surface area and porosity to significantly improve the electron transfer between the biocatalyst and the electrode, and the diffusion of mediators to the electrode. The most common materials used as electrodes in BFCs are carbon-and polymer-based. Carbonaceous materials (carbon nanotubes, graphene, carbon nanoparticles) 355,367-377 possess excellent mechanical and chemical stability, good conductivity and biocompatibility, with the ability to extract the electrons from the biocatalysts without being deactivated. Conducting polymers (polyaniline (PANI), polypyrrole (PPy), poly(ethylenimine) (PEI)) have high stability, conductivity, and surface area transforming them into attractive candidates as electrodes for BFCs. 355,[378][379][380][381] An example of carbon modified electrodes is magnesium oxide-templated mesoporous carbon (MgOC) with high surface area (B580 m 2 g À1 ) for efficient biocatalytic reactions in its mesopores and facile mass transport though its macropores. 382 The electrode was coated with a biocatalytic hydrogel composed of a conductive redox polymer d-FAD-GDH (deglycosylated flavin adenine dinucleotide-dependent glucose dehydrogenase), and a crosslinker. This MgOC modified electrode produced B30-fold higher current density (B3 mA cm À2 ) than a flat carbon electrode with the same hydrogel loading (B1 mg cm À2 ), while, after 220 days of testing, a B5% loss of the initial catalytic current was observed, indicating that immobilization of the enzyme in the mesopores of carbon support can significantly improve the enzyme's stability. 382 As these biofuel cells are bio-integrated, rather than bioinspired fuel cells, a thorough review of the electrodes employed in biofuel cells is beyond the scope of this manuscript; readers are referred to other reports in the literature. [355][356][357][358][383][384][385]

Bio-inspired batteries
Among the first bio-inspired approaches for the development of electrodes for Li-ion batteries was the utilization of viruses, such as M13. 386,387 A self-assembled layer of virus (M13)templated cobalt oxide nanowires serving as the active anode material was formed on top of polyelectrolyte multilayers acting as the battery electrolyte, and this assembly was stamped onto platinum microband current collectors. 386,387 The microbattery electrodes exhibited a similar charging/discharging curve as that of commercial cobalt oxide nanoparticles, even at a high charging rate of B255 nA, demonstrating the potential of this approach toward the design of new electrode materials. 386,387 Nowadays, the facile development of flexible and wearable electronics creates the need for batteries with high flexibility and energy density. In conventional battery designs these two characteristics are at odds with each other; 388,389 ultrathin batteries may be flexible, but have low energy density and vice versa, requiring the connection of numerous battery units in series to increase their energy density. [390][391][392] A nature-inspired approach to solve this issue is the design and manufacturing of a Li-ion battery based on the structure of the spine. The spine or vertebral column consists of 33 bones (vertebrae), which interlock with each other, providing the main support of the human body and its flexibility. To prevent friction between vertebrae, each vertebra is separated from the others by intervertebral discs.
Based on this structural characteristic of the human spine, flexible Li-ion batteries with high energy density were manufactured (Fig. 21a). 393 Each of the conventional anode/separator/ cathode/separator battery stack components were placed in different strips with multiple branches; then, each strip was wrapped around the backbone to form thick stacks of electrodes, corresponding to the vertebral column. The unwound part interconnecting these thick stacks of electrodes functions as the marrow, providing high flexibility to the battery (Fig. 21b). The spine-inspired design of the Li-ion battery demonstrated an B242 W h L À1 power density and B94% retention. At 0.5C, the discharge capacity was above 125 mA h g À1 under continuous dynamic load. Mechanical load tests revealed the flexibility of this design, since the largest strain on the interconnected joints was B0.08% compared to B1.8% and B1.1% for a prismatic and stacked pouch cell, respectively. 393 Another source of natural inspiration is cartilage. Cartilage is a tough, flexible connective tissue, serving as the center of ossification for bone growth and covering the surface of the joints, reducing the friction between adjacent bones and preventing damage (Fig. 22). 394,395 It comprises collagen fibers and a dense extracellular matrix (ECM) with a sparse distribution of cells called chondrocytes. 394 The ECM is composed of water, collagen, proteoglycans and noncollagen proteins, and retains the water within the cartilage, which is critical for its unique mechanical properties. Cartilage consists of several zones, namely, the superficial, middle, deep, and the calcified zone, distinguished by the distribution of the collagen and the orientation and shape of the chondrocytes. 394 The superficial zone contributes 10-20% to the cartilage thickness and protects the cartilage from shear, tensile, and compressive forces imposed by articulation. This zone is responsible for most of the tensile properties of cartilage and consists of collagen fibers tightly packed and aligned parallel to the articular surface. Below the superficial zone lies the middle zone, which contains thicker collagens than the superficial zone to resist compressive forces. 394 Right below the middle zone lies the deep zone, which constitutes the second layer of resistance to compressive forces. It contains thicker collagen fibers than in the previous layers perpendicular to the articular surface and chondrocytes oriented vertically to the underlying bone. Finally, the calcified zone secures the cartilage to the bone by attaching the collagen fibrils of the deep zone to subchondral bone. 394 To imitate the nanofiber networks of cartilage, aramid nanofibers were synthesized from Kevlar via layer-by-layer deposition and used as ion conducting membranes for Li-ion batteries, 396 solid electrolytes for Zn batteries, 397 and separators for redox flow batteries. 398 Ion conducting membranes are a key component of Li-ion batteries, providing high ionic mobility to Li + ions, stiffness and flexibility. The safety issues of Li-ion batteries are related to dendrite growth and anode expansion in the charged state; [399][400][401] the membranes are pierced during dendrite growth, which could lead to battery failure, short-circuiting and fire. 402 Aramid nanofibers (B1 mm length, B5-10 nm average diameter) were used as ion conducting membranes, due to their low Ohmic resistance, high mechanical flexibility, ionic conductivity and resistance to dendrite growth. 396 Investigation of dendrite growth was carried out under high current density (B10 mA cm À2 ); a copper (Cu) electrode was used to examine the expansion of Cu dendrites (growth zones of B50-100 nm and B25 nm tip diameter), 396 due to their smaller size than Li dendrites, making the suppression of Cu dendrites more challenging than those of Li. The theory of electrochemical dendrite growth 403 indicates that if the local mechanical properties of the ion conducting membrane are sufficient to prevent mechanical stress from Cu dendrites, it will also suppress Li dendrites, which are much softer.
Copper electrodes were examined via scanning electron microscopy (SEM) after a total charge of 0.006 mA h cm À2 , and dendrites with B500 nm size were formed on bare copper electrode. The size of the copper dendrites was significantly reduced to B100-200 nm after depositing a coating of aramid nanofibers on the electrode with a film thickness of B162 nm.  As the thickness of the coating was increased, the suppression of copper dendrite growth was enhanced as well, and no dendrite formation was observed for coatings with thickness of B809 nm. 396 Dendrite growth also has a detrimental effect on the energy density and cyclability of zinc (Zn) batteries, as Zn dendrites can easily traverse the inter-electrode space, piercing existing separators. Solid electrolytes can provide a solution to this issue at the expense of low mobility of divalent ions, leading to low energy density. 397 A cartilage-inspired composite of aramid nanofibers, poly-ethylene oxide (PEO) and zinc trifluoromethanesulfonate (Zn(CF 3 SO 3 ) 2 ) can serve as a solid electrolyte providing facile ion transport and excellent mechanical properties. 397 PEO and Zn(CF 3 SO 3 ) 2 act as the ion transport components of the solid electrolyte; the composition of PEO/Zn(CF 3 SO 3 ) 2 /aramid nanofibers was optimized with respect to Zn 2+ conductivity and mechanical properties with an optimal ratio of 9 : 3 : 1. These composites were thinner (B10 mm) than a commercial separator for Zn or Li batteries (B30-200 mm), the interfilament distances in the aramid nanofibers network were B10-20 and B2-4 times smaller than the average diameter of the stems (B1-2 mm) and growth points of Zn dendrites (B200 nm). Their Zn 2+ conductivity was B2.5 Â 10 À5 S cm À1 , 10-fold greater than the original Li batteries (B2.5 Â 10 À6 S cm À1 ). 397 This bio-inspired composite enabled Zn batteries to be rechargeable and reformable, due to the plasticity of Zn anodes and the reconfigurability of a cartilage-inspired fiber network. After 50 cycles at 0.2C, the battery retained B96% of its highest achievable capacity (B123 mA h g À1 ) and still exceeded 90% after 100 cycles at 0.2C (Fig. 23a). 397 The plasticity characteristic of this bio-inspired composite increased battery safety, since it was less prone to mechanical damage and could withstand elastic deformation from bending and plastic deformation. It also enabled shape modification to improve the ability of the battery to carry a load; various shapes of these bio-inspired batteries were tested as load bearing and charge storage elements in small drones (Fig. 23b and c), demonstrating their promising potential in the transportation industry. 397 The feasibility of using aramid nanofiber-based films as electrically conducting separators for non-aqueous redox flow batteries has also been explored.
An B8.5 mm film was obtained via layer-by-layer assembly of 20 layers of aramid nanofibers of B425 nm thickness each. The dense network of nanofibers comprising this multilayer structure reduced the pore size of the film, while the pore network for ionic transport remained intact ( Fig. 24a and b). 398 The average pore size was approximately 5 nm, much less than the pore size of a commercial separator, such as Celgard 2325 (B25 mm thick), with B390 nm pores on its surface ( Fig. 24c and d). The small pores of aramid nanofiber-based film impeded the mobility of vanadium ions, resulting in lower ionic conductivity (B0.1 mS cm À1 ) and permeability (B0.8 Â 10 À7 cm 2 s À1 ) than Celgard 2325 (B0.6 mS cm À1 and B7 Â 10 À7 cm 2 s À1 , respectively). 398 To further decrease the permeability of an aramid nanofiber-based separator without sacrificing its conductivity, its surface was functionalized with PDDA (poly(diallyldimethylammonium chloride)) and PSS (poly(styrene sulfonate)) polyelectrolytes. The addition of these charged PDDA/PSS layers on the surface of aramid nanofiber-based separator enabled the Donnan exclusion of the positively charged vanadium ions, 404 further reducing the permeability of the functionalized aramid nanofiber-based separator to B0.003 Â 10 À7 cm 2 s À1 , while its conductivity remained constant. 398 The low permeability of these functionalized separators translated into high coulombic efficiency (B95% compared to B55% of commercial Celgard) during commercial cycling (5 h duration) and high stability, exhibiting minimal degradation after 100 h of cycling. 398 Another feature in nature that has been used as inspiration for the improvement of the stability of electrodes in Li batteries is self-healing, an important survival feature of biological organisms increasing their life expectancy. The high-capacity electrode materials, such as silicon and sulfur, suffer from facile capacity fading and a short lifetime. In a commercial silicon anode, silicon particles are surrounded by a polymer binder, which binds them to the current collector to maintain electrical contact. Upon cycling, the stress generated by the volumetric changes throughout lithiation and delithiation of silicon particles fractures the particles and polymer layer, leading to a loss of electrical contact, and, hence, a loss of capacity. Thus, if damage on these electrodes can be repaired spontaneously, the cycle life of the negative electrodes of Li batteries will be significantly increased. 405 Silicon particles of the anode of a Li-ion battery were coated with a hydrogen-bonding-directed self-healing polymer, 406 which allows for cracks to heal autonomously and repeatedly. 407 A cycle life that is ten times longer than that of commercial silicon anodes and a high capacity (B3000 mA h g À1 ) were achieved, when silicon anodes were modified with a self-healing polymer. 405 However, the same self-healing strategy is not effective in positive electrodes (such as sulfur, oxygen, carbon dioxide), since they undergo conversion electrochemistry, i.e. multi-electron and drastic phase transfer (e.g. solid sulfur to soluble polysulfides) with rapid diffusion and uncontrolled deposition of intermediate polysulfides, [408][409][410] severely altering their structure during electrochemical reactions. [408][409][410][411][412] This uncontrolled phase transfer between solid materials (sulfur and lithium disulfide) and liquid polysulfides is the main contributor to the poor cycling stability and reversibility of positive electrodes.
Hence, a different self-healing approach was utilized for the positive anodes of Li-S batteries, inspired by the fibrinolysis reaction within blood vessels (Fig. 25). 413 The uncontrolled deposition and accumulation of inactive solid products in Li-S batteries is similar to the coagulation of thrombus, which obstructs the blood flow to healthy vessels. During fibrinolysis, thrombus is transformed into soluble fibrin fragments and plasmin solubilizes these thrombi fragments. An analogue to plasmin in Li-S batteries is polysulfides, the employed selfhealing agent responsible for the transfer of solid polysulfide compounds into solution, where they re-participate into electrochemical reactions. The cycling performance of sulfur particle anodes containing B0.3 M Li 2 S 5 in self-healing polysulfides was significantly extended to 7500 cycles at 1.2 mA cm À2 . The average coulombic efficiency was above 99%, exhibiting a very low decay rate of 0.01% per cycle. 413 Thus, novel healing agents that are smart, sustainable, and rapidly responsive hold future promises. This bio-inspired approach can also be easily implemented in other high-energy electrochemical storage and conversion systems, such as metal-O 2 batteries and fuel cells.

Bio-inspired fuel cells
In terms of fuel cells, nature-inspired or bio-inspired design is utilized in flow fields to circumvent the uneven reactant distribution issue across the catalyst layer, which leads to losses in fuel cell performance. The majority of the reports in the literature are based on the imitation of apparent features of biological structures (such as leaves, veins, etc.), without providing a clear theoretical foundation that allows to capture the underlying physical phenomena. [414][415][416][417][418] As a result, the design and channel geometries of these bio-imitating flow fields cannot be systematically reproduced or scaled-up, and are prone to disappointing fuel cell performance.
A true bio-inspired design was employed to improve the design of commercial flow fields in fuel cells, whose primary role is the distribution of reactants across the catalyst layer, electron transfer, as well as water and heat management. The human lung serves a similar role in nature; air is uniformly transported through its complex architecture to the bloodstream to oxygenate the blood cells. Its architecture comprises two regions: the dense upper region (bronchi) with 14-16 generations decreases the convective gas flow rate from the bronchial to the acinar airways located in the lower region (7-9 generations), which is dominated by diffusion driven transport, resulting in the production of constant entropy across each level in both regions and, hence, in minimal overall entropy production of the entire lung. 101,[107][108][109][110][111]419 Prior to the manufacturing of lung-inspired flow fields, modeling simulations were conducted to calculate the optimum number of generations required to ensure uniform distribution of reactants and minimal entropy production, or, in other words, the number of generations required for the convection driven flow to be equal to the diffusion driven flow (Péclet = 1). A detailed model was built in COMSOL, revealing that the ideal number of generations, N, for minimum overall entropy production is equal to 4-7; for N less than 4, gas flow is dominated by convection, whereas for N higher than 8, gas transport is driven by diffusion ( Fig. 26C and D). 109,352 The same PEFC operating conditions were used in the model and in the experimental measurements, employing flow fields with 10 cm 2 surface area, constant fuel cell temperature (70 1C), and three different RH values (50%, 75%, and 100%).
These modeling results served as the basis for the engineering of lung-inspired flow fields with N = 3, 4, and 5 generations via stereolithography, a 3D printing method creating 3D objects from successive layers of sintered steel. Flow fields (10 cm 2 ) with 4 generations exhibited a B30% increase in current and power density at 50 and 75% relative humidity (RH), compared to commercial single serpentine flow fields (Fig. 27A-D). 109 The positive effect of a fractal structure was also evident in the produced pressure drop, with the lung-inspired flow fields exhibiting B75% lower values than their commercial counterparts (less than 2 kPa and B5 kPa for lung-inspired and serpentine flow fields, respectively) for each RH tested, minimizing the parasitic losses and enhancing the fuel cell performance. The same pressure drop (less than 2 kPa) was measured when larger lung-inspired flow fields were constructed with 25 cm 2 surface area, whereas the high pressure drop of B25 kPa measured in serpentine flow fields (25 cm 2 ) was detrimental to their performance. 109 Lung-inspired flow field based PEFCs with N = 3 generations demonstrated the worst performance under all experimental conditions tested, due to the large spacing between adjacent outlets, resulting in insufficiently high reactant concentrations across the catalyst layer. 109 On the contrary, lung-inspired flow field based PEFCs with N = 5 generations exhibited a lower performance than commercial single serpentine flow field based PEFCs, since their narrow channels in their final generation were prone to flooding. At high RH conditions (100%), all lunginspired flow fields were vulnerable to flooding and fuel cell performance deteriorated. 109 Their susceptibility to flooding at high humidity conditions (100% RH) was evaluated via neutron radiography of a lunginspired PEFC (N = 4) during galvanostatic operation at various current densities (0.3, 0.5, and 0.6 A cm À2 ). 420 Neutron images revealed significant water accumulation in the interdigitated outlet channels of the fractal distributor, due to limited convective water removal as a result of narrow channels and slow gas flow in the lung-inspired flow fields. Flooding was mitigated at high current densities (0.5 and 0.6 A cm À2 ), since faster gas flow and high pressure drop enhanced the rate of water removal, resulting in an instantaneous increase in the potential of over 200 mV (at 0.6 A cm À2 operating current density). 420 Such a significant increase in fuel cell performance emphasizes the importance of an unobstructed inlet structure of the lung-inspired flow fields, as any defect in the fractal channel network, and especially in the outlet channels of the fractal distributor, render the inlet channels susceptible to clogging. 420 Fig. 27 The lung-inspired flow-field based PEFCs, demonstrate improved performance at 50 and 75% RH (N = 4) compared to conventional, serpentine flow-field based PEFCs (A and B); 10 cm 2 flow field area. When scaled (25 cm 2 flow field area), similar results are obtained for fractal flow fields at 50 and 75% RH (C and D), even though the performance of serpentine flow field based PEFCs is improved, due to an order of magnitude higher pressure drop than fractal flow field based PEFCs. Reproduced from ref. 109 with permission from the Royal Society of Chemistry.
It is evident that apart from a carefully crafted internal structure, the adoption of a water management strategy is required to ensure reliable operation of lung-inspired flow field based PEFCs. Recently, we developed a novel strategy 421 for water removal in PEFCs based on the utilization of capillary arrays laser drilled in the land of the flow fields, which allow supply or removal of water depending on the demand of the electrode. A parallel flow-field modified with capillaries exhibited B95% and B7% improvement in peak power density over the conventional parallel and serpentine flow-fields, respectively. 421 This effective water management strategy ensures reliable fuel cell operation and is currently incorporated into lung-inspired flow fields.
The complex 3D structure of lung-inspired flow fields, produced from stainless steel, faces manufacturing and cost challenges, as it necessitates the employment of the expensive and time consuming, laser sintering. Most recently, printed circuit boards 422 were used instead, as an alternative, cost-efficient material for the rapid manufacturing of lung-inspired flow fields, resulting in uniform distribution of reactants across the catalyst layer, increased performance compared to single serpentine flow fields and improved water management at high RH.

Conclusions and perspective
The design of electrocatalysts with optimized electronic and ionic mobility as well as kinetics at multiphase boundaries is gaining enormous interest in the context of renewable energy and more efficient electrochemical conversions. However, to improve their properties, it is imperative to control their architecture. Nanopores are characterized by a short diffusion length and high multiphase contact, which enlarges the electrochemically active surface area, while macropores enhance the transport properties and kinetics. Hence, the identification and practical use of systematic design principles that guide the chemical properties and structure of electrocatalytic materials to increase activity and stability is crucial.
This necessitates the development of multi-scale models that properly consider geometrical, physical and chemical phenomena across all scales. Nature can be an excellent guide to rational design, as it is full of hierarchical structures and biological catalysts that are scalable, efficient, and robust. However, the majority of research that aims to learn from nature to develop new electrocatalysts and devices is based on straightforward biomimicry or ''bio-imitation'', which mimics isolated features of biological or non-biological natural structures. This risks to lead to sub-optimal activity and stability, because the difference in context between nature and technology is not accounted for and the actual physical processes that govern the biological organism or system are neglected.
On the contrary, the nature-inspired engineering approach advocated in this review is based on maximizing the fundamental mechanistic understanding of the principles that underpin desired traits, and their context, followed by their appropriate incorporation into the design of new electrocatalysts and electrochemical devices. Nano-, meso-, and macro-scale levels are considered into the design, resulting in the engineering of robust, highly efficient, and scalable electrochemical devices and the synthesis of highly active electrocatalysts with increased surface area, number of active sites, and enhanced charge and mass transport.
At the nanoscale, the structure and function of metalloenzymes is the most popular source of inspiration currently, since they catalyze the same reactions as in electrochemical devices. Even though their synthesis procedure is tedious and time-consuming, the synthesized electrocatalysts exhibit high activity and stability, with enhanced transport properties. The advocated nature-inspired approach is also utilized in the design of electrochemical devices, where both the meso-and macroscale matter. The fractal network of the human lung serves as the basis for the design of lung-inspired flow fields for polymer electrolyte fuel cells, leading to uniform reactant distribution across the catalyst layer and improved fuel cell performance, compared to commercial single serpentine flow field based PEFCs. In batteries, the thorough study of the structure of the cartilage results in the creation of hierarchical, porous electrodes with high ionic mobility, mechanical flexibility, and resistance to dendrite growth. The investigation of the fibrinolysis reaction within blood vessels serves as the template for the creation of selfhealing anodes for Li-S batteries with high faradaic efficiencies and cycling stability, using polysulfides as the healing agent. Table 4 illustrates the potential of nature-inspired engineering to transform the design of electrocatalysts and electrochemical devices. It summarizes a number of examples on nature-inspired design discussed in this review, including advantages and outstanding challenges.
All these examples demonstrate the diversity of applications of nature-inspired engineering in the electrochemical domain, and its innate ability to provide innovative solutions to engineering challenges, leveraged by parallel advances in synthesis techniques, additive manufacturing, and computational tools. Nature-inspired chemical engineering (NICE) facilitates this process through a systematic methodology for design and innovation, translating nature-inspired concepts to computationally assisted designs, prototypes and implementations.
We are only at the beginning. Because nature contains examples much ahead of current technology in terms of material properties, scalability, and efficiency, with a need to satisfy multiple objectives all at once, it is worthwhile to thoroughly investigate the underlying properties to inspire innovation, and NICE offers a systematic methodology to accomplish this goal.