Next Article in Journal
Similar Yet Different–Structural and Functional Diversity among Arabidopsis thaliana LEA_4 Proteins
Next Article in Special Issue
DFT Study of Molecular and Electronic Structure of Ca(II) and Zn(II) Complexes with Porphyrazine and tetrakis(1,2,5-thiadiazole)porphyrazine
Previous Article in Journal
Impact of Preheating Temperature on the Separation of Whey Proteins When Combined with Chemical or Bipolar Membrane Electrochemical Acidification
Previous Article in Special Issue
Preface from the Guest Editor of Special Issue “Quantum-Chemical Modeling and Design of Chelate and Macrocyclic Metal Complexes”
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Prediction of Binding Stability of Pu(IV) and PuO2(VI) by Nitrogen Tridentate Ligands in Aqueous Solution

1
Department of Chemistry, Nuclear & WMD Protection Research Center, Korea Military Academy, Seoul 01805, Korea
2
Department of Nuclear and Energy Engineering, Jeju National University, Jeju 63243, Korea
3
Department of Architectural Engineering, Hanyang University, Seoul 04763, Korea
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2020, 21(8), 2791; https://doi.org/10.3390/ijms21082791
Submission received: 15 March 2020 / Revised: 12 April 2020 / Accepted: 14 April 2020 / Published: 17 April 2020

Abstract

:
Plutonium has potential applications in energy production in well-controlled nuclear reactors. Since nuclear power plants have great merit as environmentally friendly energy sources with a recyclable system, a recycling system for extracting Pu from spent fuels using suitable extractants has been proposed. Pu leakage is a potential environmental hazard, hence the need for chemical sensor development. Both extractants and chemical sensors involve metal–ligand interactions and to develop efficient extractants and chemical sensors, structural information about Pu ligands must be obtained by quantum calculations. Herein, six representative nitrogen tridentate ligands were introduced, and their binding stabilities were evaluated. The tridentate L6, which contains tri-pyridine chelate with benzene connectors, showed the highest binding energies for Pu(IV) and PuO2(VI) in water. Analysis based on the quantum theory of atoms in molecular analysis, including natural population analysis and electron density studies, provided insight into the bonding characteristics for each structure. We propose that differences in ionic bonding characteristics account for the Pu-ligand stability differences. These results form a basis for designing novel extractants and organic Pu sensors.

Graphical Abstract

1. Introduction

