Study on Dynamic Column Behavior and Complexation Mechanism of DBS-Modified Crown Ether-Based Silica to 90Sr

A crown ether-loaded hybrid adsorbent suitable for the separation and enrichment of strontium from high-level liquid waste was synthesized. 4′,4′(5″)-di(tert-butylcyclohexano)-18-crown-6 (DtBuCH18C6) and its modifiers dodecyl benzenesulfonic acid (DBS) and 1-dodecanol were impregnated into silica-based polymer support. The hybrid adsorbent exhibited excellent Sr(II) selectivity ability, and effective chromatographic separation and recovery of Sr(II) from simulated high-level liquid waste could be achieved with a (DtBuCH18C6 + DBS + dodec)/SiO2-P packed column. The recovery rate of Sr(II) calculated based on the mass balance was approximately 99% and over 80% for the other coexisting metal ions. An appropriate increase in the concentration of Na-DTPA eluent was favorable to improve the efficiency of the elution process because of the increased complexation capacity of [DTPA]5− to Sr(II). The developed theoretical model can simulate the dynamic breakthrough curves of the material on the basis of short column data, thereby predicting the scale-up column of the practical operation. Density functional theory calculation was used to explore the action mechanism of DBS modifiers on the Sr(II) complexation process of crown ether groups. Two Sr(II) complexation isomeric models of DtBuCH18C6 were established, and the calculation results revealed a similar complexation ability. DtBuCH18C6 could form a stable Sr(II) complexation structure with DBS coordination, which further indicated that DBS could be a ligand to promote the Sr(II) adsorption ability of crown ether materials.


Introduction
Over the last decades, the rapid development of the nuclear energy industry has increased the demand for the enhanced efficiency of high-level liquid waste (HLLW) treatment [1].As a common β-emitter nuclide within HLLW, 90 Sr has a long half-life and a high heat release rate, thereby possibly weakening the mechanical stability of glass curing products and complicating the long-term geological disposal process of HLLW [2][3][4]. 90Sr removal from HLLW offers the benefits of decreasing the cost of geological disposal, improving the utilization of the repository space, and reducing the threat of radiation to the ecosystem and humans [5].Therefore, 90 Sr decontamination of HLLW is of great significance.
Various techniques have been reported for the application of Sr(II) removal from radioactive wastewater, including adsorption [6][7][8], ion exchange [9][10][11], extraction techniques [12,13], and membrane separation [14].The extraction method possesses the advantages of high selectivity and recovery and is suitable for continuous operation under high-level radiation conditions.Therefore, it has been proposed as the primary method for large-scale Sr(II) decontamination.However, the extraction process, especially the steps of solvent extraction, consistently suffers the drawbacks of being complicated, requiring multiple steps, having a high organic phase consumption, and resulting in secondary waste generation as a result of solvent hydrolytic and radiolytic degradation [15,16].Thus, the novel technology of extraction chromatography, which is based on solid adsorbent and has the advantages of simple handling and good mechanical and radiation resistance of the material, is gaining more attention [17].
Several solid adsorbents that integrate the organic extractant with solid supports to the application of Sr(II) extraction chromatography separation have been reported [18][19][20].Considering the strong selectivity chelation of crown ether compounds with alkaline earth metal cations, solid-extraction adsorbent with crown ether as the adsorption phase has become an attractive option.4 ,4 (5 )-di(tert-butylcyclohexano)-18-crown-6 (DtBuCH18C6), which polyether cavity size matches the Sr(II) ionic radius, possesses superior hydrolytic and radiation stability and can rapidly chelate with Sr(II) inside the aqueous solution, yielding a complex with strong stability.Therefore, it has become an efficient adsorption phase for the removal of 90 Sr from HLLW [21].In our previous work, a novel crown ether solid-extraction adsorbent was successfully synthesized by impregnating DtBuCH18C6 into the macroporous silica-based polymer composite carrier (SiO 2 -P) [22].The synthesized adsorbent exhibited excellent mechanical properties and chemical and radiation resistance, which was suitable for column separation.Moreover, 1-dodecanol and dodecyl benzenesulfonic acid (DBS) were added as the molecule modifiers to further widen the adsorbent applicable acidity range, and the adsorbent containing DBS better promoted the adsorption of Sr(II) compared to the DBS-free one, particularly addressing the defect of poor adsorption capacity in low-acidity conditions [23].The fundamental adsorption performance of adsorbent with DBS modification has been investigated in the literature, demonstrating that DBS works as a counter ion.However, DBS modification significantly improves the Sr(II) adsorption ability of the adsorbent, leading to difficulties in desorption using deionized water.The dynamic Sr(II) adsorption behavior of the adsorbent applied in the column separation process remains poorly understood.Thus, an adsorbent synthesized with DBS modification needs to be further investigated.
Other aspects we would like to explore are how DtBuCH18C6 supplies the Sr(II) selective separation ability within the nitrate system and how the DBS molecule modifiers influence the complexation process.In general, the metal ion separation ability of crown ethers could be explained by the compatibility of the polyether cavity size with the cation size and the hard and soft acids and bases principle-namely, the stability of the yielded complex mainly depends on the nature of the electron donors and the cation speciation in the aqueous solution [24,25].Several experimental investigations have demonstrated that the principle of Sr(II) adsorption is mainly controlled by the interaction of Sr(II) with the oxygen in crown ethers to form Sr-O bonds [23,26].Some theoretical studies have revealed the complexation mechanism of crown ethers with Sr(II) during the extraction process [27][28][29].However, nitrate has been identified as the possible ligand that participates in the reaction of Sr(II) extraction by DtBuCH18C6 from the solutions; the same phenomenon may be observed in DBS [30,31].Thus, the extraction mechanism of DtBuCH18C6 to Sr(II) with DBS as the ligand needs to be further theoretically studied.In addition, the scaling up of column experimental studies is necessary before adsorbents can be used industrially.However, multiple expensive and time-consuming pilot-scale experiments are usually required to obtain dynamic experimental results as accurately as possible.Therefore, analytical models such as Thomas, the Clark column, and the Wolborska model have been used to evaluate the dynamic adsorption performance of materials.These models are mainly used for fitting experimental data and can hardly accurately infer the effects caused by changes in operating conditions.Therefore, the establishment of a numerical calculation model for breakthrough curve simulations is crucial for optimizing the design of dynamic adsorption processes.
In this study, a DBS-modified crown adsorbent was prepared.We constructed the Sr(II) adsorption modeling framework of DtBuCH18C6 on the basis of the optimized electrostatic potential (ESP) and calculated the complexation of DtBuCH18C6 to Sr(II) with DBS coordination to explore the possible promotion mechanism of DBS decoration to Sr(II) adsorption of a hybrid crown adsorbent, with a view toward providing a theoretical basis for the design of more efficient adsorbent materials.Further, we investigated the dynamic breakthrough property of the adsorbent, and a theoretical model was established for the breakthrough curve prediction.A mechanistic model was chosen to theoretically describe the adsorption process with a mass conservation equation, kinetic equation, and adsorption isotherm.Approaches for scale-up techniques applied to the removal of 90 Sr from HLLW were established.Four reagents were utilized for the investigation of the adsorbent Sr(II) elution capacity, and dynamic column separation of Sr(II) from the simulated HLLW was conducted with the best eluent for further performance evaluation.

