Next Article in Journal
Soft Switched Current Fed Dual Active Bridge Isolated Bidirectional Series Resonant DC-DC Converter for Energy Storage Applications
Next Article in Special Issue
Application of Cement-Based Materials as a Component of an Engineered Barrier System at Geological Disposal Facilities for Radioactive Waste—A Review
Previous Article in Journal
Machine Learning Assisted Chemical Process Parameter Mapping on Lignin Hydrogenolysis
Previous Article in Special Issue
Calculation of the Henry’s Constant and the Thickness of the Equilibrium Adsorption Layer of Radon in the Layer-by-Layer Measurement of the Sorbent Activity
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

DFT Prediction of Radiolytic Stability of Conformationally Flexible Ligands

Department of Chemistry, Lomonosov Moscow State University, Leninskie Gory 1 Bld. 3, 119991 Moscow, Russia
*
Author to whom correspondence should be addressed.
Energies 2023, 16(1), 257; https://doi.org/10.3390/en16010257
Submission received: 31 October 2022 / Revised: 20 December 2022 / Accepted: 21 December 2022 / Published: 26 December 2022
(This article belongs to the Special Issue Treatment of Radioactive Waste and Sustainability Energy)

Abstract

:
Radiolytic stability is one of the main requirements of the ligands for the reprocessing of spent nuclear fuel. The prediction of radiation stability based only on the 2D molecular structural formula allows us to accelerate and simplify the development of new ligands. Here, we used quantum chemistry to investigate the radiolytic behavior of water-soluble diglycolamides as one of the most popular ligands for spent nuclear fuel reprocessing. The accurate accounting of conformational mobility in the descriptors based on the Frontier Orbital Fukui theory allowed us to obtain a good correlation between theoretical and experimental data.

1. Introduction

Today, approximately 11% of the world’s electricity is produced with nuclear power, which is third in the world after coal and hydropower [1]. Due to the high efficiency of nuclear power plants, their total number in the world is growing every year [2]. This fact also relates to the accessibility of nuclear energy generation in various climatic zones, regardless of natural conditions. In the context of the world’s decarbonization strategy, the nuclear industry is able to reduce carbon dioxide emissions into the atmosphere.
However, a critical problem of the nuclear industry is the accumulation of radioactive waste and spent nuclear fuel (SNF). The reprocessing of SNF first allows for the return of uranium and plutonium to the nuclear fuel cycle. In addition, an important step is the extraction of minor actinides (Am, Cm, and Np), which determines the requirements for storage facilities and repositories in the long term [3]. Despite being extracted from the SNF, these metals also need to be separated from each other and from the lanthanides, which possess similar chemical behavior.
A family of diglycolamides (DGAs) is a promising class of organic ligands for various strategies of SNF reprocessing using liquid-liquid extraction. N,N,N0,N0-tetraoctyldiglycolamide (TODGA) plays a crucial role in numerous partitioning processes for advanced fuel cycles, including DIAMEX [4], 1cycle-SANEX [5], and GANEX [6]. A replacement of n-octyl fragments by a short branch with 1–2 carbon atoms makes the molecule water-soluble [7]. These DGAs are widely used in the ALSEP [8] and EXAm [9] processes. On the one hand, the DGAs are the focus of a large number of radiochemical investigations that provide the theoreticians with a sufficient amount of experimental information to be placed in the basement of new theories. On the other hand, it also raises questions about the origins of the ligands’ efficiency, selectivity, and radiation stability.
Due to the high activity of the separated elements, the radiation stability of the extractants is one of the most important requirements for the ligands. The radiolytic degradation of water-soluble DGA molecules has been investigated experimentally and theoretically in numerous studies [10,11,12,13]. The summary information on radiolytic products is presented in Table S1 (Supplementary information). In general, the key points of the DGAs radiolysis are: (a) the main way of the DGA degradation process is the breaking of the C-O bond except TMDGA; (b) the C breaking in the DGA skeleton is not observed; (c) the stability of molecules increases with the methylation of the ether bonds. The latter fact, at first glance, looks counterintuitive and does not correlate with the increase in the stability of the more branched radical, which could be assumed from the general knowledge of organic chemistry.
Theoretical modeling of the radiolysis process allows us to estimate the radiolytic stability of the molecules and suggest modifications in the ligand design process prior to experimental studies, thus reducing the doses of the number of staff. Furthermore, the search for new effective and radiolytically stable ligands may be significantly simplified and accelerated by computing the properties of the molecule based only on its structural formula before the synthesis [14].
Theoretical investigations of radiolysis in similar systems are based on the assumption that indirect radiolysis is dominant due to a low concentration of ligands in the solvent, which has also been experimentally confirmed [15,16,17]. Firstly, ionizing radiation generates radicals and ions from the solvent molecules. Next, the formed particles attack the ligand molecules. Therefore, to predict the ligand stability, we need to estimate its reactivity regarding the ions and radicals produced by solvent radiation. A set of quantum chemical approaches is used for the DGAs’ property simulation. Koubsky et al. [11,18,19] suggested using the Fukui Frontier Orbital Theory [20,21] and to calculate the electrophilic, nucleophilic, and radical Fukui functions to determine the most reactive sites in the molecule regarding the corresponding particles [22]. The authors also calculated the bond orders (using Mayer and natural population analyses) to identify the weakest bonds. The same approach was used in the work of Verlinden and colleagues [23], where they continued the discussion about less reactive hydrogen atoms near the ether bond regarding the N-alkyl substitutional groups in relation to the radical attack. Although the theoretical data in most of the cases described the experimental results qualitatively well, the theoretical quantum chemical technique was not used for the determination of all class tendencies and for the further prediction of the required modifications of the molecule.
Earlier, we showed the significance of taking into account the conformational mobility of the ligand when modeling its properties, including the quantitative estimation of the radiolitic stability [14,24]. The aim of this study is to optimize the approaches to the radiolysis theoretical investigation. We investigated a set of water-soluble DGA molecules’ radiation stability by quantum chemical approaches, using the experimental data from [10] for theory validation.

