Journal of Environmental and Public Health

Journal of Environmental and Public Health / 2021 / Article

Research Article | Open Access

Volume 2021 |Article ID 6612500 |

R. Devanathan, G. L. Balaji, R. Lakshmipathy, "Adsorption of Rare Earth Ce3+ and Pr3+ Ions by Hydrophobic Ionic Liquid", Journal of Environmental and Public Health, vol. 2021, Article ID 6612500, 8 pages, 2021.

Adsorption of Rare Earth Ce3+ and Pr3+ Ions by Hydrophobic Ionic Liquid

Academic Editor: Ping Xiang
Received15 Dec 2020
Revised31 Mar 2021
Accepted24 May 2021
Published30 May 2021


This study reports the use of hydrophobic ionic liquid (IL) based on D-galactose for the recovery of Ce (III) and Pr (III) ions from solutions. The equilibrium data were obtained by optimization of batch parameters, and various isotherms and kinetic models were utilised to predict the mechanistic process of sequestration of ions. The Arrhenius activation energies are found to be between 5–40 kJ, suggesting the physisorption process of ions onto IL. The present process is understood to be rapid and exothermic in nature according to thermodynamic experiments. The loading capacity was found to be 179.3 g L−1 and 141.5 g L−1, respectively, for Ce (III) and Pr (III) ions at pH 5 with a contact time of 30 min and dose being 0.1 g L−1. The higher uptake capacity is attributed to the presence of a highly electronegative fluorine atom in the IL. These results highlight the potential application of IL in the sequestration of Ce (III) and Pr (III) ions from any water sources.

1. Introduction

Rare earth elements (REEs) are found to have imperative applications in high-tech consumer products such as television, cell phones, and many other electronic appliances. REEs find wider applications in mundane and futuristic domains because of their chemical, electrical, optical and magnetic values [1]. A burgeoning population, modernization, and rapid industrialization have created a huge demand for the elements that are less commonly exploitable due to geochemical properties and dispersion [2]. With the ever-growing demand, the recovery of RREs from wastewater streams will be highly desirable and economical rather than depending on the mineral ore deposits.

Several techniques for the recovery and preconcentration of RREs are in practice such as coprecipitation [3], ion exchange [4], solvent extraction [5], and adsorption [6]. All the abovementioned techniques except adsorption find limitations due to the high cost of the operation, production of pollutants, and low recovery rates. Adsorption is a premier technique found to be a powerful and cost-effective method for the elimination and recovery of heavy metal ions and precious RREs. Many adsorbents such as activated carbon [6], agricultural waste and by-products [7], and organic frameworks [8] are found to be prolific in terms of cost and being ecofriendly; however, the efficiency and regeneration is always a concerning factor. Recently, the development of nanomaterials and ionic liquids has led to a new era of technological advancements in the field of wastewater treatments and recovery of REEs. Gul et al. published an interesting review on carbon-based nanoadsorbents for advanced water treatment [9]. Dhiman and Kondal reported the review of popular ZnO nanoadsorbents for the effective treatment of wastewater which highlights the role of ZnO nanoadsorbents [10]. In spite of the high efficiency exhibited by nanomaterials, the cost of production and unknown toxicity remain less prolific.

Ionic liquids are a class of salts that commonly contains short-lived species and ions and has exhibited enhanced distribution coefficients for heavy metal ions abolition [11, 12]. However, the liquid-phase ionic liquids require chelating agents to extract the cations from liquid phases [13]. Ionic liquids anchored onto nanomaterials were also explored for the extraction of cations from liquid phases [14, 15]. Hydrophobic ionic liquids have exhibited high efficiency towards cation removal and have put an end to choosing a chelating agent or support material [16]. Hydrophobic ionic liquids have exhibited greater efficiency towards the removal of cations from liquid phases [17]. Ionic liquids have been used in the removal and recovery of rare earth elements. Sol-gel materials doped with ionic liquids have been explored for the elimination of Ce (III) ions in liquid phases [18]. The application of hydrophobic ionic liquids in the removal and storing of rare earth elements are yet to be explored fully. In continuation to the earlier reports, the present investigation aims at the removal of rare earth elements such as trivalent Ce and Pr ions from liquid phases by hydrophobic PF6-anion-based ionic liquid.

2. Materials and Methods

