Supporting Information

Similar documents
1) Silicon oxide has a typical surface potential in an aqueous medium of ϕ,0

Towards the modeling of microgalvanic coupling in aluminum alloys : the choice of boundary conditions

Calculation Of The Interdiffusion Coefficient In The Cu Zn Diffusion Couple

2010 Black Engineering Building, Department of Mechanical Engineering. Iowa State University, Ames, IA, 50011

Supplementary Information for Observation of Parity-Time Symmetry in. Optically Induced Atomic Lattices

A Modulated Hydrothermal (MHT) Approach for the Facile. Synthesis of UiO-66-Type MOFs

Supporting Information

If two volatile and miscible liquids are combined to form a solution, Raoult s law is not obeyed. Use the experimental data in Table 9.

Assignment 4. Adsorption Isotherms

(Online First)A Lattice Boltzmann Scheme for Diffusion Equation in Spherical Coordinate

Problem adapted reduced models based on Reaction-Diffusion Manifolds (REDIMs)

Diffusion Mass Transfer

Amplification and Relaxation of Electron Spin Polarization in Semiconductor Devices

Three-Phase Distillation in Packed Towers: Short-Cut Modelling and Parameter Tuning

THE EFFECT OF TORSIONAL RIGIDITY BETWEEN ELEMENTS ON FREE VIBRATIONS OF A TELESCOPIC HYDRAULIC CYLINDER SUBJECTED TO EULER S LOAD

Week 9 Chapter 10 Section 1-5

between standard Gibbs free energies of formation for products and reactants, ΔG! R = ν i ΔG f,i, we

STUDY ON TWO PHASE FLOW IN MICRO CHANNEL BASED ON EXPERI- MENTS AND NUMERICAL EXAMINATIONS

Robert Eisberg Second edition CH 09 Multielectron atoms ground states and x-ray excitations

AP Physics 1 & 2 Summer Assignment

NAME and Section No. it is found that 0.6 mol of O

COMPLEX NUMBERS AND QUADRATIC EQUATIONS

EVALUATION OF THE VISCO-ELASTIC PROPERTIES IN ASPHALT RUBBER AND CONVENTIONAL MIXES

COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD

Problem Points Score Total 100

Appendix B. The Finite Difference Scheme

Be true to your work, your word, and your friend.

Electrical double layer: revisit based on boundary conditions

Airflow and Contaminant Simulation with CONTAM

modeling of equilibrium and dynamic multi-component adsorption in a two-layered fixed bed for purification of hydrogen from methane reforming products

Grand canonical Monte Carlo simulations of bulk electrolytes and calcium channels

Solutions for Euler and Navier-Stokes Equations in Powers of Time

Open Systems: Chemical Potential and Partial Molar Quantities Chemical Potential

THE IGNITION PARAMETER - A quantification of the probability of ignition

Mass Transfer Processes

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems

A particle in a state of uniform motion remain in that state of motion unless acted upon by external force.

V. Electrostatics. Lecture 25: Diffuse double layer structure

An identification algorithm of model kinetic parameters of the interfacial layer growth in fiber composites

EN40: Dynamics and Vibrations. Homework 4: Work, Energy and Linear Momentum Due Friday March 1 st

11. Dynamics in Rotating Frames of Reference

6.3.4 Modified Euler s method of integration

Spin-rotation coupling of the angularly accelerated rigid body

Tensor Smooth Length for SPH Modelling of High Speed Impact

Supplementary Information:

10.34 Numerical Methods Applied to Chemical Engineering Fall Homework #3: Systems of Nonlinear Equations and Optimization

Modeling of Dynamic Systems

B and H sensors for 3-D magnetic property testing

Design Equations. ν ij r i V R. ν ij r i. Q n components. = Q f c jf Qc j + Continuous Stirred Tank Reactor (steady-state and constant phase)

Solution Thermodynamics

Electrochemistry Thermodynamics

Neutral-Current Neutrino-Nucleus Inelastic Reactions for Core Collapse Supernovae

Number Average Molar Mass. Mass Average Molar Mass. Z-Average Molar Mass

The internal structure of natural numbers and one method for the definition of large prime numbers

CHAPTER II THEORETICAL BACKGROUND

THE VIBRATIONS OF MOLECULES II THE CARBON DIOXIDE MOLECULE Student Instructions

Supplemental Material: Causal Entropic Forces

Transfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system

Lab 2e Thermal System Response and Effective Heat Transfer Coefficient

Inductance Calculation for Conductors of Arbitrary Shape

Supplementary Notes for Chapter 9 Mixture Thermodynamics

THE MIND & THE SCIENCE OF SUCCESS - All rights reserved - - APPENDIX -

