Simulations and curve fitting

The theory for most systems involving coupled chemical reactions is rather compUcated. Analytical approximations are available only for a limited number of relatively simple processes. Semi-analytical solutions based on infinite series, integral equations, tabulated 

Several commercially available simulation packages can be used to simulate two-dimensional or three-dimensional UME and SECM problems, from relatively simple and 

Faulkner, Electrochemical Methods Fundamentals and Applications, 2nd ed., John Wiley Sons New York, 2001. 

Fundamentals of Electrochemical Analysis, 2nd ed., Harwood PWN New York, 1994. 

Girault, Analytical and Physical Electrochemistry, Marcel Dekker New York, 2004. 

---

Almost all the analysis of cyclic and linear sweep voltammograms has been done through peak currents and peak potentials. Unless digital simulation and curve-fitting by parameter adjustment is carried out, all the information contained in the rest of the wave is ignored this brings problems of accuracy and precision. Besides this, a kinetic model has to be proposed before the results can be analysed. 

Here we assume that [H+] from acidic or basic impurities that may influence the dissociation equilibrium. If the values of Ka, Kf(HAi), Qand Cs are known, we can get a(H+) by solving the quadratic equation. Figure 3.3 shows the simulated titration curves for various values of Kf(HAy)(Ca+Cs). If we compare the experimental titration curve with the simulated ones (curve-fitting method), we can get the value of K HA ) (see Fig. 6.6). The value of Kf(HAy) can also be obtained by the relation ... 

Figure 4.7. Observed and curve-fitted Si spectra of three erystalline forms of the silica polymorph tridymite. A. Room-temperature ordered monoclinic form, showing resolution of nine of the twelve Si sites, fitted to 12 pseudo-Voight lines of equal area with a Gaussian Lorentzian ratio of 0.3. B. Orthorhombic form at 142°C, fitted to six lines broadened by an ineommensurate plane-wave modulation. C. Orthorhombic form at 202°C, fitted to a single line simulated with a non-linear incommensurate modulation. Adapted from Kitchin et al. (1996), by permission of the...

D Rodbard, Y Feldman. Kinetics of two-site immunoradiometric ( sandwich ) assays. I. Mathematical models for simulation, optimization, and curve fitting. Im-munochemistry 15 71, 1978. 

The simulation may be expressed in reaction time (/) rather than cycles (m) by assigning absolute values of ni, 2. and consistent with Eqs. (3-171) and (3-172). Then these equations relate / to m. In this way rate constant estimates may be obtained by curve fitting the simulations to experimental results. 

The method is quite effective, but is not widely used now because of the ubiquity of digital computers. Zuman and Patel - 36. show circuit designs for some kinetic schemes. Williams and Bruice made good use of the analog computer in their study of the reduction of pyruvate by 1,5-dihydroflavin. In this simulation eight rate constants were evaluated variations in these parameters of 5% yielded discemibly poorer curve fits. 

Fig. 5.14 Experimental (a, black curve), fitted (a, red) and simulated (b) NIS spectrum of the Fe (Ill)-azide complex obtained at the BP86ATZVP level (J = 20 K). Bar graphs represent the corresponding intensities of the individual vibrational transitions. The blue curve represents the fitted spectrum with a background line removed (taken from [63])...

Figure 23 shows a simulation of Haure s periodic flow interruption data at time-average u = 0.86 mm/s (Curves 1 and 1 ) and u = 1.65 mm/s (Curves 2 and 2 ). Data points and the fit (dashed line) are from Fig. 19. The simulation predictions are 28 to 35% too high, which is not bad considering the 25% variation in the experimental data shown in Fig. 19. The trends, however, are not properly represented. Figure 23 predicts enhancement factors declining with r, whereas data indicate an increase at low r and...

The simulations can be made to reproduce the initial ratio of fits in Equation 9.27 using the measured 7j (ps) and fitting the distance r (nm), which is the only adjustable parameter. For canthaxan-thin and BI the experimental fits and the integrated values showed the best match in a very narrow temperature range ( 10 K) in the vicinity of the maximum enhancement in the relaxation rate. The distances obtained from the curve fits were similar to those determined from l/rM - 1/7 M0 difference, namely, 13.0 2.0 A for canthaxanthin and 10.0 2.0 A for BI. It was found for canthaxanthin, which shows no prominent peak in the relaxation rate, that the distance does not depend on l/rM -l/rM0. Using the ratio of curve fits, we can estimate the value of r for canthaxanthin as 9. 0+3.0 A in TiMCM-41 in the temperature range of 110-130K. 