The synthesis of galactose-based IL (V) was carried out as per the earlier reports [19]. The scheme of the synthesis of IL5 is represented in Scheme 1. The white solid obtained at the final step was washed with water several times under vacuum.

2.1. Batch Adsorption Studies

In order to optimize the parameters for the removal of Ce (III) and Pr (III) ions, batch mode adsorption was adopted. Variables such as H+ concentration, time of contact, IL5 dose, and preliminary rare earth metal ion concentration were optimized by varying each parameter and keeping the remaining variables unchanged. The hydronium ion concentration of the solution was honed by changing the hydronium ion concentration from pH 2–9. The contact time between the IL5 and Ce (III) and Pr (III) ions was varied from 5 to 30 min, and the IL5 dosage was considered from 10 to 50 mg. The preliminary REE’s ion amount was considered between 10 and 60 mg L−1. In addition to that mentioned above, the influence of temperature was also investigated by changing the variable from 303 to 323 K. The temperature-controlled orbital shaker was equipped to carry out all the experiments, and the test solutions were agitated at 150 rpm in a 50 ml conical flask. After each parameter optimization, the solid was recovered by centrifugation and the residual remaining was analysed with flame-AAS to quantify the residual concentration of the trivalent Ce and Pr ions. The percentage removal and maximum sequestration capacity of the IL5 were found by the following expressions:

3. Results and Discussion

The synthesized IL5 was evaluated as a potential sorbent for the sequestration of Ce (III) and Pr (III) ions from liquid phases by the batch process.

3.1. Effect of H+ Ion Concentration

The H+ ion concentration is one of the imperative parameters that must be optimized prior to any adsorption studies. The H+ ions’ concentration in the solution plays a vital role in the binding of cations onto the surface of the adsorbent. In order to fix the H+ ion concentration of the solution, investigations were performed on changing the hydronium ion concentration from 2–9, and the results are shown in Figure 1. It is observed that, with increasing pH, the sequestration capacity of the IL5 increased for both the cations and was found to be maximum at pH 5. The surge in elimination ability is because of reducing amounts of H+ ion that competes with the trivalent Ce and Pr ions for the active sites on IL5. Beyond pH 5 for Ce3+ and pH 6 for Pr3+, the removal efficiency from the aqueous solution was found to decrease, and this may be due to the surface of IL5 being positive and the formation of several hydroxides. Moving forward, pH 5 was set as the optimal H+ concentration for later investigations.

3.2. Influence of IL5 Dosage

The quantity of the sorbent required is of significance in selecting an adsorbent for any adsorption system, and any classic adsorbent displays high output at low dosage levels. In order to optimize and evaluate the efficacy of the IL5 dose, the dose was varied from 0.1 to 0.5 g L−1, and the outputs are represented in Figure 2. The results clearly depict that, at a very low dosage, the IL5 exhibited over 75% of removal efficiency. Furthermore, with the surge in dosage, the ability of removal increases, but the efficacy is substantially small and negligible. The high removal efficiency of IL5 at low dose levels indicates that IL5 is a prolific sorbent for the elimination of Ce (III) and Pr (III) ions from liquid phases.

3.3. Influence of Contact Time

Investigations were performed on varying the contact time (5 to 30 min) to fix the optimum time of contact for the sorption of trivalent Ce and Pr ions from liquid phases, and the outputs are represented in Figure 3. The adsorption of Ce (III) and Pr (III) ions was rapid within 5 min of time of contact and gradually increased with the extension of time. The swift removal efficiency noticed within 5 min of time must be due to spontaneous occupation of active sites, and slower adsorption with burgeoning time is due to exhaustion of sites on the plane of IL5. Moving forward, 30 min of time was fixed as the optimum time for other investigations.

3.4. Desorption of Ce (III) and Pr (III) Ions from IL5

Desorption and regeneration of the adsorbent are of prime importance since the desorbed cations can be recovered and used for further application. The regenerated adsorbent can be reused again in the sorption process which makes the process economical. In view of that mentioned above, the adsorbed Ce (III) and Pr (III) ions onto IL5 were desorbed by contacting with various desorbing agents such as 0.01 M HCl, 0.01 M CH3COOH, and 0.01 M NaOH, and Figure 4 represents the outcomes. It is evident that the 0.01 M HCl was superior in desorbing the Ce (III) and Pr (III) ions from IL5 followed by 0.01 M CH3COOH and NaOH. The regenerated IL5 was reused for two successive cycles.