Irregular vibrations in multi-mass discrete-continuous systems torsionally deformed

25 Reprogramming of metabolic systems by altering gene expression

DETERMINATION OF TEMPERATURE DISTRIBUTION FOR ANNULAR FINS WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY BY HPM

DEGREE REDUCTION OF BÉZIER CURVES USING CONSTRAINED CHEBYSHEV POLYNOMIALS OF THE SECOND KIND

CHEMICAL ENGINEERING

Title Natural Convection in CuSO4-H2SO4 S JOURNAL OF THE ELECTROCHEMICAL SOCI (2009), 156(9): F99-F108.

First Law: A body at rest remains at rest, a body in motion continues to move at constant velocity, unless acted upon by an external force.

,..., k N. , k 2. ,..., k i. The derivative with respect to temperature T is calculated by using the chain rule: & ( (5) dj j dt = "J j. k i.

Gasometric Determination of NaHCO 3 in a Mixture

KINETICS OF GAS HYDRATE FORMATION FROM PYROLYSIS GAS IN WATER-IN-OIL EMULSION SYSTEM

KINETICS OF PHOTOCHEMICAL DECOMPOSITION OF SOLIDS (I) Introduction. The Rate Equation

Energy, Entropy, and Availability Balances Phase Equilibria. Nonideal Thermodynamic Property Models. Selecting an Appropriate Model

Process Modeling. Improving or understanding chemical process operation is a major objective for developing a dynamic process model

Using the estimated penetrances to determine the range of the underlying genetic model in casecontrol

CHAPTER 13. Exercises. E13.1 The emitter current is given by the Shockley equation:

Normally, in one phase reservoir simulation we would deal with one of the following fluid systems:

Study on Active Micro-vibration Isolation System with Linear Motor Actuator. Gong-yu PAN, Wen-yan GU and Dong LI

CHEMICAL REACTIONS AND DIFFUSION

Mathematical Modeling to Support Gamma Radiation Angular Distribution Measurements

Perfect Fluid Cosmological Model in the Frame Work Lyra s Manifold

Army Ants Tunneling for Classical Simulations

Stability of Tip in Adhesion Process on Atomic Force Microscopy Studied by Coupling Computational Model

NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS

Publication 2006/01. Transport Equations in Incompressible. Lars Davidson

Chapter 8. Potential Energy and Conservation of Energy

This model contains two bonds per unit cell (one along the x-direction and the other along y). So we can rewrite the Hamiltonian as:

An Algorithm to Solve the Inverse Kinematics Problem of a Robotic Manipulator Based on Rotation Vectors

and Statistical Mechanics Material Properties

Lecture 12. Modeling of Turbulent Combustion

On the correction of the h-index for career length

Wilbur and Ague 4 WILBUR AND AGUE; APPENDIX DR1. Two-dimensional chemical maps as well as chemical profiles were done at 15 kv using

Relationship between fiber degradation and residence time distribution in the processing of long fiber reinforced thermoplastics

Electron-Impact Double Ionization of the H 2

On Higher Order Dynamics in Lattice-based Models Using Chapman-Enskog Method

( ) = ( ) + ( 0) ) ( )

Fundamental Considerations of Fuel Cells for Mobility Applications

The Feynman path integral

Comparison of the Population Variance Estimators. of 2-Parameter Exponential Distribution Based on. Multiple Criteria Decision Making Method

Transcription:

Supportng Informaton for Scalable, hgh performance, enzymatc cathodes based on nanomprnt lthography Dmtry Pankratov,2, Rchard Sundberg 3, Javer Sotres, Dmtry B. Suyatn 3,4, Ivan Maxmov 3, Sergey Shleev,2,5*, and Lars Montelus 3,4 Address: Bomedcal Scences, Health and Socety, Malmö Unversty, 20560 Malmö, Sweden, 2 A.N. Bach Insttute of Bochemstry, 907 Moscow, Russa, 3 Dvson of Sold State Physcs and The Nanometer Structure Consortum (nmc@lu), Lund Unversty, 2200 Lund, Sweden, 4 Neuronano Research Center, Medcal Faculty, Lund Unversty, 22 00 Lund, Sweden and 5 Kurchatov NBICS Centre, Natonal Research Centre "Kurchatov Insttute", 2382 Moscow, Russa Emal: Sergey Shleev - sergey.shleev@mah.se * Correspondng author S

Supplementary data Fgure S: Schematc llustraton of the fabrcaton process of the nanomprnted plastc electrodes. A) N stamp wth ca. 00 nm features s used to transfer the pattern nto the thermoplastc fol, B) nanomprnt patternng s performed by drect mechancal contact of the stamp and the fol at T equal to 60 C and pressure of 50 bar, C) after separaton of the stamp and the polymer, the stamp nano-pattern s replcated nto the polymer, D) 00 nm thck Au layer s thermally evaporated on the mprnted polymer. S2

