Silica gel functionalized with imidazolium group via click chemistry –new stationary phase for ion chromatography Chimica Techno Acta ARTICLE published by Ural Federal University 2021, vol. 8(4), № 20218409 eISSN 2411-1414; chimicatechnoacta.ru DOI: 10.15826/chimtech.2021.8.4.09 1 of 6 Silica gel functionalized with imidazolium group via click chemistry – new stationary phase for ion chromatography D.A. Chuprynina, I.A. Lupanova, V.V. Konshin , Dzh.N. Konshina * Kuban State University, Department of Chemistry and High Technologies, 350040 Stavropolskaya st., 149, Krasnodar, Russia * Corresponding author: jfox@list.ru This article belongs to the regular issue. © 2021, The Authors. This article is published in open access form under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/). Abstract This manuscript describes the preparation of a simple effective ion- exchange material based on silica gel, on the surface of which me- thylimidazolium bromide is fixed using a click reaction. The resulting material was used as a stationary phase for the separation and de- termination of Cl–, NO2–, NO3–, I–, and SO42– using ion exchange chromatography. The separation efficiency and retention factors for the selected anions were studied in the pH range 3.5–6.5. The pro- posed material was used for the determination of Cl–, SO42– in water and can be suggested for successful use in real water samples. Keywords click reaction ion exchange modified silica Received: 06.08.2021 Revised: 16.11.2021 Accepted: 09.12.2021 Available online: 11.12.2021 1. Introduction The development of new anion-exchange phases with im- proved efficiency and selectivity is one of the topical prob- lems of modern ion chromatography [1]. The selectivity of an anion exchanger is largely governed by the nature and structure of functional layer and the method of its attach- ment to the matrix surface. The efficiency of an anion ex- changer depends on the type of material, as well as on the morphology and the packing mode of particles. Organic pol- ymers can be considered as the most convenient and com- mon matrices for the design of anion exchangers. Neverthe- less, such materials as zirconia, alumina, and, especially, silica gel are quite often used for the preparation of or- ganomineral materials, stationary phases in chromatog- raphy and adsorbents for solid-phase extraction. The disad- vantage of silica materials is their limited stability at lower and higher pH values, especially in alkaline solutions. Sili- ca-based anion exchangers are often used in the pH range 2.0–9.5. However, compared to organic polymers, silica- based ion exchangers have the advantages of higher chro- matographic efficiency and greater mechanical stability. In addition, such materials are preferable for operation in the nonsuppressive version of ion chromatography with con- ductometric detection, since in this case it is necessary to use dilute eluents, which is possible with materials of low exchange capacity. Silica gel-based sorbents are synthesized using a con- ventional approach, which consists in the surface modifi- cation with different functional groups. Functionalization through covalent attachment of a modifier to the matrix surface has a number of advantages. First, the required amount of sorption and ion-exchange centers is governed by the structure and amount of a modifier. Second, varia- tion in the structure of a modifier can influence the ca- pacity, efficiency, and separation selectivity of an ion. Third, there are cross linking agents which enable ex- tending the working pH range to 9.2 without affecting the efficiency over the entire life cycle of a column [2]. Such surface-grafted anion exchangers have the ad- vantage of a small thickness of the ion-exchange layer that favors an increase in the rate of mass transfer upon ion exchange and thereby makes it possible to separate ions with high performance and high selectivity [3]. Various classes of organic compounds are used as the surface modifiers of stationary phase matrices [4]. In recent years, there is a growing interest in the use of ionic liquids that enable a wide variation in the nature of a cationic moiety, which influences the properties of ob- tained materials [5–7]. Among ionic liquids, imidazolium salts gained wide- spread acceptance as modifiers. Examples of their use as efficient extractants capable of forming ion-associative complexes with simple and complex anions have been described. Such complexes are readily produced and quite stable and variation of functional groups in the cat- ionic moiety of ionic liquids offer manifold possibilities to apply salts in different version of sample preparation: http://chimicatechnoacta.ru/ https://doi.org/10.15826/chimtech.2021.8.4.09 https://orcid.org/0000-0003-1864-531X https://orcid.org/0000-0002-9239-5470 http://creativecommons.org/licenses/by/4.0/ Chimica Techno Acta 2021, vol. 8(4), № 20218409 ARTICLE 2 of 6 sorption concentration, liquid-liquid microextraction, and modification of stationary phases in gas and liquid chromatography [7, 8]. Materials with attached imidazolium salts have been used successfully as stationary phases in liquid chroma- tography for separation of caffeine, theophylline, theo- bromine [9], xylose, glucose [10], ephedrine [11], and vitamins [12], as well as for separation of organic and inorganic anions [13–15]. Due to the variety of available starting reagents, not only materials for supernatant col- umns, but also solid columns with attached ionic liquids have been synthesized [16, 17]. The aim of the present work was to obtain ion-exchange materials based on silica gel having a particle size of 8–12 m with imidazolium salt covalently immobilized by click reaction and to study whether they can be used as a sta- tionary phase for ion-exchange chromatography. 2. Experimental 2.1. Reagents and instrumentation The sorbent was prepared using silica gel “Sorbfil” with a particle size of 8–12 m. IR spectra were recorded on a Shimadzu IR Prestige spectrometer in a range of 400–4200 cm–1. 13С NMR spec- tra were measured using a 400 MHz Bruker WB Avance III spectrometer operating at 9.39 Tesla equipped with a Bruker H-F/X 4 mm pencil CP/MAS probehead. 13С chemi- cal shifts were referenced to external solid TSP ((trime- thylsilyl)propionic acid sodium salt) standard. Cross- polarization technique from 1H with spinning sideband suppression (CP TOSS), contact pulse durations of 2–4 ms at a MAS rate of 10 kHz was used. The thermal stabilities of modified silica gel samples were studied on an STA 409 PC Luxx synchronous thermal analyzer (Netzsch, Gemrany) in a temperature range from 30 to 1000 °C at a heating rate of 10 оC/min in the air at- mosphere in Al2O3 ceramic crucibles. The pH value of working buffers was verified on an Expert-001 ionomer using a calibrated ESC-10608 com- bined glass electrode. Chromatographic properties of the column packed with modified silica gel were studied using a modular high- performance liquid chromatograph from Shimadzu (Kyoto, Japan) including a CTO-20A column thermostat, an LC-20AD sp mobile phase feed module, and CDD-10A vp. conductometric detector. The volume of injection loop was 20 l. Data were collected and processed using the LCsolu- tion program. 2.2. Column packing 316 Stainless steel HPLC columns (150×2 mm) (Phenom- enex) were used. The chromatographic column was packed by the suspension method under a pressure of 13 MPa. A test portion of the modified silica gel was added to a C2H5OH–CHCl3 solution (1:1, v/v). The column was packed and the sorbent was compacted on exposure to ultrasound. After packing, the column was conditioned by passing iso-propanol in a volume equal to the 20-fold vol- ume of the packed column, next – bidistilled water, and then – a working mobile phase until the background signal has become constant. 2.3. Synthesis of imidazolium-modified silica gel Acetonitrile (70 mL), 3-azidopropyl silica gel (5 g), 1-methyl-3-prop-2-yn-1-yl-1H-imidazolium (1 g), CuI (0.095 g), and N,N,N`,N`-(tetramethylethylenediamine) (750 L) were placed in a pressure flask with fluoroplastic screw cap and magnetic bar in the argon atmosphere. The resulting suspension was kept with vigorous stirring at 70 °C for 4 h. Silica gel was separated on a Schott filter, washed with acetone, water, 2 M hydrochloric acid, and again acetone and dried at 55 °C for 12 h under a residual pressure of 5 mm Hg. A portion of the resulting modified silica gel was fur- ther treated with a solution of hexamethyldisilazane in toluene for 8 h at 80 °C. 2.4. Determination of the total exchange capacity of the modified silica gel The maximum exchange capacity of the material was deter- mined by titrimetry. The modified silica gel (0.5 g) was agi- tated with 0.1 M nitric acid (20 mL) for 1 h [18]. Silica gel was filtered off and the amount of chloride ion released after the ion exchange reaction with 0.1 M HNO3 was determined by titrimetry in an aliquot portion of the filtrate. To the ali- quot portion of filtrate (5 mL), 0.05 M AgNO3 (5 mL) was added and the excess of unreacted AgNO3 was titrated with 0.05 M KSCN using a saturated solution of Fe(NH4)(SO4)2 (0.2 mL) as an indicator. Titration was terminated after a sorrel color of the solution appeared due to the formation of iron rhodanate complex. The calculated total exchange capac- ity was 0.26±0.02 mM/g. The capacity of the silanized mate- rial remained unchanged. 3. Results and discussion The test object was an organomineral material based on silica gel with acovalently attached imidazolium group obtained by the azide-alkyne cycloaddition click reaction (Scheme 1). In 13С NMR spectrum for the modified silica gel with imidazolium salt two groups of spectral signals are seen. One group contains two signals at δ 124.4 and 137.4 corre- sponding to the carbon nuclei in nitrogen-containing ring. Another group consists of spectral signals at δ 52, 44.6, 36.8, 24.1, and 17, which corresponds to the carbon atoms of aliphatic –CH2– fragments. The signal at δ 9.7 corre- sponds to the –CH3 group. Fig. 1 shows the IR spectra for the starting silica gel with covalently attached azide group and the imidazolium- bearing material obtained according to Scheme 1. Chimica Techno Acta 2021, vol. 8(4), № 20218409 ARTICLE 3 of 6 Sil-im Scheme 1 Synthesis of the modified silica gel Both spectra display a broad intense band at about 1300–1290 cm–1 corresponding to stretching vibrations of the siloxane (Si–O–Si) bond in silica. The intense ab- sorption band at 1627 cm–1 is due to bending vibrations of water absorbed on the silica gel surface. The broad intense band at 3200–3500 cm−1 corresponds to the stretching vibrations of О–H adsorbed on the water sur- face and silanol groups. The intense absorption band at 2106 cm–1 corresponds to stretching vibrations of the azido group grafted to the silica gel surface. In the spec- trum of the silica gel sample obtained after the click reaction, the stretching vibration band of the azide group disappears, which suggests a successful click re- action on the silica gel surface. The obtained batch of silica gel with a covalently im- mobilized group was divided into two portions. One por- tion was treated with a silanization reagent, hexamethyl- disilazane (Scheme 2), in order to inactivate residual silanol groups, which are additional sorption centers. The second portion of the material was used without additional treatment. The heat stability is one of the key characteristics of sorption materials, since it governs the temperature range of their possible application (Table 1). Similar temperature regions of weight loss and the presence of exothermic effect at about 350 °C can be dis- tinguished for both samples of the modified silica gel. Table 1 Thermal analysis data of silica gel samples Sil-im (h) Sil-im Temperature range, oС Weight loss,% Temperature range, oС Weight loss,% 30–154 4.7 30–165 2.8 154–85 6.1 165–325 4.4 285–566 10.7 325–460 3.6 566–950 0.2 460–900 4.1 Fig. 1 IR spectra of the modified silica gel samples Chimica Techno Acta 2021, vol. 8(4), № 20218409 ARTICLE 4 of 6 Sil-im (h) Scheme 2 Synthesis of the modified silica gel According to the literature data [19], the first section on the TG curve in a range from 80 to 170 °C is due to the evaporation of water adsorbed on the silica gel surface. Further decrease in the sample weight in a range from 170 to 900 °C corresponds to destruction of the functional or- ganic layer. The total weight loss of modified silica gel samples at 950 °C was about 16.0–21.7%. Chromatographic conditions of the modified silica gels were studied in a single-column version of ion chromatog- raphy with non-suppressed conductivity detection. In such cases, for determination of anions on stationary phases pos- sessing relatively low capacity values, diluted eluents based on aromatic acids, such as benzoic and phthalic acids, are preferred. The choice of the nature and composition of a buffer solution was caused by the fact that solutions based on phthalic acid possess high buffer capacity in the pH range recommended for stationary phases based on silica [20, 21]. A mixture of Cl–, NO2–, NO3–, I–, and SO42– anions was chosen as the model. The working parameters for chroma- tographic separation were chosen as follows: the eluent was HOOCC6H4COOK with C = 2.5 mM and pH = 4 [22]. The comparison of the efficiency of separation of the standard anion mixture by the studied materials under iden- tical conditions demonstrates a considerable decrease in the plate numbers per meter (N/m) for Sil-im(h) (Table 2). Table 2 Efficiency of anion separation (N/m) and retention time (tr, min) on the Sil-im and Sil-im (h) sorbents Anion Sil-im Sil-im (h) N tr N tr Cl– 3126 5.0 2720 11.8 NO2 – 10093 5.9 4045 14.2 NO3 – 8686 6.5 3546 17.8 I– 6273 8.7 2808 29.6 The addition of NaOH to the stationary phase results in an increase in the efficiency and higher rapidity of separation of some anions due to an increase in the eluting power of the mobile phase. Upon pH change from 4 to 6, the concentration of hydrogen phthalate, an average-strength eluting ion, in- creases and, upon pH above 6, divalent phthalate with high eluting ability becomes the main anionic form of the eluent. In addition, at pH close to 7, the residual silanol groups on the surface of the anion exchanger, which can enter into ion exchange interactions, undergo almost complete ionization resulting in an increase in the separation efficiency (Fig. 2). The most efficient separation was observed at рН = 6.5; one can note a multiple decrease in the retention, which leads to a possibility of determining anions in lower amounts. No change in the elution order was observed, which indirectly suggests a predominant ion-exchange mechanism of anion separation on the proposed stationary phase. A B Fig. 2 The effect of NaOH content in the mobile phase on the separation efficiency N (N/m) (A) and the retention K of anions on the Sil-im(h) sorbent (B) Chimica Techno Acta 2021, vol. 8(4), № 20218409 ARTICLE 5 of 6 The linearity of the method was tested using a series of inorganic anions standard solutions. Each point of the calibration plot was the average of three peak height measurements, because the baseline resolution of some anions couldn’t be achieved without further dilution of the mobile phase and significant decrease of efficiency, which could cause a problem in the analysis of the samples with complex matrices. An example of a typical chromatogram is shown in the Fig. 3. The coeffi- cient for calibration curves, the linear range and detec- tion limits (defined as a signal three times the height of the noise level) as well as the quantification limits are presented in Table 3. The comparison of the quantitative characteristics by the example of univalent inorganic anions demonstrates narrowing of the working concentration range on going from Sil-im to Sil-im(h) and, as a consequence, an increase in the limit of detection of an analyte. Fig. 3 Chromatogram obtained with the Sil-im stationary phase using anion-exchange conditions. Test mixture: 1) F–; 2) CH3COO –; 3) IO3 –; 4) Cl–; 5) NO2 –+Br–; 6) NO3 –; 7) I–; 8) SCN–; 9) SO4 2–. Chromatographic conditions: mobile phase: 2.5 mmol/L HOOCC6H4COOK with pH = 4, flow-rate: 0.3 ml/min, injection volume: 20 l and detection: non-suppressed conductivity However, the sensitivity of determination in this case increases, which is evidenced by an increase in the slope of the calibration curve. The error in the determi- nation of standard solutions of the analyzed group of anions as estimated by the added-found method differs in a regular manner: no more than 6.2% for Sil-im and no more than 4% for Sil-im(h) (Table 4). The applicability of the obtained ion-exchange mate- rial as a stationary phase in the ion-exchange chroma- tography was estimated by the example of Sil-im upon determination of inorganic anions in mineral water. Since the single-column non-suppressed ion chromatog- raphy used in this study does not possess a high sensi- tivity compared to two-column with suppressed conduc- tivity detection, although it allows one to determine simultaneously in the isocratic mode weakly and strongly retained anions, it seems interesting to analyze real objects with a high content of anions. A “Lyso- gorskaya” bottled mineral drinking water relating to a group of chloride-sulfate ones was chosen as the test object (Table 5). The content of main macrocomponents was estimated by the calibration curve. To verify the determination accuracy of Cl– and SO42–, their contents in the sample were estimated using a commercial col- umn on a DIONEX ICS-300 chromatographic system with a possibility of eluent generation and suppressed conductivity detection. The confidence intervals of ion determination ob- tained for different chromatographic determination sys- tems overlap, which suggests the accuracy of the ob- tained data. Table 3 Sensitivity factors and linear ranges of calibration curves for the studied anions Anion Sil-im(h) Sil-im Linear range, mg/L a* R2 LOD**, mg/L Linear range, mg/L a* R2 LOD**, mg/L Cl– 40–160 6227 0.9986 5.6 25–400 3494 0.9985 1.9 NO2 – 20–160 2210 0.9983 4.9 25–400 2914 0.9990 1.6 NO3 – 40–160 3859 0.9978 12.1 25–200 484 0.9982 3.8 I– 40–160 6325 0.9984 11.0 25–300 685 0.9925 11.2 *a – coefficient for calibration curves (y = ax + b) **LOD limits of detection = S∙3.3, S – standard deviation of ten independent measurements of a blank sample Table 4 Errors in the determination of the studied anions using the studied anion exchangers (n = 3, P = 0.95) Anion Sil-im(h) Sil-im Added, mg/L Found, mg/L Δ, % Added, mg/L Found, mg/L Δ, % Cl– 80 77.3±9.9 –3.4 100 94.5±14.1 –5.5 NO2 – 78.3±10.2 –2.1 94.0±14.1 –6.0 NO3 – 80.8±10.4 1.0 97.5±14.3 2.5 I– 81.5±10.6 1.9 106.2±16.2 6.2 Table 5 Assessment of the anion content in the real sample using different chromatographic systems (n = 5, P = 0.95) Anion Cl– SO4 2– Stationary phase C, mg/L Cclaimed, mg/L Sr, % C, mg/L Cclaimed, mg/L Sr, % Sil-Im 2227±387 2200–7700 1.93 7321±425 5500–9000 0.17 Seporus A-UNI (HC-1) 2845±345 4.13 7691±880 3.87 Chimica Techno Acta 2021, vol. 8(4), № 20218409 ARTICLE 6 of 6 4. Conclusions In this work, the possibility of modifying the silica gel sur- face with an imidazolium salt using a click reaction was shown. The resulting ion exchange material was used as a stationary phase in the ion exchange chromatography method to separate Cl–, NO2–, NO3–, I–, and SO42–. The study showed that the proposed material allows the de- termination of the selected anions with an error of 3.5–6.0%. Moreover, the developed material showed good stability and repeatability of the results of the determina- tion and separation of anions during the operation. Acknowledgements This publication was financially supported by the Ministry of Science and Higher Education of the Russian Federation (project no. FZEN-2020-0022). Сonflicts of Interests The authors declare that they have no competing interests. References 1. Zatirakha AV, Smolenkov AD, Shpigun OA. Preparation and chromatographic performance of polymer-based anion ex- changers for ion chromatography. Anal Chim Acta. 2016;904:33–50. doi:10.1016/j.aca.2015.11.012 2. Alpert AA, Regnier FE. Preparation of a porous microparticu- latee anion-exchange chromatography support for proteins. J Chromatogr. 1979;185:375–392. doi:10.1016/S0021-9673(00)85615-0 3. Collinson MM. Recent trends in analytical applications of or- ganically modified silicate materials. Trends Anal Chem. 2002;21(1):31–39. doi:10.1016/S0165-9936(01)00125-X 4. Jal PK, Patel S, Mishra BK. Chemical modification of silica surface by immobilization of functional groups for extractive concentration of metal ions. Talanta. 2004;62:1005–1028. doi:10.1016/j.talanta.2003.10.028 5. Wang Y, Tian M, Bi W, Row KH. Application of ionic liquids in high performance reversed-phase chromatography. Int J Mol Sci. 2009;10:2591–2610. doi:10.3390/ijms10062591 6. Pino V, Afonso AM. Surface-bonded ionic liquid stationary phases in high-performance liquid chromatography – a re- view. Anal Chim Acta. 2012;714:20–37. doi:10.1016/j.aca.2011.11.045 7. Vidal L, Riekkola ML, Canals A. Ionic liquid-modified materi- als for solid-phase extraction and separation: a review. Anal Chim Acta. 2012;715:19–41. doi:10.1016/j.aca.2011.11.050 8. Shamsi SA, Danielson ND. Utility of ionic liquids in analytical separations. J Sep Science. 2007;30:1729–1750. doi:10.1002/jssc.200700136 9. Bi W, Row KH. Comparison of different silica-based imidazo- lium stationary phases for LC in separation of alkaloids. Chromatographia. 2010;71:25–30. doi:10.1365/s10337-009-1380-9 10. Bi W, Zhou J, Row KH. Separation of xylose and glucose on different silica-confined ionic liquid stationary phases. Anal Chim Acta. 2010;677(2):162–168. doi:10.1016/j.aca.2010.08.004 11. Liu SJ, Zhou F, Xiao XH, Zhao L, Liu X, Jiang SX. Surface con- fined ionic liquid – A new stationary phase for the separation of ephedrines in high-performance liquid chromatography. Chin Chem Lett. 2004;15(9):1060–1062. 12. Qiu H, Jiang Q, Wei Z, Wang X, Liu X, Jiang S. Preparation and evaluation of a silica-based 1-alkyl-3-(propyl-3-sulfonate) im- idazolium zwitterionic stationary phase for high- performance liquid chromatography. J Chromatogr A. 2007;1163(1–2):63–69. doi:10.1016/j.chroma.2007.06.001 13. Qiu H, Jiang S, Liu X. N-Methylimidazolium anion-exchange stationary phase for high-performance liquid chromatog- raphy. J Chromatogr A. 2006;1103(2):265–270. doi:10.1016/j.chroma.2005.11.035 14. Qiu H, Jiang S, Liu X, Zhao L. Novel imidazolium stationary phase for high-performance liquid chromatography. J Chro- matogr A. 2006;1116(1–2):46–50. doi:10.1016/j.chroma.2006.03.016 15. Qiu H, Wang L, Liu X, Jiang S. Preparation and characteriza- tion of silica confined ionic liquids as chromatographic sta- tionary phases through surface radical chain-transfer reac- tion. Analyst. 2009;134:460–465. doi:10.1039/B809125K 16. Liu K, Tolley HD, Lawson JS, Lee ML. Highly crosslinked pol- ymeric monoliths with various C6 functional groups for re- versed-phase capillary liquid chromatography of small mole- cules. J Chromatogr A. 2013;1321:80–87. doi:10.1016/j.chroma.2013.10.071 17. Lin H, Ou J, Zhang Z, Dong J, Zou H. Ring-opening polymeriza- tion reaction of polyhedral oligomeric silsesquioxanes (POSSs) for preparation of well-controlled 3D skeletal hybrid monoliths. Chem Commun. 2013;49(3):231–233. doi:10.1039/c2cc36881a 18. Iamamoto MS, Gushikem Y. Adsorption of metal ions from aqueous and ethanol solutions by silica gel functionalized with pyridinium ions. J Colloid Interface Sci. 1989;129(1):162–165. doi:10.1016/0021-9797(89)90426-8 19. Protsak IS, Kuzema PO, Tertykh VA, Bolbukh YM, Kozakevich RB. Thermogravimetric analysis of silicas chemically modi- fied with products of deoligomerization of polydimethylsilox- ane. J Therm Anal Calorim. 2015;121:547–557. doi:10.1007/s10973-015-4623-z 20. Sun M, Feng J, Wang X, Duan H, Li L, Luo C. Dicationic imid- azolium ionic liquid modified silica as a novel re- versed-phase/anion-exchange mixed-mode stationary phase for high-performance liquid chromatography. J Sep Sci. 2014;37(16):2153–2159. doi:10.1002/jssc.201400176 21. Haddad PR, Jackson PE. Ion Chromatography: Principles and Applications, Elsevier: Amsterdam; 1990. 775 p. 22. Konshin VV, Chuprynina DA, Lupanova IA, Konshina DzhN. Silica gel with covalently immobilized imidazolium salt as a new stationary phase for ion chromatography. AIP Conf Proc. 2020;2280(1):050028. doi:10.1063/5.0018424 https://doi.org/10.1016/j.aca.2015.11.012 https://doi.org/10.1016/S0021-9673(00)85615-0 https://doi.org/10.1016/S0165-9936(01)00125-X https://doi.org/10.1016/j.talanta.2003.10.028 https://doi.org/10.3390/ijms10062591 https://doi.org/10.1016/j.aca.2011.11.045 https://doi.org/10.1016/j.aca.2011.11.050 https://doi.org/10.1002/jssc.200700136 https://doi.org/10.1365/s10337-009-1380-9 https://doi.org/10.1016/j.aca.2010.08.004 https://doi.org/10.1016/j.chroma.2007.06.001 https://doi.org/10.1016/j.chroma.2005.11.035 https://doi.org/10.1016/j.chroma.2006.03.016 https://doi.org/10.1039/B809125K https://doi.org/10.1016/j.chroma.2013.10.071 https://doi.org/10.1039/c2cc36881a https://doi.org/10.1016/0021-9797(89)90426-8 https://doi.org/10.1007/s10973-015-4623-z https://doi.org/10.1002/jssc.201400176 https://doi.org/10.1063/5.0018424