3.5. Adsorption Isotherms

At equilibrium, the initial concentration of the trivalent Ce and Pr ions was varied to estimate the uptake capacity of the present ionic liquid. The initial quantities of selected ions were changed from 10 to 60 mg L−1, and the uptake ability was calculated to be 179.3 and 141.5 g L−1 for Ce (III) and Pr (III) ions, respectively. The superior electronegative atom such as fluorine present in the anion part of IL is primarily responsible for supreme loading capacity exhibited by IL5. The electron cloud around the 6 F atoms attracts the positively charged Ce (III) and Pr (III) ions. In order to predict the mechanistic process of adsorption, two well-known mathematical models such as Langmuir and Freundlich adsorption isotherms were employed.

The plots of selected isotherms models as mentioned above are shown in Figures 5(a) and 5(b), and the respective constants and coefficient of correlation are tabulated in Table 1. The coefficients of correlation for Langmuir and Freundlich isotherms were observed to be close to one recommending the applicability of the models for the recovery of Ce (III) and Pr (III) ions from liquid phases by IL5. Furthermore, the relevance of the Langmuir model is backed up by the theoretical loading capacity that is the same as the experimental values. The binding of Ce (III) and Pr (III) ions onto IL5 is a multilayer adsorption process, and each layer obeys Langmuir isotherm.

IsothermParametersCe (III)Pr (III)

FreundlichKf (mg g−1)6.238.91
1/n (L mg−1)0.0210.017

LangmuirVm (mg g−1)169.2144.4
b (L mg−1)0.0230.031

3.6. Adsorption Kinetics

Kinetic investigations provide valuable information on the mechanism of removal of the analytes, and experiments were performed in line with the equilibrium studies for Ce (III) and Pr (III) ions at various preliminary REE ion concentration and temperatures. The residual concentrations of test samples obtained at preset time intervals were analysed for understanding the mechanistic aspects of adsorption. Various kinetic models helpful in predicting the mechanism of binding were adopted, and the obtained kinetic data were analysed.

The coefficients of correlation and respective constants of the models used for understanding the mechanism of binding of Ce (III) and Pr (III) ions onto IL5 are summarized in Tables 2 and 3. The coefficients of correlation at different temperatures and concentrations obtained for pseudo-second-order were found to be high compared to those of the other two models used in this study. The correlation coefficients of the pseudo-second-order kinetic model are very near to one, and the loading capacity theoretical values are very similar to experimental advising the better fit of the model. These observations suggest that the binding of Ce (III) and Pr (III) ions onto IL5 is governed by chemical reactions. However, for an adsorption process to be controlled by chemical reactions, many setups must be satisfied. These conditions are as follows [20]:(i)The constants of the rate of reaction should remain unchanged irrespective of the initial concentration of target ions(ii)The size of the adsorbent particle should not alter the rate constant(iii)The degree of shaking should not influence the rate constant

Kinetic modelParameters303 K313 K323 K
10 (mg L−1)30 (mg L−1)10 (mg L−1)30 (mg L−1)10 (mg L−1)30 (mg L−1)

Experimentalqe (mg g−1)53.5126.352.9125.851.7123.9

Pseudo-first-orderk1 (min−1)0.0230.0960.0330.1090.0370.122
qe (mg g−1)12.421.111.720.710.519.7

Pseudo-second-orderk2 (g mg−1min−1)0.0090.0070.0100.0080.0120.010
qe (mg g−1)54.1126.853.8126.253.2125.7

Elovichα (mg g−1 min−1)2.1 × 10−252.0 × 10−262.3 × 10−262.1 × 10−272.4 × 10−272.3 × 10−27
β (mg g−1)1.161.871.241.851.331.88

Kinetic modelParameters303 K313 K323 K
10 (mg L−1)30 (mg L−1)10 (mg L−1)30 (mg L−1)10 (mg L−1)30 (mg L−1)

Exploratoryqe (mg g−1)47.1107.746.7106.646.1106.1

Pseudo-first-orderk1 (min−1)0.0270.1170.0250.1130.0220.109
qe (mg g−1)11.621.311.220.610.719.8