2. Materials and Methods

The general idea of the model building is based on the concept of local reactivity, or the so-called conceptual density functional theory (DFT). According to the concept, we can predict the reactivity of the molecule or its parts based on the electron density distribution and alteration. We considered two regions of interest: the electron density condensed on atoms and between them. The first region can be described with a local reactivity index, showing its probability to perform as the center of a chemical reaction. The second one can be described as a non-integer variant of the well-known bond order. The investigation of both regions can be performed during the analysis of the Schrodinger equation solution (e.g., using the Kohn-Sham approach).
We used the local reactivity descriptors (atom indices) based on the Fukui border orbitals theory to determine the most reactive atom. The indices show the variation of local electron density as a reaction of adding or removing an electron from the whole system. Due to the large number of possible particles created during the irradiation in the solution, we preferred to use a condensed dual descriptor (CDD) as the local reactivity index [25,26]. Moreover, it was observed that the influence of the orbital relaxation was smaller in the case of the CDD calculation, and hence this descriptor is more accurate than the Fukui indices [27]. The CDD is a linear combination of condensed Fukui indices and is defined as follows:
Δ f A = 2 q N A q N + 1 A q N 1 A
The greater the reactivity of the atom is, the greater the absolute value of the CDD becomes, while the sign of the CDD shows the type of the attacking particle.
We also estimated the influence of conformation errors and steric factors on the calculated indices and bond orders.
All initial molecules were created from their 2D-molecular formulas. We performed three single-point calculations for each molecule with charges −1, 0, +1 for the condensed dual descriptor calculations by using ORCA 5.0 [28]. We used the natural population analysis for the bond order and charge calculating realized in the Janpa package [29]. The SASA (Solvent-Accessible Surface Area) approach realized in the MDTraj library [30] was used to estimate the steric accessibility of different sites in molecules. This method was based on the Golden section helix algorithm to generate the sphere with a radius equal to the sum of the atom’s and the attacking particle’s van der Waals radii. A proton was chosen as the smallest possible attacking particle. We introduced a new descriptor, SASAk, which was equal to the SASA for atom k in a molecule divided on all the surface of atom k considering the attacking particle. Since the main radiolytic products of DGAs were produced by C-C, C-O and C-N bonds, all electronic properties were calculated for non-hydrogen atoms.
We estimated molecule electronic properties in two ways—for the optimized molecule geometry (opt-approach) and for the confrontation set (conf-approach):
  • “Opt-approach”: We provided the optimization of molecules in two steps: (1) conformational analysis by the PM7 semiempirical Hamiltonian [31] in the MOPAC2016 package [32]; (2) the obtained conformation was optimized in the ORCA package on the D3-B3LYP/def2-TZVPP RIJCOSX [33,34,35,36,37] level of theory.
  • “Conf-approach”: For the conformation error estimation, we selected 20 conformations for each molecule by the following method. We generated 200 conformations for each structure using the rdkit open library. Next, the energy of each molecule considering solvation energy (where the solvent is water) was computed using the GFN2-xTB [38]. We chose the conformation with minimal energy and then we selected 19 other conformations with energy less than the minima energy adding 20 kcal/mol (to select only those conformations that thermodynamically can exist in solution) and with the maxima of the RMSD of atom coordinates (to cover the potential energy surface as widely as possible). We calculated the electronic structure of each conformation with charges −1, 0, +1 in the ORCA at the same D3-B3LYP/def2-TZVPD RIJCOSX level of theory. The MBOs, NBOs and CDD values were averaged for each molecule over all conformations, considering the probability of this conformation existing in the solution. This probability was calculated based on the Gibbs energy obtained by XTB for each conformation. It should be mentioned that despite the long description of the process, its automatization together with the absence of DFT geometry optimization led to significant time saving compared to the first one. It also excluded the human factor and solved the issue with the symmetric atoms.

