Stabilization of the Computation of Stability Constants and Species Distributions from Titration Curves

Thermodynamic equilibria and concentrations in thermodynamic equilibria are of major importance in chemistry, chemical engineering, physical chemistry, medicine etc. due to a vast spectrum of applications. E.g., concentrations in thermodynamic equilibria play a central role for the estimation of drug delivery, the estimation of produced mass of products of chemical reactions, the estimation of deposited metal during electro plating and many more. Species concentrations in thermodynamic equilibrium are determined by the system of reactions and to the reactions’ associated stability constants. In many applications the stability constants and the system of reactions need to be determined. The usual way to determine the stability constants is to evaluate titration curves. In this context, many numerical methods exist. One major task in this context is that the corresponding inverse problems tend to be unstable, i.e., the output is strongly affected by measurement errors, and can output negative stability constants or negative species concentrations. In this work an alternative model for the species distributions in thermodynamic equilibrium, based on the models used for HySS or Hyperquad, and titration curves is presented, which includes the positivity of species concentrations and stability constants intrinsically. Additionally, in this paper a stabilized numerical methodology is presented to treat the corresponding model guaranteeing the convergence of the algorithm. The numerical scheme is validated with clinical numerical examples and the model is validated with a Citric acid–Nickel electrolyte. This paper finds a stable, convergent and efficient methodology to compute stability constants from potentiometric titration curves.

[1]  H. Cesiulis,et al.  Electroreduction of Ni(II) and Co(II) from Pyrophosphate Solutions , 2010 .

[2]  A. Zaitsev,et al.  Theory for the determination of activity coefficients of strong electrolytes with regard to concentration dependence of hydration numbers , 2009 .

[3]  Hao Zhang,et al.  Key role of the resin layer thickness in the lability of complexes measured by DGT. , 2011, Environmental science & technology.

[4]  S. Roy,et al.  Application of a duplex diffusion layer model to pulse reverse plating , 2017 .

[5]  Christian Kanzow,et al.  An Augmented Lagrangian Method for Optimization Problems in Banach Spaces , 2018, SIAM J. Control. Optim..

[6]  Julian W. Vincze,et al.  The nonmonotonic concentration dependence of the mean activity coefficient of electrolytes is a result of a balance between solvation and ion-ion correlations. , 2010, The Journal of chemical physics.

[7]  M. Tolazzi,et al.  Lanthanides(III) and Silver(I) complex formation with triamines in DMSO: The effect of ligand cyclization , 2020 .

[8]  Joan Cecilia Averós,et al.  Numerical Simulation of Non-Linear Models of Reaction - Diffusion for a DGT Sensor , 2020, Algorithms.

[9]  Hans Petter Langtangen,et al.  Ordinary differential equation models , 2016 .

[10]  Lei Xu,et al.  Unraveling the complexation mechanism of actinide(iii) and lanthanide(iii) with a new tetradentate phenanthroline-derived phosphonate ligand , 2020 .

[11]  Jorge Nocedal,et al.  A trust region method based on interior point techniques for nonlinear programming , 2000, Math. Program..

[12]  A. Becke,et al.  Density-functional exchange-energy approximation with correct asymptotic behavior. , 1988, Physical review. A, General physics.

[13]  Kazufumi Ito,et al.  The augmented lagrangian method for equality and inequality constraints in hilbert spaces , 1990, Math. Program..

[14]  M. Maeder,et al.  Activity-based analysis of potentiometric pH titrations. , 2019, Analytica chimica acta.

[15]  R. Tennant.  Algebra , 1941, Nature.

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

[17]  O. Y. Zelenin.  Interaction of the Ni2+ ion with citric acid in an aqueous solution , 2007 .

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

[19]  J. Jahn.  Introduction to the Theory of Nonlinear Optimization , 1994 .

[20]  Y. Orlov,et al.  Correlations between the stability constants of metal hydroxo complexes and the solubility products of crystalline hydroxides. Series of the polarizing effect of metal cations , 2011 .

[21]  T. Mehner,et al.  On a Robust and Efficient Numerical Scheme for the Simulation of Stationary 3-Component Systems with Non-Negative Species-Concentration with an Application to the Cu Deposition from a Cu-(β-alanine)-Electrolyte , 2021, Algorithms.

[22]  N. Rawat,et al.  Experimental and theoretical approach to probe the aquatic speciation of transuranic (neptunyl) ion in presence of two omnipresent organic moieties. , 2021, Chemosphere.

[23]  Reham A. Mohamed,et al.  Complexation of chromium (III) with the antifibrinolytic drug tranexamic acid: Formation, kinetics, and molecular modeling studies , 2021 .