Pseudo-second-orderk2 (g mg−1 min−1)0.0160.0110.0170.0130.0190.015
qe (mg g−1)46.3106.945.7106.544.8105.8

Elovichα (mg g−1 min−1)1.9 × 10−262.1 × 10−262.3 × 10−262.2 × 10−262.5 × 10−262.7 × 10−26
β (mg g−1)

If any of the abovementioned conditions are not obeyed, then rate controlling is not chemical reaction kinetics even though the kinetic values follow the pseudo-second-order model [13]. Experiments were carried out at two different initial metal ion concentrations of Ce (III) and Pr (III) ions onto IL5 in order to check whether the rate of reaction is controlled by chemical reactions or not. In Tables 2 and 3, the rate constant k appears to be not inconsistent for varying initial metal ion concentration of Ce (III) and Pr (III) ions onto IL5. The inconsistency in the value of k for varying initial concentration of counter ions suggests that chemical reactions are not the rate-limiting step, even though kinetic data tend to fit well with the pseudo-second-order model. Communications observed were similar to this study for the elimination of Cd2+ ions by ionic liquid [17] and the removal of Pb2+ and Cu2+ ions by Watermelon rind [7].

The Arrhenius equation expresses temperature as a function for the second-order rate constant.

The Arrhenius activation energy calculated for the varying preliminary metal ion concentrations of Ce (III) and Pr (III) is tabulated in Table 4. Interestingly, the activation energy for both metal ions at two different initial concentrations is less than 40 kJ, suggesting that the process is physisorption. In general, chemisorption has higher activation energy between 40 and 800 kJ/mol. These observations conclude that the current sorption system is physisorption.

IonConcentration (mg L−1)Ea (kJ mol−1)

Ce (III)1015.7

Pr (III)1016.4

Furthermore, to understand the diffusion mechanism, the rate data were analysed with the intraparticle diffusion model.

If the resulting plot is linear, it suggests that the process is completely controlled by intraparticle diffusion, and the multilinear plot indicates two or more steps of a diffusion process. In the present process, the recovery of Ce (III) and Pr (III) ions by IL5 is intraparticle diffusion only and it is concluded by the liner plots obtained from the experimental results.

3.7. Thermodynamics of Adsorption

The removal of Ce (III) and Pr (III) ions by IL5 was investigated at varying temperatures in order to understand the nature of the adsorption process. It was found that the efficiency of removal and loading capacity decreased with the surge in temperature and this must be because of the weakening of the surface binding spots of IL5 that binds the Ce (III) and Pr (III) ions. Also, with the surge in temperature, the kinetic energy of the analytes increases, and as a result, the probability of counterions meeting the active sites of IL5 decreases. Furthermore, the thermodynamic parameters which reveal the nature of the adsorption process were derived and summarized in Table 5.

Metal ionTemperature (K)Sorption capacity (qe mg g−1)G° (kJ mol−1)H° (kJ mol−1)S° (J mol−1 k−1)

Ce (III)30353.5−4.56−5.22692.3

Pr (III)30347.1−4.12−4.76789.1

The variables of thermodynamics such as Gibbs free-energy change (∆G°), enthalpy (∆H°), and entropy (∆S°) derived for the present process are tabulated in Table 5. The free-energy change (∆G°) at all temperatures is found to be negative for the recovery of target analyte ions from an aqueous solution by IL5, suggesting that the process is spontaneous. However, with the rise in temperature, the negative values of ∆G° were found to decline, suggesting that the spontaneity decreases with the rise in temperature. The lowering of loading capacity also supports this observation. The exothermic nature of the present process was understood from the negative values of enthalpy (∆H°), and the positive values of entropy (∆S°) highlight that the entropy increases at the interface of Ce (III) and Pr (III) ions and IL5.

3.8. Influence of Salt Ionic Strength