Preparation of Material
The (DtBuCH18C6 + DBS + dodec)/SiO 2 -P adsorbent (Figure 1) with DBS modification was synthesized following the rotary vacuum evaporation method [23].DtBuCH18C6 (25.5 wt%) and DBS were dissolved in dichloromethane and then mixed with SiO 2 -P particles in a glass flask and then stirred gently by a rotary evaporator.The SiO 2 -P particles were synthesized by immobilizing a 20% styrene-divinylbenzene (SDB) copolymer into the SiO 2 porous.Sr(NO 3 ) 2 was of reagent grade and purchased from Sinopharm Chemical Reagent Co., Ltd, Shanghai, China.Diethylenetriaminepentaacetic acid (DPTA) and diethylenetriaminepentaacetic acid pentasodium salt solution (Na-DTPA) were purchased from Sinopharm Chemical Reagent Co., Ltd, Shanghai, China.The simulated HLLW (SHLLW) containing 3 M HNO 3 was prepared, and the main initial concentration of SHLLW is summarized in Table 1.
electrostatic potential (ESP) and calculated the complexation of DtBuCH18C6 to S with DBS coordination to explore the possible promotion mechanism of DBS decora to Sr(II) adsorption of a hybrid crown adsorbent, with a view toward providing a theo ical basis for the design of more efficient adsorbent materials.Further, we investigated dynamic breakthrough property of the adsorbent, and a theoretical model was establis for the breakthrough curve prediction.A mechanistic model was chosen to theoretic describe the adsorption process with a mass conservation equation, kinetic equation, adsorption isotherm.Approaches for scale-up techniques applied to the removal of from HLLW were established.Four reagents were utilized for the investigation of the sorbent Sr(II) elution capacity, and dynamic column separation of Sr(II) from the si lated HLLW was conducted with the best eluent for further performance evaluation.

Preparation of Material
The (DtBuCH18C6 + DBS + dodec)/SiO2-P adsorbent (Figure 1) with DBS modifi tion was synthesized following the rotary vacuum evaporation method [23].DtBuCH1 (25.5 wt%) and DBS were dissolved in dichloromethane and then mixed with SiO2-P ticles in a glass flask and then stirred gently by a rotary evaporator.The SiO2-P parti were synthesized by immobilizing a 20% styrene-divinylbenzene (SDB) copolymer the SiO2 porous.Sr(NO3)2 was of reagent grade and purchased from Sinopharm Chem Reagent Co., Ltd, Shanghai, China.Diethylenetriaminepentaacetic acid (DPTA) and di ylenetriaminepentaacetic acid pentasodium salt solution (Na-DTPA) were purcha from Sinopharm Chemical Reagent Co., Ltd, Shanghai, China.The simulated HL (SHLLW) containing 3 M HNO3 was prepared, and the main initial concentration SHLLW is summarized in Table 1.

Batch Experiments
Static batch experiments, which included the Sr(II) adsorption selectivity estim and desorption property evaluation, were performed using a thermostatic water b

Batch Experiments
Static batch experiments, which included the Sr(II) adsorption selectivity estimate and desorption property evaluation, were performed using a thermostatic water bath shaker (Tokyo PIKAKIKA Co., Ltd., Tokyo, Japan) at 25 + dodec)/SiO 2 -P was added to an aqueous solution (5 mL).The equilibrium distribution coefficient (K d , mL/g) and the elution rate (E r , %) were calculated using the following equations: where C 0 and C e (mg/L) denote the concentrations of metal ions inside the aqueous solutions before and after equilibrium adsorption, respectively.V (mL) and m (g) represent the volume of aqueous solution and mass of the adsorbent, respectively.C d (mg/L) denotes the metal equilibrium concentration of Sr(II) inside the elution solution, and Q e (mg/g) is the adsorption capacity of the adsorbent.

Dynamic Column Operation
A series of dynamic column adsorption experiments was conducted using a glass column of φ10 mm × h 100 mm that was packed with around 2 g of (DtBuCH18C6 + DBS + dodec)/SiO 2 -P adsorbent and then fed with different solutions.The breakthrough curve was determined by supplying the Sr(II) aqueous solution (5 mM Sr in 3 M HNO 3 ) at a flow of 0.3 mL/min.Eluent with three concentrations of 0.05, 0.075, and 0.1 mol/L was fed at a feeding rate of 0.5 mL/min to further investigate the dynamic Sr(II) desorption property of the adsorbent.The chromatographic column separation of Sr(II) was prepared with 3 mL SHLLW feeding, followed by 50 mL of 3 M HNO 3 and 120 mL of 0.1 M Na-DTPA solutions.A fraction collector was used for effluent collection, and concentrations of various metal cations were measured by using a sequential plasma spectrometer (ICP-AES, Shimadzu Corporation, Kyoto, Japan).

Mechanistic Model of Column Adsorption
A theoretical model based on the mass balance principle, adsorption Freundlich model, and adsorption kinetic model was developed by our group to describe dynamic column separation processes [32].The model assumed that (a) the solution flowed through with a constant velocity, (b) the density of the adsorbate did not change during the adsorption process, and (c) the temperature inside the column remained uniform and identical.
Thus, the adsorption process inside the column bed can be represented using Equations ( 3) and ( 4): The external mass coefficient K fa was employed to summarize the factors that control the adsorption rate, and the different components inside the adsorption resin with different adsorption rates were considered, as exhibited in Equation (5).
where values of k and l are positive.
Based on Equation ( 5), the kinetic model can be obtained as follows: Toxics 2023, 11, 919 The Freundlich isotherm, which can accurately describe the adsorption mechanism, was chosen to be the adsorption equilibrium model, which is given as follows: The difference method was utilized on the right side of Equation (1) for the numerical calculations.The L (m) column bed along the flow direction was divided into n stages, with each stage having a height of ∆z = L/n.The feeding liquid in each stage was simulated by the layer sequence model with the n-stage fully mixed.
The efficiency of the calculation was improved by considering the value of ∆z to be similar to the adsorbent particle diameter d P .
The final model is described as follows: and the initial and boundary conditions are defined as follows: where C (mg/L) is the concentration of Sr(II) in the solutions, t (h) is the time, u (m/h) is the flow rate, ε is the bed porosity, z (m) is the axial height along the column length, m (kg/m 3 ) is the density of the adsorbent, k and l are positive constants, K fa is the external mass coefficient, and K F is the constant of the Freundlich model.