Plutonium is very useful for energy production in a carefully controlled nuclear reactor. Currently, approximately 14% of global electricity production is derived from nuclear reactors [1]. Nonetheless, a major drawback of this strategy is the generation of a large amount of radioactive waste, usually a mixture of U, Pu, and other actinides [2]. For example, about 2000 metric tons of spent nuclear fuel are annually generated in the United States [3]. Moreover, the Oklo natural nuclear reactor in Gabon went critical around 2 billion years ago. According to Gauthier-Lafaye’s paper, the experimental data proved retention of fission products and plutonium in Gabon [4].
The nuclear power industry faces challenges related to the irradiation of nuclear fuels by used fuel recycled by a pyroprocessing system or PUREX (Pu and U Extraction) system. Liquid waste from these systems contains highly toxic actinides, including Pu, which must be effectively extracted for reuse via a closed-loop recycling program, such that it does not contaminate the environment [5,6,7,8,9]. The efficiency of such programs depends on effective recycling of the spent nuclear fuel by separating the fissile Pu, which is formed by neutron capture reactions and is intended for subsequent use in a fast neutron reactor. Leakage of even small amounts of nuclear waste is undesirable, because the radioactive materials pose environmental problems owing to their long half-lives [10]. Therefore, it is imperative to develop efficient extractants and chemical sensors, such as fluorescence sensors (which fluoresce upon binding with the target metal ion) [11,12]. Therefore, investigating the binding and structures formed between toxic actinides and organic ligands is essential for understanding various modes of complexation.
Actinide partitioning using selectively binding organophosphorus ligands like tributyl phosphate (TBP), octyl-(phenyl)-N,N-diisobutylcarbamoylmethyl phosphine oxide (CMPO), trialkyl phosphine oxide (TRPO), and diisodecyl phosphoric acid (DIDPA), along with diamide extractants like N,N′-dimethyl-N,N′-dioctylhexylethoxy malonamide (DMDOHEMA) and N,N′-dimethyl-N,N′-dibutyltetradecyl malonamide (DMDBTDMA) [13,14,15,16,17,18,19,20], has been proposed. The extractants usually have a nitrogen (soft) donor in their ligand structure. Therefore, a systematic study on various kinds of nitrogen donor ligands with Pu should be informative for further study on new types of extractants.
While numerous experimental studies on the complexation of actinides using the appropriate extractants in a PUREX system have been reported, there have been hardly any theoretical studies that provide structural information on such complexes in aqueous solution [21]. Molecular-level investigations via reliable and extensive quantum calculations such as density functional theory (DFT) can provide valuable information for the design of efficient extractants for actinide separation in the environment, as well as new organic ligands for actinide detection. Pu mainly exists in the +4 and +6 oxidation states in the environment (oxygenated system), and Pu(VI) is easily oxidized to PuO22+ in water [22].
Therefore, in this study, attempts have been made to understand the contrasting selectivity for Pu in aqueous solution using a soft donor, nitrogen, thus revealing structural information related to Pu complexation. Tridentate nitrogen donor ligands are employed as binders and extractants for lanthanides and actinides. The most commonly used N-donor ligands and those basic structures are utilized for this study (Figure 1). This research aims to clarify the differences in the binding capabilities of six tridentate soft nitrogen donor ligands with Pu (IV) and PuO2 (VI) in water, and to investigate their bonding characteristics.

2. Results and Discussion

2.1. Complexation and Bond Formation

Both Pu(IV) and PuO2(VI) form 1:1 complexes with tridentate ligands in a strongly solvating solvent like water to afford nonahydrates. The inner coordination sphere of the metal in the complexes comprises water molecules. For PuO2(VI), since the two oxygens occupy axial positions, the ligands interact with Pu in the equatorial plane, where the Pu can accommodate 4 more electron donors. Thus, the 1:1 complexes of Pu(VI) and Pu(IV) are filled with two and six water molecules, respectively [23]. Therefore, because of the lack of information about the target compound, we assessed the quality of the level of theory which was implemented and the same level of theory except the 6-31+G(d, p) for N, O, C, and H, respectively. Recently, the average Pu-Namin and Pu-Namido bond distances in Pu(IV)-triamidoamine(TREN) were measured to be 2.577 and 2.225 Å, respectively [24]. Quantum calculations on Pu(IV)-triamidoamine(TREN) without TIPS (Triisoproplysilane) were performed by using both of levels of theory. The average bond distances were found to be 2.606 (2.623) and 2.195 (2.199) Å, which is a slight overestimation (less than 1.5%) of the experimental data (structure and coordinates are presented in Supplementary Materials, Figure S2). Another study with reliable and extensive quantum calculations using a TPSSH functional with the Stuttgart small core ECP and def2-TZVP basis set on plutonium complexes has been carried out [25]. Therefore, the Pu(IV)-triamidoamine(TREN) structure was also optimized using the same level of theory, which also showed almost the same overestimation as the previous one (structure and coordinates are presented in Supplementary Materials, Figure S2). Therefore, here we performed the theoretical study at both levels of theory, both the computationally less expensive and the more expensive. With the same levels of theory, the optimized structures of each complex with water molecules are presented in Figure 2, and the average bond lengths are listed in Table 1. The Pu(IV)–N bonds are shorter than the corresponding PuO2 (VI)–N bonds, because the interaction between Pu(IV) and N is stronger than that in the other complexes in both calculations. This difference could stem from the fact that Pu(VI) forms a PuO22+ complex with weak positive charges on the central metal ion, which decrease the strength of the ionic bonding. The Pu–O (oxygen in the coordinating water molecule) bond length shows a different trend, because the oxidized structure of Pu(VI) has less steric hindrance between ligands, though the atomic charges of Pu(VI) and O are smaller than others. This hypothesis is well supported by the NPA analysis discussed in the supporting information.