It is salient to consider the influence of ionic strength of various salts on the removal and recovery of Ce (III) and Pr (III) ions from aqueous solution. In general, effluents tend to contain various ions usually due to salts and hardness of water, and in order to understand the influence of those, efforts were made to understand their role. Sorption experiments were conducted in the presence of varying salt concentrations for the removal and recovery of Ce (III) and Pr (III) ions from the aqueous solution, and the results are tabulated in Table 6. The results were found to be interesting as seen in Table 6. The presence of Na+ ions has a partial impact on the loading capacity of Ce (III) and Pr (III) ions, and this is due to Na+ ions being more electropositive than the target ions. However, the presence of Ca2+ and Mg2+ ions has not influenced the loading capacity, and this is due to the electronegativity of the Ca2+ and Mg2+ ions being the same and greater than the Ce (III) and Pr (III) ions. These results suggest that the hardness of water does not have any influence on loading capacity, but the presence of Na+ ions must be considered prior to the sorption of Ce (III) and Pr (III) ions from the aqueous solution.

Ionic strength (mol L−1)SaltSorption capacity (mg g−1)
Ce (III)Pr (III)




3.9. Adsorption Mechanism

The mechanism of removal of Ce (III) and Pr (III) ions by IL is attributed to the electrostatic attraction of the cations on the highly electronegative fluorine atoms present in the form of PF6 ions on IL. The fluorine atoms on the IL form an electron cloud over them due to their high electronegativity, and this could attract the positively charged cations electrostatically. The electrostatic attraction mechanism is supported by variation of hydronium ion concentration where the high concentration of H+ ions resulted in the low attraction of Ce (III) and Pr (III) ions, while at lower H+ ion concentration, the sorption was found to be superior. The salt ionic strength also corroborates the observations suggesting that the electrostatic attraction is the mechanism and physiosorption is taking place. Desorption and regeneration studies further support the physiosorption process and electrostatic attraction mechanism.

4. Conclusions

This study reported the potential application of carbohydrate-based hydrophobic ionic liquid with PF6 ion as a prolific sequester for the recovery of Ce (III) and Pr (III) ions from aqueous solution. Adsorption batch investigations were performed for optimizing the imperative variables such as pH, time of contact, the dosage of IL, and preliminary Ce (III) and Pr (III) ion concentrations. Adsorption isotherm investigations suggested that the present process is multilayer adsorption and each layer obeys Langmuir isotherm. Profound kinetic studies suggested that the process is physiosorption following a pseudo-second-order kinetic model. The physiosorption is also assisted by the Arrhenius activation energy calculated from various temperatures. The thermodynamic studies highlight the spontaneous and exothermic process of the present system. The results conclude that the present ionic liquid is a potential precursor for the recovery of Ce (III) and Pr (III) ions from water sources.

Data Availability

The data of this study are available in the corresponding author’s Google Sites page.

Conflicts of Interest