Density Functional Theory Calculation
The optimized structures and interaction energies of DtBuCH18C6 coordinated with Sr(II) in the presence and absence of the DBS modifier or nitrate ions were calculated using the m06-2x density functional as implemented in the Gaussian 16 program package [33,34].During the process of complex geometry optimization and the frequency calculation, the SDD basic set was used for Sr(II), and the rest of the atoms were performed at the def2-svp basic set level.The component of "quantitative analysis of molecular surface" in Multiwfn 3.8 was used to calculate the electrostatic potential of the optimized models in conjunction with VMD 1.9.3 to depict the electrostatic potential (ESP) pictures.Single-point energy calculation was performed with def2-tzvp in conjunction with a popular recessive solvent (COSMO, water) model for all the elements, and D3 dispersion correction was added to describe weak interactions.The correction of the zero-point energy and the basis set superposition error correction (BSSE) were involved in the calculation of the interaction energies, and the expression is as follows [35]: where ∆E represents the interaction free energy of the ion-crown ether complex.E complex , E DtBuCH18C6 , and E Sr represent the energies of the ion-crown ether complex, DtBuCH18C6, and Sr(II) ion, respectively.E BSSE represents the BSSE energy.

Influence of Coexisting Ions on Sr(II) Separation
The SHLLW with an acidity of 3 M HNO 3 was utilized for evaluating the adsorption selectivity of the (DtBuCH18C6 + DBS + dodec)/SiO 2 -P adsorbent, and the results are summarized in Figure 2. The separation factor (SF) was determined as K d , Sr /K d , M n+ .The K d of Sr(II) was higher than that of the other coexisting ions by up to 54.29 mL/g, thus indicating the excellent Sr(II) separation property of the adsorbent.As a result of the selective adsorption capacity of crown ether being mainly controlled by the matches of the polyether cavity size with metal cations, the removal of the adsorbent to Ba(II) and Na(I) ions could be explained with the similar size of them to Sr(II), indicating the effective competition between Ba(II) and Na(I) in the Sr(II) separation process of the adsorbent [36,37].The calculated K d of the adsorbent to Zr(IV) was 11.53 cm 3 /g.Considering the coexistence of Mo(V) inside the HLLW, the decreased concentration of Zr(IV) could be attributed to the coprecipitation of Zr(IV) with Mo(V) to form the Zr(MoO 4 ) 2 •2H 2 O precipitate [38,39].As a result, the adsorbent exhibited significant Sr(II) separation property from the simulated HLLW, indicating that (DtBuCH18C6 + DBS + dodec)/SiO 2 -P possessed excellent Sr(II) adsorption selectivity ability.

Influence of Coexisting Ions on Sr(II) Separation
The SHLLW with an acidity of 3 M HNO3 was utilized for evaluating the adsorption selectivity of the (DtBuCH18C6 + DBS + dodec)/SiO2-P adsorbent, and the results are sum marized in Figure 2. The separation factor (SF) was determined as Kd,Sr/Kd,M n+ .The Kd o Sr(II) was higher than that of the other coexisting ions by up to 54.29 mL/g, thus indicating the excellent Sr(II) separation property of the adsorbent.As a result of the selective ad sorption capacity of crown ether being mainly controlled by the matches of the polyethe cavity size with metal cations, the removal of the adsorbent to Ba(II) and Na(I) ions could be explained with the similar size of them to Sr(II), indicating the effective competition between Ba(II) and Na(I) in the Sr(II) separation process of the adsorbent [36,37].The cal culated Kd of the adsorbent to Zr(IV) was 11.53 cm 3 /g.Considering the coexistence o Mo(V) inside the HLLW, the decreased concentration of Zr(IV) could be attributed to the coprecipitation of Zr(IV) with Mo(V) to form the Zr(MoO4)2•2H2O precipitate [38,39].As a result, the adsorbent exhibited significant Sr(II) separation property from the simulated HLLW, indicating that (DtBuCH18C6 + DBS + dodec)/SiO2-P possessed excellent Sr(II) ad sorption selectivity ability.

Elution Behavior of Adsorbent with Different Eluents
The Sr(II) elution rates of various eluents to (DtBuCH18C6 + DBS + dodec)/SiO2-P are shown in Figure 3. Compared to the deionized water, nitric acid, and DTPA+NaOH mix ture, solutions of Na-DTPA exhibited the best recovery rate, up to 76.03%.According to our previous investigation, the adsorbed Sr(II) onto DtBuCH18C6-loaded resin could be eluted with deionized water [22].However, H2O showed an invalid elution ability with regard to the DBS-modified adsorbent.In addition, 1 mol/L HNO3 recovered only 45.90% of the adsorbed Sr(II) onto the (DtBuCH18C6 + DBS + Dodec)/SiO2-P adsorbent, and the low Sr(II) recovery rate would result in the poor reuse performance of the material.The elution behaviors of the deionized water and nitric acid could be explained as follows: the DBS modifier could offer the coordination anion for the Sr(II) complexation o DtBuCH18C6, thus leading to the low desorption efficiency of the deionized water and nitric acid to the DBS-modified adsorbent [30].

Elution Behavior of Adsorbent with Different Eluents
The Sr(II) elution rates of various eluents to (DtBuCH18C6 + DBS + dodec)/SiO 2 -P are shown in Figure 3. Compared to the deionized water, nitric acid, and DTPA+NaOH mixture, solutions of Na-DTPA exhibited the best recovery rate, up to 76.03%.According to our previous investigation, the adsorbed Sr(II) onto DtBuCH18C6-loaded resin could be eluted with deionized water [22].However, H 2 O showed an invalid elution ability with regard to the DBS-modified adsorbent.In addition, 1 mol/L HNO 3 recovered only 45.90% of the adsorbed Sr(II) onto the (DtBuCH18C6 + DBS + Dodec)/SiO 2 -P adsorbent, and the low Sr(II) recovery rate would result in the poor reuse performance of the material.The elution behaviors of the deionized water and nitric acid could be explained as follows: the DBS modifier could offer the coordination anion for the Sr(II) complexation of DtBuCH18C6, thus leading to the low desorption efficiency of the deionized water and nitric acid to the DBS-modified adsorbent [30].As a commonly used metal chelate, DTPA could coordinate with Sr(II) in the solu tions to form an annular structure, thus occupying the Sr(II) reaction sites and preventing the complexation reactions of Sr(II) with adsorbents [40].The main coordination group o DTPA to metal ions is [DTPA] 5− , which suggests that the agents of DTPA and Na-DTPA can be used for Sr(II) elution of the (DtBuCH18C6 + DBS + dodec)/SiO2-P adsorbent.Fo the mixed solution of DTPA with NaOH, NaOH was added to increase the dissolution o DTPA and adjust the main chelate ion concentrations inside the solution.However, the elution efficiency of Na-DTPA was higher than that of the DTPA+NaOH mixture.This phenomenon could be ascribed to the pH value of the DTPA+NaOH mixture solution which was 3.6, lower than the best Sr(II) complexation condition of pH = 12.The ionic concentration of [DTPA] 5− inside the Na-DTPA solution was higher than DTPA + NaOH which could offer more chelate sites for Sr(II) complexation.As a result, Na-DTPA exhib ited excellent Sr(II) elution performance, which could be used as the eluent for the dy namic desorption property investigation of the (DtBuCH18C6 + DBS + dodec)/SiO2-P ad sorbent.