3. Results and Discussion

The general structure of DGAs investigated in this study with marked atoms is presented in Figure 1. First, we compared the experimental results of the main degradation products of the molecules with the theoretical ones predicted based on the bond orders obtained by different methods. We also compared each molecule using the CDD and SASA calculations with experimental values of the radiolysis constant.
Concerning the applied side of the issue, the calculation of Mayer bond orders is much faster compared to the natural ones, thus, these calculations are the most common, and the construction of the computational approach based on MBOs would be time-saving. Calculated Mayer and Natural bond orders for opt- and conf- approaches are presented in Figure 2 (and Table S1, supplementary information). First, the MBOs obtained for the optimized geometries presented abnormally small values in all cases (0.56 for TMGA, 0.49 for TEDGA, 0.21 for MeTEDGA and 0.19 for Me2TEDGA). Such bond orders rather indicated the absence of covalent bonds between atoms, which was impossible in our case. Nonetheless, the comparison of the MBOs revealed the lowest bond orders (and predicted the leakage) of the ether bond for the TMDGA molecule and the elimination of the methyl group (by C3-C18 bond breaking). Since these processes were not observed in the experiments, we concluded that the approach was inapplicable.
Moreover, we noted the significant difference between the symmetric bonds in some cases, which also could not be explained by anything other than the imperfection of the proposed computational method. Though a flexible molecule could be stable in non-symmetric conformations, the overall properties of the molecule should be a result of the properties of the whole set of stable 3D structures. The conformationally corrected property, in its turn, should generally reproduce the symmetry of the structure.
The MBOs determined by using the conf-approach did not show extremely low values. Furthermore, the bond orders were very similar for the symmetrical atoms in the molecules. However, again, the most reactive sites were not obtained correctly. In the case of TMDGA, the method predicted the C2-C3 bond breaking, which was not observed experimentally. The removal of the methyl group was again calculated for the Me2TEDGA molecule. Thus, we finally excluded this approach from further investigation.
Unfortunately, the NBOs simulated for the opt-approach still predicted the removal of the methyl group for the MeTEDGA as well as C2-C3 and C5-C6 bond breaking for the Me2TEDGA molecules. However, the conf-approach showed results closest to the experimental ones. We predicted ether bond breaking in all cases with good values of bond orders. Meanwhile, the C3-O4 and C5-O4 bonds were strong in the TMDGA substances. At this step, we may conclude that the analysis of bond leakage (even when using the accurate computational approach and considering the conformational mobility) cannot be used solely as a predictive tool for the estimation of the radiolitic stability. It can be explained by comparing the probabilities of the irradiation interaction with the ligand in the solution. Considering the volumes of the ligand molecule and the solvent one, as well as the typical concentrations in the radiochemical processes, the probabilities ratio can be estimated as nearly 102–103, and the solvent irradiation appeared to be the most probable one. The computational approach should preferably consider the probable reactions with the products of the solvent radiolysis rather than the direct ligand degradation. Thus, we decided to try other descriptors for the DGAs radiolytic behavior explanation using the NPA analysis and conf-approach.
The proposed extended approach was based on the following assumptions; We considered the radiolitic stability as a reciprocal of the reactivity of the atoms, i.e., the probable targets of the solvent fragments. We calculated the CDD and SASAk values for each molecule in each conformation and then averaged CDD × SASAk values, considering the probability of encountering this conformation in the water solution. The results are presented in Table 1.
First, it is important to note that the calculated CDD × SASAk values for the C3 atom in the MeTEDGA molecule as well as the C3 and C5 positions in the Me2TEDGA structure were equal to zero. This fact meant that these molecular sites were not sterically accessible for the attacking particle. Thus, the DGA’s methylation led to the ether bonds’ protection, and therefore the stability of molecules increased in the row TEDGA < MeTEDGA < Me2TEDGA.
Moreover, the CDD × SASAk values were decreased for the C2 and C6 atomic positions. The attack of these sites leads to the breaking of bonds between the nitrogen and carbonyl atoms in the structure. While the products of N1-C2 and C6-N7 bond breaking were presented in the experimental data of radiolytic degradation, the decrease in C2 and C6 CDD × SASAk values in the rows TMDGA, TEDGA, MeTEDGA, Me2TEDGA was by accompanied by the increase in radiation stability in this row.
The CDD × SASAk values were small for C10–C17 for all molecules with ethyl groups. This fact allowed us to make a suggestion about the same tendency for all carbon atoms in the long N-alkyl branches in lipophilic DGAs (such as TODGA, TEHDGA, etc.). Therefore, the main radiolytic degradation processes were taking place in positions 1–7 for all the DGAs molecules except TMDGA. This result was confirmed by the results of gamma irradiation experiments for lipophilic DGAs, where no products of the N-alkyne chain rupture were observed [39,40,41]. Moreover, the increase in radiation stability in the case of methylation in the C3 and C5 positions was validated by the TODGA and its methylating derivative, gamma-radiolysis [42].
Next, we applied the conformational reactivity descriptor (CRD) presented in our previous work [12]. The calculated CRD descriptor appeared to be linearly correlated with the experimental values of the dose constants (Table S3, Figure 3). While the TMDGA molecule has the lowest radiolytic stability, this ligand has the highest value of CRD. The radiolysis dose constants obtained by experimental investigation decrease in the row TMDGA > TEDGA > MeTEDGA > Me2TEDGA, and the CRD has the same behavior in this series. The coefficient of linear determination (R2) between the CRD and dose constant values is equal to 0.82. This fact validates the possibility of using the CRD descriptor for the estimation of the organic molecule’s radiolytic stability.
We may also note the decrease in model accuracy compared to our previous work, where the R2 value was 0.97. Though the average flexibility of the molecules was the same, here, we used automated subsampling together with semiempirical quantum chemistry in the first stage instead of the computationally expensive DFT calculations. Thus, we may consider the current computational scheme as much more helpful for the ligand design, as it reduces the required computational time several times.