2.2. Comparison of Stability

The reactions below denote the complex formation reaction for each Pu oxidation state.
Lx + Pu(IV) + 6H2O→PuLx(H2O)6
Lx + PuO22+(VI) + 2H2O→PuO2Lx(H2O)6
(x = 1, 2, 3, 4, 5, 6)
To quantitatively estimate the strength of stability for the complexes in aqueous solution, all possible spin multiplicities are considered, because most of the unpaired electrons reside in the 5f orbitals and lead to different multiplicities depending on the ligand binding. The lowest energy was observed for the triplet state in the case of the Pu(IV) structures, and for the quintet state in the case of the PuO2(VI) structures (see Table S4 in Supplementary Materials for details of the calculated energy for each spin multiplicity). Since the reactants with Pu (PuO2) and the water molecules are the same in each reaction, the energy difference from the lowest energy was used for calculating the stabilization energy for each functional: Ediff = E [PuL6(H2O)6] − (E [PuLx(H2O)6] − E [Lx]) and Ediff’ = E [PuL6(H2O)6] − (E [PuO2(H2O)2] − E [Lx]), (x = 1, 2, 3, 4, 5, 6). According to the analysis, L6 formed the most stable complex with both Pu(IV) and PuO22+(VI) species in both levels of theory (Table 2).
After applying two different functionals and basis sets, of all the optimized structures the L6 complex was the most stable, and its stability trend for the series of studied ligands is very similar to all except the L2 and L3 cases, indicating that complexation with the pyridine nitrogen on rigid scaffold leads to high stability and Pu complexation is strongly dependent on the ligand structures (Figure 3). It is possible that the ligand complexation properties are almost the same in both Pu(IV) and Pu(VI), based on QTAIM analysis of the bonding properties. L6 may be a better extractant for Pu(IV) than for PuO2(VI), as it showed a much bigger energy difference between Pu(IV) and L6 in aqueous solution. In addition, L6 would selectively extract Pu(IV) ions in a solution containing various other metal ions and ligand-like materials. Its stability of complexation shows its tendency to increase stability on more rigid structures, which implies a nitrogen-doped graphene structure could be a key structure toward extracting plutonium in water. Nitrogen doping technology has been developed in various fields [26], and it may be a good start for application of graphene-based materials in nuclear chemistry in the future.

2.3. Electronic Structure Analysis

Wiberg bond indices are a direct indicator of covalency [27]. Generally, Pu binding with a primary amine results in a larger Wiberg index, indicating more strongly covalent characteristics, which is the same trend as those of the other indicators described below. One can thus conclude that stronger ionic bonding is related to the stability of the complex formed with a secondary amine such as pyridine. The NPA charge values (Table S1) show that all the Pu atoms act as electron receptors, which carry a charge of less than 4 + (Pu(IV)) or 2 + (PuO2(VI)). Interestingly, the change of charges on Pu(IV) after forming complexes, from 4 to 1.43~1.54, is greater than the change of charge on PuO2(VI) cases, from 2 to 1.26~1.30. This is closely related to the bonding properties, which can be attributed to more covalent bonding forms in the Pu(IV) structures. This is strongly supported by the following topology analysis. Furthermore, the populations of the 7s, 7p, 5f, and 6d atomic orbitals of both structures differ from the canonical numbers (most electrons reside on 5f orbital), indicating that electron transfer from the ligand to Pu makes a major contribution to the ionic bond (see Supplementary Materials, Table S3).
Topology analysis was carried out to investigate each bond more closely. A BCP (3, −1) existed between two binding atoms (see Figure S1 and Table S2 in Supplementary Materials for BCPs). Generally, the electron density at the BCP is used as a standard for evaluating the bonding characteristics. A larger ρ(r) and smaller ∇2ρ(r) indicate a more covalent bond, while a smaller ρ(r) and larger ∇2ρ(r) indicate an ionic bond. This information clarified that Pu(IV)–N bonds had a stronger covalent nature than PuO2 (VI)–N bonds. In addition, the −V(r)/G(r) ratios (G(r): local kinetic energy density, V(r): local potential energy) are shown in Table 3. A − V(r)/G(r) ratio of less than 1 indicates a typical ionic bond, while a ratio greater than 1 represents a covalent bond. Accordingly, Pu(IV) showed stronger covalent bonding characteristics than PuO2(VI). There was no significant difference in the Pu–O bond between the two structures. Even though the two structures did not show exactly the same trend of ligand stability, we may conclude that stronger ionic bonding results in better stability of the complexes. Moreover, these studied bonding properties can be harnessed to extract specific oxidation states of plutonium in water and, in future research, this kind of investigation on bonding properties can be utilized in the development of new types of extractants by comparing different functional groups in ligands.

