A novel multipurpose Excel tool for equilibrium speciation based on Newton-Raphson method and on a hybrid genetic algorithm.

A new algorithm for simulation of chemical equilibria is developed, based on classical Newton-Raphson method applied to mass balance. This tool, named EST (Equilibrium Speciation Tool), is improved by using a robust Genetic Algorithm. In addition, EST works by using Excel spreadsheets and therefore offers the innovation of a great simplicity and versatility. In fact, it allows the users to simulate, or to obtain from experimental data, desired chemical-physical parameters as well as to interact with other available or freely created Excel tools. The reliability of this utility is here proved by comparison with some published data by other authors, concerning both complicated homogeneous and heterogeneous equilibria. In addiction its flexibility is tested computing thermodynamic parameters by using experimental calorimetric data referred to the complex formation of cobalt(II) with a macrocyclic ligand. A brief review and comparison of the relative robustness and quickness of main numerical methods are also reported.

[1]  Patrick Siarry,et al.  Genetic and Nelder-Mead algorithms hybridized for a more accurate global optimization of continuous multiminima functions , 2003, Eur. J. Oper. Res..

[2]  G. H. Nancollas,et al.  EQUIL. General computational method for the calculation of solution equilibriums , 1972 .

[3]  D. D. Perrin,et al.  Computer calculation of equilibrium concentrations in mixtures of metal ions and complexing species. , 1967, Talanta.

[4]  E. Joseph Billo,et al.  Excel for Chemists: A Comprehensive Guide , 2001 .

[5]  J. Allison,et al.  MINTEQA2/PRODEFA2, a geochemical assessment model for environmental systems: Version 3. 0 user's manual , 1991 .

[6]  P. Polese,et al.  SOLVERSTAT: a new utility for multipurpose analysis. An application to the investigation of dioxygenated Co(II) complex formation in dimethylsulfoxide solution. , 2003, Talanta.

[7]  Lutgarde M. C. Buydens,et al.  Evolutionary optimisation : a tutorial , 1998 .

[8]  P. Vaughan.  Evaluation of numerical techniques applied to soil solution speciation including cation exchange , 2002 .

[9]  A Pappa-Louisi,et al.  On the use of genetic algorithms for response surface modeling in high-performance liquid chromatography and their combination with the Microsoft Solver. , 2002, Journal of chromatography. A.

[10]  D. King.  A General Approach for Calculating Speciation and Poising Capacity of Redox Systems with Multiple Oxidation States: Application to Redox Titrations and the Generation of p∊-pH Diagrams , 2002 .

[11]  R. D. Levie.  A simple expression for the redox titration curve , 1992 .

[12]  G. Ginzburg.  Calculation of all equilibrium concentrations in a system of competing complexation. , 1976, Talanta.

[13]  David N. Blauch.  Java Classes for Managing Chemical Information and Solving Generalized Equilibrium Problems , 2002, J. Chem. Inf. Comput. Sci..

[14]  M. Meloun,et al.  Multiparametric curve fitting-XI: POLET computer program for estimation of formation constants and stoichiometric indices from normalized potentiometric data. , 1986, Talanta.

[15]  John Walkenbach,et al.  Excel 2003 Bible , 2001 .

[16]  Modifications of the classical genetic algorithm for non-linear fitting applied to response surface modeling in HPLC , 2001 .

[17]  L. G. Sillén,et al.  High-speed computers as a supplement to graphical methods--V. HALTAFALL, a general program for calculating the composition of equilibrium mixtures. , 1967, Talanta.

[18]  P. Gans,et al.  SUPERQUAD: an improved general program for computation of formation constants from potentiometric data , 1985 .

[19]  A. Ivaska,et al.  Determination of component concentrations in mixtures of weak and strong acids and bases by linear algebraic methods. , 1980, Talanta.

[20]  Graeme Puxty,et al.  Application of a genetic algorithm: near optimal estimation of the rate and equilibrium constants of complex reaction mechanisms , 2004 .

