Gouy-Chapman model (1910) The double layer is not as compact as in Helmholtz rigid layer.

Similar documents
Chapter 2. Electrode/electrolyte interface: ----Structure and properties

V. Electrostatics. Lecture 25: Diffuse double layer structure

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

1. Mean-Field Theory. 2. Bjerrum length

Electrical double layer: revisit based on boundary conditions

Electrochemical Equilibrium Electromotive Force

1st Year Thermodynamic Lectures Dr Mark R. Wormald BIBLIOGRAPHY

Appendix II Summary of Important Equations

ˆ (0.10 m) E ( N m /C ) 36 ˆj ( j C m)

Lecture 3 Electrostatic effects November 6, 2009

Week3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity

PES 1120 Spring 2014, Spendier Lecture 6/Page 1

Review of Classical Thermodynamics

Electrostatic Potential from Transmembrane Currents

Introduction to Statistical Methods

Mass Transfer Processes

Differentiating Gaussian Processes

Chapter 18, Part 1. Fundamentals of Atmospheric Modeling

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.

Week 9 Chapter 10 Section 1-5

Physics 607 Exam 1. ( ) = 1, Γ( z +1) = zγ( z) x n e x2 dx = 1. e x2

Supplementary Notes for Chapter 9 Mixture Thermodynamics

Solution Thermodynamics

Answers Problem Set 2 Chem 314A Williamsen Spring 2000

Open Systems: Chemical Potential and Partial Molar Quantities Chemical Potential

The equation of motion of a dynamical system is given by a set of differential equations. That is (1)

Grand canonical Monte Carlo simulations of bulk electrolytes and calcium channels

Measurement of Radiation: Exposure. Purpose. Quantitative description of radiation

Frequency dependence of the permittivity

Understanding Electrophoretic Displays: Transient Current Characteristics of Dispersed Charges in a Non-Polar Medium

Physics 114 Exam 2 Spring Name:

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE

Snce h( q^; q) = hq ~ and h( p^ ; p) = hp, one can wrte ~ h hq hp = hq ~hp ~ (7) the uncertanty relaton for an arbtrary state. The states that mnmze t

NUMERICAL DIFFERENTIATION

2 Finite difference basics

Physics 114 Exam 2 Fall 2014 Solutions. Name:

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

Interconnect Modeling

Copyright 2017 by Taylor Enterprises, Inc., All Rights Reserved. Adjusted Control Limits for U Charts. Dr. Wayne A. Taylor

PHY688, Statistical Mechanics

( ) 1/ 2. ( P SO2 )( P O2 ) 1/ 2.

CHAPTER 14 GENERAL PERTURBATION THEORY

Homework 4. 1 Electromagnetic surface waves (55 pts.) Nano Optics, Fall Semester 2015 Photonics Laboratory, ETH Zürich

Physics 181. Particle Systems

Inductance Calculation for Conductors of Arbitrary Shape

Fundamental Considerations of Fuel Cells for Mobility Applications

STATISTICAL MECHANICS

Canonical transformations

COMPUTER AIDED ANALYSIS OF THE MOS CAPACITOR IN LOW FREQUENCY CONDITIONS

Problem Points Score Total 100

2016 Wiley. Study Session 2: Ethical and Professional Standards Application

Chapter 13: Multiple Regression

Linear Approximation with Regularization and Moving Least Squares

Mass transfer in multi-component mixtures

Module 3: Element Properties Lecture 1: Natural Coordinates

Introduction to Vapor/Liquid Equilibrium, part 2. Raoult s Law:

Supplementary Information. Fundamental Growth Principles of Colloidal Metal Nanoparticles

Gravitational Acceleration: A case of constant acceleration (approx. 2 hr.) (6/7/11)

4.2 Chemical Driving Force

Numerical Investigation of Electroosmotic Flow. in Convergent/Divergent Micronozzle

Kernel Methods and SVMs Extension

MAGNETISM MAGNETIC DIPOLES

Econ107 Applied Econometrics Topic 3: Classical Model (Studenmund, Chapter 4)

and Statistical Mechanics Material Properties

Lecture 7: Boltzmann distribution & Thermodynamics of mixing

1. Mathematical models of the chromatographic process

This chapter illustrates the idea that all properties of the homogeneous electron gas (HEG) can be calculated from electron density.

Lecture 8 Modal Analysis

ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM

MEASUREMENT OF MOMENT OF INERTIA

Complex Numbers. x = B B 2 4AC 2A. or x = x = 2 ± 4 4 (1) (5) 2 (1)

Chapter 2: Electric Energy and Capacitance

PART I: MULTIPLE CHOICE (32 questions, each multiple choice question has a 2-point value, 64 points total).

Global Sensitivity. Tuesday 20 th February, 2018

Numerical Heat and Mass Transfer

Temperature. Chapter Heat Engine

Solutions of execrises

Psychology 282 Lecture #24 Outline Regression Diagnostics: Outliers

Thermal-Fluids I. Chapter 18 Transient heat conduction. Dr. Primal Fernando Ph: (850)