3. Computational Details

The computational level of B3LYP functional that is mainly considered in our study has recently been demonstrated to be reliable for actinide complexes [28,29,30]. The B3LYP functional and TPSSH were also utilized for geometry optimization in the gas phase, with tight energy and geometry convergence criteria, and without symmetry constraints. The ECP60MWB relativistic effective core potential and its associated basis set, developed by the Stuttgart–Cologne group, was chosen for describing plutonium, while 6-31G(d) for B3LYP functional and def2-TZVP for TPSSH functional were used for N, O, C, and H, respectively. These basis sets have been successfully utilized for describing actinide-organic complexes including plutonium in recent reports [31,32,33,34]. Spin-orbit effects were not considered, because the quasi-relativistic effective core potential takes into account the spin-orbit coupling effect and produces reliable electronic structures for studies on actinide complexes. To compare all structures under the same conditions, the initial ligand structures were required to be non-protonated. After optimization by B3LYP functional, the local minimum structure was confirmed by checking all positive frequencies. The theoretical calculations were performed using the Gaussian 09 software package [35]. Bader’s atoms-in-molecule (AIM) model, natural population analysis (NPA) of charge, and topological analysis of electron density can provide information about chemical bonds through bond critical point (BCP) positions. These methods were adopted here to understand the nature of bonding by estimating the bond strengths and electron density in the bonds. Analysis based on the quantum theory of atoms in molecules (QTAIM) was carried out using the Multiwfn program [36].

4. Conclusions

In this work, representative nitrogen tridentate ligands were introduced and their binding properties were investigated by using two different levels of theory to obtain information that will aid the design of plutonium extractants and metal ion sensors in aqueous solution. Pu(IV) and PuO22+(VI) with six tridentate ligands in solution were simulated. Interestingly, L6, which is a tri-pyridine chelate with benzene connectors, showed the highest binding energy for both Pu(IV) and PuO2(VI) ions in water. QTAIM analysis and NPA studies revealed the nature of the Pu–N and Pu–O bonds, and showed that highly ionic bonds account for the enhanced bonding strength and, consequently, the stability of these complexes. L6-like structures might be promising candidates for strong extractants and sensors in water, and L6-based structures may be used as platforms for designing novel extractants and Pu sensors. More importantly, since studies on radioactive actinides have safety limitations, alternative and safe theoretical studies can provide more insight on the nature of actinide chemistry before doing experiments, which will open up new opportunities to discover new chelating agents and sensors in the future.

Supplementary Materials

The following are available online at https://www.mdpi.com/1422-0067/21/8/2791/s1.

Author Contributions

Conceptualization, K.J. and S.M.W.; Validation, S.M.W. and S.B.; Formal analysis, H.J.J.; Investigation, K.J.; Data curation, H.J.J., K.J. and S.B.; Writing—Original draft preparation, K.J.; Writing—Review and editing, all authors. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by NRF research fund, grant number NRF-2020R1C1C1007888 and KISTI Supercomputing center, grant number KSC-2018-CHA-0015.

