{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2025,10,12]],"date-time":"2025-10-12T01:57:07Z","timestamp":1760234227673,"version":"build-2065373602"},"reference-count":34,"publisher":"MDPI AG","issue":"5","license":[{"start":{"date-parts":[[2021,4,27]],"date-time":"2021-04-27T00:00:00Z","timestamp":1619481600000},"content-version":"vor","delay-in-days":0,"URL":"https:\/\/creativecommons.org\/licenses\/by\/4.0\/"}],"funder":[{"DOI":"10.13039\/100009117","name":"Technische Universit\u00e4t Chemnitz","doi-asserted-by":"publisher","award":["ID 29_2021"],"award-info":[{"award-number":["ID 29_2021"]}],"id":[{"id":"10.13039\/100009117","id-type":"DOI","asserted-by":"publisher"}]}],"content-domain":{"domain":[],"crossmark-restriction":false},"short-container-title":["Computation"],"abstract":"<jats:p>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\u2019 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\u2013Nickel electrolyte. This paper finds a stable, convergent and efficient methodology to compute stability constants from potentiometric titration curves.<\/jats:p>","DOI":"10.3390\/computation9050055","type":"journal-article","created":{"date-parts":[[2021,4,27]],"date-time":"2021-04-27T06:19:11Z","timestamp":1619504351000},"page":"55","update-policy":"https:\/\/doi.org\/10.3390\/mdpi_crossmark_policy","source":"Crossref","is-referenced-by-count":6,"title":["Stabilization of the Computation of Stability Constants and Species Distributions from Titration Curves"],"prefix":"10.3390","volume":"9","author":[{"ORCID":"https:\/\/orcid.org\/0000-0002-9229-8793","authenticated-orcid":false,"given":"Stephan Daniel","family":"Schwoebel","sequence":"first","affiliation":[{"name":"Materials and Surface Engineering Group, Institute of Materials Science and Engineering, Chemnitz University of Technology, D-09107 Chemnitz, Germany"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"ORCID":"https:\/\/orcid.org\/0000-0001-5792-8625","authenticated-orcid":false,"given":"Dominik","family":"H\u00f6hlich","sequence":"additional","affiliation":[{"name":"Materials and Surface Engineering Group, Institute of Materials Science and Engineering, Chemnitz University of Technology, D-09107 Chemnitz, Germany"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Thomas","family":"Mehner","sequence":"additional","affiliation":[{"name":"Materials and Surface Engineering Group, Institute of Materials Science and Engineering, Chemnitz University of Technology, D-09107 Chemnitz, Germany"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"ORCID":"https:\/\/orcid.org\/0000-0002-2390-9159","authenticated-orcid":false,"given":"Thomas","family":"Lampke","sequence":"additional","affiliation":[{"name":"Materials and Surface Engineering Group, Institute of Materials Science and Engineering, Chemnitz University of Technology, D-09107 Chemnitz, Germany"}],"role":[{"role":"author","vocabulary":"crossref"}]}],"member":"1968","published-online":{"date-parts":[[2021,4,27]]},"reference":[{"doi-asserted-by":"crossref","unstructured":"Dumpala, R., Srivastava, A., and Rawat, N. (2021). Experimental and theoretical approach to probe the aquatic speciation of transuranic (neptunyl) ion in presence of two omnipresent organic moieties. Chemosphere, 273.","key":"ref_1","DOI":"10.1016\/j.chemosphere.2021.129745"},{"doi-asserted-by":"crossref","unstructured":"Mahmoud, S., Taha, M., Mohamed, R., Khaled, E., and Abdel-khalek, A. (2021). Complexation of chromium (III) with the antifibrinolytic drug tranexamic acid: Formation, kinetics, and molecular modeling studies. J. Mol. Liq., 329.","key":"ref_2","DOI":"10.1016\/j.molliq.2021.115513"},{"key":"ref_3","first-page":"52","article-title":"Electroreduction of Ni (II) and Co (II) from Pyrophosphate Solutions","volume":"16","author":"Cesiulis","year":"2010","journal-title":"Medziagotyra"},{"doi-asserted-by":"crossref","unstructured":"Schwoebel, S.D., Mehner, T., and Lampke, T. (2021). 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-(\u03b2-alanine)-Electrolyte. Algorithms, 14.","key":"ref_4","DOI":"10.3390\/a14040113"},{"doi-asserted-by":"crossref","unstructured":"Aver\u00f3s, J., Llorens, J., and Uribe-Kaffure, R. (2020). Numerical simulation of non-linear models of reaction-diffusion for a DGT sensor. Algorithms, 13.","key":"ref_5","DOI":"10.3390\/a13040098"},{"key":"ref_6","doi-asserted-by":"crossref","first-page":"4869","DOI":"10.1021\/es200609v","article-title":"Key Role of the Resin Layer Thickness in the Lability of Complexes Measured by DGT","volume":"45","author":"Mongin","year":"2011","journal-title":"Environ. Sci. Technol."},{"key":"ref_7","doi-asserted-by":"crossref","first-page":"46","DOI":"10.1080\/00202967.2016.1214354","article-title":"Application of a duplex diffusion layer model to pulse reverse plating","volume":"95","author":"Green","year":"2017","journal-title":"Trans. IMF"},{"key":"ref_8","doi-asserted-by":"crossref","first-page":"1726","DOI":"10.1039\/D0QI00200C","article-title":"Unraveling the complexation mechanism of actinide(iii) and lanthanide(iii) with a new tetradentate phenanthroline-derived phosphonate ligand","volume":"7","author":"Xu","year":"2020","journal-title":"Inorg. Chem. Front."},{"key":"ref_9","doi-asserted-by":"crossref","first-page":"119392","DOI":"10.1016\/j.ica.2019.119392","article-title":"Lanthanides(III) and Silver(I) complex formation with triamines in DMSO: The effect of ligand cyclization","volume":"503","author":"Melchior","year":"2020","journal-title":"Inorganica Chim. Acta"},{"key":"ref_10","doi-asserted-by":"crossref","first-page":"3098","DOI":"10.1103\/PhysRevA.38.3098","article-title":"Density-functional exchange-energy approximation with correct asymptotic behavior","volume":"38","author":"Becke","year":"1988","journal-title":"Phys. Rev. A"},{"key":"ref_11","doi-asserted-by":"crossref","first-page":"311","DOI":"10.1016\/S0010-8545(98)00260-4","article-title":"Hyperquad simulation and speciation (HySS): A utility program for the investigation of equilibria involving soluble and partially soluble species","volume":"184","author":"Alderighi","year":"1999","journal-title":"Coord. Chem. Rev."},{"key":"ref_12","doi-asserted-by":"crossref","first-page":"1739","DOI":"10.1016\/0039-9140(96)01958-3","article-title":"Investigation of equilibria in solution. Determination of equilibrium constants with the HYPERQUAD suite of programs","volume":"43","author":"Gans","year":"1996","journal-title":"Talanta"},{"unstructured":"Hinze, P., Pimau, R., Ulbrich, M., and Ulbrich, S. (2009). Optimization with PDE Constraints, Springer. [1st ed.]. Mathematical Modelling: Theory and Applications.","key":"ref_13"},{"doi-asserted-by":"crossref","unstructured":"Richter, M. (2015). Inverse Probleme-Grundlagen, Theorie und Anwendungsbeispiele, Springer. [1st ed.].","key":"ref_14","DOI":"10.1007\/978-3-662-45811-2"},{"doi-asserted-by":"crossref","unstructured":"Langtangen, H.P., and Pedersen, G.K. (2016). Ordinary differential equation models. Scaling of Differential Equations, Springer International Publishing.","key":"ref_15","DOI":"10.1007\/978-3-319-32726-6"},{"doi-asserted-by":"crossref","unstructured":"M\u00e4kel\u00e4, M., and Neittaanm\u00e4ki, P. (1992). Nonsmooth Optimization-Analysis and Algorithms with Applications to Optimal Control, World Scientific Publishing Co. Pte. Ltd.. [1st ed.].","key":"ref_16","DOI":"10.1142\/1493"},{"doi-asserted-by":"crossref","unstructured":"Lange, K. (2013). Optimization, Springer. [2nd ed.]. Springer Texts in Statistics.","key":"ref_17","DOI":"10.1007\/978-1-4614-5838-8"},{"doi-asserted-by":"crossref","unstructured":"Engquist, B. (2015). Homotopy Methods. Encyclopedia of Applied and Computational Mathematics, Springer.","key":"ref_18","DOI":"10.1007\/978-3-540-70529-1"},{"doi-asserted-by":"crossref","unstructured":"Floudas, C.A., and Pardalos, P.M. (2009). Globally convergent homotopy methodsGlobally Convergent Homotopy Methods. Encyclopedia of Optimization, Springer.","key":"ref_19","DOI":"10.1007\/978-0-387-74759-0"},{"unstructured":"Martell, A.M. (1992). Determination and Use of Stability Constants, VCH Puplishers Inc.. [2nd ed.].","key":"ref_20"},{"key":"ref_21","first-page":"33","article-title":"Chapter 3-Basic Thermodynamics","volume":"Volume 1","author":"Saunders","year":"1998","journal-title":"CALPHAD: Calculation of Phase Diagrams"},{"key":"ref_22","doi-asserted-by":"crossref","first-page":"155","DOI":"10.1016\/j.jct.2008.10.003","article-title":"Theory for the determination of activity coefficients of strong electrolytes with regard to concentration dependence of hydration numbers","volume":"41","author":"Zaitsev","year":"2009","journal-title":"J. Chem. Thermodyn."},{"key":"ref_23","doi-asserted-by":"crossref","first-page":"154507","DOI":"10.1063\/1.3489418","article-title":"The nonmonotonic concentration dependence of the mean activity coefficient of electrolytes is a result of a balance between solvation and ion-ion correlations","volume":"133","author":"Vincze","year":"2010","journal-title":"J. Chem. Phys."},{"key":"ref_24","doi-asserted-by":"crossref","first-page":"49","DOI":"10.1016\/j.aca.2019.05.002","article-title":"Activity-based analysis of potentiometric pH titrations","volume":"1075","author":"Karimvand","year":"2019","journal-title":"Anal. Chim. Acta"},{"unstructured":"Bosch, S. (2009). Algebra, Springer.","key":"ref_25"},{"key":"ref_26","doi-asserted-by":"crossref","first-page":"341","DOI":"10.1007\/BF01585750","article-title":"The augmented lagrangian method for equality and inequality constraints in hilbert spaces","volume":"46","author":"Ito","year":"1990","journal-title":"Math. Program."},{"key":"ref_27","doi-asserted-by":"crossref","first-page":"272","DOI":"10.1137\/16M1107103","article-title":"An Augmented Lagrangian Method for Optimization Problems in Banach Spaces","volume":"56","author":"Kanzow","year":"2018","journal-title":"SIAM J. Control. Optim."},{"unstructured":"Jahn, J. (2007). Introduction to the Theory of Nonlinear Optimization, Springer.","key":"ref_28"},{"unstructured":"Nocedal, J., and Wright, S. (2006). Numerical Optimization, Springer.","key":"ref_29"},{"doi-asserted-by":"crossref","unstructured":"Polak, E. (1997). Unconstrained Optimization. Optimization: Algorithms and Consistent Approximations, Springer.","key":"ref_30","DOI":"10.1007\/978-1-4612-0663-7_1"},{"doi-asserted-by":"crossref","unstructured":"Ulbrich, M., and Ulbrich, S. (2012). Restringierte Optimierung. Nichtlineare Optimierung, Springer.","key":"ref_31","DOI":"10.1007\/978-3-0346-0654-7"},{"key":"ref_32","doi-asserted-by":"crossref","first-page":"149","DOI":"10.1007\/PL00011391","article-title":"A trust region method based on interior point techniques for nonlinear programming","volume":"89","author":"Byrd","year":"2000","journal-title":"Math. Program."},{"key":"ref_33","doi-asserted-by":"crossref","first-page":"346","DOI":"10.1134\/S1070328407050065","article-title":"Interaction of the Ni2+ ion with citric acid in an aqueous solution","volume":"38","author":"Zelenin","year":"2007","journal-title":"Russ. J. Coord. Chem."},{"key":"ref_34","doi-asserted-by":"crossref","first-page":"975","DOI":"10.1134\/S0036023611060192","article-title":"Correlations between the stability constants of metal hydroxo complexes and the solubility products of crystalline hydroxides. Series of the polarizing effect of metal cations","volume":"56","author":"Orlov","year":"2011","journal-title":"Russ. J. Inorg. Chem."}],"container-title":["Computation"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/www.mdpi.com\/2079-3197\/9\/5\/55\/pdf","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2025,10,11]],"date-time":"2025-10-11T05:53:18Z","timestamp":1760161998000},"score":1,"resource":{"primary":{"URL":"https:\/\/www.mdpi.com\/2079-3197\/9\/5\/55"}},"subtitle":[],"short-title":[],"issued":{"date-parts":[[2021,4,27]]},"references-count":34,"journal-issue":{"issue":"5","published-online":{"date-parts":[[2021,5]]}},"alternative-id":["computation9050055"],"URL":"https:\/\/doi.org\/10.3390\/computation9050055","relation":{},"ISSN":["2079-3197"],"issn-type":[{"type":"electronic","value":"2079-3197"}],"subject":[],"published":{"date-parts":[[2021,4,27]]}}}