Lecture Notes on Linear Regression

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

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

Gasometric Determination of NaHCO 3 in a Mixture

II. Equilibrium Thermodynamics. Lecture 8: The Nernst Equation

Math1110 (Spring 2009) Prelim 3 - Solutions

Comparative Studies of Law of Conservation of Energy. and Law Clusters of Conservation of Generalized Energy

Relaxation in water /spin ice models

8.022 (E&M) Lecture 4

Chapter 8. Potential Energy and Conservation of Energy

V T for n & P = constant

Important Instructions to the Examiners:

Correlation and Regression. Correlation 9.1. Correlation. Chapter 9

Joint Statistical Meetings - Biopharmaceutical Section

π e ax2 dx = x 2 e ax2 dx or x 3 e ax2 dx = 1 x 4 e ax2 dx = 3 π 8a 5/2 (a) We are considering the Maxwell velocity distribution function: 2πτ/m

Kinematics of Fluids. Lecture 16. (Refer the text book CONTINUUM MECHANICS by GEORGE E. MASE, Schaum s Outlines) 17/02/2017

Osmotic pressure and protein binding

STATISTICAL MECHANICAL ENSEMBLES 1 MICROSCOPIC AND MACROSCOPIC VARIABLES PHASE SPACE ENSEMBLES. CHE 524 A. Panagiotopoulos 1

PHYS 705: Classical Mechanics. Newtonian Mechanics

STAT 511 FINAL EXAM NAME Spring 2001

3. Be able to derive the chemical equilibrium constants from statistical mechanics.

Transcription:

CHE465/865, 6-3, Lecture 1, 7 nd Sep., 6 Gouy-Chapman model (191) The double layer s not as compact as n Helmholtz rgd layer. Consder thermal motons of ons: Tendency to ncrease the entropy and make the layer dffuse (whle electrostatc attracton tends to keep ons close to nterface and make the layer more rgd). etal Soluton Dstrbuton of ons: determned by electrostatcs and stat. mechancs Assumptons: Ions: pont ons, hghly moble Solvent: delectrc contnuum etal: perfect conductor σ x Dstrbuton of ons wth dstance x from nterface: apply Boltzmann statstcs zf n = n exp n : partcle densty (#partcles/cm3) Regons of large ϕ ( x) S ( ϕ ( x) ϕ ) RT : reduced densty of postve ons, enhanced densty of anons (compared to bulk).

Charge densty nvolvng all onc speces: ( ) ρ x = n z e wth e = 1.6 1 zf = n z e exp S ( ϕ ( x) ϕ ) RT Galvan potental n soluton obeys Posson s equaton: d ϕ dx ( x) ρ ( x) = Resultng expressons s the Posson-Boltzmann equaton: d ϕ dx ( x) 1 z Fϕ = n z q x = x εε εε S exp, wher e ϕ ( ) ϕ ( ) ϕ RT Note: Ths equaton s precsely equal to the Debye-Hückel theory of onc nteractons n dlute solutons, whch determnes the charge dstrbuton around a central on. -19 C Next: Let s consder a z-z-electrolyte (*) ρ ( x) ( ) ϕ ( ) zfϕ x zf x = n zq exp exp RT RT

Debye-Hückel approxmaton Let s further assume that the potental varaton s small, zf ϕ RT ( x) << 1.e. potental varatons are ϕ << 5 mv at room-t. Ths wellknown approxmaton corresponds to the so-called Debye-Hückel approxmaton. If we use ths approxmaton n the charge dstrbuton (*) of the z-z electrolyte and then nsert ths charge dstrbuton nto the Posson-Boltzmann equaton, then we arrve at the lnear Posson-Boltzmann equaton: d ϕ dx = κ ϕ Here, κ s the nverse Debye-length, ( zf ) 1/ c 1 κ = = εε RT L D, c : bulk electrolyte concentraton (far away from nterface) κ or ts nverse L D are mportant characterstcs of the electrolyte. Debye-length L D : quanttatve measure of wdth of space charge regon wthn whch ϕ ( x) decreases from ϕ to tremendous mportance n electrochemstry and bology. S ϕ

The followng effects are mportant to remember: Larger electrolyte concentraton c smaller.a.w: the double layer becomes less dffuse Hgher temperature larger L D.a.w.: the double layer becomes more dffuse L D Potental and charge dstrbuton n the electrolyte soluton are gven n Debye-Hückel approxmaton by: σ ϕ = εε κ ( x) exp ( κ x) ( x) = exp ( x) ρ σ κ κ.e. they are exponental functons of the dstance from the nterface,.e. the excess surface charge densty σ on the metal s balanced by an exponentally decayng space-charge layer n soluton. The double layer capacty n ths approxmaton s gven by C d,dh εε = εε κ = L.e. t s gven by the plate capactor formula wth the Debye length as the effectve plate separaton. D

Table: Debye length at varous electrolyte concentratons (1-1 electrolyte) c / mol l -1 1-4 1-3 1-1 -1 L D / Å 34 96 3.4 9.6 Hgher concentratons: steeper potental drop n soluton narrower space charge regon, smaller L D hgher double layer capactes, C d,dh In practce: The Debye-Hückel (DH) )approxmaton s not vald at large electrolyte concentraton. It works well as long as the potental varaton does not exceed ϕ 8 mv.

