Carboxymethyl-beta-cyclodextrin (CMCD) has been proposed for remediation of metal-contaminated sediments. This research presents stability constants for CMCD-lead complexes, and demonstrates a rigorous methodology for estimating stability constants for metal-complexing agents. The conditional stability constant for the lead-CMCD aqueous complex was determined to be 10(5.18) with the 95% confidence interval ranging from 10(5.14) to 10(5.22). The best fit for experimental data was made by assuming a reaction between divalent CMCD(2-) and Pb(2+) and using the WATEQ activity coefficient formulation. The optimized value was derived from experimental data with the geochemical model PHREEQC coupled to UCODE_2005, a parameter optimization program. Like FITEQL, UCODE has a built-in option to optimize parameter values by minimizing the weighted sum of squared residuals (WSSR). However, our approach not only allows rapid, automatic optimization of the stability constant, but also allows determination of uncertainties in estimated parameter values and statistical analysis to assess the appropriateness of the conceptual model. The automation of the process allows testing of multiple conceptual models and the final values produced are internally consistent with the PHREEQC database. In this case five different conceptual models to describe the metal complexation and protonation reactions of CMCD were considered.