Stability constants for Zirconium(IV) complexes with EDTA, CDTA, and DTPA in perchloric acid solutions
Not Available
SEARCH · Search NASA
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Not Available
Stability constants provide insight into ion complexation in water. While computational studies have been shown to model the energy of the complexation successfully using a thermodynamic cycle approach, it does not extend to calculating the stability constants for 1:1 lanthanide to ligand complexes in solution. Using B3LYP and 6-31+G* Pople basis with small core effective core potential (ECP) on the lanthanum ion, and a solvent model based on the full solute electron density (SMD) solvation model we computed and compared with previously published stability constants of the ligands: acetate, acetohydroximate, acetylacetonate, methanoate, tropolonate, hydroxide, catecholate, malonate, oxalate, phthalate, and sulfate. The best R 2 values for the thermodynamic cycle can only be determined by separating the mono and divalent ions to achieve an R 2 value of 0.86 and 0.74 for mono and divalent ions, respectively. We show that by optimizing the lanthanide-ligand structures in implicit solvent, we achieve an improved correlation between experimental and computed stability constants of R 2 value of 0.89 for the combined mono and divalent ions.
Calculating stability constants for lanthanide--nitrate complexes in aqueous solution is challenging due to the complex free-energy landscapes of the participating species. In this work, we compare cluster-continuum solvation and condensed-phase approaches using universal machine learning interatomic potentials (MLIPs) for determining the stability constants of lanthanide--nitrate complexes in aqueous solutions. Within the cluster--continuum solvation framework at the B3LYP level of theory, reactions involving lanthanide coordination numbers of both 8 and 9 are found to be relevant. After an empirical linear free-energy correction, the cluster--continuum results fall on the same order-of-magnitude scale as the experimental stability constants. By contrast, MACE-MP0 MLIP underestimates lanthanide hydration numbers and gives PMF-derived stability constants with large deviations from experiment, whereas MACE-MATPES-R2SCAN improves both hydration structure and the stability-constant scale but still does not quantitatively reproduce the detailed lanthanide trend. Across both approaches, nitrate binding is best viewed as a labile coordination motif rather than a fixed mono- or bidentate structure, with hydration-shell structure influencing which configurations are favored. Overall, the cluster--continuum calculations provide a practical semi-quantitative baseline for the experimental stability-constant scale, while the explicit-solvent MLIP benchmarks show clear progress from MACE-MP0 to MACE-MATPES-r2SCAN but also highlight the need for lanthanide-targeted training, fine-tuning, or improved long-range and polarization treatments to obtain predictive thermodynamics for complex aqueous lanthanide chemistry.
Stability constants of simple reactions involving addition of the NO 3 − ion to hydrated metal complexes, [M(H 2 O) x ] n + are calculated with a computational workflow developed using cloud computing resources.
Ligand selectivity to specific lanthanide (Ln) ions is key to the separation of rare earth elements from each other. Ligand selectivity can be quantified with relative stability constants (measured experimentally) or relative binding energies (calculated computationally). The relative stability constants of EDTA (ethylenediaminetetraacetic acid) with La 3+ , Eu 3+ , Gd 3+ , and Lu 3+ were predicted from relative binding energies, which were quantified using electronic structure calculations with relativistic effects and based on the molecular structures of Ln–EDTA complexes in solution from density functional theory molecular dynamics simulations. The protonation state of an EDTA amine group was varied to study pH ~7 and ~11 conditions. Further, simulations at 25 °C and 90 °C were performed to elucidate how structures of Ln–EDTA complexes varying with temperature are related to complex stabilities at different pH conditions. Relative stability trends are predicted from computation for varying Ln 3+ ions (La, Eu, Gd, Lu) with a single ligand (EDTA at pH ~11), as well as for a single Ln 3+ ion (La) with varying ligands (EDTA at pH ~7 and ~11). Changing the protonation state of an EDTA amine site significantly changes the solution structure of the Ln–EDTA complex resulting in a reduction of the complex stability. As a result, increased Ln–ligand complex stability is correlated to reduced structural variations in solution upon an increase in temperature.
Uranyl ion, UO 2 2+ , and its aqueous complexes with organic and inorganic ligands can be the dominant species for uranium transport on the Earth surface or in a nuclear waste disposal system if an oxidizing condition is present. As an important biodegradation product, oxalate, C 2 O 4 2– , is ubiquitous in natural environments and is known for its ability to complex with the uranyl ion. Oxalate can also form solid phases with uranyl ion in certain environments thus limiting uranium migration. Therefore, the determination of stability constants for aqueous and solid uranyl oxalate complexes is important not only to the understanding of uranium mobility in natural environments, but also to the performance assessment of nuclear waste disposal. In this work, we developed a thermodynamic model for the UO 2 2+ –Na + –H + –Cl – –ClO 4 – –C 2 O 4 2– –NO 3 – –H 2 O system to ionic strength up to ~11 mol•kg –1 . We constrained the stability constants for UO 2 C 2 O 4 (aq) and UO 2 (C 2 O 4 ) 2 2– at infinite dilution based on our evaluation of the literature data over a wide range of ionic strengths up to ~11 mol•kg –1 . We also obtained the solubility constants at infinite dilution for solid uranyl oxalates, UO 2 C 2 O 4 •3H 2 O, based on the solubility data over a wide range of ionic strengths. The developed model will enable for the accurate stability assessment of oxalate complexes affecting uranium mobility under a wide range of conditions including those in deep geological repositories.
Correction for ‘The solution structures and relative stability constants of lanthanide–EDTA complexes predicted from computation’ by Ravi D. O’Brien et al. , Phys. Chem. Chem. Phys. , 2022, 24 , 10263–10271, https://doi.org/10.1039/D2CP01081J.
The new LOGKPREDICT program integrates HostDesigner molecular design software with the machine learning (ML) program Chemprop. By supplying HostDesigner with predicted log K values, LOGKPREDICT enhances the computer-aided molecular design process by ranking ligands directly by metal–ligand binding strength. Harnessing reliable experimental data from a historic National Institute of Standards and Technology (NIST) database and data from the International Union of Pure and Applied Chemistry (IUPAC), we train message passing neural net algorithms. The multi-metal NIST-based ML model has a root mean square error (RMSE) of 0.629 ± 0.044 (R 2 of 0.960 ± 0.006), while two versions of lanthanide-only IUPAC-based ML models have, respectively, RMSE of 0.764 ± 0.073 (R 2 of 0.976 ± 0.005) and 0.757 ± 0.071 (R 2 of 0.959 ± 0.007). For relative log K predictions on an out-of-sample set of six ligands, demonstrating metal ion selectivity, the RMSE value reaches a commendably low 0.25. Here we showcase the use of LOGKPREDICT in identifying ligands with high selectivity for lanthanides in aqueous solutions, a finding supported by recent experimental evidence. We also predict new ligands yet to be verified experimentally. Therefore, our ML models implemented through LOGKPREDICT and interfaced with the ligand design software HostDesigner pave the way for designing new ligands with predetermined selectivity for competing metal ions in an aqueous solution.
Abstract not provided.
The trivalent f-element coordination chemistry of a novel aminopolycarboxylate complexant was investigated. The novel reagent is an octadentate complexant that resembles diethylenetriamine-N,N,N',N",N"-pentaacetic acid (DTPA), but a single N-acetate pendant arm was substituted with a N-2-pyrazinylmethyl functional group. Thermodynamic studies of ligand protonation and trivalent lanthanide, americium and curium, complexation by N-2-pyrazinylmethyldiethylenetriamine-N,N',N",N"-tetraacetic acid (DTTA-PzM) emphasize the strong electron withdrawing influence of the N-2-pyrazinylmethyl group. Particularly, DTTA-PzM is more acidic compared to a N-2-pyridinylmethyl-substituted structural equivalent, DTTA-PyM, with a substantial lowering of pK7, corresponding to the protonation of a second aliphatic amine site. The participation of the pyrizyl nitrogen in the metal ion coordination sphere is observed from the fluorescence lifetime decay measurements of metal hydration and the interpretation of the stability constants for ML– and MHL(aq) complexes. The overall conditional stability constants for the trivalent f-element complexation by DTTA-PzM complexes decrease, relative to DTTA-PyM, as expected based on the lower basicity of pyrazine in water relative to pyridine. Replacement of the N-2-pyridinylmethyl group with N-2-pyrazinylmethyl, while enhancing the total acidity of DTTA-PzM, also reduces its softness, as manifested by a small lowering of β 101 Am/Nd and liquid–liquid separation of trivalent lanthanides from trivalent americium. Despite this, the 4f/5f separation is doubled when DTTA-PzM replaces DTPA as an aqueous complexant in solvent extraction.
Technetium-99 is a high yield (~6% fission yield) fission product and long-lived (2.13×10 5 year half-life) component of nuclear waste that will be disposed of in a geological repository. Some 99 Tc has been released into the environment due to nuclear fuel and weapon production activities at sites such as Hanford, WA. Strongly complexing ligands such as ethylenediamine- N , N , N ', N '-tetraacetic acid (EDTA) are known to increase Tc(IV) solubility and mobility in environmental systems and an accurate quantification of the complexation of Tc(IV) with EDTA is important for predicting its behavior in a geological repository. A liquid–liquid extraction system utilizing 0.2 M TOPO in dodecane was used to measure the stability constants of Tc(IV)-EDTA in 0.50 m NaNO 3 at variable temperatures (14.0±0.1, 25.0±0.1, and 32.0±0.1°C). The acid dependence of the apparent stability constants in the p C H range of 2.00–2.70 indicated the formation of TcO(EDTA) 2- (log β 101 =17.9±0.3, 25.0±0.1°C) and a protonated complex TcO(H)(EDTA) - (log β 111 =20.5±0.1, 25.0±0.1°C). The associated thermodynamic parameters Δ r G 101 =-101.7±0.4 kJ·mol -1 , Δ r H 101 =-47±9 kJ·mol -1 , Δ r S 101 =179±36 J·mol -1 ·K -1 , Δ r G 111 =-117.2±0.3 kJ·mol -1 , Δ r H 111 =-23±5 kJ·mol -1 , and Δ r S 111 =315±63 J·mol -1 ·K -1 (0.50 m NaNO 3 , 25.0±0.1°C) were determined by van’t Hoff analysis. The formation of each Tc(IV)–EDTA complex is exothermic and present favorable entropy terms.
The design of metal-concerned solvent extraction systems frequently leverages thermodynamically derived differences in selectivity. An alternative approach, leveraging kinetic control, has been considered much less seriously. Our recent manuscript describing DOTA (1,4,7,10-tetraazacyclododecane-1,4,7,10-tetraacetic acid) complexation kinetics across the lanthanide series shows the observed rate constant steadily increases across the series, with some non-monotonic behavior observed at terbium and thulium. This contrasts the thermodynamic stability constant trend, where lanthanide-DOTA stability constants initially increase and then plateau as a function of ionic radii after samarium. To leverage the kinetic differences of DOTA with the lanthanides across the series, kinetically based separations must be utilized. Since DOTA has very slow complexation kinetics, a separations system must expedite DOTA-metal complexation to allow a separation approximating practical application. Here in this report, a DOTA-based solvent extraction system, where DOTA is the aqueous holdback reagent and bis-2,4,4-trimethylpentylphosphinic acid (Cyanex 272) is the organic phase extractant, is demonstrated and compares the separations chemistry of Nd, Eu and Am. The slowness of DOTA complexation was addressed by heating the system. Results showed, in general, separation between metals is better during early phase contact, and diminishes under longer contact times. Under all conditions, separations are better than would be predicted based on a thermodynamic basis. This report suggests that while slowly complexing ligands classically used for biological applications may not be appropriate for thermodynamically designed metal separations, their use for kinetically based systems may be appropriate and enable a new design basis for f-element separations.
Development of a defensible source-term model (STM), usually a thermodynamical model for radionuclide solubility calculations, is critical to a performance assessment (PA) of a geologic repository for nuclear waste disposal. Such a model is generally subjected to rigorous regulatory scrutiny. In this article, we highlight key guiding principles for STM model development and validation in nuclear waste management. We illustrate these principles by closely examining three recently developed thermodynamic models with the Pitzer formulism for aqueous H + —Nd 3+ —NO 3 - (—oxalate) systems in a reverse alphabetical order of the authors: the XW model developed by Xiong and Wang, the OWC model developed by Oakes et al., and the GLC model developed by Guignot et al., among which the XW model deals with trace activity coefficients for Nd(III), while the OWC and GLC models are for concentrated Nd(NO 3 ) 3 electrolyte solutions. The principles highlighted include the following: (1) Principle 1. Validation against independent experimental data: A model should be validated against experimental data or field observations that have not been used in the original model parameterization. We tested the XW model against multiple independent experimental data sets including electromotive force (EMF), solubility, water vapor, and water activity measurements. The results show that the XW model is accurate and valid for its intended use for predicting trace activity coefficients and therefore Nd solubility in repository environments. (2) Principle 2. Testing for relevant and sensitive variables: Solution pH is such a variable for an STM and easily acquirable. All three models are checked for their ability to predict pH conditions in Nd(NO 3 ) 3 electrolyte solutions. The OWC model fails to provide a reasonable estimate for solution pH conditions, thus casting serious doubt on its validity for a source-term calculation. In contrast, both the XW and GLC models predict close-to-neutral pH values, in agreement with experimental measurements. (3) Principle 3. Honoring physical constraints: Upon close examination, it is found that the Nd(III)-NO 3 association schema in the OWC model suffers from two shortcomings. Firstly, its second stepwise stability constant for Nd(NO 3 ) 2+ (log K 2 ) is much higher than the first stepwise stability constant for NdNO 3 2+ (log K 1 ), thus violating the general rule of (log K 2 –log K 1 ) < 0, or $\frac{K1}{K2}$>1. Secondly, the OWC model predicts abnormally high activity coefficients for Nd(NO 3 ) 2 + (up to ~900) as the concentration increases. (4) Principle 4. Minimizing degrees of freedom for model fitting: The OWC model with nine fitted parameters is compared with the GLC model with five fitted parameters, as both models apply to the concentrated region for Nd(NO 3 ) 3 electrolyte solutions. The latter appears superior to the former because the latter can fit osmotic coefficient data equally well with fewer model parameters. The work presented here thus illustrates the salient points of geochemical model development, selection, and validation in nuclear waste management.
Targeted α therapy (TAT) of soft-tissue cancers using the α particle-emitting radionuclide 223 Ra holds great potential because of its favorable nuclear properties, adequate availability, and established clinical use for treating metastatic prostate cancer of the bone. Despite these advantages, the use of 223 Ra has been largely overshadowed by other α emitters due to its challenging chelation chemistry. A key criterion that needs to be met for a radionuclide to be used in TAT is its stable attachment to a targeting vector via a bifunctional chelator. The low charge density of Ra 2+ arising from its large ionic radius weakens its electrostatic binding interactions with chelators, leading to insufficient complex stability in vivo. In this study, we synthesized and evaluated macropa-XL as a novel chelator for 223 Ra. It bears a large 21-crown-7 macrocyclic core and two picolinate pendent groups, which we hypothesized would effectively saturate the large coordination sphere of the Ra 2+ ion. The structural chemistry of macropa-XL was first established with the nonradioactive Ba 2+ ion using X-ray diffraction and X-ray absorption spectroscopy, which revealed the formation of an 11-coordinate complex in a rare anti pendent-arm configuration. Subsequently, the stability constant of the [Ra(macropa-XL)] complex was determined via competitive cation exchange with 223 Ra and 224 Ra radiotracers and compared with that of macropa, the current state-of-the-art chelator for Ra 2+ . A moderate log K ML value of 8.12 was measured for [Ra(macropa-XL)], which is approximately 1.5 log K units lower than the stability constant of [Ra(macropa)]. This relative decrease in Ra 2+ complex stability for macropa-XL versus macropa was further probed using density functional theory calculations. Additionally, macropa-XL was radiolabeled with 223 Ra, and the kinetic stability of the resulting complex was evaluated in human serum. Furthermore, although macropa-XL could effectively bind 223 Ra under mild conditions, the complex appeared to be unstable to transchelation. Collectively, this study sheds additional light on the chelation chemistry of the exotic Ra 2+ ion and contributes to the small, but growing, number of chelator development efforts for 223 Ra-based TAT.
Separation of An(III) and Ln(III) ions will benefit the recycling of used nuclear fuel (UNF). For this purpose, many ligands have been tested over the years, and several separation processes have been successfully demonstrated on the laboratory scale. Current research aims at the development of new ligands that are built only with carbon, hydrogen, oxygen, and nitrogen (CHON), as they can be incinerated completely without secondary waste production. Here, we tested a new class of water-soluble ligands, the bistriazinyl-octa-carboxylic acids. One member, in particular, 2,6-bis-[5,6-di(3,4-dicarboxyphenyl)-1,2,4- triazin-3-yl]-pyridine (BTPOA), was found to be suitable for the selective separation of Am(III) and Cm(III) ions from Ln(III) ions and may act as a CHON alternative to its sulfonated analogue (SO 3 -Ph-BTP). BTPOA exhibited good extraction results and a high selectivity for Am(III) over Eu(III) ions. This ligand’s complexation of metal ions was further studied using potentiometric spectroscopy, as well as time-resolved laser-induced spectroscopy with Cm(III) and Eu(III) in aqueous HClO 4 and HNO 3 media. Conditional stability constants of each formed species were determined. In the HClO 4 system, Cm(III) formed three species (1:1, 1:2, and 1:3) through the stepwise addition of a single BTPOA molecule. On the other hand, in HNO 3 , Cm(III) formed two 1:2 complexes and one 1:3 complex, while the stepwise formation of three species was observed for Eu(III). The stability constants are comparable to the values for SO 3 -Ph-BTP. The radiolytic behavior of BTPOA was also investigated using electron pulse irradiation measurements to determine absolute rate coefficients (k) under ambient temperature conditions for the reaction of BTPOA with typical UNF reprocessing radical radiolysis products the hydrated electron (e aq − , k = (1.60 ± 0.02) × 10 10 M −1 s −1 ), the hydrogen atom (H • , k = (2.17 ± 0.03) × 10 9 M −1 s −1 ), and hydroxyl ( • OH, k = (6.95 ± 0.06) × 10 9 M −1 s −1 ) and nitrate (NO 3 • , k = (0.37 ± 0.02) × 10 7 M −1 s−1) radicals. These rate coefficients indicate that the radiolytic longevity of BTPOA should increase with HNO 3 concentration, owing to the consumption of e aq − /H • , by nitrate anions, and the replacement of • OH by the less reactive NO 3 • .
Multi-Scale modeling of mercury (Hg) geochemical speciation and reactions has been performed by integrating atomistic quantum chemical calculations with continuum scale speciation models. Major progress has been made in the improvement of quantum chemical models to calculate critical thermodynamic data for Hg complexes in aquatic environments. Rapid and reliable quantum chemical approaches have been developed for calculating acid dissociation constants (pK a ) and stability constants (log K), with calculated mean unsigned errors of 0.5 and 1.5 log units, respectively for ligand molecules and Hg complexes. At the continuum scale, systematic analysis of uncertainty propagation in mercury (Hg) speciation modeling has been conducted and was used to identify environmental conditions under which thermodynamic constant uncertainties are significant and recommended to be accounted for. The integrated framework for multi-scale modeling of mercury geochemistry is open to the research community through the web-based multiscale modeling aqueous speciation resource, AQUA-MER. The improved quantum chemical approaches for thermodynamic constant calculations are accessible through AQUA-MER and can be used to provide the missing constants in the continuum scale speciation calculations. In addition to low molecular mass Hg complex speciation, modeling natural aquatic environments also involve the transport of high molecular weight dissolved organic matter (DOM) in reactive flows simultaneously with equilibrium and kinetic reactions. To this end, atomistic MD simulations were performed to capture the details of aggregation, mechanisms and distribution of functional groups in DOM at the molecular level. The elemental composition and calculated bulk properties of the DOM models are in close agreement with experimental measurements. A travel-time based reactive transport model in the hyporheic zone of stream corridors was established for the multicomponent Hg-DOM-S system and implemented through PFLOTRAN.
Actinide oxalates are chemical compounds important to nuclear industry, ranging from actinide separation in waste reprocessing, to production of specialty actinides, and to disposal of high level nuclear waste (HLW) and spent nuclear fuel (SNF). In this study, the solubility constants for Pr 2 (C 2 O 4 ) 3 •10H 2 O and Nd 2 (C 2 O 4 ) 3 •10H 2 O by performing solubility experiments in HNO 3 and mixtures of HNO 3 and H 2 C 2 O 4 at 23.0 ± 0.2 °C have been determined. The targeted starting materials, Pr 2 (C 2 O 4 ) 3 •10H 2 O and Nd 2 (C 2 O 4 ) 3 •10H 2 O, were successfully synthesized at room temperature using PrCl 3 , NdCl 3 and oxalic acid as the source metrials. Then, we utilized the targeted solubility-controlling phases to conduct solubility measurements. There was no phase change over the entire periods of experiments, demonstrating that Pr 2 (C 2 O 4 ) 3 •10H 2 O and Nd 2 (C 2 O 4 ) 3 •10H 2 O were the solubility-controlling phases in our respective experiments. Based on our experimental data, we have developed a thermodynamic model for Pr 2 (C 2 O 4 ) 3 •10H 2 O and Nd 2 (C 2 O 4 ) 3 •10H 2 O in the mixtures of HNO 3 and H 2 C 2 O 4 to high ionic strengths. The model for Pr 2 (C 2 O 4 ) 3 •10H 2 O reproduces well the reported experimental data for Pu 2 (C 2 O 4 ) 3 •10H 2 O, which are not utilized for the model development, demonstrating that Pr(III) is an excellent analog for Pu(III). Similarly, the model for Nd 2 (C 2 O 4 ) 3 •10H 2 O reproduces the solubility of Am 2 (C 2 O 4 ) 3 •10H 2 O and Cm 2 (C 2 O 4 ) 3 •10H 2 O. The Pitzer model was used for the calculation of activity coefficients. Based on the published, well established model for dissociation constants for oxalic acid and stability constants for actinide-oxalate complexes [i.e., AmC 2 O 4 + , and Am(C 2 O 4 ) 2 - ] to high ionic strengths, we have obtained the solubility constants (log 10 K 0 ) for the following reactions at 25 °C, Pr 2 (C 2 O 4 ) 3 •10H 2 O ⇌ 2Pr 3+ + 3C 2 O 4 2- + 10H 2 O(l). Nd 2 (C 2 O 4 ) 3 •10H 2 O ⇌ 2Nd 3+ + 3C 2 O 4 2- + 10H 2 O(l). to be -30.82 ± 0.30 (2σ), and - 31.14 ± 0.35 (2σ), respectively. These values for can be directly applied to Pu 2 (C 2 O 4 ) 3 •10H 2 O, Am 2 (C 2 O 4 ) 3 •10H 2 O and Cm 2 (C 2 O 4 ) 3 •10H 2 O. The model established for actinide oxalates by this study provides the needed knowledge with regard to solubilities of actinide/REE oxalates at various ionic strengths, and is expected to find applications in many fields, including the geological disposal of nuclear waste and the mobility of REE under the surface conditions, as Pr 2 (C 2 O 4 ) 3 •10H 2 O and Nd 2 (C 2 O 4 ) 3 •10H 2 O can be regarded as the pure Pr and Nd end-members of deveroite, a recently discovered natural REE oxalate with the following stoichiometry, (Ce 1.01 Nd 0.33 La 0.32 Pr 0.11 Y 0.11 Sm 0.01 Pb 0.04 U 0.03 Th 0.01 Ca 0.04 ) 2.01 (C 2 O 4 ) 2.99 •9.99H 2 O. Regarding its importance in the geological disposal of nuclear waste, Am 2 (C 2 O 4 ) 3 •10H 2 O/Pu 2 (C 2 O 4 ) 3 •10H 2 O/Cm 2 (C 2 O 4 ) 3 •10H 2 O can be the source-term phase for actinides, as demonstrated by the instance in the disposal in clay/shale formations. This is exemplified by the stability of Am 2 (C 2 O 4 ) 3 •10H 2 O in comparison with Am(OH) 3 (am), Am(OH) 3 (s) and AmCO 3 (OH)(s) under the relevant geological repository conditions.