Conflicts of Interest

The authors declare no conflict of interest

References

  1. Bratt, D. Re-igniting the Atom: The Political Consequences of the Global Nuclear Revival. Whitehead J. Dipl. Int’l Rel. 2010, 11, 59. [Google Scholar]
  2. Hajima, R.; Hayakawa, T.; Kikuzawa, N.; Minehara, E. Proposal of nondestructive radionuclide assay using a high-flux gamma-ray source and nuclear resonance fluorescence. J. Nucl. Sci. Technol. 2008, 45, 441–451. [Google Scholar] [CrossRef]
  3. Mueller, M. 5 Fast Facts About Nuclear Energy | Department of Energy. Available online: https://www.energy.gov/ne/articles/5-fast-facts-about-spent-nuclear-fuel (accessed on 12 April 2020).
  4. Gauthier-Lafaye, F.; Holliger, P.; Blanc, P.L. Natural fission reactors in the Franceville basin, Gabon: A review of the conditions and results of a “critical event” in a geologic system. Geochim. Cosmochim. Acta 1996, 60, 4831–4852. [Google Scholar] [CrossRef]
  5. Series, P. The Actinide and Transactinide Elements; Springer: Dordrecht, The Netherlands, 1940; ISBN 9781402035555. [Google Scholar]
  6. Shirai, O.; Iizuka, M.; Iwai, T.; Suzuki, Y.; Arai, Y. Electrode reaction of plutonium at liquid cadmium in LiCl-KCl eutectic melts. J. Electroanal. Chem. 2000, 490, 31–36. [Google Scholar] [CrossRef]
  7. Serp, J.; Konings, R.J.M.; Malmbeck, R.; Rebizant, J.; Scheppler, C.; Glatz, J.P. Electrochemical behaviour of plutonium ion in LiCl—KCl eutectic melts. J. Electroanal. Chem. 2004, 561, 143–148. [Google Scholar] [CrossRef]
  8. Wolf, S.F.; Bates, J.K.; Buck, E.C.; Dietz, N.L.; Fortner, J.A.; Brown, N.R. Physical and chemical characterization of actinides in soil from Johnston Atoll. Environ. Sci. Technol. 1997, 31, 467–471. [Google Scholar] [CrossRef]
  9. Sengupta, A.; Mohapatra, P.K.; Iqbal, M.; Huskens, J.; Verboom, W. A highly efficient solvent system containing functionalized diglycolamides and an ionic liquid for americium recovery from radioactive wastes. Dalton Trans. 2012, 41, 6970. [Google Scholar] [CrossRef]
  10. Slater, J. Nuclear Wastelands: A Global Guide to Nuclear Weapons Production and its Health and Environmental Effects. BMJ 1995, 311, 1100–1101. [Google Scholar] [CrossRef]
  11. Wen, J.; Dong, L.; Tian, J.; Jiang, T.; Yang, Y.Q.; Huang, Z.; Yu, X.Q.; Hu, C.W.; Hu, S.; Yang, T.Z.; et al. Fluorescent BINOL-based sensor for thorium recognition and a density functional theory investigation. J. Hazard. Mater. 2013, 263, 638–642. [Google Scholar] [CrossRef]
  12. Sessler, J.L.; Melfi, P.J.; Seidel, D.; Gorden, A.E.V.; Ford, D.K.; Palmer, P.D.; Tait, C.D. Hexaphyrin(1.0.1.0.0.0). A new colorimetric actinide sensor. Tetrahedron 2004, 60, 11089–11097. [Google Scholar] [CrossRef]
  13. Reilly, S.D.; Gaunt, A.J.; Scott, B.L.; Modolo, G.; Iqbal, M.; Verboom, W.; Sarsfield, M.J. Plutonium(iv) complexation by diglycolamide ligands—Coordination chemistry insight into TODGA-based actinide separations. Chem. Commun. 2012, 48, 9732. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Gelis, A.V.; Lumetta, G.J. Actinide lanthanide separation process—ALSEP. Ind. Eng. Chem. Res. 2014, 53, 1624–1631. [Google Scholar] [CrossRef]
  15. Brown, J.; McLachlan, F.; Sarsfield, M.; Taylor, R.; Modolo, G.; Wilden, A. Plutonium Loading of Prospective Grouped Actinide Extraction (GANEX) Solvent Systems based on Diglycolamide Extractants. Solvent Extr. Ion Exch. 2012, 30, 127–141. [Google Scholar] [CrossRef]
  16. Ansari, S.A.; Pathak, P.; Mohapatra, P.K.; Manchanda, V.K. Chemistry of diglycolamides: Promising extractants for actinide partitioning. Chem. Rev. 2012, 112, 1751–1772. [Google Scholar] [CrossRef]
  17. Serrano-Purroy, D.; Baron, P.; Christiansen, B.; Glatz, J.P.; Madic, C.; Malmbeck, R.; Modolo, G. First demonstration of a centrifugal solvent extraction process for minor actinides from a concentrated spent fuel solution. Sep. Purif. Technol. 2005, 45, 157–162. [Google Scholar] [CrossRef]
  18. Madic, C.; Ouvrier, N. EUROPART: EUROpean research program for the PARTitioning of minor actinides from high active wastes arising from the reprocessing of spent nuclear fuels. Radiochimica Acta 2008, 96, 183–185. [Google Scholar] [CrossRef]
  19. Musikas, C. Solvent extraction for the chemical separations of the 5f elements. Inorg. Chim. Acta 1987, 140, 197–206. [Google Scholar] [CrossRef]
  20. Ozawa, M.; Koma, Y.; Nomura, K.; Tanaka, Y. Separation of actinides and fission products in high-level liquid wastes by the improved TRUEX process. J. Alloys Compd. 1998, 271–273, 538–543. [Google Scholar] [CrossRef]
  21. Pahan, S.; Boda, A.; Ali, S.M. Density functional theoretical analysis of structure, bonding, interaction and thermodynamic selectivity of hexavalent uranium (UO22+) and tetravalent plutonium (Pu4+) ion complexes of tetramethyl diglycolamide (TMDGA). Theor. Chem. Acc. 2015, 134, 41. [Google Scholar] [CrossRef]
  22. Rai, D.; Moore, D.A.; Rosso, K.M.; Felmy, A.R.; Bolton, H. Environmental mobility of Pu(IV) in the presence of ethylenediaminetetraacetic acid: Myth or reality? J. Solut. Chem. 2008, 37, 957–986. [Google Scholar] [CrossRef]
  23. Zeng, J.; Yang, X.; Liao, J.; Liu, N.; Yang, Y.; Chai, Z.; Wang, D. A computational study on the complexation of Np(V) with N,N,N′,N′-tetramethyl-3-oxa-glutaramide (TMOGA) and its carboxylate analogs. Phys. Chem. Chem. Phys. 2014, 16, 16536–16546. [Google Scholar] [CrossRef]
  24. Brown, J.L.; Gaunt, A.; King, D.M.; Liddle, S.; Reilly, S.D.; Scott, B.; Wooles, A. Neptunium and plutonium complexes with a sterically encumbered triamidoamine (TREN) scaffold. Chem. Commun. 2016, 52, 5428–5431. [Google Scholar] [CrossRef] [Green Version]
  25. Šulka, M.; Cantrel, L.; Vallet, V. Theoretical study of plutonium(IV) complexes formed within the PUREX process: A proposal of a plutonium surrogate in fire conditions. J. Phys. Chem. A 2014, 118, 10073–10080. [Google Scholar] [CrossRef]
  26. Wang, H.; Maiyalagan, T.; Wang, X. Review on recent progress in nitrogen-doped graphene: Synthesis, characterization, and its potential applications. ACS Catal. 2012, 2, 781–794. [Google Scholar] [CrossRef]
  27. Mayer, I. Bond order and valence indices: A personal account. J. Comput. Chem. 2007, 28, 204–221. [Google Scholar] [CrossRef]
  28. Wang, D.; van Gunsteren, W.F.; Chai, Z. Recent advances in computational actinoid chemistry. Chem. Soc. Rev. 2012, 41, 5836. [Google Scholar] [CrossRef]
  29. Wu, Q.Y.; Wang, C.Z.; Lan, J.H.; Xiao, C.L.; Wang, X.K.; Zhao, Y.L.; Chai, Z.F.; Shi, W.Q. Theoretical investigation on multiple bonds in terminal actinide nitride complexes. Inorg. Chem. 2014, 53, 9607–9614. [Google Scholar] [CrossRef]
  30. Gendron, F.; Páez-Hernández, D.; Notter, F.P.; Pritchard, B.; Bolvin, H.; Autschbach, J. Magnetic properties and electronic structure of neptunyl(VI) complexes: Wavefunctions, orbitals, and crystal-field models. Chem. A Eur. J. 2014, 20, 7994–8011. [Google Scholar] [CrossRef]
  31. Zhang, X.; Wu, Q.; Lan, J.; Yuan, L.; Xu, C.; Chai, Z.; Shi, W. Highly selective extraction of Pu (IV)and Am (III)by N,N′′-diethyl-N,N′′-ditolyl-2,9-diamide-1,10-phenanthroline ligand: An experimental and theoretical study. Sep. Purif. Technol. 2019, 223, 274–281. [Google Scholar] [CrossRef]
  32. Jian, T.; Yu, X.; Dan, D.; Albrecht-Schmitt, T.E.; Autschbach, J.; Gibson, J.K. Gas-Phase Complexes of Americium and Lanthanides with a Bis-Triazinyl Pyridine: Reactivity and Bonding of Archetypes for F-Element Separations. J. Phys. Chem. A 2020. [Google Scholar] [CrossRef]
  33. Wang, C.; Wu, Q.-Y.; Wang, C.-Z.; Lan, J.-H.; Nie, C.-M.; Chai, Z.-F.; Shi, W.-Q. Theoretical insights into selective separation of trivalent actinide and lanthanide by ester and amide ligands based on phenanthroline skeleton. Dalton Trans. 2020, 49, 4093–4099. [Google Scholar] [CrossRef]
  34. Wu, Q.Y.; Lan, J.H.; Wang, C.Z.; Zhao, Y.L.; Chai, Z.F.; Shi, W.Q. Understanding the interactions of neptunium and plutonium ions with graphene oxide: Scalar-relativistic DFT investigations. J. Phys. Chem. A 2014, 118, 10273–10280. [Google Scholar] [CrossRef]
  35. Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; Barone, V.; Mennucci, B.; Petersson, G.A.; et al. Gaussian 09; Gaussian, Inc.: Wallingford, CT, USA, 2009; pp. 2–3. [Google Scholar]
  36. Lu, T.; Chen, F. Multiwfn: A multifunctional wavefunction analyzer. J. Comput. Chem. 2012, 33, 580–592. [Google Scholar] [CrossRef]