Fgure S2: Typcal CVs of Au (black curve) and NIL/Au (red curve) electrodes. Condtons: 0.5 M H 2 SO 4 ; 00 mv s scan rate; second cycle. S3

Addtonal AFM studes Fgure S3: Atomc force mcroscopy mages and correspondng heght profles of nanostructured Au electrodes. The mages show 2D and 3D representatons of a representatve area of the samples. Scan area: 500 nm 500 nm. Colour heght scale: 0 nm 245.7 nm. The heght profles (postons hghlghted n red n the topographcal mages) show the depth of the nano-cavtes along the two dfferent drectons defned by the prmtve vectors of the hexagonal lattce defned by the nanocavtes. S4

Addtonal boelectrochemcal studes Fgure S4: CVs of BOx-modfed Au electrodes n ar-saturated (black curve) and O 2 -saturated (blue curve) buffers. Condtons: PBS; scan rate: 20 mv s ; second cycle. Fgure S5: CVs of BOx-modfed NIL/Au electrodes n ar-saturated (black curve) and O 2 -saturated (blue curve) buffer. Condtons: PBS; scan rate: 20 mv s ; second cycle. S5

Fgure S6: CVs of Au electrode (black curve) and Au electrodes modfed wth enzyme solutons of dfferent concentratons, μg ml : 0.4 (blue curve), 4 (green curve), 40 (navy curve), and 400 (red curve). Condtons: O 2 -saturated PBS; 20 mv s scan rate; second cycle. Fgure S7: CVs of a non-bomodfed NIL/Au electrode (black curve) and NIL/Au electrodes modfed wth enzyme solutons of dfferent concentratons, μg ml : 0.4 (blue curve), 4 (green curve), 40 (navy curve), and 400 (red curve). Condtons: O 2 -saturated PBS; 20 mv s scan rate; second cycle. S6

Ellpsometry studes The tme evoluton of the adsorbed amount, thckness, and refractve ndex of BOx layer formed on a planar Au surface were montored by means of null-ellpsometry (Fgure S8). For ths, the BOx layer was ntally formed by mmersng the Au electrodes n a BOx buffer soluton at a concentraton of 40 μg ml for 60 mn (adsorpton was also montored by means of ellpsometry, correspondng to the frst 60 mnutes of the plots n Fgure S8). Then, the ambent soluton was replaced wth proten-free buffer and the layer montored for 5 more than hours. Ellpsometry data clearly showed that the adsorbed amount of BOx on the Au surface was approxmately constant wth a value of ca. 2.6 mg m 2 2 geom (4.4 pmol cm geom or 2.6 pmol cm 2 real ) durng the whole process. a Fgure S8: a) Adsorbed amount, b) Thckness and c) Refractve ndex correspondng to the formaton (ntal 60 mn), rnsng (followng 5 mn) and stablzaton (followng 5 hours) of a BOx layer on a planar Au surface montored by null-ellpsometry. b c S7

Theoretcal bass of boelectrochemcal nvestgatons and modellng For proper mathematcal elaboraton of boelectrochemcal data, where a boelectrocatalytc process has a mxed knetcs regme, possble dffuson lmtatons should be excluded. In general, the rate of a boelectrocatalytc process can be descrbed by a mxed knetcs equaton (Equaton S []). ET cat ss (S) where, s the observed current, ET s the lmtng current of the heterogeneous electron transfer (ET; Equaton S2; vde nfra), cat s the lmtng current of the bocatalytc process (Equaton S3), and s the lmtng dffuson current (Equaton S4 s s at close to nfnty). ET nfa real n' F( E E k Г 0 exp RT 0 ') (S2) n ' s the number of electrons n the slow electrochemcal step, F s the Faraday constant, k 0 s the standard heterogeneous ET constant, Г s the surface concentraton of the enzyme, α s the charge transfer coeffcent, E s the electrode potental, E 0 s the equlbrum potental of the electrode process, R s the gas constant, and Т s the temperature n K. If ntramolecular ET (IET) s not a lmtng step n the enzymatc process of O 2 boelectroreducton (our prevous studes of BOx showed that ths was the case at ph 7.4 [2]), the boelectrocatalytc current ( cat ) can be expressed as the electrochemcal form of the Mchaels Menten equaton: S8