Breakthrough Curves
Given its advantages of high separation efficiency and large volumes of batch treat ment capacity, the application of extraction chromatography to the region of HLLW de contamination has attracted increasing attention.However, the successful development o a chromatography method of Sr(II) separation from the HLLW needs the ion dynamic behavior inside the column bed to be described accurately.Thus, the dynamic adsorption performance of the (DtBuCH18C6 + DBS + Dodec)/SiO2-P adsorbent to Sr(II) at a flow rate of 0.3 mL/min was investigated.The results were simulated through the developed theo retical model.
Figure 4 presents the breakthrough curve of (DtBuCH18C6 + DBS + dodec)/SiO2-P as the function of the feeding volume.The plots exhibit a similar trend to the typical S shaped line, and the values of C/Co increase rapidly to 0.8 and then slowly reach 1.This finding demonstrates the excellent dynamic adsorption property and high efficiency o the adsorbent.The breakthrough point of 5% in column experiments is often used as a momentous indicator for evaluating the adsorption performance of adsorbents.The breakthrough capacity (B.Cap.) and total capacity (T.Cap.) were acquired when the Sr(II concentration in the effluent reached 5% (C/Co = 5%) and 100% (C/Co = 100%) of the initia pumping concentration, respectively.Table 2 lists  As a commonly used metal chelate, DTPA could coordinate with Sr(II) in the solutions to form an annular structure, thus occupying the Sr(II) reaction sites and preventing the complexation reactions of Sr(II) with adsorbents [40].The main coordination group of DTPA to metal ions is [DTPA] 5− , which suggests that the agents of DTPA and Na-DTPA can be used for Sr(II) elution of the (DtBuCH18C6 + DBS + dodec)/SiO 2 -P adsorbent.For the mixed solution of DTPA with NaOH, NaOH was added to increase the dissolution of DTPA and adjust the main chelate ion concentrations inside the solution.However, the elution efficiency of Na-DTPA was higher than that of the DTPA+NaOH mixture.This phenomenon could be ascribed to the pH value of the DTPA+NaOH mixture solution, which was 3.6, lower than the best Sr(II) complexation condition of pH = 12.The ionic concentration of [DTPA] 5− inside the Na-DTPA solution was higher than DTPA + NaOH, which could offer more chelate sites for Sr(II) complexation.As a result, Na-DTPA exhibited excellent Sr(II) elution performance, which could be used as the eluent for the dynamic desorption property investigation of the (DtBuCH18C6 + DBS + dodec)/SiO 2 -P adsorbent.

Breakthrough Curves
Given its advantages of high separation efficiency and large volumes of batch treatment capacity, the application of extraction chromatography to the region of HLLW decontamination has attracted increasing attention.However, the successful development of a chromatography method of Sr(II) separation from the HLLW needs the ion dynamic behavior inside the column bed to be described accurately.Thus, the dynamic adsorption performance of the (DtBuCH18C6 + DBS + Dodec)/SiO 2 -P adsorbent to Sr(II) at a flow rate of 0.3 mL/min was investigated.The results were simulated through the developed theoretical model.
Figure 4 presents the breakthrough curve of (DtBuCH18C6 + DBS + dodec)/SiO 2 -P as the function of the feeding volume.The plots exhibit a similar trend to the typical S-shaped line, and the values of C/C o increase rapidly to 0.8 and then slowly reach 1.This finding demonstrates the excellent dynamic adsorption property and high efficiency of the adsorbent.The breakthrough point of 5% in column experiments is often used as a momentous indicator for evaluating the adsorption performance of adsorbents.The breakthrough capacity (B.Cap.) and total capacity (T.Cap.) were acquired when the Sr(II) concentration in the effluent reached 5% (C/C o = 5%) and 100% (C/C o = 100%) of the initial pumping concentration, respectively.As shown in Figure 5, the 5% breakthrough period of 0.3 mL/min was 3.05 h on the dynamic adsorption results of (DtBuCH18C6 + DBS + dodec)/SiO2-P under t sidered flow rate used in designing the adsorption column theoretical model, the in coefficient Kfa of the theoretical model was calculated using Equation (15), and the tion isotherm parameters were obtained through Freundlich model fitting [23].Th lated parameters of the built final column theoretical model are listed in Table 3.
where Cn (t = 0) (g/cm 3 ) represents the Sr(II) concentration of the effluent at the be phase of the column experiment, and Δz represents the minute height of the height division (m).

Simulation of Breakthrough Curves
As shown in Figure 5, the 5% breakthrough period of 0.3 mL/min was 3.05 h.Based on the dynamic adsorption results of (DtBuCH18C6 + DBS + dodec)/SiO 2 -P under the considered flow rate used in designing the adsorption column theoretical model, the involved coefficient K fa of the theoretical model was calculated using Equation (15), and the adsorption isotherm parameters were obtained through Freundlich model fitting [23].The calculated parameters of the built final column theoretical model are listed in Table 3.
where C n (t = 0) (g/cm 3 ) represents the Sr(II) concentration of the effluent at the beginning phase of the column experiment, and ∆z represents the minute height of the column height division (m).
In accordance with the Runge-Kutta numerical methods, in conjunction with the column adsorption experiment conditions, the equations of the developed theoretical model were resolved.The comparison of the model fitting curves with the experimental data is shown in Figure 5.The theoretical model curves were consistent with the experimental breakthrough curves, suggesting that the numerical model could feasibly predict the dynamic adsorption process.The close values of the 5% and 100% breakthrough periods obtained through the model fitting and the experiment data further proved the appropriate model development and parameter calculations (Table 4).The similar results of prediction and experimental data at 0.5 mL/min in our previous work also proved the feasibility of a reasonable mechanical model and appropriate calculations of the parameters.where Cn (t = 0) (g/cm 3 ) represents the Sr(II) concentration of the effluent at the be phase of the column experiment, and Δz represents the minute height of the height division (m).In the actual separation process, the increasing length of the adsorption column in simulation experiments would lead to problems such as a long reaction cycle time, lower experimental efficiency, and higher costs.This problem is especially serious in radioactive conditions, thereby limiting the development of a large number of experiments.The length of the column also affects the loading amount of the adsorbent and the column pressure.With the use of the constants from the above simulation, the breakthrough curves of different bed heights were calculated at a flow rate of 0.3 m 3 /min (Figure 6).Given the high affinity of crown ether for Ba(II), with the assumption that Ba was totally co-adsorbed as Sr in HLLW, the initial concentration of the feed solution was calculated as 2500 g/m 3 (PWR, 235U:4.5%,burnup: 4.5 × 10 4 MW/t, 5 years cooling down) [23].As seen in Figure 5, the breakthrough curve was a regular shift with a column length change.The 5% breakthrough time increased from 0.43 to 4.6 h in column lengths ranging from 5 to 40 cm.The longer column contributed to a longer breakthrough time and a larger processing volume for one cycle.This result is similar to the results of Fujie et al. [41]-that is, a standard breakthrough curve would appear to move parallel to the time axial under constant flow velocity and sufficiently long column condition.Under specific experimental conditions, the appropriate design of the adsorption column length is of great significance for improving adsorbent utilization and economy.The column was proposed to have a length of around 40 cm, according to our previous study, so that the requirements of efficient processing and a relatively large processing volume could be optimized.The effect of the column's height/diameter ratio and heat transfer should be clarified in detail in our subsequent studies.

Dynamic Desorption Performance
In accordance with the Sr(II) desorption performance of the considered four elutes to the (DtBuCH18C6 + DBS + dodec)/SiO2-P adsorbent in the above section, the Na-DTPA solution was utilized to further investigate the dynamic Sr(II) elution ability of the test column.The dynamic ionic elution properties of varying concentrations of the Na-DTPA solution to the adsorption column are described in Figure 7.As the Na-DTPA concentration increases, the Sr(II) elution plots of the column tend to present with a sharper and stronger elution peak.At concentrations of 0.05, 0.075, and 0.1 mol/L, volumes of 94.7, 70.2, and 53.4 mL of the eluent were required to realize the complete elution of Sr(II).The elution rates achieved by the three Na-DTPA concentrations were 82.2%, 86.4%, and 91.2%, respectively.The experiment results indicated that increasing the concentration of Na-DTPA could be beneficial to enhance the efficiency of the elution process because of the improved complexion of Sr(II) with [DTPA] 5− .Thus, the 0.1 mol/L Na-DTPA solution could be used in the design of the column separation process of HLLW to attain Sr(II) separation and recovery.

Dynamic Desorption Performance
In accordance with the Sr(II) desorption performance of the considered four elutes to the (DtBuCH18C6 + DBS + dodec)/SiO 2 -P adsorbent in the above section, the Na-DTPA solution was utilized to further investigate the dynamic Sr(II) elution ability of the test column.The dynamic ionic elution properties of varying concentrations of the Na-DTPA solution to the adsorption column are described in Figure 7.

Dynamic Desorption Performance
In accordance with the Sr(II) desorption performance of the considered four elutes to the (DtBuCH18C6 + DBS + dodec)/SiO2-P adsorbent in the above section, the Na-DTPA solution was utilized to further investigate the dynamic Sr(II) elution ability of the test column.The dynamic ionic elution properties of varying concentrations of the Na-DTPA solution to the adsorption column are described in Figure 7.As the Na-DTPA concentration increases, the Sr(II) elution plots of the column tend to present with a sharper and stronger elution peak.At concentrations of 0.05, 0.075, and 0.1 mol/L, volumes of 94.7, 70.2, and 53.4 mL of the eluent were required to realize the complete elution of Sr(II).The elution rates achieved by the three Na-DTPA concentrations were 82.2%, 86.4%, and 91.2%, respectively.The experiment results indicated that increasing the concentration of Na-DTPA could be beneficial to enhance the efficiency of the elution process because of the improved complexion of Sr(II) with [DTPA] 5− .Thus, the 0.1 mol/L Na-DTPA solution could be used in the design of the column separation process of HLLW to attain Sr(II) separation and recovery.As the Na-DTPA concentration increases, the Sr(II) elution plots of the column tend to present with a sharper and stronger elution peak.At concentrations of 0.05, 0.075, and 0.1 mol/L, volumes of 94.7, 70.2, and 53.4 mL of the eluent were required to realize the complete elution of Sr(II).The elution rates achieved by the three Na-DTPA concentrations were 82.2%, 86.4%, and 91.2%, respectively.The experiment results indicated that increasing the concentration of Na-DTPA could be beneficial to enhance the efficiency of the elution process because of the improved complexion of Sr(II) with [DTPA] 5− .Thus, the 0.1 mol/L Na-DTPA solution could be used in the design of the column separation process of HLLW to attain Sr(II) separation and recovery.

Column Separation of Simulated HLLW
Given the complex composition inside HLLW, the effect of competitive ions on the dynamic Sr(II) separation property of the material needs to be investigated.Therefore, SHLLW was used for the dynamic column separation study of the synthesized material, and the results are shown in Figure 8.

Column Separation of Simulated HLLW
Given the complex composition inside HLLW, the effect of competitive ions on the dynamic Sr(II) separation property of the material needs to be investigated.Therefore, SHLLW was used for the dynamic column separation study of the synthesized material, and the results are shown in Figure 8. Aside from the elements of Sr(II) and Ba(II), all rested impurity nuclides were eluted from the adsorption column after 3 mol/L nitric acid solution feeding.The similar chemical properties of Ba (II) and Sr (II) led to the simultaneous adsorption of Ba (II) on (DtBuCH18C6 + DBS + dodec)/SiO2-P [36].Flowing of a 20 mL nitric acid solution could successfully achieve the separation of Sr(II) with most coexisting cations in the SHLLW, indicating the excellent Sr(II) dynamic selective adsorption capacity of (DtBuCH18C6 + DBS + dodec)/SiO2-P.With the better elution efficiency of the Na-DTPA solution to Sr(II) under low acidity, the deionized water was retained for use after the feeding of HNO3 to decrease the H + concentration of the filled layers.Sr(II) remained anchored in the column after the offering of water proved that the modification of DBS resulted in a wider acidity range of DtBuCH18C6-loaded SiO2-P material, and the adsorbent exhibited good Sr(II) adsorption performance at low acidity.The recovery rate of each element is shown in Table 5.All nuclides could achieve a recovery of more than 80%, and the rate of Sr(II) exceeded 99%, which revealed that Sr(II) could be successfully separated from HLLW through the (DtBuCH18C6 + DBS + Dodec)/SiO2-P adsorption column.
However, some Ba(II) impurities existed in the eluent during the Sr(II) elution process of the packed adsorption column by using the Na-DTPA solution, revealing that the synthesized adsorbent possessed a certain Ba(II) complexation capacity.This condition could be reasonably ascribed to the similar chemical properties of Ba(II) and Sr(II) and that the DtBuCH18C6 adsorbent has a poor ability to separate Ba(II) and Sr(II) elements.Aside from the elements of Sr(II) and Ba(II), all rested impurity nuclides were eluted from the adsorption column after 3 mol/L nitric acid solution feeding.The similar chemical properties of Ba (II) and Sr (II) led to the simultaneous adsorption of Ba (II) on (Dt-BuCH18C6 + DBS + dodec)/SiO 2 -P [36].Flowing of a 20 mL nitric acid solution could successfully achieve the separation of Sr(II) with most coexisting cations in the SHLLW, indicating the excellent Sr(II) dynamic selective adsorption capacity of (DtBuCH18C6 + DBS + dodec)/SiO 2 -P.With the better elution efficiency of the Na-DTPA solution to Sr(II) under low acidity, the deionized water was retained for use after the feeding of HNO 3 to decrease the H + concentration of the filled layers.Sr(II) remained anchored in the column after the offering of water proved that the modification of DBS resulted in a wider acidity range of DtBuCH18C6-loaded SiO 2 -P material, and the adsorbent exhibited good Sr(II) adsorption performance at low acidity.The recovery rate of each element is shown in Table 5.All nuclides could achieve a recovery of more than 80%, and the rate of Sr(II) exceeded 99%, which revealed that Sr(II) could be successfully separated from HLLW through the (DtBuCH18C6 + DBS + Dodec)/SiO 2 -P adsorption column.However, some Ba(II) impurities existed in the eluent during the Sr(II) elution process of the packed adsorption column by using the Na-DTPA solution, revealing that the synthesized adsorbent possessed a certain Ba(II) complexation capacity.This condition could be reasonably ascribed to the similar chemical properties of Ba(II) and Sr(II) and that the DtBuCH18C6 adsorbent has a poor ability to separate Ba(II) and Sr(II) elements.

Density Functional Theory Computational Results
Two geometric models of 4,4 -di-tert-butylcyclohexane-18-crown-6 and 4,5 -di-tertbutylcyclohexane-18-crown-6 structures were established by using GaussView 6.0 software to derive a further understanding of the complexation properties of DtBuCH18C6 to metal ions [42] (Figure 9).The ESP diagrams of the ground state obtained after the optimization of the built models are presented in Figure 10, and the red and blue regions in the picture indicate the positive and negative ESP areas of DtBuCH18C6, respectively.The centers of the DtBuCH18C6 structures exhibited a relatively low ESP, revealing the distribution of the main adsorption sites.

Density Functional Theory Computational Results
Two geometric models of 4,4′-di-tert-butylcyclohexane-18-crown-6 and 4,5′-di-t butylcyclohexane-18-crown-6 structures were established by using GaussView 6.0 s ware to derive a further understanding of the complexation properties of DtBuCH18C metal ions [42] (Figure 9).The ESP diagrams of the ground state obtained after the o mization of the built models are presented in Figure 10, and the red and blue region the picture indicate the positive and negative ESP areas of DtBuCH18C6, respectively.centers of the DtBuCH18C6 structures exhibited a relatively low ESP, revealing the dis bution of the main adsorption sites.

Density Functional Theory Computational Results
Two geometric models of 4,4′-di-tert-butylcyclohexane-18-crown-6 and 4,5′-di-tertbutylcyclohexane-18-crown-6 structures were established by using GaussView 6.0 software to derive a further understanding of the complexation properties of DtBuCH18C6 to metal ions [42] (Figure 9).The ESP diagrams of the ground state obtained after the optimization of the built models are presented in Figure 10, and the red and blue regions in the picture indicate the positive and negative ESP areas of DtBuCH18C6, respectively.The centers of the DtBuCH18C6 structures exhibited a relatively low ESP, revealing the distribution of the main adsorption sites.On the basis of the ESP results of the optimized DtBuCH18C6 structure, DtBuCH18C6-Sr(II) complexation models were developed, and the Sr(II) atom was set in the center of the crown ether ring.Figure 11 displays the geometric optimized DtBuCH18C6-Sr(II) complexation models.The frequency calculation results showed that no virtual frequency existed, indicating the complexation models after optimization were the most stable structures.Sr(II) remained in the center of the DtBuCH18C6 ring and close to the oxygen atoms inside the crown ether ring, and the center-to-center distances of the Sr-O bonds varied from 2.63 Å to 2.80 Å, indicating that DtBuCH18C6 was mainly dependent on the adsorption sites of the oxygen atoms to realize the complexation ability of Sr(II).In addition, both optimized complexation structures were almost identical, presenting with unidirectional bending of the crown ether structure.This finding reveals the similar Sr(II) complexation capacity of the two DtBuCH18C6 structures.the center of the crown ether ring.Figure 11 displays the geometric optimized DtBuCH18C6-Sr(II) complexation models.The frequency calculation results showed that no virtual frequency existed, indicating the complexation models after optimization were the most stable structures.Sr(II) remained in the center of the DtBuCH18C6 ring and close to the oxygen atoms inside the crown ether ring, and the center-to-center distances of the Sr-O bonds varied from 2.63 Å to 2.80 Å, indicating that DtBuCH18C6 was mainly dependent on the adsorption sites of the oxygen atoms to realize the complexation ability of Sr(II).In addition, both optimized complexation structures were almost identical, presenting with unidirectional bending of the crown ether structure.This finding reveals the similar Sr(II) complexation capacity of the two DtBuCH18C6 structures.

Mulliken Charge Distribution
The overall coordination numbers of Sr(II) inside both complexes were calculated by using the method introduced in Reference [43], and the results were 6-coordinate for the two structures.Localized orbital bonding analysis was utilized to determine the valence state of the atoms, and the Sr valence states obtained inside both complexes were +2 at a judgment threshold of 50%.These two calculation results further suggest interaction between Sr(II) and O atoms.Table 6 shows the atoms' Mulliken charge distribution of DtBuCH18C6-Sr(II) complexes in the water and gas phases to quantify the charge change of both complex structures with atomic displacement.The Sr atom has a positive charge, and all O atoms have negative charges with similar values, thereby revealing the similar transfer charges number from the oxygen ligand of the crown ether ring to the center atom.Moreover, the interaction forces of Sr(II) with O atoms in different positions were close.The net charge of Sr(II) inside the water phase is greater than that in the gas phase, indicating a greater transfer charge number from the oxygen ligand to Sr(II) in the gas phase.

Mulliken Charge Distribution
The overall coordination numbers of Sr(II) inside both complexes were calculated by using the method introduced in Reference [43], and the results were 6-coordinate for the two structures.Localized orbital bonding analysis was utilized to determine the valence state of the atoms, and the Sr valence states obtained inside both complexes were +2 at a judgment threshold of 50%.These two calculation results further suggest interaction between Sr(II) and O atoms.Table 6 shows the atoms' Mulliken charge distribution of DtBuCH18C6-Sr(II) complexes in the water and gas phases to quantify the charge change of both complex structures with atomic displacement.The Sr atom has a positive charge, and all O atoms have negative charges with similar values, thereby revealing the similar transfer charges number from the oxygen ligand of the crown ether ring to the center atom.Moreover, the interaction forces of Sr(II) with O atoms in different positions were close.The net charge of Sr(II) inside the water phase is greater than that in the gas phase, indicating a greater transfer charge number from the oxygen ligand to Sr(II) in the gas phase.

Binding Energy Analysis
The quantified calculated energies of Sr(II) adsorption by DtBuCH18C6 are listed in Table 7.The results showed that the total energy of the system decreased after the formation of the DtBuCH18C6-Sr(II) complexes, indicating that the yield process of the complexes was energetically favorable [35].Moreover, the binding energy is closely related to the stability and binding ability of the complex.The crown ether structure of the 4 ,4 isomer processed a lower adsorption energy compared to the 4 ,5 type, indicating that the stability of the complex that formed after the adsorption of Sr(II) by 4 ,4 -Di-tert-butyldicyclohexyl-18crown-6 was stronger.However, the interaction free energies were similar, further proving that the Sr(II) complexation ability of the two crown ether structures was similar.During the Sr(II) extraction process of DtBuCH18C6, the nitrate ions and DDBS inside the aqueous solution could be the ligand ions that participate in the complexation reaction.Therefore, two complexation models of DtBuCH18C6 to Sr(II) with the coordination of different ions were developed to investigate the effect of ligand ions on the complexation process.The ESP graphics of benzenesulfonic acid (BSA) and DBS were calculated (Figure 12), including the HOMO-LUMO gap and the atomic dipole moment corrected Hirshfeld (ADCH) atomic charges, to simplify the calculation process (Table 8).The results showed that the gap values, ADCH atomic charges, and ESP images of both ligands were relatively close, indicating that using BSA to replace DBS as the ligand for the complexation calculation is reasonable.
The quantified calculated energies of Sr(II) adsorption by DtBuCH18C6 are liste Table 7.The results showed that the total energy of the system decreased after the mation of the DtBuCH18C6-Sr(II) complexes, indicating that the yield process of the c plexes was energetically favorable [35].Moreover, the binding energy is closely relate the stability and binding ability of the complex.The crown ether structure of the 4′,4″ mer processed a lower adsorption energy compared to the 4′,5″ type, indicating that stability of the complex that formed after the adsorption of Sr(II) by 4′,4″-Di-tert-buty cyclohexyl-18-crown-6 was stronger.However, the interaction free energies were sim further proving that the Sr(II) complexation ability of the two crown ether structures similar.During the Sr(II) extraction process of DtBuCH18C6, the nitrate ions and DDBS side the aqueous solution could be the ligand ions that participate in the complexa reaction.Therefore, two complexation models of DtBuCH18C6 to Sr(II) with the coo nation of different ions were developed to investigate the effect of ligand ions on the c plexation process.The ESP graphics of benzenesulfonic acid (BSA) and DBS were ca lated (Figure 12), including the HOMO-LUMO gap and the atomic dipole moment rected Hirshfeld (ADCH) atomic charges, to simplify the calculation process (Table 8).results showed that the gap values, ADCH atomic charges, and ESP images of both ands were relatively close, indicating that using BSA to replace DBS as the ligand for complexation calculation is reasonable.On the basis of the ESP mapping images, the ionic coordination models of DtBuCH18C6-Sr(II) with BSA or the nitrate ligands (Figure 13) were established, and the optimized structures of the two coordination models are shown in Figure 14.The frequency calculation results showed no virtual frequency, indicating that the optimized model possessed the most stable structure.Therefore, both DBS and nitrate ions can be complexed with DtBuCH18C6-Sr(II) through a bidentate coordination mode.With the exception of the misattributed C atoms, the total coordination number of Sr(II) inside each coordination model was calculated to be 14, including 10 O atoms and 2 S or N atoms, which is consistent with the optimized structure.
Toxics 2023, 11, x FOR PEER REVIEW On the basis of the ESP mapping images, the ionic coordination mo DtBuCH18C6-Sr(II) with BSA or the nitrate ligands (Figure 13) were established, optimized structures of the two coordination models are shown in Figure 14.T quency calculation results showed no virtual frequency, indicating that the op model possessed the most stable structure.Therefore, both DBS and nitrate ions complexed with DtBuCH18C6-Sr(II) through a bidentate coordination mode.With ception of the misattributed C atoms, the total coordination number of Sr(II) insi coordination model was calculated to be 14, including 10 O atoms and 2 S or N which is consistent with the optimized structure.

Interaction Force and Binding Energy Analysis
Figure 15 shows the reduced density gradient (RDG) isosurface of the cente DtBuCH18C6-Sr(II) complexes with BSA or the nitrate ligand coordination.The indicated that weak interaction attractions existed between Sr(II) and the O atom ionic ligand, and a weak repulsion existed between Sr(II) with the N or S atoms.Th van der Waals force interactions were suggested to occur between the benzene ring and the crown ether.On the basis of the ESP mapping images, the ionic coordination mo DtBuCH18C6-Sr(II) with BSA or the nitrate ligands (Figure 13) were established, optimized structures of the two coordination models are shown in Figure 14.T quency calculation results showed no virtual frequency, indicating that the op model possessed the most stable structure.Therefore, both DBS and nitrate ions complexed with DtBuCH18C6-Sr(II) through a bidentate coordination mode.With ception of the misattributed C atoms, the total coordination number of Sr(II) insi coordination model was calculated to be 14, including 10 O atoms and 2 S or N which is consistent with the optimized structure.

Interaction Force and Binding Energy Analysis
Figure 15 shows the reduced density gradient (RDG) isosurface of the cente DtBuCH18C6-Sr(II) complexes with BSA or the nitrate ligand coordination.The indicated that weak interaction attractions existed between Sr(II) and the O atom ionic ligand, and a weak repulsion existed between Sr(II) with the N or S atoms.Th van der Waals force interactions were suggested to occur between the benzene ring and the crown ether.

Interaction Force and Binding Energy Analysis
Figure 15 shows the reduced density gradient (RDG) isosurface of the center of the DtBuCH18C6-Sr(II) complexes with BSA or the nitrate ligand coordination.The results indicated that weak interaction attractions existed between Sr(II) and the O atoms of the ionic ligand, and a weak repulsion existed between Sr(II) with the N or S atoms.The direct van der Waals force interactions were suggested to occur between the benzene ring of DBS and the crown ether.Thermodynamic calculations for both complexes were performed at the def2-s sis group level, and the results are listed in Table 9.The negative reaction energie cated that both coordination reactions were energetically favorable for the forma stable complex structures.Therefore, the promotion mechanism of DBS could be a ligand that participated in the complexation process of DtBuCH18C6 on Sr(II), t enhancing the complexation adsorption ability of the crown ether groups on Sr(II aqueous solution.Furthermore, DBS, as an organic BSA, exhibited certain hydroph [44,45], and the complexation behavior of DBS with Sr(II) could explain the weak d tion ability of deionized water and nitric acid to (DtBuCH18C6 + DBS + Dodec) adsorbent after Sr(II) adsorption.The interaction free energy of the crown ether nated with nitrate showed a greater value than that with BSA, revealing that the i tion of the DtBuCH18C6-Sr(II) complex with the nitrate ion was stronger than that o

Conclusions
Static and dynamic selective removal performance studies of a DBS-modified based DtBuCH18C6 adsorbent for Sr(II) in simulated HLLW revealed that the synth material had excellent Sr(II) separation ability and could be used for Sr(II) decont tion in HLLW.(DtBuCH18C6 + DBS + dodec)/SiO2-P achieved 88.2% adsorption c utilization during the dynamic test at a flow rate of 0.3 cm 3 /min.A mathematical a tion dynamics model was built: a standard breakthrough curve would appear to parallel to the time axial under a constant flow velocity and sufficiently long colum dition.The close values of the 5% and 100% breakthrough periods from model fitti the experiment data proved the feasibility of a reasonable mechanical model and priate calculations of the parameters.Chromatographic separation of Sr(II) fro SHLLW was accomplished by using stepwise elution techniques.The Sr(II) ions w Thermodynamic calculations for both complexes were performed at the def2-svp basis group level, and the results are listed in Table 9.The negative reaction energies indicated that both coordination reactions were energetically favorable for the formation of stable complex structures.Therefore, the promotion mechanism of DBS could be DBS as a ligand that participated in the complexation process of DtBuCH18C6 on Sr(II), thereby enhancing the complexation adsorption ability of the crown ether groups on Sr(II) in the aqueous solution.Furthermore, DBS, as an organic BSA, exhibited certain hydrophobicity [44,45], and the complexation behavior of DBS with Sr(II) could explain the weak desorption ability of deionized water and nitric acid to (DtBuCH18C6 + DBS + Dodec)/SiO 2 -P adsorbent after Sr(II) adsorption.The interaction free energy of the crown ether coordinated with nitrate showed a greater value than that with BSA, revealing that the interaction of the DtBuCH18C6-Sr(II) complex with the nitrate ion was stronger than that of DBS.

Conclusions
Static and dynamic selective removal performance studies of a DBS-modified silicabased DtBuCH18C6 adsorbent for Sr(II) in simulated HLLW revealed that the synthesized material had excellent Sr(II) separation ability and could be used for Sr(II) decontamination in HLLW.(DtBuCH18C6 + DBS + dodec)/SiO 2 -P achieved 88.2% adsorption column utilization during the dynamic test at a flow rate of 0.3 cm 3 /min.A mathematical adsorption dynamics model was built: a standard breakthrough curve would appear to move parallel to the time axial under a constant flow velocity and sufficiently long column condition.The close values of the 5% and 100% breakthrough periods from model fitting and the experiment data proved the feasibility of a reasonable mechanical model and appropriate calculations of the parameters.Chromatographic separation of Sr(II) from the SHLLW was accomplished by using stepwise elution techniques.The Sr(II) ions were effectively eluted with 0.1 mol/L Na-DTPA, and the dynamic recovery was up to 99%.The optimized structures and interaction energies of DtBuCH18C6 coordinated with Sr(II) and the DBS modifier were investigated by using DFT calculations.The adsorption sites of oxygen atoms played a major role in the complexation process, and the center distance of the Sr-O bond in the optimized complexation model was 2.63 Å to 2.80 Å.The crown ether could

Figure 3 .
Figure 3.Effect of various solutions on the elution rate of Sr(II).

Figure 3 .
Figure 3.Effect of various solutions on the elution rate of Sr(II).

3. 3 . 1 .
DtBuCH18C6-Sr(II) Complexation Model On the basis of the ESP results of the optimized DtBuCH18C6 structu

3. 3 . 1 .
DtBuCH18C6-Sr(II) Complexation Model On the basis of the ESP results of the optimized DtBuCH18C6 structure, DtBuCH18C6-Sr(II) complexation models were developed, and the Sr(II) atom was set in
Mo 1.2 g/L Al 5.3 g/L Fe 4.8 g/L Ba 1.1 g/L
Table 2 lists the column parameters, including the breakthrough volume, dynamic adsorption amount, and calculated utilization rate of the adsorption column.A relatively high column utilization rate (5% B. Cap./T.Cap.) of 88.17% was obtained, indicating the effective removal of Sr(II) by the (DtBuCH18C6 + DBS + dodec)/SiO 2 -P packed column.Toxics 2023, 11, x FOR PEER REVIEW 88.17% was obtained, indicating the effective removal of Sr(II) by the (DtBuCH18C + dodec)/SiO2-P packed column.

Table 2 .
Dynamic binding capacity and the column utilization.

Table 2 .
Dynamic binding capacity and the column utilization.

Table 3 .
Calculated mass coefficient K fa and parameters of the Freundlich model.

Table 4 .
Comparison of the predicted breakthrough curve and actual experimental result.

Table 5 .
Recovery rate of each ion inside the SHLLW.

Table 5 .
Recovery rate of each ion inside the SHLLW.

Table 8 .
HOMO-LUMO gap and ADCH atomic charges of the DBS and BSA ligands.

Table 8 .
HOMO-LUMO gap and ADCH atomic charges of the DBS and BSA ligands.

Table 9 .
Calculated energies for DtBuCH18C6-Sr(II) coordinated with DBS or the nitrate lig

Table 9 .
Calculated energies for DtBuCH18C6-Sr(II) coordinated with DBS or the nitrate ligands.