4. Conclusions

In summary, we showed that the previously proposed approach using CDD × SASA and CRD as descriptors of reactivity allows us to obtain a qualitative and quantitative estimation of the organic ligand radiolytic stability. In addition, our method made it possible to critically reduce the human factor in determining the molecule’s electronic properties by starting automated modeling from a 2D structural formula. It should be noted that, despite the long description of the process of the conf-approach calculation, its automation, combined with the lack of optimization of the DFT geometry, leads to significant time savings compared to the traditional quantum chemistry way of geometry optimization.
We also showed that the methods using only optimized geometry for the molecule electronic structure calculation do not allow to simulate the properties of conformationally flexible molecules accurately. In particular, the bond orders calculated in this way do not correlate with the experimentally detected radiolysis products. Regarding the DGA synthesis, the fact that the radiolytic stability grows during the methylation of the ether bonds becomes the general recommendation for the new ether-containing ligands development. It is the steric protection of the C-O bond that leads to a decrease in the probability of its breaking despite the seemingly more stable radicals formed. We could also propose an increase in radiolytic stability after increasing the lengths of the alkyl chains. Though it is important to control the solubility of DGA in the corresponding solvent, these recommendations remain the same for water- and organic-soluble molecules.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/en16010257/s1, Table S1: Summary of radiolysis products of water-soluble DGAs; Table S2: Calculated Mayer and natural bond orders for TMDGA, TEDGA, MeTEDGA and Me2TEDGA molecules by opt- and conf- approaches; Table S3: Dose constants for the investigated molecules and CRD value (0.85 quantile of CDD × SASAk).

