Consec_ABC_SVD_v7

Consec_ABC_SVD_v7 - Irreversible Two-step Consecutive...

Info iconThis preview shows pages 1–4. Sign up to view the full content.

View Full Document Right Arrow Icon

Info iconThis preview has intentionally blurred sections. Sign up to view the full version.

View Full DocumentRight Arrow Icon

Info iconThis preview has intentionally blurred sections. Sign up to view the full version.

View Full DocumentRight Arrow Icon
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: Irreversible Two-step Consecutive Process: Analysis of Noisy UV-Vis Absorbance Data Dr. Kalju Kahn UC Santa Barbara, 2004-2010 This Notebook illustrates how to generate multiwavelength absorbance data containing random noise for the irreversible process A fi k1 B fi k2 C. The lack of isosbestic point in the three-component system is illustrated. First, we generate concentration profiles for A, B, and C for 40 time points with 0.2 second spacing. The concentrations follow consecutive first-order kinetics. In[1]:= Remove @ "Global` * " D twostep = DSolve @ 8 concA @ time D == - k1 concA @ time D , concA @ D == A0, concB @ time D == k1 concA @ time D- k2 concB @ time D , concB @ D == 0, concC @ time D == k2 concB @ time D , concC @ D == < , 8 concA @ time D , concB @ time D , concC @ time D< , time D Flatten; In[3]:= cA = concA @ time D . twostep Simplify cB = concB @ time D . twostep Simplify cC = concC @ time D . twostep Simplify Out[3]= A0 - k1 time Out[4]=- A0 I - k1 time- - k2 time M k1 k1- k2 Out[5]= A0 I k1- - k2 time k1 + I- 1 + - k1 time M k2 M k1- k2 Data Import and Setup We will use a previously generated noisy absorbance data matrix that describes multi-wavelength kinetics merasured from 200 to 700 nm over 8 seconds with 0.2 second time interval. The data was generated with k 1 = 0.5 and k 2 = 1.0. Our goal here is to obtain these rate coefficient values via analysis of the absorbance data. In[6]:= timeval = Table @ i, 8 i, 40 <D 5 N; wavelength = Table @ x, 8 x, 200, 700 <D ; matSignal = Import @ "dataABC.dat", "Table" D ; Dimensions @ matSignal D Out[9]= 8 501, 40 < In[10]:= ListPlot3D @ matSignal, DataRange fi 88 0, 8 < , 8 200, 700 < , 8 0, 1.2 << , AxesLabel fi 8 "Time, s", " , nm", "A" < , Mesh fi False D Out[10]= Action Spectra In[11]:= wavelength = Table @ x, 8 x, 200, 700 <D ; signal1 = Transpose @ 8 wavelength , matSignal @@ All, 1 DD < D ; signal6 = Transpose @ 8 wavelength , matSignal @@ All, 6 DD < D ; signal12 = Transpose @ 8 wavelength , matSignal @@ All, 12 DD < D ; signal18 = Transpose @ 8 wavelength , matSignal @@ All, 18 DD < D ; signal24 = Transpose @ 8 wavelength , matSignal @@ All, 24 DD < D ; In[17]:= Needs @ "PlotLegends` " D ListPlot @8 signal1, signal6, signal12, signal18, signal24 < , PlotRange fi 88 200, 700 < , 8 0, 1.2 << , AxesLabel fi 8 " , nm", "A" < , PlotLegend fi 8 "t = 1.2", "t = 2.4", "t = 3.6", "t = 4.8", "t = 6.0" < , LegendPosition fi 8 0.4, 0.10 < , LegendSize fi 8 0.45, 0.45 < , PlotLabel fi "Action Spectra", LabelStyle fi 8 Medium, FontFamily fi "Helvetica" < , Background fi ColorData @ "Atoms", "He" DD Out[18]= 200 300 400 500 600 700 , nm 0.0 0.2 0.4 0.6 0.8 1.0 1.2 A Action Spectra t = 6.0 t = 4.8 t = 3.6 t = 2.4 t = 1.2 2 Consec_ABC_SVD_v7.nb The action spectra do not show an isosbestic point, suggesting that the interconversion of the reactant to the product involves at least one intermediate. We can use singular value decomposition to find the number of species contributing to the observed signalleast one intermediate....
View Full Document

This note was uploaded on 01/09/2011 for the course CHEM 111 taught by Professor Kahn during the Fall '08 term at UCSB.

Page1 / 15

Consec_ABC_SVD_v7 - Irreversible Two-step Consecutive...

This preview shows document pages 1 - 4. Sign up to view the full document.

View Full Document Right Arrow Icon
Ask a homework question - tutors are online