[21]  T. J. Wolery,et al.  EQ3/6, a software package for geochemical modeling of aqueous systems: Package overview and installation guide (Version 7.0) , 1992 .

[22]  C. D. Gelatt,et al.  Optimization by Simulated Annealing , 1983, Science.

[23]  M. F. Cardoso,et al.  The simplex-simulated annealing approach to continuous non-linear optimization , 1996 .

[24]  Milan Meloun,et al.  The effect of influential data, model and method on the precision of univariate calibration. , 2002, Talanta.

[25]  A. Holstad.  Numerical solution of nonlinear equations in chemical speciation calculations , 1999 .

[26]  David L. Parkhurst,et al.  Phreeqe--A Computer Program for Geochemical Calculations , 1980 .

[27]  A. Bianchi,et al.  The Use of Calculated Species Distribution Diagrams to Analyze Thermodynamic Selectivity , 1999 .

[28]  CALCULATION OF CHEMICAL EQUILIBRIUM IN COMPLEX SYSTEMS: SYSTEM RESTRICTIONS , 2000 .

[29]  C. Rigano,et al.  IONIC STRENGTH DEPENDENCE OF FORMATION CONSTANTS. XVII: THE CALCULATION OFEQUILIBRIUM CONCENTRATIONS AND FORMATION CONSTANTS , 1993 .

[30]  S. Rajeev,et al.  Discrete Optimization of Structures Using Genetic Algorithms , 1992 .

[31]  N. Metropolis,et al.  Equation of State Calculations by Fast Computing Machines , 1953, Resonance.

[32]  S. Gómez,et al.  Multiphase equilibria calculation by direct minimization of Gibbs free energy with a global optimization method , 2002 .

[33]  Gabriela A. Ibañez,et al.  EPSILON: A Versatile Microcomputer Program for the Spectrophotometric Data Analysis of Metal-ligand Equilibria , 1998, Comput. Chem..

[34]  Olivier Raguin,et al.  Equilibrium expert: an add-in to Microsoft Excel for multiple binding equilibrium simulations and parameter estimations. , 2002, Analytical biochemistry.

[35]  Robert de Levie,et al.  A GENERAL SIMULATOR FOR ACID-BASE TITRATIONS , 1999 .

[36]  William H. Press,et al.  The Art of Scientific Computing Second Edition , 1998 .

[37]  Malcolm C. Harrison,et al.  Data-structures and programming , 1972 .

[38]  Bruce D. McCullough,et al.  Assessing the Reliability of Statistical Software: Part I , 1998 .

[39]  P. Brassard,et al.  A feasible set for chemical speciation problems , 2000 .

[40]  Gianpiero Colonna,et al.  A hierarchical approach for fast and accurate equilibrium calculation , 2004, Comput. Phys. Commun..

[41]  Craig M. Bethke,et al.  Geochemical reaction modeling , 1996 .

[42]  J. Grime.  Analytical solution calorimetry , 1985 .

[43]  John H. Holland,et al.  Adaptation in Natural and Artificial Systems: An Introductory Analysis with Applications to Biology, Control, and Artificial Intelligence , 1992 .

[44]  William H. Press,et al.  Numerical recipes in C. The art of scientific computing , 1987 .

[45]  R. Missen,et al.  Computer calculation of ionic equilibria using species- or reaction-related thermodynamic data , 1987 .

[46]  A simplex search for conditions giving maximal and/or minimal concentrations of species in distribution plots , 1988 .

[47]  J. Barbosa,et al.  PKPOT, a program for the potentiometric study of ionic equilibria in aqueous and non-aqueous media , 1995 .

[48]  A. Shimizu,et al.  Fortran Editions of Haltafall and Letagrop. , 1970 .

[49]  David R. Williams,et al.  Computer simulation of metal-ion equilibria in biofluids: models for the low-molecular-weight complex distribution of calcium(II), magnesium(II), manganese(II), iron(III), copper(II), zinc(II), and lead(II) ions in human blood plasma , 1977 .