Author Contributions

Conceptualization, A.S. and V.P.; methodology, A.M.; software, M.Y.; investigation, A.S.; writing—original draft preparation, A.S.; writing—review and editing, A.M. and V.P. All authors have read and agreed to the published version of the manuscript.

Funding

The research was supported by RSF (grant 19-73-20115).

Data Availability Statement

Not applicable.

Acknowledgments

The research was carried out using the equipment of the shared research facilities of HPC computing resources at Lomonosov Moscow State University.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Mathew, M.D. Nuclear energy: A pathway towards mitigation of global warming. Prog. Nucl. Energy 2022, 143, 104080. [Google Scholar] [CrossRef]
  2. Alam, F.; Sarkar, R.; Chowdhury, H. Nuclear power in emerging economies and human resouce Development: A review. Energy Procedia 2019, 160, 3–10. [Google Scholar] [CrossRef]
  3. Taylor, R.F. Reprocessing and Recycling of Spent Nuclear Fuel; Woodhead Publishing: Cambridge, UK, 2015; ISBN 978-1-78242-212-9. [Google Scholar]
  4. Magnusson, D.; Christiansen, B.; Glatz, J.P.; Malmbeck, R.; Modolo, G.; Serrano-Purroy, D.; Sorel, C. Demonstration of a TODGA based extraction process for the partitioning of minor actinides from a PUREX raffinate part III: Centrifugal contactor run using genuine fuel solution. Solvent Extr. Ion Exch. 2009, 27, 26–35. [Google Scholar] [CrossRef]
  5. Wilden, A.; Schreinemachers, C.; Sypula, M.; Modolo, G. Direct Selective Extraction of Actinides (III) from PUREX Raffinate using a Mixture of CyMe4BTBP and TODGA as 1-cycle SANEX Solvent. Solvent Extr. Ion Exch. 2011, 29, 190–212. [Google Scholar] [CrossRef]
  6. McLachlan, F.; Wilden, A.; Brown, J.; Modolo, G.; Taylor, R.; Sarsfield, M. 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]
  7. Sasaki, Y.; Suzuki, H.; Sugo, Y.; Kimura, T.; Choppin, G.R. New water-soluble organic ligands for actinide cations complexation. Chem. Lett. 2006, 35, 256–257. [Google Scholar] [CrossRef]
  8. Gelis, A.V.; Lumetta, G.J. Actinide lanthanide separation process—ALSEP. Ind. Eng. Chem. Res. 2014, 53, 1624–1631. [Google Scholar] [CrossRef]
  9. Vanel, V.; Bollesteros, M.-J.; Marie, C.; Montuir, M.; Pacary, V.; Antégnard, F.; Costenoble, S.; Boyer-Deslys, V. Consolidation of the EXAm Process: Towards the Reprocessing of a Concentrated PUREX Raffinate. Procedia Chem. 2016, 21, 190–197. [Google Scholar] [CrossRef]
  10. Wilden, A.; Mincher, B.J.; Mezyk, S.P.; Twight, L.; Rosciolo-Johnson, K.M.; Zarzana, C.A.; Case, M.E.; Hupert, M.; Stärk, A.; Modolo, G. Radiolytic and hydrolytic degradation of the hydrophilic diglycolamides. Solvent Extr. Ion Exch. 2018, 36, 347–359. [Google Scholar] [CrossRef]
  11. Koubský, T.; Luštinec, J. Application of quantum mechanical simulations for studying the radiolytic stability of prospective extractants in the nuclear fuel cycle. J. Radioanal. Nucl. Chem. 2018, 318, 2407–2413. [Google Scholar] [CrossRef]
  12. Roscioli-Johnson, K.M.; Zarzana, C.A.; Groenewold, G.S.; Mincher, B.J.; Wilden, A.; Schmidt, H.; Modolo, G.; Santiago-Schübel, B. A Study of the γ-Radiolysis of N,N-Didodecyl-N′,N′-Dioctyldiglycolamide Using UHPLC-ESI-MS Analysis. Solvent Extr. Ion Exch. 2016, 34, 439–453. [Google Scholar] [CrossRef]
  13. Ravi, J.; Selvan, B.R.; Venkatesan, K.A.; Antony, M.P.; Srinivasan, T.G.; Vasudeva Rao, P.R. Radiolytic stability of di-2-ethylhexyl-dioctyl diglycolamide. Radiochim. Acta 2014, 102, 451–457. [Google Scholar] [CrossRef]
  14. Smirnova, A.; Mitrofanov, A.; Matveev, P.; Baygildiev, T.; Petrov, V. A search of a quantitative quantum-chemical approach for radiolytic stability prediction. Phys. Chem. Chem. Phys. 2020, 22, 14992–14997. [Google Scholar] [CrossRef] [PubMed]
  15. Sharma, S.; Ghosh, S.K.; Naik, D.B.; Dhami, P.S.; Sharma, J.N. Comparative studies on radiolytic degradation of deuterium labeled and unlabeled tributyl phosphates. J. Radioanal. Nucl. Chem. 2014, 302, 583–591. [Google Scholar] [CrossRef]
  16. Mincher, B.J.; Elias, G.; Martin, L.R.; Mezyk, S.P. Radiation chemistry and the nuclear fuel cycle. J. Radioanal. Nucl. Chem. 2009, 282, 645–649. [Google Scholar] [CrossRef]
  17. Galán, H.; Núñez, A.; Espartero, A.G.; Sedano, R.; Durana, A.; de Mendoza, J. Radiolytic Stability of TODGA: Characterization of Degraded Samples under Different Experimental Conditions. Procedia Chem. 2012, 7, 195–201. [Google Scholar] [CrossRef] [Green Version]
  18. Koubský, T.; Kalvoda, L. Application of ab-initio molecular electronic structure calculations of radiolytic and hydrolytic stabilities of prospective extractants. J. Radioanal. Nucl. Chem. 2015, 304, 227–235. [Google Scholar] [CrossRef]
  19. Koubský, T.; Fojtíková, J.; Kalvoda, L. Radical degradation stability of ether linkage in N,N,N’,N’-tetraoctyldiglycolamide and related organic extractants: A density functional study. Prog. Nucl. Energy 2017, 94, 208–215. [Google Scholar] [CrossRef]
  20. Stulberg, J. The physician quality reporting initiative—A gateway to pay for performance: What every health care professional should know. Qual. Manag. Health Care 2008, 17, 2–8. [Google Scholar] [CrossRef]
  21. Parr, R.G.; Yang, W. Density Functional Approach to the Frontier-Electron Theory of Chemical Reactivity. Am. Chem. Soc. 1984, 106, 4049–4050. [Google Scholar] [CrossRef]
  22. Yang, W.; Parr, R.G. Hardness, softness, and the fukui function in the electronic theory of metals and catalysis. Proc. Natl. Acad. Sci. USA 1985, 82, 6723–6726. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Verlinden, B.; Van Hecke, K.; Wilden, A.; Hupert, M.; Santiago-sch, B.; Egberink, R.J.M.; Verboom, W.; Kowalski, P.M.; Modolo, G.; Verwerft, M. Gamma radiolytic stability of the novel modified diglycolamide 2,20-oxybis(N,N- didecylpropanamide) (mTDDGA) for groupedactinide extraction. RSC Adv. 2022, 12, 12416–12426. [Google Scholar] [CrossRef] [PubMed]
  24. Mitrofanov, A.; Andreadi, N.; Matveev, P.; Zakirova, G.; Borisova, N.; Kalmykov, S.; Petrov, V. An(III)/Ln(III) solvent extraction: Theoretical and experimental investigation of the role of ligand conformational mobility. J. Mol. Liq. 2021, 325, 115098. [Google Scholar] [CrossRef]
  25. Morell, C.; Grand, A.; Toro-Labbé, A. New dual descriptor for chemical reactivity. J. Phys. Chem. A 2005, 109, 205–212. [Google Scholar] [CrossRef]
  26. Morell, C.; Grand, A.; Toro-Labbé, A. Theoretical support for using the Δf(r) descriptor. Chem. Phys. Lett. 2006, 425, 342–346. [Google Scholar] [CrossRef]
  27. Martínez-Araya, J.I. Why is the dual descriptor a more accurate local reactivity descriptor than Fukui functions? J. Math. Chem. 2014, 53, 451–465. [Google Scholar] [CrossRef]
  28. Neese, F.; Wennmohs, F.; Becker, U.; Riplinger, C. The ORCA quantum chemistry program package. J. Chem. Phys. 2020, 152, 224108. [Google Scholar] [CrossRef]
  29. Nikolaienko, T.Y.; Bulavin, L.A.; Hovorun, D.M. JANPA: An open source cross-platform implementation of the Natural Population Analysis on the Java platform. Comput. Theor. Chem. 2014, 1050, 15–22. [Google Scholar] [CrossRef]
  30. McGibbon, R.T.; Beauchamp, K.A.; Harrigan, M.P.; Klein, C.; Swails, J.M.; Hernández, C.X.; Schwantes, C.R.; Wang, L.P.; Lane, T.J.; Pande, V.S. MDTraj: A Modern Open Library for the Analysis of Molecular Dynamics Trajectories. Biophys. J. 2015, 109, 1528–1532. [Google Scholar] [CrossRef] [Green Version]
  31. Stewart, J.J.P. Optimization of parameters for semiempirical methods VI: More modifications to the NDDO approximations and re-optimization of parameters. J. Mol. Model. 2013, 19, 1–32. [Google Scholar] [CrossRef]
  32. Stewart, J.J.P. MOPAC2016. 2016. [Google Scholar]
  33. Stephens, P.J.; Devlin, F.J.; Chabalowski, C.F.; Frisch, M.J. Ab Initio Calculation of Vibrational Absorption and Circular Dichroism Spectra Using Density Functional Force Fields. J. Phys. Chem. 1994, 98, 11623–11627. [Google Scholar] [CrossRef]
  34. Weigend, F.; Ahlrichs, R. Balanced basis sets of split valence, triple zeta valence and quadruple zeta valence quality for H to Rn: Design and assessment of accuracy. Phys. Chem. Chem. Phys. 2005, 7, 3297–3305. [Google Scholar] [CrossRef] [PubMed]
  35. Weigend, F. Accurate Coulomb-fitting basis sets for H to Rn. Phys. Chem. Chem. Phys. 2006, 8, 1057. [Google Scholar] [CrossRef] [PubMed]
  36. Grimme, S.; Antony, J.; Ehrlich, S.; Kreig, H. A consistent and accurate ab initio parametrization of density functional dispersion correction (DFT-D) for the 94 elements H-Pu. Chem. Phys. 2010, 132, 154104. [Google Scholar] [CrossRef] [Green Version]
  37. Neese, F.; Wennmohs, F.; Hansen, A.; Becker, U. Efficient, approximate and parallel Hartree–Fock and hybrid DFT calculations. A ‘chain-of-spheres’ algorithm for the Hartree–Fock exchange. Chem. Phys. 2009, 356, 98–109. [Google Scholar] [CrossRef]
  38. Bannwarth, C.; Ehlert, S.; Grimme, S. GFN2-xTB—An Accurate and Broadly Parametrized Self-Consistent Tight-Binding Quantum Chemical Method with Multipole Electrostatics and Density-Dependent Dispersion Contributions. J. Chem. Theory Comput. 2019, 15, 1652–1671. [Google Scholar] [CrossRef] [Green Version]
  39. Zarzana, C.A.; Groenewold, G.S.; Mincher, B.J.; Mezyk, S.P.; Wilden, A.; Schmidt, H.; Modolo, G.; Wishart, J.F.; Cook, A.R. A Comparison of the γ-Radiolysis of TODGA and T(EH)DGA Using UHPLC-ESI-MS Analysis. Solvent Extr. Ion Exch. 2015, 33, 431–447. [Google Scholar] [CrossRef]
  40. Deepika, P.; Sabharwal, K.N.; Srinivasan, T.G.; Vasudeva Rao, P.R. Studies on the use of N,N,N′,N′-Tetra(2-Ethylhexyl) Diglycolamide (TEHDGA) for Actinide Partitioning II: Investigation on Radiolytic Stability. Solvent Extr. Ion Exch. 2011, 29, 230–246. [Google Scholar] [CrossRef]
  41. Sharma, J.N.; Ruhela, R.; Singh, K.K.; Kumar, M.; Janardhanan, C.; Achutan, P.V.; Manohar, S.; Wattal, P.K.; Suri, A.K. Studies on hydrolysis and radiolysis of tetra(2-ethylhexyl)diglycolamide (TEHDGA)/isodecyl alcohol/n-dodecane solvent system. Radiochim. Acta 2010, 98, 485–491. [Google Scholar] [CrossRef]
  42. Galán, H.; Zarzana, C.A.; Wilden, A.; Núñez, A.; Schmidt, H.; Egberink, R.J.M.; Leoncini, A.; Cobos, J.; Verboom, W.; Modolo, G.; et al. Gamma-radiolytic stability of new methylated TODGA derivatives for minor actinide recycling. Dalt. Trans. 2015, 44, 18049–18056. [Google Scholar] [CrossRef]