The authors declare that they have no conflicts of interest.


  1. G. B. Haxel, J. B. Hedrick, and G. J. Orris, Rare Earth Elements: Critical Resources For High Technology, U.S. Geological Survey, Reston VA, USA, 2002.
  2. K. Veronese, Rare: The High-Stakes Race to Satisfy Our Need for the Scarcest Metals on Earth, Prometheus Books, Buffalo, NY, USA, 2015.
  3. N. Freslon, G. Bayon, D. Birot, C. Bollinger, and J. A. Barrat, “Determination of rare earth elements and other trace elements (Y, Mn, Co, Cr) in seawater using Tm addition and Mg(OH)2 co-precipitation,” Talanta, vol. 85, no. 1, pp. 582–587, 2011. View at: Publisher Site | Google Scholar
  4. V. Hatje, K. W. Bruland, and A. R. Flegal, “Determination of rare earth elements after pre-concentration using NOBIAS-chelate PA-1resin: method development and application in the San Francisco Bay plume,” Marine Chemistry, vol. 160, pp. 34–41, 2014. View at: Publisher Site | Google Scholar
  5. F. Xie, T. A. Zhang, D. Dreisinger, and F. Doyle, “A critical review on solvent extraction of rare earths from aqueous solutions,” Minerals Engineering, vol. 56, pp. 10–28, 2014. View at: Publisher Site | Google Scholar
  6. D. S. R. Murty and G. Chakrapani, “Preconcentration of rare earth elements on activated carbon and its application to groundwater and sea-water analysis,” Journal of Analytical Atomic Spectrometry, vol. 11, no. 9, pp. 815–820, 1996. View at: Publisher Site | Google Scholar
  7. R. Lakshmipathy and N. C. Sarada, “Metal ion free watermelon (Citrullus lanatus) rind as adsorbent for the removal of lead and copper ions from aqueous solution,” Desalination and Water Treatment, vol. 57, no. 33, pp. 15362–15372, 2016. View at: Publisher Site | Google Scholar
  8. X. Liu, H. Pang, X. Liu et al., “Orderly porous covalent organic frameworks-based materials: superior adsorbents for pollutants removal from aqueous solutions,” The Innovation, vol. 2, no. 1, Article ID 100076, 2021. View at: Publisher Site | Google Scholar
  9. A. Gul, N. G. Khaligh, and N. M. Julkapli, “Surface modification of carbon-based nanoadsorbents for the advanced wastewater treatment,” Journal of Molecular Structure, vol. 1235, Article ID 130148, 2021. View at: Publisher Site | Google Scholar
  10. V. Dhiman and N. Kondal, “ZnO nanoadsorbents: a potent material for removal of heavy metal ions from wastewater,” Colloid and Interface Science Communications, vol. 41, Article ID 100380, 2021. View at: Publisher Site | Google Scholar
  11. U. Domanska and A. Rekawek, “Extraction of metal ions from aqueous solutions using imidazolium basedionic liquids,” Journal of Solution Chemistry, vol. 38, pp. 739–751, 2009. View at: Publisher Site | Google Scholar
  12. G.-T. Wei, Z. Yang, and C.-J. Chen, “Room temperature ionic liquid as a novel medium for liquid/liquid extraction of metal ions,” Analytica Chimica Acta, vol. 488, no. 2, pp. 183–192, 2003. View at: Publisher Site | Google Scholar
  13. H. C. J. Camiel, N. A. Macías-Ruvalcaba, M. Aguilar-Martínez, and M. N. Kobrak, “Metal extraction to ionic liquids: the relationship between structure, mechanism and application,” International Reviews in Physical Chemistry, vol. 34, pp. 591–622, 2015. View at: Publisher Site | Google Scholar
  14. M. E. Mahmoud and H. M. Al-Bishri, “Supported hydrophobic ionic liquid on nano-silica for adsorption of lead,” Chemical Engineering Journal, vol. 166, no. 1, pp. 157–167, 2011. View at: Publisher Site | Google Scholar
  15. M. E. Mahmoud, “Surface loaded 1-methyl-3-ethylimidazolium bis (trifluoromethylsulfonyl) imide [EMIM+Tf2N−] hydrophobic ionic liquid on nano-silica sorbents for removal of lead from water samples,” Desalination, vol. 266, no. 1-3, pp. 119–127, 2011. View at: Publisher Site | Google Scholar
  16. N. Kozonoi and Y. Ikeda, “Extraction mechanism of metal ion from aqueous solution to the hydrophobic ionic liquid, 1-Butyl-3-methylimidazolium nonafluorobutanesulfonate,” Monatshefte für Chemie-Chemical Monthly, vol. 138, no. 11, pp. 1145–1151, 2007. View at: Publisher Site | Google Scholar
  17. R. Jayachandra, S. R. Reddy, and R. Lakshmipathy, “D-galactose based hydrophobic ionic liquid: a new adsorbent for the removal of Cd2+ ions from aqueous solution,” Environmental Progress & Sustainable Energy, vol. 38, no. S1, pp. S139–S145, 2019. View at: Publisher Site | Google Scholar
  18. L. Zhang, E. Chen, and X. Yu, “Preparation of sol-gel materials doped with ionic liquid and extractant for cerium (III) extraction,” Journal of Radioanalytical and Nuclear Chemistry, vol. 298, no. 2, pp. 1055–1063, 2013. View at: Publisher Site | Google Scholar
  19. R. Jayachandra, R. Lakshmipathy, and S. R. Reddy, “Hydrophobic d-galactose based ionic liquid for the sequestration of Pb2+ ions from aqueous solution,” Journal of Molecular Liquids, vol. 219, pp. 1172–1178, 2016. View at: Publisher Site | Google Scholar
  20. C. E. Harland, Ion Exchange: Theory and Practice, The Royal Society of Chemistry, Londom, UK, 1994.

Copyright © 2021 R. Devanathan et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Related articles

No related content is available yet for this article.
 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

No related content is available yet for this article.

Article of the Year Award: Outstanding research contributions of 2021, as selected by our Chief Editors. Read the winning articles.