[50]  B. Darvell,et al.  A rapid algorithm for solution of the equations of multiple equilibrium systems-RAMESES. , 1988, Talanta.

[51]  J. Madariaga,et al.  SPECA: a program for the calculation of thermodynamic equilibrium constants from spectrophotometric data. , 1994, Talanta.

[52]  John A. Nelder,et al.  A Simplex Method for Function Minimization , 1965, Comput. J..

[53]  Lee F. Brown,et al.  Modeling Precipitation From Concentrated Solutions with the EQ3/6Chemical Speciation Codes , 1998, Comput. Chem..

[54]  P. Polese,et al.  Mixed nitrogen/oxygen ligand affinities for bipositive metal ions and dioxygen binding to cobalt(II) complexes. , 2004, Dalton transactions.

[55]  Liwu Li,et al.  Java: Data Structures and Programming , 2002, Springer Berlin Heidelberg.

[56]  G. Eriksson.  An algorithm for the computation of aqueous multi-component, multiphase equilibria , 1979 .

[57]  William R. Smith,et al.  Chemical Reaction Equilibrium Analysis: Theory and Algorithms , 1982 .

[58]  S Walsh,et al.  Non-linear curve fitting using Microsoft Excel solver. , 1995, Talanta.

[59]  M. Haugg,et al.  LAKE : A PROGRAM SYSTEM FOR EQUILIBRIUM ANALYTICAL TREATMENT OF MULTIMETHOD DATA, ESPECIALLY COMBINED POTENTIOMETRIC AND NUCLEAR MAGNETIC RESONANCE DA TA , 1996 .

[60]  Claudimir Lucio do Lago,et al.  An equilibrium simulator for multiphase equilibria based on the extent of reaction and Newton-Raphson method with globally convergent strategy (SEQEx2) , 1996 .

[61]  Antonio Brunetti,et al.  A Fast and precise genetic algorithm for a non-linear fitting problem , 2000 .

[62]  C. Rigano,et al.  The calculation of equilibrium concentrations in large multimetal/multiligand systems , 1986 .

[63]  P. Gans,et al.  Hyperquad simulation and speciation (HySS): a utility program for the investigation of equilibria involving soluble and partially soluble species , 1999 .

[64]  William L. Goffe.  SIMANN: A Global Optimization Algorithm using Simulated Annealing , 1996 .

[65]  G. R. Hext,et al.  Sequential Application of Simplex Designs in Optimisation and Evolutionary Operation , 1962 .

[66]  J. M. Casas,et al.  Aqueous speciation of sulfuric acid–cupric sulfate solutions , 2000 .

[67]  P. Gans,et al.  Investigation of equilibria in solution. Determination of equilibrium constants with the HYPERQUAD suite of programs. , 1996, Talanta.

[68]  D. L. Parkhurst,et al.  User's guide to PHREEQC (Version 2)-a computer program for speciation, batch-reaction, one-dimensional transport, and inverse geochemical calculations , 1999 .

[69]  A. Martell,et al.  Program PKAS: a novel algorithm for the computation of successive protonation constants , 1982 .

[70]  Y. K. Rao.  IECMQ: a computer program for calculating equilibrium-state in aqueous systems , 2000 .

[71]  Dmitrii A. Kulik,et al.  Gibbs energy minimization approach to modeling sorption equilibria at the mineral-water interface: Thermodynamic relations for multi-site-surface complexation , 2002 .

[72]  A. Safavi,et al.  Evaluation of Formation Constants, Molar Absorptivities of Metal Complexes, and Protonation Constants of Acids by Nonlinear Curve Fitting Using Microsoft Excel Solver and User-Defined Function , 1999 .

[73]  D. Leggett.  Machine computation of equilibrium concentrations-some practical considerations. , 1977, Talanta.