Figure 1. The structural formula of the investigated DGAs molecules.
Figure 1. The structural formula of the investigated DGAs molecules.
Energies 16 00257 g001
Figure 2. The breaking bonds from the experimental data (black lines) and the values of bond orders obtained by different calculation methods (color lines).
Figure 2. The breaking bonds from the experimental data (black lines) and the values of bond orders obtained by different calculation methods (color lines).
Energies 16 00257 g002
Figure 3. Correlation between the experimental radiolysis constants and the calculated CRD values.
Figure 3. Correlation between the experimental radiolysis constants and the calculated CRD values.
Energies 16 00257 g003
Table 1. CDD × SASAk values ×10−3 obtained for each non-hydrogen atom in the DGAs molecules.
Table 1. CDD × SASAk values ×10−3 obtained for each non-hydrogen atom in the DGAs molecules.
AtomN1C2C3O4C5C6N7O8O9C10
TMDGA−1.564.970.590.410.505.07−1.53−7.83−7.640.99
TEDGA−0.0081.66−0.060.01−0.051.81−0.001−6.39−6.450.20
MeTEDGA−0.311.1408.1−0.0060.79−0.02−1.26−8.460.15
Me2TEDGA−0.010.2503.100.27−0.45−1.45−5.580.16
AtomC11C12C13C14C15C16C17C18C19-
TMDGA1.091.001.12-------
TEDGA0.110.010.240.18−0.03−0.030.17---
MeTEDGA0.250.160.190.16−0.04−0.060.140.07--
Me2TEDGA0.230.290.320.01−0.04−0.150.340.180.32-
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Smirnova, A.; Yablonskiy, M.; Petrov, V.; Mitrofanov, A. DFT Prediction of Radiolytic Stability of Conformationally Flexible Ligands. Energies 2023, 16, 257. https://doi.org/10.3390/en16010257

AMA Style

Smirnova A, Yablonskiy M, Petrov V, Mitrofanov A. DFT Prediction of Radiolytic Stability of Conformationally Flexible Ligands. Energies. 2023; 16(1):257. https://doi.org/10.3390/en16010257

Chicago/Turabian Style

Smirnova, Anastasiia, Maksim Yablonskiy, Vladimir Petrov, and Artem Mitrofanov. 2023. "DFT Prediction of Radiolytic Stability of Conformationally Flexible Ligands" Energies 16, no. 1: 257. https://doi.org/10.3390/en16010257

APA Style

Smirnova, A., Yablonskiy, M., Petrov, V., & Mitrofanov, A. (2023). DFT Prediction of Radiolytic Stability of Conformationally Flexible Ligands. Energies, 16(1), 257. https://doi.org/10.3390/en16010257

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