General Case: Nonlnear Posson-Boltzmann-equaton For a z-z-electrolyte wth concentraton c the Posson- Boltzmann equaton has an explct soluton. Ths wll be added as an appendx. Here only the man results wll be gven. Potental dstrbuton n soluton etal Soluton ϕ (nterface) ϕ, = ϕ - ϕ S ϕ S (bulk) σ x The relaton between potental and charge densty of the dffuse layer s gven by dϕ 1/ zfϕ, = = ( 8RT c ) snh RT σ εε εε dx x= Dfferentatng ths expresson gves the dfferental capacty n the Gouy-Chapman model: σ εε zfϕ, d,gc = = cosh, where ϕ, = pzc ϕ, LD RT C E E

The varaton of ths capactance wth electrode potental E s thus gven by C d,gc Ths model predcts correctly the mnmum n capacty at the pzc. For large E, t predcts an unlmted rse of the capacty. E-E pzc Why would the DL capacty rse unlmtedly n ths approach? Whch assumpton s responsble for ths unphyscal behavour? It s assumed that ons are pont charges, whch could approach the electrode surface arbtrarly close wth ncreasng E. Ths leads to a very small charge separaton (consder the plate capactor as an analogue!). The capacty rses wthout lmt. What s the soluton out of ths capacty catastrophe??? Ions have a fnte sze! Arbtrarly close approach s not possble. Ions are stopped from approachng the electrode at dstances that correspond to ther rad. A further refnement of the double layer models takes ths effect nto account.

The Stern odel (194) accounts for fnte sze of ons combnes the Helmholtz and Gouy- Chapmann models two parts of double layer: (a) compact layer ( rgd layer ) of ons at dstance of closes approach (b) dffuse layer. The compact layer, x < x H, s chargefree (lnear varaton of potental)! etal σ compact layer x H dffuse layer Soluton athematcally: two capactors n seres wth total capacty 1 1 1 C C C x εε H = + = + L εε cosh z d dff,h dff,gc, Far from E pzc (large ϕ, = E Epzc ): D Fϕ RT Cdff,H << Cdff,GC Cd Cdff,H Helmholtz, rgd Close to E pzc (small ϕ, = E Epzc ): Cdff,H >> Cdff,GC Cd Cdff,GC Gouy-Chapman, dffuse Remember: the smaller guy always wns (well, not always, but here)!

Effect of electrolyte concentraton: smaller electrolyte concentraton smaller κ (L D larger), double s more dffuse! C dff,gc becomes more mportant The Stern model reproduces gross features of real systems. Potental dstrbuton n nterfacal layer and capacty varaton wth E are shown n the fgures below. What happens upon ncreasng the electrolyte concentraton? Potental dstrbuton Double layer capactance charge free moble ons

Further refnement: Grahame model (1947) Some ons (usually anons) loose hydraton shell smaller rad! Ions of smaller rad could approach the electrode closer. Dstngush three dfferent regons! Inner Helmholtz plane (IHP): through centers of small, partally solvated ons Outer Helmholtz plane (OHP): through centers of fully solvated ons Outsde OHP: dffuse layer

Ths s the best model so far. It shows good correspondence to the expermental data. Why are the potental dstrbutons (n (b)) dfferent for postve and negatve E? Ths s due to a dfference between anons and catons. Anons have a less rgd solvaton shell. They become more easly desolvated. Thus, desolvated anons could form the IHP. Ths happens at postve E. Catons keep ther solvaton shell. They, thus, cannot approach the electrode closer than to the OHP.

Appendx: Electrocapllary measurements Only applcable to lqud electrodes (metals), based on measurement of surface tenson (Lppmann s method). The method was developed specfcally for mercury electrodes. Ths s a null-pont technque,.e. the measurement of a physcal property s performed when the system has reached an equlbrum state. Such measurements are rather accurate. Balanced forces: Surface tenson of mercury n the capllary counterbalances the force of gravty π r γ cosθ = π r ρ hg c c Hg

r c s the capllary radus, γ s the surface tenson, θ s the contact angle h s the heght of the capllary column of Hg. In the measurement, the contact angle θ s measured wth a mcroscope. Above relaton then provdes values of the surface tenson γ. As the expermental parameter, the electrode potental E can be vared. What happens upon changng E? When the system s balanced (reachng a new equlbrum mnmum of Gbbs free energy!), the change n surface energy s exactly balanced by the electrcal work, A γ = QdE s d where A s s the surface area and Q s the excess charge on the surface. Therefore, the charge densty on the Hg electrode s gven by Q dγ σ = =. As de It s, thus, possble to measure the surface charge as a functon of E.