To use this method, the sample is dissolved in a system containing two phases (e.g., water and octanol) such that the solution is at least about 5 x 10-4 M. The solution is acidified (or basified) and titrated with base (or acid) under controlled conditions. The shape of the ensuing titration curve is compared with the shape of a simulated curve, which is created in silico. The estimated p0Ka values (together with other variables used to construct the simulated curve such as substance concentration factor, CO2 content of the solution and acidity error) are allowed to vary systematically until the simulated curve fits as closely as possible to the experimental curve. The p0Ka values required to achieve the best fit are assumed to be the correct measured p0Ka values. This computerized calculation technique is called refinement , and is described elsewhere [14, 15]. 

B (right) Observed chemical shift (<50bs—<5l) as a function of total ligand concentration, LT. The points are values calculated from the exact simulation, and the curve is the best fit to Eq. (14), assuming no exchange contribution the value of Kj estimated from the fit is 6.2xl(T6 M. From Ref. [13] with permission. 

The relationship between the temperature difference, AT, and the input power is shown in Fig. 4.5 for microhotplate simulations and measurements. The simulated values are plotted together with the mean value of the experimental data for a set of three hotplates of the same wafer. The experimental curve was fitted with a second-order polynomial according to Eq. (3.24). As a result of the curve fit, the thermal resistance at room temperature, tjo, is 5.8 °C/mW with a standard deviation of 0.2 °C/mW, which is mainly due to variations in the etching process. 

The comparison of simulation and measurement data of an uncoated membrane is shown in Fig. 4.7. The temperature curves, T to T4, were measured with the on-membrane temperature sensors. The graphs of the simulated temperatures are denoted Si to S4. The temperature discrepancy between simulation and experiment was less than 5% for all sensors. The general shape of the temperature distribution was correctly modeled within measurement accuracy. It has to be noted that no additional fitting parameters were used for these simulations. 

In a detailed investigation of the kinetic behaviour of bases generated from (fluoren-9-ylidene)methane derivatives the problem has been overcome by computer simulation of current-time transients expected for the extended mechanism (including reaction 5), The program used for the comparison of simulated and experimental curves allows both kp and k to vary independently until the RMS deviation between the two i/t curves is minimised. The equilibrium constant for reproportionation (kf/kj) is calculable from values of Ep, (1) and Ep (2). It is important to realise that there may be any number of pairs of values of k and k which can give a good fit between experimental and simulated i/t curves. 

Digital simulation software, which is now commercially available, is useful in analyzing cyclic voltammograms of complicated electrode reactions [67]. If we assume a possible reaction mechanism and can get simulated CV curves that fit the experimental CV curves, we can confirm the reaction mechanism and obtain thermodynamic and kinetic parameters concerning the electron transfer and chemical processes. By the development of simulation softwares, cyclic voltammetry has become a very powerful technique. On the contrary, without a simulation software, cyclic voltammetry is not as convenient.14)... 

Fig. 28. X-ray diffraction rocking curve of (Ga.Mn)As/GaAs superlaltices with 20 periods (solid line). The nominal thicknesses of the GaAs and (Ga,Mn)As layers and the Mn composition x are 11.4 nm, 12.1 nm, and 0.0S4, respectively. The dashed line shows the simulated rocking curve. The fit to the experimental curve is obtained with GaAs and (Ga,Mn)As thicknesses and x of 11.14 nm, 11.79 nm, and 0.056, respectively (Shen...

Figure 1. Adsorption isotherms for Xe in A1P04-31 at T = 100, 200, and 300 K (circles, squares, and diamonds respectively) as computed with GCMC simulations. Solid curves are fits to the data using the LUD isotherm.

The rate constants were deduced from the experimental curves by curve fitting and simulation using the AnaCin Software (ref.3), assuming all the reactions to be first order in the organic reactant. The calculated reaction rate constants (in min ) depend on the weight of catalyst and are then referred to lg catalyst. 

The application of this rate law to the simulation of electrochemical behavior requires two dimensionless input parameters ktf and KC. When these are supplied, three-dimensional chronoamperometric or chronocoulometric working surfaces [34] are generated. These working surfaces both indicate first-order behavior when KC is large and second-order behavior when KC is small. Intermediate values of KC produce the variable reaction orders between one and two that are observed experimentally when the bulk olefin concentration is varied. Appropriate curve fitting of the experimental i(t,C) data to the simulation results in the evaluation of k and K details appear in the referenced work. 

---