Open Collections

UBC Theses and Dissertations

UBC Theses Logo

UBC Theses and Dissertations

NMR study of PF₆̄ doped polypyrrole : a potential candidate for artificial muscles Tso, Chien-Hsin 2006

Your browser doesn't seem to have a PDF viewer, please download the PDF to view this item.

Notice for Google Chrome users:
If you are having trouble viewing or searching the PDF with Google Chrome, please download it here instead.

Item Metadata


831-ubc_2006-0685.pdf [ 4.1MB ]
JSON: 831-1.0085182.json
JSON-LD: 831-1.0085182-ld.json
RDF/XML (Pretty): 831-1.0085182-rdf.xml
RDF/JSON: 831-1.0085182-rdf.json
Turtle: 831-1.0085182-turtle.txt
N-Triples: 831-1.0085182-rdf-ntriples.txt
Original Record: 831-1.0085182-source.json
Full Text

Full Text

N M R Study of P F ~ Doped Polypyrrole: A Potential Candidate for Art i f ic ial Muscles by Chien-Hsin Tso B.Sc, Tsing Hua University, Taiwan, 2001 A THESIS SUBMITTED IN PARTIAL F U L F I L M E N T OF T H E REQUIREMENTS FOR T H E D E G R E E OF Master of Science in The Faculty of Graduate Studies (Physics) The University Of British Columbia August 2, 2006 © Chien-Hsin Tso2006 11 Abstract Conducting polymer actuators, such as polypyrrole, are promising materials for novel applications because they are more muscle-like and potentially less expensive than com-mon motors. A variety of N M R techniques are employed to examine the ion content and dynamics of PFg doped polypyrrole films in different oxidation states. The ion content decreases linearly with decreasing electrochemical potential, directly confirming the ion insertion mechanism of polymer actuation. With known ion content and deposition cur-rent, a doping level (dopant ion/pyrrole) of 0.26 was determined for the as-grown film. A Ti relaxation study reveals that the rotational correlation time of PFg ions in the oxidized film (10.76 ps) is similar to that in the solvent (8.08 ps), suggesting the ions are located inside solvent pockets rather than at stable sites in the polymer matrix. ID Nuclear Overhauser Effect (NOE) difference experiments confirm the solvated ion envi-ronment in oxidized films, but reveal a drastic decrease in the NOE enhancement factor in reduced films, implying that polypyrrole undergoes a significant structural change when reduced. This change leads to a much less solvated environment as experienced by the dopant ions in the reduced state. Translational motion of the PFg ions in the oxidized films at two orientations is probed via self-diffusion coefficient measurements made using pulsed-field gradient (PFG) NMR. The D values obtained at different diffusion times range from 3 x 1 0 - 8 cm 2/s to 5 x 10~9 cm 2/s. The echo attenuation obtained by P F G . N M R experiments not only gives direct information on ionic diffusion but also serves "as a potential tool to explore the pore morphology of doped polypyrrole films. I l l Contents Abstract i i Contents iii List of Tables . . v List of Figures vi Acknowledgements viii 1 Introduction 1 1.1 Introductory Remarks . 1 1.2 An Introduction to Nuclear Magnetic Resonance (NMR) 2 1.3 Principles of NMR . . . • 2 1.4 Effects of Pulses 6 1.5 Ti, T 2 Relaxation Times '. . 7 1.6 N M R Spectroscopy: FID Signal and the Fourier Transform 8 1.7 Internal Spin interactions 9 1.7.1 Chemical shift < 9 1.7.2 J-couplings 10 1.7.3 Dipole-dipole couplings 11 2 A n Introduction to PFg Doped Polypyrrole 12 2.1 PFg Doped Polypyrrole Belongs to the Class of Conducting Polymer . . 12 2.2 PFg Doped Polypyrrole as Muscle-Like Actuators 14 2.3 Comparison between Conducting Polymer Actuators and Other Actuation Technologies 15 2.4 Synthesis of PFg Doped Polypyrrole 16 2.4.1 Electrodeposition 16 2.4.2 Electrochemical Reduction 18 2.4.3 Conductivity Measurement 18 3 PFg Ion Content in the PPy Films as a Function of Oxidation States 20 3.1 Experimental Method . . : 20 3.1.1 ID 3 1 P Spectroscopy 20 Contents iv 3.2 Results 20 3.3 Discussion 22 4 Rotational Motions of PFg Ions in the P P y Fi lms 24 4.1 Relaxation Theory 24 4.2 Experimental Method 28 4.2.1 Relaxation Measurements 28 4.3 Results 29 4.4 Discussion 31 5 Solvent Accessibility of PFg Ions in the P P y Fi lms 33 5.1 Nuclear Overhauser Effect (NOE) 33 5.2 Experimental Method 34 5.2.1 ID NOE Difference Spectroscopy . 34 5.3 Results 34 5.4 Discussion 35 6 Translational Motions of PFg Ions in the P P y Fi lms 37 6.1 Pulsed-Field Gradient (PFG) NMR 37 6.1.1 Self-Diffusion in Restricted Geometries . 40 6.2 Experimental Method , 42 6.2.1 Pulsed-Field Gradient N M R Measurement 42 6.3 Results 48 6.4 Discussion 49 7 Conclusion 53 References 55 V List of Tables 1.1 Properties of most common nuclei for N M R studies 6 3.1 Doping levels as a function of oxidation state 23 4.1 Ti relaxation times and the rotational correlation times r c of two solvent samples ~ 30 4.2 3 1 P Ti relaxation times as a function of temperature of PPy sample Si~S3. 31 4.3 T i relaxation times and the rotational correlation times T c of two PPy samples: St and S2 at room temperature 31 4.4 T 2 relaxation times of three PPy samples: Si~S3 31 5.1 NOE enhancement factor rj of the five samples 35 6.1 Spectra width Au from one-D image gradient calibration at gradient am-plitude x (0-1) and the corresponding gradient strength. 46 6.2 Self-diffusion coefficients of solvent samples 48 6.3 Diffusion coefficients of the oxidized S1 and the less oxidized S2 PPy sample. 50 vi List of Figures 1.1 Energy levels of spin 1/2 nucleus in the presence of the magnetic field. . . 3 1.2 Spin precession 4 1.3 Schematic descriptions of N M R phenomenon 5 1.4 The magnetization is tipped down by 6 by a RF pulse with duration tp. . 7 1.5 Internal spin interactions for spin 1/2 10 1.6 Line splitting for a two spin 1/2 system (AX) with J-coupiing constant JAX- 10 2.1 Chemical structure of some common conducing polymers 13 2.2 Conductivities of polypyrrole relative to other materials 13 2.3 Redox reaction of PPy in the presence of PFg anions 14 2.4 Volume changes in PFg doped polypyrrole when it is subjected to a redox reaction 15 2.5 Chemical structures of TEAP,PC and pyrrole monomer 17 2.6 Polypyrrole deposition setup and synthesis conditions 17 2.7 Time response of reduction current during potentiostatic discharging of PPy films 18 3.1 ID 3 1 P spectrum of doped PPy films along with KH2PO4 standard solution. 21 3.2 Number of PFg ions per unit mass of PPy films at different oxidation states. 22 4.1 Modulation of the local dipolar field by molecular rotations 25 4.2 The spin-lattice relaxation time as a function of rotational correlation time for dipolar relaxations 28 4.3 Inversion-recovery pulse sequence for Ti measurement. 28 4.4 Spin echo pulse sequence for T 2 measurement 29 4.5 The effect of temperature on Tx depends on the the location of r c with respect, to the 7\ minimum 30 4.6 3 1 P Ti relaxation times as functions of temperature of PPy sample SrS3. 32 5.1 Pulse sequences for steady state NOE difference spectroscopy 35 5.2 NOE enhancement factor rj of the five samples. 36 6.1 Pulse field gradient spin-echo sequence (PFGSE) . 37 6.2 Features of molecule motions when undergoing free diffusion and restricted diffusion 40 6.3 Expected behaviors of D vs. A in free diffusion and restricted diffusions. 41 List of Figures vii 6.4 A gradient coil set designed for higher coil constant and better gradient profile 43 6.5 Signal attenuation of XE intensity in H 2 0 for gradient strength calibration. 44 6.6 One-D image gradient calibration 45 6.7 Pulse sequences for diffusion measurements 47 6.8 Echo attenuation with respect to gradient strength g of 0.05M PC sample. 48 6.9 Diffusion coefficients measured by P F G N M R experiment 49 6.10 The sheets of PPy thin films were stacked either parallel or perpendicular to B field 49 6.11 Signal attenuation 1(g) from S[ oxidized PPy sample at A=50 ms and 3s. 51 V l l l Acknowledgements I gratefully thank my advisor Dr. Carl Michal, from whom I really have learned a lot in this research project. I also thank my fellow grad students in the NMR lab: Xiang Li , Phil Eles and Tom Depew for being there whenever I've needed help. Many thanks go to Mya Warren and Dr. John Madden for sharing their experience in the synthesis and the electrochemical experiments of polypyrrole. 1 Chapter 1 Introduction 1 . 1 I n t r o d u c t o r y R e m a r k s Materials scientists have long sought to develop lightweight materials that grow and shrink significantly in length or volume when subjected to electric stimulation. Such substances could serve as drivers for novel motion-generating devices, (generally called actuators) possible replacements for the ubiquitous electric motor, which is often too bulky and heavy for smaller-scale applications. Electroactive polymers are a promis-ing alternative to materials commonly used for existing actuators because they display sufficient physical response to electrical excitation to power new classes of actuators. Polypyrrole (PPy) is one of the most studied electroactive polymers due to its en-vironmental stability, ease of fabrication and high conductivity when doped. [1] Because PPy films expand and contract in response to applied electrochemical potentials, PPy has been proposed for use as a soft actuator material. [2] Although many research groups have had success in developing real applications of PPy actuators [3, 4, 5], the relationships between microscopic structural changes and macroscopic movements are still not fully understood, partly due to the difficulties in investigating the insoluble and structurally disordered conductive PPy films [6]. 1 H and 1 3 C solid-state N M R have been applied to PPy [7, 8, 9], however, low spectral resolution from overlapping broad lines limits the information available from such studies. Cross-polarization/magic angle spinning (CP-MAS) can improve the resolution somewhat,[10, 11] but requires films to be mechanically ground into fine powders, raising concerns about changes to the material structure. Dopant ions inside PPy films are believed to play an important role in both the macroscopic and microscopic behavior of actuators as actuation is generally attributed to dopant ion insertion/expulsion into or out of the bulk polymer films. [12] There appears to be little research reported to date in the literature directly investigating the ions in PPy films. In this thesis we have employed N M R techniques to study PFg ions inside polypyrrole as a function of oxidation state, with the aims of understanding the microscopic details of the actuation, as well as elucidating the structure of the material. The fundamentals of NMR as well as the operation of Fourier spectroscopy are briefly reviewed in the introduction chapter, followed by an introduction to the material under study, PFg doped polypyrrole. Chapter 2 presents quantitative measurements of PFg doping levels as a function of oxidation state by ID 3 1 P NMR, which provides direct evidence of the ion insertion actuation mechanism. Chapter 3 explores the rotational mobility of PFg ions by measuring their T\ relaxation time constant. Chapter 1. Introduction 2 In chapter 4, nuclear Overhauser effect (NOE) experiments are performed to reveal the solvent accessibility of PFg ions in the polymer films. Chapters 3 and 4 provide a picture of the environment in which the ions reside. Chapter 5 deals with measuring translational motion of dopant ions inside polypyrrole by pulsed-field gradient (PFG) NMR. The diffusion coefficient of ions, a critical quantity in predicting actuator behavior, can be readily obtained without the need for parameter fitting via some intermediate modeling process. Lastly, findings presented in the preceding chapters are summarized in the conclusion chapter. 1 . 2 A n I n t r o d u c t i o n t o N u c l e a r M a g n e t i c R e s o n a n c e ( N M R ) Nuclear magnetic resonance spectroscopy is one of the most versatile and informative tools in many areas of science today. Often NMR readily provides difficult to acquire or otherwise unobtainable information. There are two essential features that have rendered its power. First, it uses built-in atomic labels directly from the sample. Second, it is selective and noninvasive. In addition, new NMR techniques have continued to prolif-erate, and commercially available instrumentation has reached the state of the art in terms of sensitivity, flexibility, and computational power. As examples of the diversity of the techniques of N M R consider magnetic resonance imaging, which has become a major clinical diagnostic tool, and the use of multi-dimensional solution state N M R for solving the structures of proteins. 1 . 3 P r i n c i p l e s o f N M R Each nucleus carries an intrinsic nuclear spin / , (which may be 0) and has its own nuclear magnetic moment ~jt according to Equation 1.1, where 7 is the gyromagnetic ratio. jt = j 7 (1.1) When the nuclear magnetic moment is placed in an external magnetic field ~§, the different spin states acquire different magnetic potential energies. In the presence of the static magnetic field B, an NMR sample, which normally involves a huge number of spins, produces a small amount of spin polarization (net magnetization) in equilibrium. A radio frequency (RF) signal of proper frequency can induce a transition between spin states, leading to a nuclear magnetic resonance. If the radio frequency is switched off, then the spins will relax back to equilibrium producing a measurable amount of RF signal at the resonant frequency associated with the spin flip. Take a spin 1/2 nucleus as an example: a spin 1/2 nucleus has two spin states, |f) spin up and \[) spin down. The Hamiltonian for a magnetic moment in the presence of Chapter 1. Introduction 3 m= -1/2 m=l/2, -1/2 B A E m= -1/2 No magnetic field Magnetic field Figure 1.1: Energy levels of spin 1/2 nucleus in the presence of the magnetic field. magnetic magnetic field ~B* = B'z is given in Equation 1.2: H=---ft-~S. (1.2) For positive 7 , spins parallel to the magnetic field have lower energy than spins that are anti-parallel to the magnetic field. By plugging Equation 1.1 into Equation 1.2, we find: H = - 7 T • ~B* = -jBIz (1.3) where H and Iz are quantum mechanical operators. Evaluating the energies of the spin states, we find: E = (ip I H I ip) = (tp I —-yBIz I ip) = (ip I —jBmh I (p) — mhuiL, (1-4) in which E is the potential energy of the spin 1/2 in the magnetic field, | <p) is its spin state, U>L is the Larmor frequency defined as U>L = —7 -6 , and m is the z-component quantum number, which can be 1/2 or —1/2 for a spin 1/2 nucleus (see Figure 1.1). Nuclei with different 7 values lead to different observation frequencies for different nuclei on a particular spectrometer. For an ensemble of spins 1/2 in the static ~~§ field, the distribution of nuclei in the two energy states ±1/2HUL in equilibrium is given by the Boltzmann distribution, with the lower energy state more populated than the upper: AT p-Eh/kT !Hl - 1 - p-(Eh-Ei)/h!r _ —AE/kT _ -hwL/kT / , r\ N ~ e-Et/kT ~ e ~ e - e ' ^ ' where Nh and Ni are the numbers of spins in the higher and lower energy states respec-tively. It is the small excess number of spins in the lower energy state (around 1 0 - 4 Chapter 1. Introduction 4 (a) B (b) AA M B A oo, M Figure 1.2: The behavior of M (a) at equilibrium, and (b) after perturbation after a pulse. or less) that gives the system a net magnetization M", and it is the net magnetization that determines the NMR signal. In thermal equilibrium, the net magnetization will be aligned with the static ~~§ field (Figure 1.2 (a)) and will stay there unless the system is disturbed in some way. If we somehow misalign the induced magnetization from the applied field (e.g. by applying a radio frequency pulse, as described in Section 1.4), then the Zeeman Hamiltonian (Equation 1.3) will cause the system to time evolve according to: | ip(t)) = exp (—ijrt) \ (p), and the expectation values of Ix>y will be functions of cosujit and sinujLt. Physically this means that the net magnetization precesses around the z-axis. Once the magnetization has been moved away from the positive z-direction, the static magnetic field causes it to rotate at its Larmor frequency (Figure 1.2 (b)). This motion is known as Larmor precession. A simple N M R experiment then consists of two steps, first a radio frequency (RF) pulse that is oscillating at an appropriate frequency is applied to disturb the system, transitions between the energy levels will be induced, ie, the nuclear spin system will "resonate." After the perturbation, the system will relax back to thermal equilibrium. If we were to wind a coil of wire around the axis perpendicular to the B* field, the component of the precessing magnetization in the x-y plane would induce an oscillating current in the coil which then generates a signal known as a free induction decay (FID). The resulting FID contains a wealth of information about the spin system, and the Fourier transform of the FID is the desired frequency domain N M R spectrum ready for subsequent analysis. The N M R phenomenon is shown schematically in Figure 1.3. We can also write down the classical description of the equation of motion of the net macroscopic magnetization, the Bloch Equations: ^ - = jA4x~t1 (1.6) dt  K ' The evolution of the spin magnetization in principle can be fully described by the quantum mechanical formalism if the spin Hamiltonian is written down. However the classical treatment of the spin system can often be useful to provide an intuitive under-standing, and for simple NMR experiments the classical treatment provides a complete Chapter 1. Introduction 5 (a) z (b) z (c) z X x (d) z . ( e ) (f) A B Volts Intensity Figure 1.3: Schematic descriptions of N M R phenomenon, (a) In equilibrium, the net magnetization from an NMR sample is along the static field, (b) The N M R experiment consists of disturbing this equilibrium, generally with a burst of radio frequency energy, (c) The spins dephase according to their chemical shifts and spin-spin interactions, (d) The spin system strives to re-establish its equilibrium via spin-spin and spin lattice processes, (e) Meanwhile, the induced voltage is monitored in the x-y plane by a resonant radio frequency coil. This signal is the free induction decay (FID), (f) The FID is then Fourier transformed to give the desired frequency-domain spectrum. Chapter 1. Introduction 6 Table 1.1: Properties of most common nuclei for N M R studies. Nucleus Spin Quantum Number (/) Natural Abundance (%) 7 ( l O W / T s ) 1 H 1/2 99.98 26.7520 2 H 1 0.0156 4.1067 1 3 C 1/2 ' 1.108 6.7265 1 5 N 1/2 .0.365 -2.7108 1 9 F 1/2 100 25.167 3 1 p 1/2 100 10.829 and accurate picture of the phenomenon. Table 1.3 lists the properties of the most common nuclei in NMR studies. 1 . 4 Effects of Pulses It has been pointed out earlier that if the magnetization can be moved away from its equilibrium position, then precession occurs and a signalmay be detected. The excitation is produced by a second magnetic field B\, which oscillates at the appropriate radio frequency. This field is created by a current in a coil wound perpendicular to £?. The coil is often the same coil in which the signal is received, and by convention we define the coil axis as the x axis. If the excitation frequency Ui is equal to the Lamor frequency UL, U)I=U)L = jB. (1.7) The field Bi interacts with the magnetization A? to produce a torque which moves the magnetization towards the x-y plane. The magnetization will precess around the RF field at an angular frequency 7 B 1 . At the same time, it is precessing about the magnetic field at a frequency jB. This rather complicated motion, simultaneous precession about two axes, is referred to as a nutation. The final position of the magnetization will depend upon the length of time for which the radio frequency is applied. (See Figure 1.4) The B\ field is applied as a pulse of duration tp, which usually lasts for a few microseconds. The tip angle 6 which the magnetization is tipped from the z axis is easily calculated as: 0 = lBxtp (1.8) We can determine the relationship between the pulse width and the tip angle by varying the duration and examining the resulting signal intensity. The signal will be a maximum for a tip angle of ir/2, a null at ir, a negative maximum at 37r/2, and return to a null at 2n with the magnetization back along the z axis. Chapter 1. Introduction 7 z M ( t p ) M s y X Figure 1.4: The magnetization is tipped down by 9 by a RF pulse with duration tp. The nuclear spin relaxation times (Ti and T 2) refer to characteristic time scales for spin magnetization to regain equilibrium after an RF pulse perturbation. Relaxation times probe molecular dynamics because molecular motions will generate fluctuating magnetic fields at the site of the nuclear spins, causing transitions in the spin states and thus bringing the system back to thermal equilibrium. The spin relaxation times Ti and T 2 are essential parameters in the NMR experiments because not only do they provide valuable information as to motion and dynamics, but they also set time scales over which various aspects of the spin system stay coherent after being disturbed from equilibrium. Ti is the spin-lattice or longitudinal relaxation time, and is a measure of how fast the magnetization relaxes back along the z-axis. T 2 is the spin-spin or transverse relaxation time, and is a measure of how fast the magnetization dephases in the transverse(x-y) plane. Phenomenologically, the time dependence of macroscopic magnetization, having com-ponents Mx, My, and Mz, during the relaxation is described as 1 . 5 T i , T2 Relaxation Times Mz(t) Mx{t) My(t) My(0) sm(u0t)e T 2 M0 + (Mz(0)-M0)e-^ Mx(0)cos(uj0t)e~% (1.9) (1.10) (1.11) (1.12) where M0 is the equilibrium magnetization. The Bloch equations including relaxation Chapter 1. Introduction 8 can then be written as: dA4 7 A ? x ~B* - R(tf - M y ) (1.13) R 1 . 6 N M R S p e c t r o s c o p y : F I D S i g n a l a n d t h e F o u r i e r T r a n s f o r m Once an RF pulse with spectrometer reference frequency w r e / is applied to disturb the spin system from equilibrium, the magnetization relaxes towards equilibrium, and an oscillating current (FID) is induced in the coil. This signal is amplified, filtered, and captured by analogue-to-digital converters to allow for processing and analysis on a com-puter. Because the signals oscillate as quickly as several hundred MHz, the quadrature receiver combines the N M R signal, which oscillates at the Larmor frequency UJQ, with the reference signal to generate a new signal with the relative Larmor frequency Q.Q. The range of f20 values is typically in the kHz range, reducing the speed requirements of the analogue-to-digital converters, and reducing the amount of data that must be collected and processed. Depending on the exact value of coref, there may be some signals with oscillation frequencies larger than toref, and some with oscillation frequencies smaller than u>ref due to chemical and other resonance shifts. In order to separate frequencies above uref from those below, the receiver supplies two output signals, which for exponentially decaying resonances, are of the form: The two signals can be combined into a complex signal s(t), s(t) = s^(£) + is sit) and positive and negative values of Qo can then be distinguished. Note that a more general form of s(t) is s(t) = X^ 5'(0> because normally the NMR signal is a superposition of signals with different relative Larmor frequencies in the system. The time domain signal s(t) can then be converted into a frequency spectrum by Fourier transformation, ^0 — ^0 ~ ^ref (1.14) sA(t) oc c o s ( £ V ) e " t / T 2 sB(t) oc sin(f2 0t)e" t / : r 2 (1.15) (1.16) if we substitute s(i) into Equation 1.16, we obtain: St(tt) = a J L(Q;Q ; ) (1.17) Chapter 1. Introduction 9 where L(Q; Q;) is caled the complex Lorentzian, and a; is the amplitude. The real part of L(Q; Cli) is called the absorption Lorentzian Re{L(f2; fij)}: R e { L ( f i ; f i ' ) } = ( i 7 w r | ^ <18> and the imaginary part is called the dispersion Lorentzian Im{L(f2; Qi)}: •i-Wn=ft)}- ( J ^ ^ , (!•") The full-width-at-half-height (FWHH) of the absorption Lorentzian is just ^ Hz, with the important implication that a fast signal decay leads to a broad peak and a slow signal decay leads to a narrow peak. 1 . 7 I n t e r n a l S p i n i n t e r a c t i o n s The external spin interactions involve the interaction with the magnetic field generated by the external static B field' and an oscillating Bi field which is perpendicular to it. The internal spin interactions involve interaction with the magnetic and electric fields generated by the sample itself. The internal spin interactions are the interactions we are interested in because they encode information (eg. distances, angles, bonding) about the spin system. For spin 1/2 nuclei, the relevant internal spin interactions are chemical shift, J-couplings, and dipole-dipole couplings. A l l of these interactions result in small perturbations of the resonance away from the Larmor frequency. The effect of each internal spin interaction is described below. 1.7.1 Chemical shift Nuclei of the same type can achieve the resonance condition at different frequencies because the local field each nucleus experiences is affected by its electronic environment. These bond-dependent shifts are known as chemical shifts. When a static field B* is applied, the ~S field will induce electron currents that produce small magnetic fields opposed to B*. (See Figure 1.5 (a)) As a result, the local magnetic field felt by each nucleus is given by: Blocal = B(l - o) (1.20) where a is the screening or shielding constant and o depends on the details of the electron wavefunctions in the vicinity of the nucleus. A chemical shift can be expressed as a frequency differences in Hertz with respect to a reference signal from a standard solution, but since it will depend on the ~~§ field, chemical shift is typically reported in ppm (parts per million): . . Chemical shift from reference 6 Chemical shift WP™) = Spectrometer frequency X ° ( } Chapter 1. Introduction 10 (a) shielding Figure 1.5: Internal spin interactions for spin 1/2. (a) Chemical shift (b) J-couplings (c) Direct dipole-dipole couplings. f A f x J A X J A X - * • Figure 1.6: Line splitting for a two spin 1/2 system (AX) with J-coupling constant JAX-1.7.2 J-couplings The J-coupling between nuclear spins arises from the influence of the bonding electrons on the magnetic fields created by the nuclear spins.(See Figure 1.5 (b)) While the chemical shift indicates the local electronic environment, the J-coupling provides a direct spectral manifestation of the chemical bond. If two spins 1/2, A and X , are coupled to each other with a coupling constant J Hz, and if J is larger than the linewidths, then the peaks of A and X will each be split into two peaks separated by J Hz. (See Figure 1.6) The peak of A in system A X n will be split to 2nl + 1 lines if there are n spins of X interacting with the nucleus A . The value of the coupling constant depend on the molecular geometry and the number of bonds between coupled nuclei. For protons in organic materials, three-bond J-couplings are typically around 7 Hz. One-bond J-couplings between *H and 1 3 C are typically around 135 Hz, while one-bond J-couplings between directly-bonded 1 3 C spins are usually around 50 Hz. In solid state, J-coupling typically does not result in a clear separation of line splitting because typical linewidths obtained in the solid state (100Hz - 10's of kHz) are often greater than the J-coupling. Chapter 1. Introduction 11 1.7.3 Dipole-dipole couplings Dipole-dipole couplings represent the direct magnetic interactions of nuclear spins with each other. Since each nuclear spin represents a magnetic moment, it generates a mag-netic field Biocai, looping around in the surrounding space according to the direction of the spin magnetic moment. A second nuclear spin interacts with this magnetic field. (See Figure 1.5) The local field Biocai at a nucleus I, generated by a nucleus S (in the high field approximation where only the secular part of the dipole field is retained) is given = * < W / " - ' > (1.22) RIS where ps is the magnetic moment of S, rjs is the internuclear distance, and 9 is the angle between the internuclear vector and the ^ field. Thus the dipolar Hamiltonian is given by: HD = ~^7j f i 2 (3cos 2 .g / s - 1)7,5, (1.23) where po is the permeability constant. The constant term in front of Equation 1.23 corresponds to the magnitude of the dipolar interaction which is: dI8 = i p ^ ' (1-24) Equation 1.24 provides the through-space coupling in units of Hz. For example, two  1H spins separated by a distance of 2 A experience an interaction djs = —15 kHz. The interaction goes down according to the inverse cube of internuclear distance, and scales linearly with the gyromagnetic ratio of each interacting spins. Note that the dipolar Hamiltonian in Equation 1.23 is orientation-dependent, and since the 6 changes as the molecule rotates the consequence is that in an isotropic liquid where molecules can rotate freely, the intramolecular dipolar-dipole couplings will average to zero. Chapter 2 12 A n Introduction to P F 6 Doped Polypyrrole 2 . 1 P F g D o p e d P o l y p y r r o l e B e l o n g s t o t h e C l a s s o f C o n d u c t i n g P o l y m e r Conducting polymers are conjugated polymers, namely organic compounds that have con-tinuous overlapping IT orbitals, resulting in delocalized bonds along the polymer backbone that create a band structure similar to that found in semiconductors. Chemical struc-tures of some common conducting polymers including polypyrrole are shown in Figure 2.1. When electrons or holes are introduced to the polymer chains, a process called doping, those positive charge carriers (holes) and negative charge carriers (electrons) enable electric conduction. The typical range of the electrical conductivity of PPy rel-ative to other materials is shown in Figure 2.2. The electronic states of the polymers can be reversibly controlled by doping (oxidation-reduction processes), and the oxidation state is related to tunable changes in properties: structural, mechanical or optical, which leads to the use of conducting polymers in a host of applications including actuators, electrochromatic devices, and supercapacitors. In electrochemical doping, the conducting polymer is in electrical contact with an electrode in an electrochemical cell. Electrons are added to or removed from the polymer via the electrode, thereby changing the oxidation state. The oxidation process leads to extra holes or electrons on the polymer chains (called p-doping or n-doping), which are balanced by the flux of ions to or from the electrolyte to maintain the electroneutrality. These ions are referred to as dopants or counterions. PFg doped polypyrrole in this study is prepared electro chemically and is a p-doped conjugated polymer when oxidized. PFg dopant ions are incorporated inside the polymer to balance the charge. (See Figure 2.3) The doping level of the oxidized state (as-grown state, see experimental section) reaches one dopant ion for every three to four monomers. Electrical conductivity of polypyrrole in the oxidized state is typically 3 ~ 4 x 104 S/m compared to ~ 1 x 101 S/m in the fully reduced state. Chapter 2. An Introduction to PF6 Doped Polypyrrole 13 P o l y p y r r o l e ,Ox R e v CTs/cm)V l t y ' ° 6 ^ 10"' IP' 10"'° 10'3 ^ s Figure 2 .2: Conductivities of polypyrrole relative to other materials. Chapter 2. An Introduction to P F g Doped Polypyrrole 14 PPy + P F 6 oxidation reduction PPy + (PF 6 ) + e" + (PF 6 ) 0 n + n e Figure 2.3: Redox reaction of PPy in the presence of P F 6 anions. 2 . 2 P F 6 D o p e d P o l y p y r r o l e a s M u s c l e - L i k e A c t u a t o r s As noted above, many properties of technological interest related to conducting polymers result from oxidation-reduction (redox) processes. When a neutral conducting polymer is electrochemically oxidized by applying sufficient potential in an electrolyte, positive charges (holes) are generated along the backbone and solvated counterions (dopant ions) are forced to enter the polymer to maintain electroneutrality. This promotes the open-ing of the polymeric structure and a significant increase in free volume. (See Figure 2.4) Opposite processes occur during reduction: electrons are injected into the polymer, positive charges are eliminated, and the dopant ions and solvent molecules are expelled into the solution. The redox properties of conducting polymers that are responsible for reversible and controllable changes in volume open new possibilities for the development of actuators often referred to as artificial muscles. PFg doped polypyrrole is one of the most promising polymer actuators for artificial muscles because it can attain moderate strain (~ 2%) at low working potential. It is believed that the actuation is due to the dimensional changes produced by electrochem-ically inserting PFg ions into the polymer electrode. The cations in normal cases are selected to be big enough to prevent then from entering the polymer, so that the volume changes are attributed to dopant ions only. Chapter 2. An Introduction to PF6 Doped Polypyrrole 15 Figure 2.4: Volume changes in P F 6 doped polypyrrole when it is subjected to a redox reaction. 2 . 3 C o m p a r i s o n b e t w e e n C o n d u c t i n g P o l y m e r A c t u a t o r s a n d O t h e r A c t u a t i o n T e c h n o l o g i e s Nature's skeletal muscles demonstrate the most impressive actuator that engineers have been trying to mimic. They provide billions of work cycles involving contraction of more than 20%, generate stresses of ~ 0.35 MPa, and contract at 50% per second with maximum efficiency of ~ 40% [13]. State-of-the-art conducting polymer actuators operate at voltages of a few volts and can generate high strains (~ 26%), high strain rates (~ 12%/s ), and stresses up to 100 times greater than those of mammalian muscles (~ 34 MPa) [13]. However, the drawbacks related to conducting polymers are: (1) the need for encapsulation: they generally need to be wet and so must be sealed within flexible coatings; (2) low energy conversion efficiency: typically < 1%; (3) relatively short cycle life: the actuator performance usually degrades after several thousand cycles [14, 15]. While ionic electroactive polymers (EAPs), including conducting polymers work on the basis of electrochemistry - the mobility or diffusion of charged ions; there are other types of EAPs that are driven by electric fields such as dielectric elastomers . Actuation of dielectric elastomers is mainly caused by Maxwell stress, which results from the attraction between charges on opposite electrodes and the repulsion between like charges. These systems can generate strains of 120%, stresses of 3.4 MPa, and a peak strain rate of 34000%/s for 12% strain [13]. The major shortcoming of dielectric elastomers is the requirement of high activation voltage - typically 1 — 5 kV, which limits their use in the mobile devices. Some other common materials for actuation besides EAPs are piezoelectric ceramics and shape memory alloys. In piezoelectric materials, mechanical stress causes crystals to electrically polarize and vice versa. When a voltage is applied across a piezoelectric material, it will deform. The efficiency of piezoelectric materials are very high, up to Chapter 2. An Introduction to PF6 Doped Polypyrrole 16 90%, and the response is fast. The disadvantage is the small strain of less than 1%, yet they have found use in producing motions on the order of nm in the control of scanning tunneling microscopes. Shape memory alloys (SMA) actuate through phase changes in crystal structures. A temperature difference of only about 10°C is necessary to initiate this phase change. These materials can generate strains of up to 8%, stress of up to 200 MPa but require the inefficient conversion from thermal energy to mechanical energy. [13, 16] There are many other types actuators around which are reviewed in Madden et al. [13]. A l l of these actuators have their advantages and disadvantages, and efforts are being made to improve their performance. Because polymer actuator can generally be small, light and cheap in addition to its low operating voltages and high stress, many future actuator applications are expected to be based on polymer materials. 2 . 4 S y n t h e s i s o f P F g D o p e d P o l y p y r r o l e The conditions during synthesis of PPy determine the physical properties of the resulting polymer to a great extent. In this work we follow the methods of Yamaura et al. [17] for synthesizing PPy, and the subsequent electrochemical experiments are done by methods developed in collaborator Dr. John Madden's lab. 2.4.1 Electrodeposition The oxidative polymerization of polypyrrole (PPy) was carried out galvanostatically at a current density of 0.125 mA/cm 2 at —30°C from a solution of 0.6 M purified pyrrole, 0.05 M tetraethylammonium hexafluorophosphate (C^Hs^N+PFg (TEAP) and 1% dis-tilled water by volume in propylene carbonate C4H 6 03 (PC). The chemical structures of TEAP, PC, and pyrrole monomers are shown in Figure 2.5. The electrolyte was de-oxygenated with nitrogen before use. Polypyrrole was deposited onto a polished glassy carbon substrate for a duration of 8 hours, resulting in nominal film thickness of approx-imately 16 pm ~ 26 pm and density of ~ 1.5 g/cm 3. A copper sheet is used as counter electrode. Nonconducting Kapton tape is spirally wound to mask the glassy carbon cru-cible's surface, so that strips of polypyrrole films with a width of 1 cm can be produced. The as-grown film from the above experimental conditions is in its oxidized state with an equilibrium electrochemical potential ~ 0.4 V vs. Ag/AgC104, and ~ 0.5 V vs Ag/AgCl in different depositions. Figure 2.6 shows the schematic deposition setup and the syn-thesis conditions. Note that in the ion content experiment, all potentials are reported against Ag/AgC104- However we found that the potential of the Ag/AgC10 4 reference electrode drifts on a time scale of weeks, thus in all other experiments the potentials are reported against Ag/AgCl . Chapter 2. An Introduction to P F 6 Doped Polypyrrole 17 Pyrrole Molecular Formula C 4 H 5 N Molecular Weight 67.09 Density 0.967 g/ml TEAP Molecular Formula q F g N P Molecular Weight 275.22 C H 3 Molecular Formula C 4 H 6 0 3 I / Molecular Weight 102.09 0 O Figure 2.5: Chemical structures of TEAP,PC and pyrrole monomer. •> Galvanostatic depostion at 0.125 mA/cm Deposition electrolyte: 0.05M TEAP, 0.6M pyrrole, and 1% water in PC Working electrode: Glassy Carbon .> Counter electrode: Copper sheet -> Polypyrrole thin film Figure 2.6: Polypyrrole deposition setup and synthesis conditions. Chapter 2. An Introduction to PFe Doped Polypyrrole 18 0.02 r -0.1 I -0.12 1 ' ' 1 1 0 500 1000 1500 2000 Time (s) Figure 2.7: Time response of reduction current during potentiostatic discharging of PPy films. 2.4.2 Electrochemical Reduction Subsequent reductions of the as-grown film were done potentiostatically in a three-electrode cell at several potentials in PC containing 0.05 M TEAP. For the ion content study, samples S2, S 3 , and S4 were reduced to potentials of E = 0.2 V, —0.2 V, and —0.4 V vs Ag/AgClC>4 respectively. For the other studies, sample S2-S'4 were reduced to E — 0.25 V, -0.2 V, -0.4 V vs Ag/AgCl . Different oxidation states of polymer films refer to their open circuit potentials vs a reference electrode. As-grown P F g doped films are reduced directly on the glass carbon crucible to achieve the best electrical contact between the working electrode and the polymer. Carbon fiber paper is used as counter electrode to maximize the electrode surface area. When each potentiostatic discharg-ing had reached its equilibrium (See Figure 2.7), a strip of polypyrrole film with area 1 cmx4 cm was peeled off the glassy carbon surface and rinsed thoroughly with propylene carbonate. For the ion content experiment, the films were left in ambient environment for two days and flattened into sheets to make them easier to handle. For all the other measurements, the films were soaked in PC containing 0.05M T E A P for two days and rinsed with pure PC immediately before use. 2.4.3 Conductivity Measurement The electronic properties of conducting polymer films were characterized through con-ductivity measurements. It is advantageous to have high conductivity since it reduces the IR drop (potential difference) along the polymer. The conductivity is measured us-ing a four-point probe method to avoid contact resistance. Four thin copper wires are Chapter 2. An Introduction to PF6 Doped Polypyrrole 19 forced to contact the two ends of the polymer strip and are directly connected to an HP multimeter for four-point measurement. Conductivity o is then calculated based on: R - l - l - L - (2.1) where I is the length of the film between the center two probes, w is the width and t is the thickness of the film. Since conductivity is achieved by percolation of charge carriers through ordered re-gions, it is dependent upon microscopic polymer morphology, and therefore is extremely sensitive to various experimental parameters, for example, the electrode preparation, de-position current density, and temperature. The maximal conductivity of our as-grown film is obtained to be 323 S/cm. 20 Chapter 3 PFg - Ion Content in the PPy Films as a Function of Oxidation States 3 . 1 E x p e r i m e n t a l M e t h o d 3.1.1 ID 3 1 P Spectroscopy Our ID 3 1 P NMR experiments were performed on a homebuilt N M R spectrometer work-ing at 8.4 T with a solid state probe tuned for phosphorus at 147 MHz. Sheets of PPy films were placed neatly inside the sample tube and oriented so that the sheets' normal was always perpendicular to the RF coil axis to minimize the RF loading of the probe's resonant circuit, which manifests itself as a lowering of the quality factor (Q) of the probe. The ID N M R phosphorous spectra of PPy films were obtained by single 7r/2 pulse excitation and were referenced to an external potassium phosphate ( K H 2 P O 4 ) standard solution in order to obtain quantitative information of the film samples. The 3 1 P peak in K H 2 P O 4 is well separated from the 3 1 P peak in T E A P allowing correct intensity integration. The time between each acquisition was always >5 times longer than the longest Ti to ensure full relaxation. 3 . 2 R e s u l t s Four pieces of PFg doped PPy (Si to 5 4 ) were prepared in four different oxidation states at E — 0.43 V, 0.2 V, -0.2 V, -0.4 V vs Ag/AgC10 4 . Each piece, along with a capillary of the standard solution KH2PO4 was placed in an NMR tube. Figure 3.1 shows a typical ID 3 1 P spectrum of the sample. The sharp peak near 1 ppm corresponds to the phosphorus of K H 2 P 0 4 . The broad peak is due to PFg ions inside the PPy film. The N M R line from PFg ions in the PC solvent showed the expected J-splitting pattern, but when in the films, an additional broadening masks the J-splittings. The lineshape of ions in the films can be reproduced from the J-split line observed in the solvent by convoluting with a 600 Hz F W H M Gaussian function. Relaxation measurements, described in Chapter 4, reveal that the ions tumble rapidly, suggesting that this 600 Hz broadening is likely due either to couplings to other spins inside the films or to magnetic field inhomogeneity rather than intramolecular P-F dipolar couplings. Because the number of K H 2 P O 4 molecules in the standard is known, the number of Chapter 3. PF6 Ion Content in the PPy Films as a Function of Oxidation States 21 5 0 0 - 5 0 - 1 0 0 - 1 5 0 - 2 0 0 Chemical shift (ppm) - 2 5 0 Figure 3.1: The broad peak is the PPy peak and the sharp peak is KH2PO4 peak. P F 6 ions in the films can be obtained from a comparison of the integrated intensity of the peaks according to Integrated area of PPy peak Number of 3 1 P in P F J ions Number of 3 1 P in K H 2 P 0 4 ~ Integrated area of KH2PO4 peak (3-1) The number of dopant ions divided by the film's initial mass (mass of the film before reduction) is plotted against the electrochemical potential E in Figure 3.2, and the four points can be fit into a straight line. The result provides direct evidence for ion insertion mechanism of actuation since the number of dopant ions increases linearly with increasing oxidation potential. With known total deposition charges spent, sample weight, and number of PFg ions in each sample, the number of pyrroles and solvent molecules in each sample can be extracted from Equations 3.2 and 3.3, where Q is the total deposition charge, Np, Ni, Nsoi are the number of pyrroles, PFg ions, and solvent molecules, and m is the total mass. Q 2NP + Ni (3.2) m = Npmp + Ni'nii + Nsol maol (3-3) Chapter 3. P F 6 Ion Content in the PPy Films as a Function of Oxidation States 22 too C O O O c3 .2 "3 a PL, 4) 1 7e+17 6e+17 5e+17 4e+17 3e+17 2e+17 le+17 0 linear fit s -0.6 -0.4 -0.2 0 0.2 E(V) vs Ag/Ag + 0.4 0.6 Figure 3.2: Number of P F 6 ions per unit mass of PPy films at different oxidation states. The doping level (expressed in terms of ion to pyrrole ratio, Ni/Np), as well as solvent molecule to pyrrole ratio (Nsoi/Np) as functions of oxidation state are shown in Table 3.1. Conductivity information of each sample is also included in Table 3.1 for reference. 3 . 3 D i s c u s s i o n Previous investigations regarding doping levels of doped conducting polymer have relied on elemental analysis and X-ray photoelectron (XPS) spectroscopy to measure dopant concentration [18, 19], but both techniques are destructive and insensitive to chemical environment, and XPS also suffers from the fact that only a thin layer of the polymer is probed. ID 3 1 P NMR is a noninvasive tool used here to measure PFg ion content, and because 3 1 P is only contained in PFg ions, selective information can be obtained without any ambiguity. To the best of our knowledge, systematic measurements on doping level vs. oxidation state have not been carried out previously. Our experiments have shown that the ion content increases in a linear fashion with the electrochemical potential of the polymer,, directly confirming the ion intercalation mechanism of actuation. With known ion content, the number of pyrroles and solvent molecules per unit mass can be easily extracted. The obtained doping level for the as-grown (S-L) sample (Ni/Np = 1/3.83 = 0.26) is consistent with previously determined values which is 0.25 [20, 21] and 0.27 [22] in the literature. This suggests that all the Chapter 3. PF6 Ion Content in the PPy Films as a Function of Oxidation States 23 Table 3.1: Doping levels of samples S1-S4 in terms of ion to pyrrole ratio (Ni/NP)and solvent molecule to pyrrole ratio (Nsoi/NP ). Conductivities are determined by a four point probe method. PFg doped PPy samples E(V) vs Ag/Ag+ Conductivity (S / cm) Ni/Np Nsol/NP Si (as-grown) 0.43 250 1/3.83 1/0.67 s2 0.2 189 1/6.23 1/0.92 S3 -0.2 13 1/10.63 1/0.63 Si -0.4 0.23 1/21.26 1/0.74 ions are detected, and examining ion content via ID 3 1 P N M R spectroscopy seems to be a reliable and straightforward method. The ion concentration inside the oxidized film is calculated to be at least 2 M (number of moles of PFg ions divided by the film volume). Note that we assume a uniform deposition in this ion content analysis and it is believed that the cation (tetraethylammonium, N(C2H 5)^") is unable to enter the films because it is much larger than PFg ion [22, 23]. In addition, doped polypyrrole resembles a capacitor in a way that the amount of charge can be stored in relation to the strength of the applied potential ( Q oc V), but charge does not accumulate on two conductors separated by a dielectric. Instead the charge accumulates inside the polymer film, at interfaces between the polymer and the enclosed solvent. Because of the porous nature of the films, they have a high surface area and thus a high capacitance. 24 Chapter 4 Rotational Motions of P F ^ Ions in the PPy Films Because the relaxation of 3 1 P in PFg ions is due to the dipolar couplings from 1 9 F that are modulated by the rotation of the ions, the rotational motions of PFg ions in the PPy films can be accessed by measuring the 3 1 P Ti relaxation time. Before describing the experimental results, we begin by presenting a brief introduction to relaxation theory and the analytical relations between the rotational dynamics and the Ti relaxation. The most important relaxation mechanism for spins :l/2 is often the through-space dipolar coupling. There are intramolecular dipolar interactions inside a molecule and the intermolecular dipolar interaction among spins of different molecules. In the interior of a molecule, the variation of dipolar coupling of spins arises almost uniquely from the rotation of the molecule, the variation of the distance between the spins due to vibrations being negligible. (See Figure 4.1) We will discuss the case of intramolecular interaction only as it is the case of 3 1 P relaxation in PFg ions. The general approach to relaxation theory we will pursue involves the spin density matrix. We shall describe the framework of this semi-classical relaxation theory below and show how the Ti and rotational correlation time r c are linked together in the spin-1/2 case. For details of the derivation we refer to Abragam [24]. We consider a spin system of Hamiltonian H = H0 + Hi(i), where H0 is the main static Hamiltonian, and Hi(t) represents all of the fluctuating interactions responsible for relaxation. The density matrix operator a contains all of the physically significant information about the ensemble of spins. We can write the equation of the evolution of the density matrix (evolution of the spin system): We can isolate the effect of the fluctuating Hamiltonian by switching to the Heisenberg picture, in which it is the observables that vary with time. In the Heisenberg picture, we consider a(t) = exp(iH0t)a exp(—iH0t). Equation 4.1 becomes: 4 . 1 R e l a x a t i o n T h e o r y da ~dt •i[{HQ + Hyl<x]. (4.1) d a ~dt = -^(0 ,5] , (4.2) Chapter 4. Rotational Motions of PF6 Ions in the PPy Films 25 Local Dipolar Field Molecular Tumbling Figure 4.1: Modulation of the local dipolar field by molecular rotations. By integrating Equation 4.2 and taking its average over all identical molecules of the sample (represented by the bar), we obtain: o{t) = ?(0) - i f [H\(t'),a(t')]dt! (4.3) Jo An approximate solution to Equation 4.3 can be found by iteration: o{t) = o{Q)-i [H^MVW - dt' df[W1(t),[H1(t''),o(0)}} + .-. (4.4) Jo Jo Jo The first integral on the right hand side vanishes because the dipolar coupling averages to zero in solution, so the leading non-zero term becomes: o(t) - o(0) = - f dt' f dt''[%(t),[%{t"),d{0)]\ (4.5) Jo Jo To deal with this integration, we assume the fluctuating Hamiltonian H\(t) takes the form: Hx{t) = Y,VaFa(t) (4.6) a where the Va is the spin operators and the Fa are random functions. The autocorrelation function of the random function Fa is defined as: Fa(t)F;(l') = GaP(\t-t'\) (4.7) Chapter 4. Rotational Motions of PF6 Ions in the PPy Films 26 Equation 4.5 can be written in the following form: a f dt' f dt"Ga0{\t'-t/'^e^1'-^"^ (4.8) Jo Jo The integral in Equation 4.8 is exactly the Fourier transform of the autocorrelation function Gap, which can be represented by a spectral density function Ja/3(cj): roo Jap(w)= / G a P ( \ t - t ' \ ) e M - i u \ t - t ' \ ) d ( t - t ' ) (4.9) Jo It can be shown that the dominant contributions to Equation 4.8 are the terms where a. — P, therefore from now on we retain only the relevant terms, and will denote the function Jaa by Ja. With the assumption that o(t)—o(0) is small, we have approximately: a{t) - o(0) ^ do t dt Equation 4.8 becomes : do dt = (? - °eq)]]JaM. (4.11) This is the master equation for the evolution of the density matrix under the effect of a random perturbation. The expectation value of the physical observable Q in an ensemble can be obtained from (Q) = Tr(Qo), therefore the master equation for the relaxation of a physical variable is the following: ^ = - J > K ) { ( [ [ Q , Va] ,V2]> - <[[Q, Va], V^])eq} (4.12) a In the case of dipolar relaxation between two spins of 1/2, the dipolar Hamiltonian takes the form: H1,2 = - ^ { t . t - 3 I l z I 2 z } (4.13) where r is the distance between I\ and I2. In order to perform relaxation calculations, we must express E\2 as a sum of products of spin operators and random field operators. By expressing the I\, I2 in spherical coordinates, we obtain the proper representation for the dipolar Hamiltonian where 9 and (p are random functions of time: / / i , 2 = -^53F*((?^)rm(/ 1 >/ 2) . (4.14) Chapter 4. Rotational Motions of PFe Ions in the PPy Films 27 with F0 = y | ( 3 c o s 2 0 - l ) (4.15) F±i = sin 9 cos 9 exp(±i(p) F±2 = 3/2sin2<9cos6>exp(±i0) To =• -j=(3hzI2z-t-t) T±i = T l / 2 ( / i ± / 2 , + hzh±) T ± 2 = l / 2 / 1 ± / 2 ± In principle, once we have the Hamiltonian in the appropriate representation, then we can evaluate the evolution of the desired observable Iz using the master equation (Equation 4.12). We now show the results for longitudinal relaxation (Ti) in the case of two spins 1/2, by substituting Iz, the longitudinal spin operator, in place of Q in Equation 4.12 and simplifying. We find: d(Iu) dt = -Pl({hz) - (hzU) - Pcross{{hz) - (hz)eq) (4.16) with Pi = 1 1 7i2722fi2 { 3 1 6 \ 10 r-6 T° \ 1 + o;?rc2 + 1 + ( W l - u;2)2rc2 + 1 + ( ( W l + w 2) 2r c 2 ) { ' \ _ 1 7 ^ % f 6 1 \ (Am P C T 0 S S ~ 10 T-6 T c \ l + ( ( a ; 1 + u ; 2 ) 2 r 2 1 + ( W l _ ^ 2 J where pi corresponds to the longitudinal relaxation rate 1/T 1 ; and to\ and u2 are the Larmor frequencies of spin 1 and spin 2 ( i i , i 2 ) . r c is a time scale that characterizes the decay time of the autocorrelation function G(\ t — t |) and is referred to the rotational correlation time in the case of two spin 1/2 dipolar relaxation. Through Equation 4.17, we are able to probe the rotational correlation time r c of the nucleus by measuring Ti . If we plot the longitudinal time constant Ti against the rotational correlation time r c , typically we get the curve shown in Figure 4.2. Molecules in gases have the lowest correlation times at room temperature, followed by small molecules in non-viscous liquids. Notice that there also exists a cross-relaxation between the longitudinal magnetization of the spin i i with that of i 2 in Equation 4.17. This coupling is the basis of the Overhauser effect, described in Chapter 5. Chapter 4. Rotational Motions of P F 6 Ions in the PPy Films 28 Figure 4.2: The spin-lattice relaxation time as a function of rotational correlation time for dipolar relaxations. pi pi/2 Figure 4.3: Inversion-recovery pulse sequence for Ti measurement. 4 . 2 E x p e r i m e n t a l M e t h o d 4.2.1 Relaxation Measurements 3 1 P relaxation measurements were performed on a homebuilt N M R spectrometer working at 8.4T with a variable temperature probe tuned for phosphorus at 147MHz. A commonly used pulse sequence for measuring Ti is called inversion-recovery. It is given by 7r — r — | — AQ.(See Figure 4.3) The first pulse in the sequence is a IT pulse which tips the equilibrium magnetization Mz down to —Mz. The magnetization relaxes back towards thermal equilibrium during the interval r, and the process is monitored by the second pulse, which projects the net magnetization into the x-y plane. The signal intensity is therefore r dependent and can be fit with the function J(r) = IQ(\ — 2e T i ) to extract Ti . T 2 relaxation can arise from numerous mechanisms and often cannot be correlated as easily to a single motional parameter. It does however give some measure of the motion of the spins. The faster the motion and the more homogeneous the spin's environment, Chapter 4. Rotational Motions of PF6 Ions in the PPy Films 29 pi/2 pi 111 « 1 / 2 > ft A Figure 4.4: Spin echo pulse sequence for T 2 measurement. the longer the T 2 . The T 2 measurement is routinely done by the spin echo sequence | — | — 7r — | — y4Q.(See Figure 4.4) The first | pulse tips the magnetization down to the x-y plane, and spins will start to dephase during the | period. Then a TT pulse refocuses the magnetization producing an echo signal. The experiment is arrayed with different^r, and the corresponding peak intensity is an exponential decay with r via I(T) = I§e T 2 . 4 . 3 R e s u l t s In the case of the PFg ion, we assume that the dominant mechanism for phosphorus relaxation is the P — F intramolecular dipolar interaction from the isotropic rotation of ions. Neglecting the contributions from distant fluorines (not directly bonded), the Ti relaxation equation that describes the interactions of the 6 neighboring fluorines on phosphorus is given by Equation 4.19 if the motions of the fluorines are uncorrelated with that of the phosphorus. In the PFg ion, the motions are clearly correlated, but Equation 4.19is a good approxamation nevertheless [24]. 1 = m \ 2 1 YpTFh-Ti n F U J 10 r 6 3 Jprl 1 + (up - WF)2T* 1 + (ujP + upfrl (4.19) where np is the number of directly bonded 1 9 F to 3 1 P . We examined the 3 1 P Ti relaxation times on three wet PPy samples 5 , 1-5 , 3 along with two PC solvent samples containing 0.05 M and 2 M TEAP-respectively. Because in solution state the molecules are tumbling quickly, they are in a fast motion regime where OJ2T^ <C 1, the relationship between rotational correlation time r c and Ti can be reduced to a simple form: i / / A ) \ 2 ipiW rc can then be extracted easily from Equation 4.20. Results from the two solvent samples at room temperature are shown in Table 4.1. Chapter 4. Rotational Motions of PF6 Ions in the PPy Films 30 Table 4.1: Ti relaxation times and the rotational correlation times r c of two solvent samples. Sample 3 1 P Ti (s) Tc (PS) 0.05M PC 2M PC 4.10 4.10 8.08 8.08 T, Figure 4.5: The effect of temperature on Ti depends on the the location of r c with respect to the Ti minimum. As for probing rotational motions of PFg in the films, it is insufficient to only measure Ti at room temperature because one Ti corresponds to two rotational correlation times according to Figure 4.2. Therefore, Ti of three PPy samples are measured as functions of temperature, because the effect of temperature on Tx depends on the the location of r c with respect to the Ti minimum. (See Figure 4.5) For systems having long correlation times, warming the sample (reducing the correlation time) reduces Ti ; for systems with short correlation times, warming the sample generally increases T\. The spin-lattice relaxation time constants as functions of temperature for three PPy samples that are in different oxidation states (Sj-S 3) are shown in Table 4.2 and plotted in Figure 4.6. Ti of the oxidized sample Sx and the less oxidized sample S 2 increase with temper-ature, suggesting that PFg ions are in the short correlation time limit. The reduced sample S 3 has a nearly constant Ti over much of the temperature range. This plateau in Ti vs. temperature suggest that there may be additional relaxation processes. Since the PFg ions of samples Si and S 2 are in the fast motion regime and the 3 1 P Ti relaxation is due to the rotational motions of the PFg ion, we are able to obtain r c of PFg ions in S[ and S'2 samples at room temperatures through Equation 4.20. Results are shown in Table 4.3. We have also measured T 2 of both 1 9 F and 3 1 P in the film samples. Results are shown Chapter 4. Rotational Motions of PF6 Ions in the PPy Films 31 Table 4.2: 3 1 P Ti relaxation times as a function of temperature of PPy sample S[-S'3. Temperature (°c) 3 1 P Ti of S[ (s) 3 1 P Ti of S2 (s) 3 1 P T i o f S 3 ( s ) -30 2.02 1.80 ' 1.32 -15 (-20 for S2) 2.22 1.78 1.03 0 2.70 2.27 1.72 22 3.08 2.96 1.72 40 3.36 3.37 . 1.72 Table 4.3: Ti relaxation times and the rotational correlation times T c of two PPy samples: S[ and S2 at room temperature. Sample 3 1 P Ti (s) Tc (PS) s[ . s2 3.08 2.96 10.76 11.19 in Table 4.4. 4.4 D i s c u s s i o n The rotational correlation times of PFg ions observed in the as-grown sample S[ (10.76 ps) and the less oxidized sample S'2 (11.19 ps) are not very different from the rotational correlation time of ions found in the solvent (8.08 ps). This suggests that ions are contained in liquid pockets containing solvent molecules instead of at stable positions within the polymer matrix. Therefore, we believe the dopant ions reside in pores of the amorphous region of polypyrrole. The reduction of r c in the reduced sample S3 is not so straight-forward as multiple relaxation mechanisms are possibly involved. However we suspect the rotational tumbling is further slowed down such that in this case the ions are on the slow motion side of the T l minimum (See Figure 4.2). The fact that each relaxation measurement is well fit by a single Ti relaxation constant indicates that either the motion of ions in the entire PPy sample is homogenous, or that Table 4.4: T2 relaxation times of three PPy samples: S[-S'3. Sample 3 1 P T 2 (ms) 1 9 F T 2 (ms) s[ 1.88 0.79 S2 1.78 0.58 s'3 0.19 0.10 Chapter 4. Rotational Motions of PF6 Ions in the PPy Films 32 PFg ions located in different pools communicate rapidly compared to T\. It is worth mentioning that the rotational correlation time of molecules traditionally can be derived from the Debye equation: T < - T ¥ V ( 4 - 2 1 ) with known solvent viscosity r]vis and estimated Stokes radius rs. The calculated r c of PFg ions in solvent from Debye equation at room temperature is 7.71 ps by using rs .= 2.08 A. (1.579+0.5, bond length of P-F plus the atomic radius of F) [25] which is in agreement with the r c (8.08 ps) determined from the relaxation study. However the conditions of applying Debye equation to calculate rc of ions in films are not well defined, therefore the use of Debye equation is not suitable for the studying rotational correlation times here. Discussions regarding T2 measurements are integrated into the discussion of Chapter 6. 33 Chapter 5 Solvent Accessibility of PF^~ Ions in the PPy Films 5 . 1 N u c l e a r O v e r h a u s e r E f f e c t ( N O E ) Nuclear spin cross-relaxation in liquids is caused by mutual spin flips (transitions be-tween spin states) of dipolar-coupled spins which are induced by motional processes. Cross-relaxation leads to a transfer of magnetization between the spins, known as the nuclear Overhauser effect(NOE). The NOE is extremely important in the determination of molecular structure because it provides constraints on the distances between various nuclei. These constraints allow the complete determination of the 3D structure of large biomolecules. From the above relaxation section, we've already arrived at the equation that describes how I\z of the two-spin system evolves with time. The most striking manifestation of the Overhauser effect is the steady state value of (I\z) when the spin I2z is saturated, so that (I2z) = 0. To satisfy the steady state condition, we have = 0, then Equation 4.17 becomes: d ^ I u ) = 0 = -Pl((Ilz) - (Ilz)eq) + pCToss(hz)eq (5-1) dt (hz) = (hz)eq + —(hzU (5-2) Pi I T \ l i e rcross (hz) eq I where rj is the NOE enhancement factor. Due to the cross relaxation, the signal of spin 1 is enhanced by a factor of 1 + r\ when saturating spin 2. The NOE enhancement factor depends on the correlation time, hence the motion of the molecules. In the limit of very rapid molecular motion (small r c), becomes: ( / l 2 ) = ( / l z ) e 9 { l + ^ } = ( / l , ) e g [ l + ^ } (5-3) with the equilibrium magnetization being proportional to the gyromagnetic ratios. Here we see a rather disappointing result that the NOE enhancement factor rj is predicted to Chapter 5. Solvent Accessibility of PF6 Ions in the PPy Films 34 be independent of r, however this applies only for the completely artificial circumstance of two dipoles totally isolated from all other sources of relaxation. In a reality, two- spin system may show distance and correlation time dependent NOE enhancement: ncc^ ' (5.4) when other competing sources of relaxation are present. This is the basis for using the NOE to compare internuclear distances in different molecules. Cross-relaxation and nuclear Overhauser effects can be studied by ID Nuclear Over-hauser Difference Spectroscopy where NOE enhancement is obtained by selectively sat-urating one spin while observing the intensity changes in other spins. 5 . 2 E x p e r i m e n t a l M e t h o d 5.2.1 ID NOE Difference Spectroscopy In measuring the NOE enhancement between solvent proton a H and 3 1 P in PFg ions, we saturate the solvent 1 H resonance with a long train of weak 7r/2 pulses and then compare the intensity of 3 1 P with its equilibrium value. This can be simply done by integrating the 3 1 P peak with and without 1 H radiation, and then taking the difference. The actual pulse sequences used are (See Figure 5.1): (1) [ f ( :H) - r ^ — " _ T 2 - f ( 3 1P) - AQ. (2) [f OH) — T l j P ° w e r - o n — r 2 — f ( 3 1P) — AQ where n = 50 is the number of repeats and T\ — 2 ms, r 2 = 248 ms, which are arbitrarily chosen as long as they are <C Ti of : H . Note that the pulses tuned to the solvent protons are weak enough to significantly affect only the *H in PC but not any protons attached to a more rigid polymer backbone. Sequence (1) is equivalent to a simple | pulse on 3 1 P , and Sequence (2) is the same as (1) except for the additional lH radiation. The NOE difference can be found by subtracting the perturbed and unperturbed spectra. Identical experimental conditions are ensured in both sequences in order to obtain the correct difference intensity. 5 . 3 R e s u l t s ID NOE difference experiments are done on three wet PPy samples S'rS'3 at E == 0.53 V, 0.25 V, -0.2 V vs Ag/AgCl and on two PC samples containing 0.05 M TEAP, 2 M T E A P respectively in order to perform a comparative study among them. : H radiation on protons of PC (C4H 603) is expected to affect the population distribution of 1 9 F of PFg ions if they are close enough in space to relax each other. The fluctuating fields from the 1 9 F will then in turn relax their neighboring 3 1 P. If the relaxation processes involve relevant transition pathways which are correlation time r c and the distance r dependent, a positive 3 1 P signal increase will be observed due to J H saturation. In order to compare the NOE results unambiguously, two solvent samples are used as control experiments, Chapter 5. Solvent Accessibility of P F 6 Ions in the PPy Films 35 r P i / 2 Figure 5.1: Pulse sequences for steady state NOE difference spectroscopy. Table 5.1: NOE enhancement factor 77 of the five samples. Sample NOE 77 0.05M PC 2M PC ' s[ S'2 s'3 0.107 ±0.02 0.163 ±0.004 0.163 ±0.006 0.17 ±0.02 -0.035 ±0 .03 and identical experimental conditions are applied. The NOE enhancement factor 77 of 3 1 P is calculated from Equation 5.5. (Perturbed — Unperturbed)31 P peak intensity . . 01 (^ D.Oj Unperturbed P peak intensity Results from the 5 samples are shown in Table 5.1 and Figure 5.2. 5.4 D i s c u s s i o n Samples S[ and S'2 exhibit similar 77 values to those of the solvent samples, implying that ions in the oxidized films are surrounded by solvent molecules in an environment similar to the free solvent. The much smaller (consistent with zero) NOE enhancement Chapter 5. Solvent Accessibility of P F 6 Ions in the PPy Films 36 of the reduced sample S'3 seems to indicate that PFg ions in the reduced sample are more tightly coupled to the PPy backbone, and are not contact with the solvent.(see Figure 5.2). This could suggest that between the oxidation state 0.25V7 and —0.2V, the films may undergo some significant structural changes, causing ions in S 3 to reside in a much less solvated environment than ions in Si and S 2 . 0.05M PC Figure 5.2: NOE enhancement factor rj of the five samples. 37 Chapter 6 Translational Motions of PF^~ Ions in the PPy Films 6 . 1 P u l s e d - F i e l d G r a d i e n t ( P F G ) N M R The use of pulsed-field gradient N M R to study translational self-diffusion originating from Brownian motions was first suggested by McCall, Douglass and Anderson[26] and first demonstrated by Stejskal and Tanner in 1965[27]. The Stejskal and Tanner pulse sequence (pulsed gradient spin echo sequence) involves a spin echo plus two magnetic field gradient pulses:| - T(PFG) - IT - T(PFG) (Figure 6.1). By imposing two gradient pulses separated by A with amplitude g and duration S in a spin echo sequence, spins are labeled spatially as their local Larmor frequency becomes u>(r) = jB + jg • f. In the absence of a gradient, there is no loss of phase coherence in the NMR signal following the | pulse. At the time of the first gradient pulse, there is a position-dependent phase shift, which may be said to record the spin's position along g. Following the gradient pulse, this phase shift persists until it is inverted by the IT pulse. A second phase change occurs at the time of the second gradient pulse. For nuclei that have not moved between the gradient pulses (during A , the interpulse separation), the phase shift of the second gradient pulse completely reverses the phase shift of the first gradient pulse, and the magnetization refocuses perfectly. Net motion results in incomplete refocusing at the end of the spin echo sequence, leading to an attenuated echo signal. Figure 6.1: Pulse field gradient spin-echo sequence (PFGSE). Chapter 6. Translational Motions of PF6 Ions in the PPy Films 38 Here we'll show how one can correlate the echo attenuation with diffusion in the case of free diffusion using the macroscopic approach, starting from the Bloch equations (Equation 1.13) containing diffusion terms: ^ = 7 M x B - f i ( M - M „ ) + flv2M (6.1) where D is the spin diffusion coefficient and the B field now is ~$ = (0,0,5 + ^ ( i ) ~(?(t) is 0 except during, the gradient pulses. If we define the transverse magnetization as m = Mx + iMy, and plug it into Equation 6.1, we obtain: dm m 2 — = -iLU0m, - xj(g • r). - — + DV m, (6.2) Ot ±2 In the absence of diffusion (D = 0), m relaxes exponentially with a time constant T 2 , and may be expressed as m = ipe-iu>ot-t/T2 ^ w n e r e ^ j s t h e precessing magnetization unaffected by the transverse relaxation. If we substitute rn, into Equation 6.2, we have: dip ' — = -iy(g • r)lp (6.3) The simplified equation can be solved in two parts. Between the | pulse (/; = 0) and the 7r pulse (t = r) ip is given by ip = A exp(—ijf • F) (6.4) where 7ty)= f g(i)dt (6.5) Jo A corresponds to the value of I/J immediately following the | pulse. During the period from the f to 7r pulse each spin acquires a position dependent phase — jf • F. Because the effect of the TX pulse is to set back the phase of ip by twice the amount that it has advanced, in the period following the TT pulse we have: ip = Aexpi-i-yr • (F - 2/)) (6.6) where / = F(r) . We may use a single expression to represent the behavior of ip from the | pulse to the echo: ip{r,t) = Aexp{-iyr-(F-2H(t-T)f)) (6.7) where H is the unit step function H(t — r) = 0 when t < T, H(t — r) = l when t > r. We can now consider the solution of 'ip in the presence of diffusion by allowing A to be a function of t. By substituting tp(r, t) into: ^ = -ij(g • r)iP + DV2vb (6.8) Chapter 6. Translational Motions of P F 6 Ions in the PPy Films 39 We find. d A ® =_J2D[F_2H(t-T)f]2A(t) ' (6.9) dt Integrating Equation 6.9 from t = 0 to t = 2T, the echo attenuation E — from the Stejskal and Tanner pulse sequence can be obtained. ln[^ ] = ln[E}= j*—?DF*dt + jT'-j2D[F - 2jfdt = - 7 D F2dt-±f (T Fdt + Afr) (6.10) 0 The calculation of the echo attenuation (Equation 6.10) is quite straightforward but rather tedious. For example, using the time dependence of g(t) in Figure 6.1, the integral of F(t) from t '= 0 to t = 2r is calculated as follows, F(t) = / Odt-r gdt+ 0dt+ . gdt J0 Jti Jti+S Jti+A = g(t + 6 - ti - A) (6.11) And according to the definition / = F(r) , f is / = F(T)= fg(t)dt Jo = / 0dt+ gdt+ 0dt = g5 (6.12) Jo Jti Jti+5 After performing the symbolic algebra analysis, we obtain the relationship between signal attenuation and the diffusion coefficient for free diffusion (Equation 6.13). l n [ ^ ] = H E { 9 > A ) ] = - ^ W ( A - | ) . (6.13) The term | accounts for the finite width of the gradient pulse. Rearranging Equation 6.13, we find: A(2T) = A(0) exp[- 7V JD^ 2(A - |)] (6.14) From the Bloch equations including the diffusion of magnetization, we can derive the necessary relationship analytically as shown in Equation 6.13. However, in the case of restricted diffusion the macroscopic approach becomes impossible to formulate, and a ' propagator formalism may be used to describe the result of the P F G N M R experiment. Chapter 6. Translational Motions of P F 6 Ions in the PPy Films 40 Time scale Free diffusion Restricted diffusion Figure 6.2: Features of molecule motions when undergoing free diffusion and restricted diffusion. 6.1.1 Self-Diffusion in Restricted Geometries In the PFG experiment, we probe the particle's motion by taking a measurement at time t = ti and a second measurement at time t = t\ + A . The key point is that in the P F G experiment the echo attenuation gives information on the displacement along the gradient axis that has occurred during the period A, and this displacement can be related to the diffusion coefficient. If a particle is diffusing within a restricted environment, the displacement along the z-axis will be a function of A , the diffusion coefficient, and the size and shape of the restricting geometry. Consequently, if the boundary effects are not taken into account and we analyze the data using the equation for free diffusion, we will get an apparent diffusion coefficient Dapp, that is in general smaller- than the actual diffusion coefficient. Considering molecules diffusing inside a single pore of diameter a (characteristic length of the geometry), the observed motion will be strongly dependent on the time scale (see Figure 6.2). For short diffusion times A <§C a2jD the diffusion will appear to be unrestricted because molecules have not diffused far enough to feel the effect of boundary. The diffusion coefficient determined will be the same as that observed for freely diffusing molecules. For long diffusion times A 3> a2/D, because maximum mean displacement is limited by the pore diameter a, the root mean displacement ((n — r 0) 2} becomes time-independent (where we have defined r 0 as the starting position of a given molecule, and r\ its position a time A later), therefore the diffusion coefficient will be proportional to 1/A. At intermediate timescales (A « a2j'D) the behavior may be quite difficult to characterize. Because at this time scale, Dart of the molecule's diffusion is Chapter 6. Translational Motions of PF6 Ions in the PPy Films 41 (a) Free diffusion (b) Restricted diffusion (c) Restricted diffusion (isolated pores) • (connected pores) A A • A Figure 6.3: Expected behaviors of D vs. A in free diffusion and restricted diffusions. restricted by the barrier and the mean squared displacement will not scale linearly with A , thus the measured diffusion coefficient (Dapp) will appear to be time dependent. For diffusion in connected pores, one can imagine the long time behavior will be different from that in isolated diffusion because the maximum mean displacement is no longer limited by the pore diameter. Therefore the observed diffusion coefficient at long times will not decrease as 1/A. Figure 6.3 shows the expected behaviors of D vs. A in free diffusion, diffusion in isolated pores, and diffusion in connected pores. To understand restricted diffusion' in P F G experiments, it is convenient to use a propagator formalism approach. One needs to consider the quantity P(ro | f i , A) , the conditional probability of a molecule moving from position r 0 to rx in the time A . P ( r 0 | r l 5 A) is often termed the diffusion propagator and it also obeys Fick's second law of diffusion [28], a p ( r „ j n , A ) = D v 2 f , ( r n | r i | A ) ( 6 l 5 ) In the case of unrestricted Brownian motion, the solution of P(r0 \ ri, A) is a Gaussian distribution over molecular displacement, with a distribution width which expands with time: P(r0 | r i > A) = ( 4 7 r D i ) - 3 / 2 e x p ( - ( n ' j f ^ j (6.16) The echo attenuation measured by the P F G experiment is just the Fourier transform of ) Chapter 6. Translational Motions of PF6 Ions in the PPy Films 42 the diffusion propagator: E(g, A) = f [ p(r0,0)P(r0 | ru A) x exp (ijSg • (n - r 0)) dr0drx (6.17) where p(r0, 0) is the starting molecular density (spin density). If we substitute Equation 6.16 into Equation 6.17, we retain the expected echo attenuation equation as in Equation In the case of restricted diffusion, the solution of Equation 6.15 which is the diffu-sion propagator will deviate from its classic Gaussian form because different boundary conditions are imposed. We have to substitute the relevant p(ro,0) and P(r0 \ ri, A) in order to evaluate the corresponding echo attenuation equation. By understanding the the influence of restrictions on the form of the propagator, and hence the echo atten-uation function in a P F G experiment, one can extract information not only about the motions of the molecules but also the geometry of the boundaries and thus about the pore morphology of the surrounding medium. However P(ro | r\, A) tends to be difficult to obtain for nontrivial geometries. 6 . 2 E x p e r i m e n t a l M e t h o d 6.2.1 Pulsed-Field Gradient N M R Measurement A l l PFG NMR experiments were performed on a homebuilt N M R spectrometer working at 8.4 T with a diffusion probe tuned for fluorine 1 9 F at 342 MHz. Measurements of signal intensity A(2T) are made with several values of gradient strength g, then fit to Equation 6.14 for given A , 7, 5 to extract the diffusion coefficient D. The diffusion coefficient can then be obtained as a function of diffusion time A by performing the same measurement at different A's. Gradient C o i l Design The gradient coil design of the gradient probe used here is inspired by Cory's design [29] which was developed for solid state studies where diffusion coefficients are small. A Maxwell pair of coils is generally employed for producing gradients along the z direction in superconducting magnets because it is simple to construct, and good gradient profiles can be obtained with very few turns. A Maxwell pair consists of a set of circular loops with currents of equal amplitudes but opposite directions. This configuration generates magnetic field gradients symmetric about the mid-plane and about the central axis. The gradient profile and the gradient coil constant (defined as the gradient strength per unit current: T/m-A) in principle can be calculated via the Biot-Savart Law. In order to obtain larger coil constant and better gradient linearity, a gradient set is constructed from two building blocks. Figure 6.4 shows the configuration where 12 copper wires are wound on the inner side that is close to the sample, and 84 wires are wound on the outer side. The predicted coil constant is 0.27 T/m-A. 6.13. Chapter 6. Translational Motions of P F 6 Ions in the PPy Films 43 3.5 cm 0.5 cm 1 " C Z J D " Figure 6.4: A gradient coil set designed for higher coil constant and better gradient profile. The winding direction of upper two sets is opposite with the lower two sets. The total dimension is 3.5 cm x 4.5 cm. Gradient Calibration The gradient strength needs to be carefully calibrated to obtain an accurate diffusion coefficient D from a fit to Equation 6.13. Two methods were used on the home-built diffusion probe: (1) back-calculation of the coil constant from a diffusion experiment on H 2 0 using D = 2.3 x 1CT5 cm 2/s for pure H 2 0 at 25°C [30]; (2) a spin-echo experiment with PFG turned on between | and ir pulse during the acquisition, which gives a ID image of the sample and the effective gradient across it. The back-calculation method relies on the correct D value of the standard compound. A l l recent results in the literature are consistent with the value of D = 2.3 x 10~5 cm 2/s for pure H 2 0 at 25°C, however since the N M R probe temperature inside the magnet may be as low as 18 — 21°C, we estimate D to be 2.03 x I O - 5 cm 2/s at 20°C according to [31]. The 1 H signal intensity of H 2 0 is obtained as a function of unknown gradient strength g (g = gmax% where x is gradient amplitude from 0 to 1), and is fit to Equation 6.14 for g m a x , given that A,j,5, and D (2.03 x 10~5 cm2/s) are known. (See Figure 6.5) The instrumental maximum gradient strength using the back calculation method is determined to be 480 G/cm. A cylinder sample filled with pure water was also used for ID image calibration. Fourier transforms return images of the sample cross-sections with respect to the gradient direction. For a cylinder sample with its long axis perpendicular to the gradient direction (z), we obtained a half circular image. (See Figure 6.6) The width of the spectrum Au is related to the sample radius r and the gradient strength g via: Au = ^ (6.18) Table 6.1 shows the results of Au at different gradient amplitude x (0-1) and the cor-Chapter 6. Translational Motions of P F 6 Ions in the PPy Films 44 1 r — i i i i I 1 ' fit - -\ \ \ \ -\ \ -\ \ \ \ -\ \ \ \ \ -\ N N. -^ _ 1 1 1 1 1 -++ 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 Gradient amplitude x Figure 6.5: Signal attenuation of 1 H intensity in H 2 O for gradient strength calibration. responding g values calculated from Equation 6.18. From the fit of the data g(x), we obtained the gradient strength as g = 485.53 • x + 1.22 (G/cm), yielding a maximum g of 488 G/cm. We note that the maximum gradient strength based on the predicted coil constant of 27 G/cm-A and the maximum gradient current supplied 15.2 Amp (<?max = 27 x 15.2 = 410.2 G/cm) is in reasonable agreement with both of the empirical 5max values. Since the ID image calibration is the more direct method, along with the fact that we don't have a firm empirical value for D of water at the sample temperature, we will use the ID calibration result. Pulse Sequences for Diffusion Measurements The first and the most simple pulse sequence for measuring diffusion is the Stejskal and Tanner sequence, also known as the pulsed field gradient spin-echo (PFGSE) sequence: ^ — T(PFG)—7T — T(PFG) (See Figure 6.1). There have been many other pulse sequences developed since then to specifically meet different experimental needs, for example, to eliminate eddy current effects, and to deal with imperfection of gradient pulses etc. One commonly employed pulse sequence is the stimulated echo pulse sequence (PFG-STE): f - TI(PFG) - \ - r 2 - § - n(PFG) - AQ (see Figure 6.7(a)). This sequence was designed to avoid T 2 relaxation effects by storing the magnetization along the z axis during r 2 , so that the relaxation depends primarily on Ti , which is usually much longer than X 2 for macromolecules. Longer diffusion times (A) can then be used. Chapter 6. Translational Motions of PF6 Ions in the PPy Films 45 1 Figure 6.6: One-D image gradient calibration. A longitudinal encode-decode (LED) pulse sequence was proposed by Gibbs and John-son [32] to ease eddy current effects by delaying the acquisition until the effects have dis-sipated. Eddy currents are inevitably generated in the surrounding conducting surfaces around the gradient coils due to the rapid rise of the gradient pulses, causing amplitude or phase distortion of the echo. These effects often lead to a spurious extra echo attenuation which can be easily misinterpreted as a higher diffusion coefficient. Many such sequences are based on STE, as in the following: § -Tl(PFG)-§ - r 2 - f -rx(PFG)-\-Te-\-AQ where Te is the settling time.(See Figure 6.7(b)) However the LED P F G sequence does not solve the problem of the eddy current tail from the first gradient pulse extending to the second transverse evolution period, therefore another method involving a self-compensating bipolar pulse pair (BPP) is employed [33, 34]. Each regular gradient pulse is replaced by a pair of pulses having different polarities and separated by IT pulses. The eddy currents caused by two gradient pulses with identical areas but different polarities will tend to cancel each other. Finally, the pulse sequence used in this work to determine the diffusion coefficient is a stimulated echo LED sequence incorporating bipolar gradient pulses (BPP-LED pulse sequence) shown in Figure 6.7(c). This sequence was tested extensively on H 2 0 and a PC sample containing 0.05 M T E A P to ensure that reliable D at relevant A time scales is obtained (where the eddy current effects are not present). Similar to the simple PFGSE sequence, the signal attenuation  Aj?f.) due to diffusion is given by the following equation: A(2r) = ,4(0) exp [ - 7 2 </ -W(A Tg/2)\ (6.19) Chapter 6. Translational Motions of PF6 Ions in the PPy Films 46 Table 6.1: Spectra width Ais from one-D image gradient calibration at gradient amplitude x (0-1) and the corresponding gradient strength. Gradient amplitude x (kHz) Gradient strength g (G/cm) 0.1 84 49.31 0.2 168 98.63 0.3 252 147.94 0.4 334 196.10 0.5 418 245.42 0.6 498 292.44 0.7 580 340.60 0.8 664 389.91 0.9 748 439.23 1.0 832 488.54 where Tg is the time separation between a bipolar pair. A 16 step phase cycle for the BPP-LED pulse sequence is adopted from Johnson et al. [35]. Chapter 6. Translational Motions of PF6 Ions in the PPy Films 47 (a) pi/2 pi/2 pi/2 Tl 5 i A (b) pi/2 pi/2 pi/2 pi/2 pi/2 Ti n T e r - i 1 1 pi/2 pi pi/2 pi/2 pi pi/2 pi/2 T e 5/2 H H T 1 5/2 H H Tg 1 V • 1 1 •4 : • T -4 ; • A Figure 6.7: (a) Stimulated echo(PFGSTE) pulse sequence. (b)A longitudinal eddy-current delay(LED) pulse sequence, (c) A bipolar pulse pair(BPP) LED pulse sequence. Chapter 6. Translational Motions of PFg Ions in the PPy Films 48 1 1 1 1 1 1 1 * ~ — «~ fit Kg) * s % \ \ \ \ V \ X % X X . X X • I 1 1 1 1 ~ 1 1 I 1 I I I l _ I I I I I 0 50 100 150 200 250 300 350 400 450 500 Gradient strength g (G/cm) Figure 6.8: Echo attenuation with respect to gradient strength g of 0.05M PC sample. 6 . 3 R e s u l t s The diffusion coefficient D of two wet PPy samples S[ and S'2 at E = 0.53 V, 0.25 V vs Ag/AgCl and two PC samples containing 0.05 M , 2 M T E A P are measured by a stimulated echo LED sequence incorporating bipolar gradient pulses. Solvent samples again were used for control experiments. We optimized the three experimental parameters S, Tg, Te under the constraint of Ti , T2 relaxation, and then kept them fixed in all diffusion measurements. The numbers are '5 = 300 ps, Tg = 100 /is, Te = 100 ms. Figure 6.8 shows the P F G signal intensities 1(g) at A = 50 ms of sample 0.05 M PC, and the diffusion coefficient of PFg is determined through fitting to Equation 6.19. The same fitting process was repeated for all samples at different diffusion times (A's). Self-diffusion coefficients of the 0.05 M and 2 M PC samples measured over different observation times (A) are shown in Table 6.2 and plotted in Figure 6.9 (a). Table 6.2: Self-diffusion coefficients of solvent samples. Sample D(A) (cm2/s) 0.05 PC D(50 ms) D(100 ms) D(500 ms) 3.50 x 10~6 ± 1.7% 3.43 x 10~6 ± 1.4% 3.48 x 10~6 ± 2,1% 2M PC £>(50 ms) L>(500 ms) L>(1000 ms) 1.99 x 1 0 " 6 ± 3.6% 1.59 x 10~6 ± 2.5% 1.73 x 10"6 ± 1.1% Diffusion coefficients were measured from the most oxidized sample S[ (see Table 6.3, Chapter 6. Translational Motions of P F 6 Ions in the PPy Films 49 (a) xlO 0.2 0.4 0.6 0.8 1 1.2 Diffusion time A (s) (b) xlO 0.5 1 1.5 2 2.5 Diffusion time A (s) Figure 6.9: (a) D of 0.05M,2M PC samples and (b)D of the oxidized PPy sample S[ over different observation times A . NMR tube Polypyrrole films Horizontal orientation Vertical orientation Figure 6.10: The sheets of PPy thin films were stacked either parallel or perpendicular to B field. Figure 6.9 (b)) with films oriented in two directions, once with film's normal perpendicular to the static B field axis (S[ V) and once perpendicular (^H) , so that diffusion in both planes can be sampled. (See Figure 6.10) For sample S'2, partly because there are fewer ions in the less oxidized films (smaller signal to noise ratio), data points suffer from intrinsic noise, yielding Dapp with large errors. The Dapp obtained under these circumstance may lack of credibility; nevertheless experimental values are given below: Dapp=(3.51 x 10 - 9 ± 77%) cm 2/s at A = l s, and (6.76 x 10~9 ± 63%) cm 2/s at A=1.5 s. Diffusion measurements on sample S3 are not feasible due to its short T 2 (97ps), signals basically decay away before the acquisition. 6 . 4 D i s c u s s i o n When we use the P F G method to measure diffusion in free solution, the length of the observation time we choose (A) should not affect the value of the diffusion coefficient Chapter 6. Translational Motions of PF6 Ions in the PPy Films 50 Table 6.3: Diffusion coefficients of the oxidized S[ and the less oxidized S'2 PPy sample. Films are oriented parallel (H) and perpendicular (V) to the magnetic field. Sample D(A) (cm2/s) 7J(50 ms) D(100 ms) D(1000 ms) 2.57 x lO" 8 ±29.9% 1.24 x 10"8 ± 32.1% 4.47 x 10" 9 ± 14.2% £>(1500 ms) 4.66 x lO" 9 ± 12.3% D(50 ms) £>(500 ms) D(1000ms) 4.57 x lO" 8 ±47.8% 1.04 x 10"8 ± 21.1% 4.98 x 10" 9 ± 25.4% D(1500 ms) D(2000 ms) D(3000 ms) 4.31 x 10"9 ± 24.8% 4.76 x 10"9 ± 27.5% 6.35 x 10" 9 ± 11.2% 5 2 V £>(1000 ms) D(1500 ms) 3.51 x l O " 9 ± 77% 6.76 x 10~9 ± 63% extracted from the measurement. This is what we observed in the solvent samples; the diffusion coefficients of both samples are essentially the same from A=50 ms to 1 s. From the solvent measurements, we. are confident that the D values beginning from 50 ms are not affected by the disturbance of the eddy current effect. Measurements made at A < 50 ms may bear the influence of the eddy current effect (even though the pulse sequence used is devoted to overcome this problem), and also suffer from the fact that very little echo attenuation is produced, if the diffusion coefficient is small, according to Equation 6.19. This results in a less accurate determination of D values. Figure 6.11 shows the signal attenuation data from S1 oxidized PPy sample at A=50 ms and 3 s. One can see that due to the slow translational motions of PFg ions in the films and short diffusion time A , only little attenuation is achieved. We can however try implementing stronger gradients, thereby increasing the degree of signal attenuation. The fact that D<xpp over a range of A's is not distinctly different between the two oriented Sl samples, suggests that the films are homogeneous in morphology, at least on the length scales probed here (~ 1.4/im). If there were large oriented crystalline regions inside the films, then we should have observed dissimilar behaviors in terms of diffusion coefficients as a function of A . The Dapp values measured from PFG experiments can be compared to. the D values from the numerical fitting of impedance spectroscopy data. Impedance analysis based on diffusion models yields values of diffusion coefficients as large as 10~7 cm 2/s for CIO4 ions in P P y / P C [36]and 10~8 cm 2/s for PFg ions in P P y / P C [37]which are believed rather high for diffusion of ionic species in solids. These high values have been explained in terms of solvent ingress, effectively assisting the ionic diffusion progress [38, 39]. Our D value for PFg in PPy films (~ 5 x 10~9 cm2/s)' is two order of magnitude slower than that typically obtained from the impedance data, and we believe reflects ionic motions Chapter 6. Translational Motions of PF6 Ions in the PPy Films 51 (a) 2.1 3 i=r 1.9 (3 S, Oxidized PPy A =50ms fit 1(g) 0 50 100 150 200 250 300 350 400 450 500 0 50 100 150 200 250 300 350 400 450 500 Gradient strength g (G/cm) Gradient strength g (G/cm) Figure 6.11: Signal attenuation 1(g) from S[ oxidized PPy sample at A=50 ms and 3s. in the polymer/electrolyte matrix in steady state. As seen in Figure 6.9 (b) the A dependence of Dapp of 5 : films appears to be an indication of restricted diffusion. The diffusion coefficient begins to drop quickly from 50 ms, implying that at least some of ions are already "feeling" the boundaries. We can conjecture an upper limit of pore size a by plugging D(50 ms) and A== 50 ms into D = The upper limit of the characteristic length scale (pore size) a in S[ is determined to be ~ 0.4 pm. At longer times, diffusing molecules probe the connectivity of the pore space. The diffusion coefficient reaches a plateau at 2-3 s instead of decreasing with inverse time, suggesting that the pores are interconnected rather than being isolated voids inside the polymer. The data from sample S'2 provides a higher bound of the diffusion coefficient in that oxidation state (E = 0.25 V vs Ag/AgCl), D must be < 1 0 - 8 cm 2/s. Because if the diffusion coefficient of diffusion for the ions was greater, we would have observed a greater echo attenuation. The short 1 9 F T 2 of the reduced S3 sample in Table 4.4 that impeded us from mea-suring the diffusion coefficient actually serves as an index as far as how 1 9 F spins are being relaxed. Generally speaking the stronger the spin coupling to the environment, the shorter the T 2 (this also reflects on the peak linewidth). We have found that the T 2 values of 1 9 F and 3 1 P in the reduced film sample S'3 are much shorter than in the two oxidized samples (S[ and S'2). This dramatic difference in T2 would be consistent with a structure in the reduced state in which the dopant ions interact more strongly, with the polymer backbone, whereas at high doping levels (e.g. the oxidized state) the interactions with the solvent dominate. This picture is consistent with the results from the Ti relaxation study and ID NOE difference experiments that suggest the existence, of an abrupt structural change between oxidation states 0.25V and — 0.2V vs Ag/AgCl . Chapter 6. Translational Motions of PF6 Ions in the PPy Films 52 The same conclusion has also been suggested from an impedance study on PFg doped PPy in various states of oxidation as the impedance characteristics of the oxidized and reduced PPy are distinctively different [40]. Chapter 7 Conclusion 53 We have demonstrated that NMR spectroscopy is a versatile tool for the study of dopant ions in polymer systems. Advantages of NMR over other techniques lie in its selectivity, which is employed to determine the ion content, and its ability to detect motions in different time scales, which enable us to probe the rotational and translational motions of PFg inside the PPy films. By investigating these quantities systematically as a function of oxidization state, it is possible to relate the microscopic structural changes to macroscopic material behaviors. The ion content decreases linearly with decreasing electrochemical potential (E), di-rectly confirming the ion insertion mechanism of polymer actuation. With known ion con-tent, sample mass and deposition current, the pyrrole/ion and pyrrole/solvent ratios can be determined. The obtained doping level for the as-grown (Si) sample (Ni/Np = 1/3.83) is consistent with the literature [20, 21, 22]. This suggests that all the ions are detected, and examining ion content via ID 3 1 P NMR spectroscopy seems to be a reliable and straightforward method. PFg doped polypyrrole, as fabricated using the described electrodeposition method, are found to be porous in nature. 3 1 P relaxation and ID NOE difference experiments showed that PFg ions in the oxidized films (S[) and less oxidized films (S2) are located inside pores containing propylene carbonate solvent molecules. As the PPy films are reduced, PFg ions that remained inside the films became less splvated, -possibly because the pores are literally smaller. However we believe that there are still many solvent molecules inside the reduced PPy films according to the solvent molecule to pyrrole ratio (Nsoi/Np) in Table 3.1. Those results were based on dry samples, and are likely underestimates of the true solvent content of wet films. P F G NMR experiments provide direct measurements regarding the self-diffusion co-efficient of PFg ions inside polypyrrole films. The D values of PFg in PPy films from the P F G study (~ 10~9cm2/s) appear to reflect ionic motions in the polymer/electrolyte matrix in steady state without any external perturbations, and can be used in PPy actu-ator modeling, and thus in predicting its behaviour. The time-scale dependence of Dapp of the oxidized films (5X) at short times had let us conjecture an upper limit of pore size a of ~ 0.4/-4 m, and the values at longer times enables us to probe the connectivity of the pore space. Data suggests that the pores are interconnected rather than being isolated voids inside the polymer. Suggested experiments for the future include: (1) further investigate the Ti vs. tem-perature behavior of the reduced PPy sample as other relaxation mechanisms (rather than ionic rotations) may provide information on the polymeric environment that ions Chapter 7. Conclusion 54 experience in the reduced state. (2) improve the instrumental limits of the PFG exper-iments, for example, to increase the gradient strength (reduce the gradient duration) to overcome the short 1 9 F T 2 problem so that extensive P F G experiments can be performed in the less oxidized and reduced films as well. A shielded magnetic gradient coil system in the N M R probe head would prevent the gradient pulse from affecting the region outside the gradient coils, thereby reducing the generation of eddy currents. 55 References [1] T.A. Skotheim, R.L. Elsenbaumer, and J.R. Reynolds eds. The Handbook of Con-ducting Polymers. Mercel-Dekker, New York, second edition, 1977. [2] P.I. Lazarev ed. Molecular Electronics. Kluwer Academic Publishers, Boston, 1991. [3] E.J. Smela. Micromech. Microeng., 9:1-18, 1999. [4] J.D.W. Madden, B. Schmid, M . Hechinger, S.R. Lafontaine, P G . A . Madden, F.S. Hover, R. Kimball, and I. Hunter. IEEE Journal of Oceanic Engineering, 24:199-201, 2004. [5] S. Hara, T. Zama, W. Takashima, and K. Kaneto. Synthetic Metals, 149:199-201, 2005. [6] Y . Nogani, J .P Pouget, and T. Ishiguro. Synthetic Metals, 62:257-263, 1994. [7] X . G . L i , L .X . Wang, M.R. Huang, Y.Q. Lu, M.F. Zhu, A. Menner, and J. Springer. Polymer, 42:6095-6105, 2001. [8] O. Cakmak, M . Basturkmen, and D. Kisakurek. Polymer, 45:5451-5458, 2004. [9] F. Devreux, G. Bidan, A.A. Syed, and C.J. Tsintavis. Physique, 46:1595-1601, 1985. [10] B. Wehrle, H.H. Limbach, J. Mortensen, and J. Heinze. Synthetic Metals, 38:293-298, 1990. [11] T.C. Clarke, J.C. Scott, and G.B. Street. IBM J Res Develope, 27:313-320, 1983. [12] O. Inganas and Q. Pei. Solid state ionics, 60:161-166, 1993. [13] J.D.W. Madderi, N.A. Vandesteeg, P A . Anquetil, P G . A . Madden, A. Takshi, R.Z. Pytel, S.R. Lafontaine, P.A. Wieringa, and I.W. Hunter. IEEE Journal of Oceanic Engineering, 29:706-728, 2004. -[14] T.F. Otero and J .M. Sansinera. Bioelectrochem. Bioenerg., 38:411-414, 1995. [15] G .M. Spinks, B. X i , D. Zhou, V. Truong, and G.G. Wallace. Synthetic Metals, 140: 273-280, 2004. [16] R. H. Baughman. Nature, 308:63-65, 2005. References 56 [17] M . Yamaura, T. Hagiwara, and K. Iwata. Synthetic Metals, 26:209-224, 1988. [18] H. Okuzaki, T. Kondo, and T. Kunigi. Polymer, 40:995-1000, 1999. [19] E. Naudin, P. Dabo, D. Guay, and D. Belanger. Synthetic Metals, 132:71-79, 2002. [20] M . Salmon, A.F . Diaz, A . J . Logan, M . Krounby, and J. Bargon. Mol. Crys. Liq. Crys., 83:1297-1308, 1982. [21] H. Zhao, W.E. Price, C O . Too, and G.G. Wallace. J. Membr. Sci., 119:199-212, 1996. [22] M . Yamaura, K. Sato, T. Hagiwara, and K. Iwata. Synthetic Metals, 48:337-354, 1992. [23] P G . A . Madden, J.D.W. Madden, P A . Anquetil, N.A. Vandesteeg, A. Takshi, R.Z. Pytel, and I.W. Hunter. IEEE Journal of Oceanic Engineering, 29:696-705, 2004. [24] A. Abragam. The Principles of Nuclear Magnetism. Oxford University Press, Lon-don, 1961. [25] W. F. Kuhs and G. Mattern. Acta Cryst, B50:741-746, 1994. [26] D. W. McCall, D.C. Douglass, and E.W. Anderson. Ber. Busenges. Phys. Chem., 67:336, 1963. [27] E.O. Stejskal and J.E. Tanner. J. Chem. Phys., 42:228-292, 1965. [28] E.O. Stejskal. J. Chem. Phys., 43:3597-3603, 1965. [29] W. Zhang and D.G. Cory. J. Magn. Reson, 132:144-149, 1998. [30] H. Weingartner. Z. Phys. Chem., 132:129-149, 1982. [31] F. Franks ed. Water - A comprehensive Treatise. Plenum Press, London, 1972. [32] S.J. Gibbs and C.S. Jonson Jr. J. Magn. Reson, 93:395-402, 1991. [33] R.F. Karlicek and I.J. Lowe. J. Magn. Reson, 37:75-91, 1980. [34] R .M. Cotts, M.J.R. Hoch, T. Sun, and J.T. Marker. J. Magn. Reson, 83:252-266, 1989. [35] D. Wu, A. Chen, and C.S. Johnson Jr. J. Magn. Reson A, 115:260-264, 1995. [36] K. West, M.A. Careem, and S. Skaarup. Solid State Ionics, 60:153-159, 1993. [37] J.D. Madden. Doctoral Dissertation. Massachusetts Institute of Technology, Boston, 2000. References 57 [38] B.W. Johnson, D.C. Read, R Christensen, A. Hamnett, and R.D. Armstrong. Jour-nal of Electroanalytical Chemistry, 364:103-109, 1994. [39] J. Mostany and B. R. Scharifker. Synthetic Metals, 87:179-185, 1997. [40] M.R. Warren and J.D. Madden. Synthetic Metals, 156:724-730, 2006. 


Citation Scheme:


Citations by CSL (citeproc-js)

Usage Statistics



Customize your widget with the following options, then copy and paste the code below into the HTML of your page to embed this item in your website.
                            <div id="ubcOpenCollectionsWidgetDisplay">
                            <script id="ubcOpenCollectionsWidget"
                            async >
IIIF logo Our image viewer uses the IIIF 2.0 standard. To load this item in other compatible viewers, use this url:


Related Items