Figure 1. Structures of the six tridentate nitrogen ligands.
Figure 1. Structures of the six tridentate nitrogen ligands.
Ijms 21 02791 g001
Figure 2. Optimized structures of Pu (IV) and PuO2(VI) with the suggested nitrogen tridentate ligands in aqueous solution (structural chemical formulas and coordinates of all structures are listed in SI).
Figure 2. Optimized structures of Pu (IV) and PuO2(VI) with the suggested nitrogen tridentate ligands in aqueous solution (structural chemical formulas and coordinates of all structures are listed in SI).
Ijms 21 02791 g002
Figure 3. Relative energy of each structure compared with the most stable complex (Pu–L6) using different functionals: (a) B3LYP, (b) TPSSH.
Figure 3. Relative energy of each structure compared with the most stable complex (Pu–L6) using different functionals: (a) B3LYP, (b) TPSSH.
Ijms 21 02791 g003
Table 1. Average bond length (Å) of Pu–tridentate structures. Each number shows the distance between Pu and N or O in the optimized structure.
Table 1. Average bond length (Å) of Pu–tridentate structures. Each number shows the distance between Pu and N or O in the optimized structure.
LigandB3LYPTPSSH
Pu(IV)PuO2(VI)Pu(IV)PuO2(VI)
Pu–NPu–OPu–NPu–OPu–NPu–OPu–NPu–O
L12.4882.5092.5362.4822.5012.5292.5222.473
L22.5442.5042.5582.4762.5612.5092.5512.467
L32.5092.5132.54732.4832.5742.5302.5362.478
L42.4572.5112.5162.4982.5602.5622.5022.489
L52.4982.5292.5172.5122.5702.5712.4992.533
L62.5182.5352.5242.5082.6002.5812.5062.535
Table 2. Calculated energy difference compared with L6 ligand energy (kcal/mol) for each structure using both functionals.
Table 2. Calculated energy difference compared with L6 ligand energy (kcal/mol) for each structure using both functionals.
LigandB3LYPTPSSH
Pu(IV) (Ediff)PuO2(VI) (Ediff’)Pu(IV) (Ediff)PuO2(VI) (Ediff’)
L157.25113.72454.77523.182
L298.95517.83597.51127.587
L381.10120.50678.19430.178
L431.2285.28839.21914.959
L513.3622.56411.5462.240
L60000
Table 3. Wiberg indices and −V(r)/G(r) ratios used for analyzing bond characteristics.
Table 3. Wiberg indices and −V(r)/G(r) ratios used for analyzing bond characteristics.
LigandPu(IV)PuO2(VI)
Wiberg Index−V(r)/G(r)Wiberg Index−V(r)/G(r)
Pu–NPu–OPu–NPu–OPu–NPu–OPu–NPu–O
L10.4140.2781.1210.9190.3690.3170.9290.894
L20.7940.7271.1700.9850.7160.6081.0530.958
L30.7740.5741.0480.9500.3550.3170.9281.275
L40.4200.2751.1430.9290.3490.3130.9210.898
L50.7340.5591.0030.9470.3480.3070.9230.897
L60.3890.2791.0490.9170.3450.3110.9211.650