cat app nfarealkcat С cat oxygen (S3) C K oxygen M where app k cat s the apparent rate constant for the boelectrocatalytc process and K M s the Mchaels constant. The Levch equaton (Equaton S4) defnng steady-state current lmted by the transport of substrate molecules at a rotatng electrode can be used to estmate mass-transfer lmtatons: 2 / 3 / 6.62nFD C A ss geom 0 oxygen / 2 (S4) here, s the steady-state dffuson current, n s the total number of electrons, D s the s s dffuson coeffcent of the substrate (9.7 0 6 cm 2 s for O 2 ), C s the bulk oxygen concentraton of oxygen, O 2 (2.5 0 7 mol cm 3 and.2 0 6 mol cm 3 of O 2 (ar and O 2 saturated solutons, respectvely [3]), s the knematc vscosty of the soluton (0.0 cm 2 s at 25 C, a typcal value for aqueous solutons), and s the angular frequency (n rad s ). For a boelectrocatalytc process strctly lmted by O 2 dffuson, drect dependence between С oxygen and j max should be obtaned,.e., j max should decrease by a factor of ca. 5 when the O 2 concentraton was decreased from.2 mm to 0.25 mm. For electrocatalytc currents, whch are lmted by reacton knetcs, j max can be smply expressed as the electrochemcal form of the Mchaels Menten equaton (Equaton S3). When the O 2 concentraton was decreased from.2 mm down to 0.25 mm n our studes,.e., when ar-saturated buffer was used nstead of oxygen, the maxmal current S9

densty (j max ), whch corresponds to boelectrocatalytc O 2 reducton, decreased by a factor of ca. 3, whch follows Equaton S3, suggestng just a mnor O 2 dffuson lmtaton. In the present work, recorded CVs were analysed usng the knetc scheme recently elaborated by Clment et al. [4]. The equatons defnng the currents are as follows: DET FA k ( P ) k )] (S5) [ 2P F 0' k k0 exp E ET RT (S6) F 0' k2 k0 exp E ET RT (S7) cat O 2 O kcat nfa (S8) K M 2 where k and k 2 are potental dependent DET rate constants defned by the Butler Volmer formalsm (Equatons S2, S6, and S7), whereas P represents the fracton of adsorbed BOx molecules wth the reduced T copper centre. The maxmal boelectrocatalytc current of O 2 electroreducton, as well as current dependence on the appled potental, can be descrbed by the followng summarsed equaton: j F exp RT j max 0' k5 F 0' E E exp E E T T k0 RT (S) Usng ths equaton, modellng studes of obtaned boelectrocatalytc sgnals were performed (Fgure S9) and basc parameters of boelectrocatalytc reducton of O 2 were calculated (Table S), assumng K M and Г are equal to 0.2 mm and 2.6 pmol cm 2 real, S0

respectvely. For the calculatons, the K M value was taken from our prevous studes of the enzyme mmoblsed on bare polycrystallne Au surfaces [5], whereas the Г value was obtaned from the ellpsometry studes (vde supra). Fgure S9: Expermental versus modelled voltammograms of Au (trangles) and NIL/Au (crcles) bocathodes. (ponts expermental data, lnes modelled curves). Condtons: O 2 -saturated PBS; 20 mv s scan rate; second cycle; an enzyme soluton wth a concentraton equal to 40 μg ml was used for the bomodfcaton. S

Table S: Calculated bo-electrocatalytc parameters of Au and NIL/Au boelectrodes. Parameter Au NIL/Au Maxmal bo-electrocatalytc current densty 8 ± 3 58 ± 6 (j cat ; A cm 2 ; ar-saturated buffer) Maxmal bo-electrocatalytc current densty 48 ± 5 86 ± (j cat ; A cm 2 ; O 2 -saturated buffer) k app cat (s ) 30 ± 4 39 ± 5 k 0 (s ) 27 ± 3 27 ± 3 S2

References. Sucheta, A.; Cammack, R.; Wener, J.; Armstrong, F. A. Bochemstry 993, 32, 5455 5465. do:0.02/b0007a023 2. Shleev, S.; Andoralov, V.; Falk, M.; Remann, C. T.; Ruzgas, T.; Srnec, M.; Ryde, U.; Rulíšek, L. Electroanalyss 202, 24, 524 540. do:0.002/elan.2020088 3. Truesdale, G. A.; Downng, A. L. Nature 954, 73, 236. do:0.038/73236a0 4. Clment, V.; Zhang, J.; Frs, E. P.; Østergaard, L. H.; Ulstrup, J. J. Phys. Chem. C 202, 6, 232 243. do:0.02/jp2086285 5. Pankratov, D.; Sotres, J.; Barrantes, A.; Arnebrant, T.; Shleev, S. Langmur 204, 30, 2943 295. do:0.02/la402432q S3