Share and Cite

MDPI and ACS Style

Jeong, K.; Jeong, H.J.; Woo, S.M.; Bae, S. Prediction of Binding Stability of Pu(IV) and PuO2(VI) by Nitrogen Tridentate Ligands in Aqueous Solution. Int. J. Mol. Sci. 2020, 21, 2791. https://doi.org/10.3390/ijms21082791

AMA Style

Jeong K, Jeong HJ, Woo SM, Bae S. Prediction of Binding Stability of Pu(IV) and PuO2(VI) by Nitrogen Tridentate Ligands in Aqueous Solution. International Journal of Molecular Sciences. 2020; 21(8):2791. https://doi.org/10.3390/ijms21082791

Chicago/Turabian Style

Jeong, Keunhong, Hye Jin Jeong, Seung Min Woo, and Sungchul Bae. 2020. "Prediction of Binding Stability of Pu(IV) and PuO2(VI) by Nitrogen Tridentate Ligands in Aqueous Solution" International Journal of Molecular Sciences 21, no. 8: 2791. https://doi.org/10.3390/ijms21082791

APA Style

Jeong, K., Jeong, H. J., Woo, S. M., & Bae, S. (2020). Prediction of Binding Stability of Pu(IV) and PuO2(VI) by Nitrogen Tridentate Ligands in Aqueous Solution. International Journal of Molecular Sciences, 21(8), 2791. https://doi.org/10.3390/ijms21082791

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop