Plasma-assisted preparation of nano-(ZrC, ZrO2)@carbon composites from Zr-loaded sulfonated styrene–divinylbenzene copolymers

We have developed a simple method to prepare nano-(ZrC0.93, ZrO2-polymorphs)@carbon composites with graphite/amorphous carbon content and adjustable Zr/C ratio based on using a multistep tube furnace and plasma-assisted heat treatment of zirconium-loaded sulfonated styrene–divinylbenzene (STY-DVB) copolymers. Pre-pyrolysis of zirconium-loaded sulfonated STY-DVB ion exchangers with 2 and 8 mass % DVB at temperatures between 1000 and 1400 °C for 2 h produced nano-ZrO2@C intermediates with particle sizes of ~ 30–60 nm with no ZrC formation. Plasma processing of nano-ZrO2@C resulted in nano-(ZrC0.93, ZrO2)@C composites with 11% (under a He atmosphere) (C/Zr = 73) or 13% (under a H2 atmosphere) (C/Zr = 58) ZrC0.93 content. Three polymorphs of the zirconium dioxide (tetragonal, monoclinic and cubic, between 18 and 27 nm) were found in the products. The amounts of tetragonal and monoclinic ones are comparable to that of ZrC0.93. The average particle size of ZrC0.93 prepared in this way was found to be 21–23 nm. The BET surface area of the nano-(ZrC0.93, ZrO2)@C(graphite) composites prepared in He and H2 was over 250 and 300 m2/g, respectively. We developed a reproducible and easy method to prepare nano-(ZrC, ZrO2)@C products by setting the DVB content, sulfonation degree, Zr loading and the thermal treatment conditions, which have an influence on the ZrC and graphite/amorphous carbon content of nano-ZrO2@C intermediates. The zirconium-loaded sulfonated styrene–divinylbenzene (STY-DVB) copolymers (2 and 8 mass% DVB) or their thermal decomposition was characterized with vibrational spectroscopy, thermal analysis and DSC or powder XRD, BET, XPS and HRTEM methods, respectively.


Introduction
Zirconium carbide-carbon (ZrC@C) composites have been considered strategic materials in emitters and coatings for nuclear reactor fuels due to their high melting point, good mechanical properties, small neutron absorption cross section, high sorption capacity and radiolytic stability toward radioactive isotopes [1][2][3][4]. Zirconium dioxide and its composites with carbon are widely used in various branches of science and industry, e.g., in pharmaceutical industry [5][6][7], catalysis and nanotechnology [8,9] and chemical and environmental engineering [10][11][12][13]. Synthesis of various polymorphs of ZrO 2 for various purposes is an intensively studied area of zirconium chemistry [14].
There are two main routes of methods to synthesize ZrC@C composites: (a) the solid-state reactions between elementary carbon and Zr, ZrH 2 or ZrO 2 , [15,16] and (b) the reactions of zirconium or its compounds with organic carbon sources. The zirconium source and carbon compounds may be Zr and methyl halides [17,18]. Zirconium acetate or zirconium oxychloride can also be reacted with sucrose [19,20]. Sacks et al. used zirconium n-propoxide/n-propanol (as a zirconia source) and phenolic resin or glycerol as a carbon source, [21], and Martiz et al. used a plasma-assisted synthesis route from analogous starting materials [22]. Scales et al. synthesized porous, mechanically stable bead-like ZrC@C composites via carbothermal reduction of polyacrylonitrile-ZrCl 4 composites [3] and quaternary ammonium salt group containing styrene-divinylbenzene copolymers loaded with anionic zirconium complexes [4].
Zirconium carbide ceramic materials formed in direct reactions of Zr or Zr compounds and carbon or carbon precursors, even in the case of the most reactive ZrH 2 , always contains some carbon vacancy [16]. The reaction of carbondeficient ZrC x structures with the free carbon results in the formation of non-stoichiometric Zr 1-x C (carbon-rich) materials and increases their sinterability [23].
The mixed composites have great new utilization possibilities-they combine the advantageous effects of nano-ZrC and nano-ZrO 2 on a carbon matrix, e.g., the conductivity of ZrO x C y composites in the molten state leads to emerging new electrolytic technologies to prepare zirconium.
Earlier we developed a simple way to prepare beadshaped stable carbon composites containing nanosized metal oxides from metal-ion-loaded sulfonated or iminodiacetatefunctionalized styrene-divinylbenzene (STY-DVB) copolymers [24][25][26]. Scales et al. prepared ZrC from quaternary ammonium salt-functionalized anion exchangers loaded with anionic Zr-complexes [3]. These results prompted us to develop the preparation of nano-ZrO 2 @C intermediates from Zr-loaded sulfonated STY-DVB copolymers so that we could study their direct transformation in RF thermal plasma into nano-ZrC@C materials.
We tested zirconium loading onto the sulfonated STY-DVB resins from aqueous zirconyl salt (chloride, sulfate and nitrate) solutions, to decide which was optimal. The degree of sulfonation and Zr loading can control the Zr/carbon ratio, whereas the DVB content has an influence on the mechanical properties of the resulting composites. The sulfonated polymer serves as a carbon source, solidification agent and ensures the homogeneous distribution of zirconium in the Zr-compound(s)@C composites produced in a simple hot N 2 tube-furnace pyrolysis process.
The labels of the prepared materials are generated as X-DVB-Zr-compound-T-t, where X stands for the DVB content in the initial sample, T corresponds to the temperature of pyrolysis (1000, 1200, 1400 °C), whereas t represents reaction time (2 or 8 h). The labels of the prepared samples are summarized in Table 1.
In this paper, we present the preparation and properties of nano-ZrO 2 @C and nano-(ZrC, ZrO 2 )@C composites prepared with a tube furnace and by subsequent RF thermal plasma treatment of Zr-loaded sulfonated STY-DVB resins. Also discussed in detail are the characterization of initial X-DVB and X-DVB-Zr polymers, including determining the sulfonation degree and the distribution of sulfonic acid groups on STY and DVB rings, and the influence on zirconium loading of anions in the zirconyl salts used.

Preparation of samples
The sulfonated 2-DVB and 8-DVB resins (Varion KS-2 and KS-8) were soaked in water for 24 h, and the swelled beads were filled into a 20-cm-long column. Zr loading was performed with the use of a zirconyl chloride solution containing 0.5 mass % Zr, in 3 M HCl until reaching the maximum possible loading with zirconium was reached. The 2-DVB-Zr and 8-DVB-Zr samples swelled in water formed during Zr loading were conditioned for 3 h at 90 °C and 2 h at 110 °C in air at atmospheric pressure in a programmable furnace, with a heating rate of 10 °C min −1 . The beads were placed into a porcelain crucible and left to cool in a desiccator filled with CaO.

RF plasma processing
Due to the high amount of water contained in the resins and the short residence time of the particles in RF plasma, the 2-DVB-Zr and 8-DVB-Zr samples could not be sufficiently carbonized directly. Therefore, a ZrO 2 @C sample prepared by pyrolysis at 1000 °C for 2 hours 1 was exposed to in-flight thermal plasma processing in an inert (He) and a reducing (H 2 ) atmosphere. The experimental setup consisted of a reactor chamber with a radio-frequency inductively coupled plasma torch (PL-35, TEKNA Ltd, Canada) mounted on the top, an RF generator (LEPEL, USA) with a maximum power of 30 kW at 4-5 MHz, a cyclone, a filter unit and a vacuum pump.
The experimental system was presented in detail in [27,28]. The solid samples were delivered into the plasma via a powder feeder (Praxair) with helium as carrier gas, with a flow rate of 5 L⋅min −1 . The feeding probe was placed coaxially into the middle of the induction coil. The feeding rate of the powder was 3 g⋅min −1 . The plasma and sheath gas were mixtures of argon and helium with a flow rate of 11 L⋅min −1 Ar and 6 L⋅min −1 He for the plasma gas as well as 35 L⋅min −1 Ar and 5 L⋅min −1 He for the sheath gas. The same test was repeated with hydrogen added into the sheath gas with a flow rate of 5 L⋅min −1 .
The samples prepared from resins with 2 mass % DVB content were not as stable mechanically as samples prepared from resins containing 8 mass % DVB. The size of nano-ZrO 2 (~ 60 nm) particles prepared at 1000, 1200 or 1400 °C for 2 h did not change considerably with increasing temperature. The sample prepared at 1400 °C in 8 h contained smaller (~ 20 nm) and probably more reactive ZrO 2 particles than the samples prepared at 1000 or 1200 °C, but this sample contained an unknown cubic phase (a = 0.76 nm); therefore, we did not consider it further for the preparation of ZrC@C composites. Table E1 shows the most important operating conditions for the particular tests. The tests were performed under a pressure of 70 kPa, and the power was set to 25 kW.

Elemental analysis
CHS analysis: Samples were dried in a Heraeus vacuum oven (Thermo Electron, Germany) in glass vials at 140 °C for 19 h, filled with nitrogen, tightly capped and sealed with plastic film to maintain dryness until further analysis with a Carlo Erba 1106 instrument.

Specific surface area measurements
We calculated the specific surface area (SSA) using the Brunauer-Emmett-Teller (BET) equation based on nitrogen adsorption measurements at 77 K (Autosorb 1C, Quantachrome, seven-point isotherm). Desorption of the samples was performed under vacuum for 24 h at 100 °C. The average size of the particles was calculated with the equation: SSA [m 2 g −1 ] = 6/ρd, where ρ is absolute particle density [g cm −3 ] and d is particle diameter [µm], assuming spherical particles.

Transmission Electron Microscopy
Transmission electron microscopy (TEM) was carried out with a 200 keV Talos Thermo Scientific transmission electron microscope. Grains were crushed in ethanol and deposited onto copper grids covered by Lacey carbon. We obtained bright-field TEM (BFTEM), high-resolution TEM (HRTEM) and high-angle annular dark-field (HAADF) scanning TEM images as well as selected-area electron diffraction (SAED) patterns. In TEM mode, the current was ~ 600 pA and the SAED patterns were obtained with a camera length of 520 mm. We performed elemental mapping in STEM mode in order to study the chemical composition and heterogeneity of the grains. In STEM mode, the current was ~ 200 pA. Beam convergence was 10.5 mrad, camera length was 125 mm, and dwell time was 10 µs. We used a "Super-X" detector system with four silicon drift detectors built into the microscope column. Fast Fourier transforms (FFTs) were calculated from the HRTEM images with Gatan Digital Micrograph 3.6.1 software. The details of TEM spectroscopy are given in [29].

Thermal studies
We performed simultaneous thermogravimetric analysis (TGA), differential scanning calorimetry (DSC) and mass spectrometric evolved gas analysis (MS-EGA) using a SET-ARAM Labsys Evo and a Pfeiffer Vacuum OmniStar instrument under high-purity helium (6.0) with a heating rate of 10 °C min −1 from room temperature to 1600 °C. The details of thermal measurements are given in [30,31].

Scanning Electron Microscopy
We used a scanning electron microscope (SEM, Zeiss type DSM 982 GEMINI with a field emission tungsten cathode) in establishing the morphology of the resins. Prior to the analysis, all samples were coated with an Au layer of 2-3 nm thickness. Backscattered electron images of ground and polished samples were acquired; the accelerating voltage was 10 kV.

Powder X-ray diffraction
Powder XRD patterns were acquired with a Philips Bragg-Brentano parafocusing goniometer with the use of Cu weighted Kα radiation 1.5406/1.5444 Å in the range of 4-70º 2θ, a step size of 0.02º and an interval time of 1 s. Due to problems of determining the background, no goodnessof-fit parameters are given. The unit cell dimensions of the phases in all the patterns are well determined. The mass% of each phase should be considered an estimate rather than an accurate value. The details of the measurements are given in [31].

Vibrational spectroscopy
A Bruker Alpha FT-IR (Bruker, Ettingen, Germany) instrument and a Biorad-Digilab FTS-30-FIR far-IR instrument with a Smart iTR (attenuated total reflectance) accessory were used to collect the analytical and far-IR spectra, respectively. Spectra were made from 16 scans and acquired with a resolution of 4 cm −1 [32][33][34].
For Raman spectroscopy, a Horiba Jobin-Yvon LabRAMtype micro-spectrometer was used with an external laser source of 532 nm Nd-YAG. The laser beam was focused on an objective of 20X (numerical aperture = 0.4). A confocal hole of 1000 µm was used in the confocal system, and a grating monochromator of 1800 mm −1 was used for light dispersion. The detected wavenumber was scanned with a resolution of 3 cm −1 with an accumulation time of 10 s er point for a spectral range of 100-4000 cm −1 .

X-ray photoelectron spectroscopy
X-ray photoelectron spectra of the RF-plasma-treated samples were recorded with a Kratos XSAM 800 instrument (fixed analyzer transmission mode, 40 eV pass energy, Mg K α1,2 (1253.6 eV) excitation, < 1·10 −7 Pa analysis chamber pressure). Survey spectra were recorded for both samples in the 100-1300 eV kinetic energy range. High-resolution spectra of the characteristic photoelectron lines of the Zr, O, S and C were recorded with 0.1 eV steps and a dwell time of 1 s. Quantitative analysis (after removal of the Shirleytype background) was performed with the XPS MultiQuant program with experimentally determined photoionization cross-section data [35] and asymmetry parameters [36] Results and discussion

Preparation and properties of Zr-loaded ion-exchangers (2-DVB-Zr and 8-DVB-Zr)
Two styrene-divinylbenzene copolymer-based sulfonated ion-exchangers with 2 and 8 mass % divinylbenzene content (2-DVB and 8-DVB) were loaded with zirconyl chloride (ZrOCl 2 .8H 2 O) in a 3 M HCl solution. With the elemental analysis data of 2-DVB and 8-DVB resins (CHS analysis, ESI Table 2) and the zirconium content of the eluate formed  (21) during Zr loading, and the amounts of resin and zirconyl salt used, the Zr/sulfonic acid ratio (loading rate), water content and the Zr/C ratio in the starting X-DVB-Zr samples can be calculated in the following way: The DVB content in the samples determines the distribution of carbon content between the aromatic rings and aliphatic parts (-CH-CH 2 -groups). With the analyzed carbon content, the number of styrene and DVB units in a weighted amount of sample can be determined.
The carbon content and DVB/STY ratio determine the number of hydrogen atoms that belong to the polymer chain (5 hydrogen atoms on every styrene ring, 4 hydrogen atoms on DVB rings and 3 or 6 hydrogen atoms belong to the aliphatic parts of the unsulfonated styrene and DVB rings, respectively).
The number of these hydrogen atoms does not change on sulfonation because formally only an SO 3 group is inserted into a C-H aromatic bond. All the other hydrogen content determined by the analysis must belong to the swelling water in the resin samples.
The sulfur is exclusively present in the form of -SO 3 H groups. Thus, the sulfur content allows the calculation of the amount of sulfonic acid amount in the sample directly. In other words, the average sulfonation degree of the rings located in the cross-linked resin network can be calculated from the ratio of S/C and the ratio of the STY and DVB rings. These values were found to be 1.69 and 1.63 for the sulfonated styrene/DVB samples containing 2 and 8 mass % DVB, respectively.
With the amount of sulfonic acid groups and the bound zirconium content (it was determined indirectly from the zirconium content of the eluate [37]), the loading rate (Zr/ SO 3 H group), the Zr/C ratio of X-DVB-Zr samples can be easily calculated. The latter value is of paramount importance in predicting the forming of ZrC.
The Zr/C atomic ratios calculated from the measurements were found to be 0.452 and 0.329 for the 2-DVB-Zr and 8-DVB-Zr samples, respectively. In principle, if the carbon loss during carbonization (degradation into volatile organic materials or CO/CO 2 ) is less than 58 mass % for 2-DVB-Zr samples or 66% for 8-DVB-Zr samples, then the residual carbon content in the solid phase is larger than is required to form stoichiometric ZrC. Depending on the amount of volatile carbon compounds (it depends on pressure, temperature, time and heating method) and the carbon loss during pyrolysis, the residual Zr/C ratio in nano-ZrO 2 @C and nano-ZrC@C composites can be increased and only ZrC products with a carbon vacancy (ZrC 0.63−0.98 [23]) may be created.
The anion of the zirconium salt used during zirconium loading onto the sulfonated polymers has special importance. As is well-known, the zirconium (IV) ions hydrolyze and polymerize in their aqueous solutions, even under highly acidic conditions. [38] The exclusive form of zirconium(IV) in solutions above pH > 0 is generally a zirconyl ion (ZrO 2+ ), but the counter ion has remarkable effect on the nature of species containing Zr. Zirconyl nitrate (ZrO(NO 3 [39] and turns into a non-ionic species. According to this, only 2 mass % of the initial zirconyl nitrate content was adsorbed with sulfonated cationic (Amberlite IR 100) ion exchangers [39].
Zirconium sulfate also strongly hydrolyzes in aqueous solutions, with the formation of basic zirconium sulfates, including zirconyl sulfate. The bond between sulfate ion and zirconium is strong. Since the affinity of zirconium to sulfate is higher than to water, this interaction remains stable in aqueous solutions as well. During the electrolysis of aqueous Zr(SO 4 ) 2 solutions, hydrogen evolution was observed on the cathode, and the species in the solution containing Zr was assigned as H 2 [ZrO(SO 4 ) 2 .3H 2 O] [39,40]. This complex contains the zirconium in an anionic species, and it is not supposed to bind onto cation exchangers.
Zirconyl  [38] in the spectra of the Zr-loaded resins. The non-compensated charges of ZrO 2+ ions are neutralized with neighboring "free" sulfonate and OH groups. We could not reach a saturation degree higher than ~ 80% even with an elongated sorption time on a dynamic column. It may be attributed to the differences between the steric hindrances of p-and o-positions of sulfonic acid groups on the STY and DVB rings, and probably only the easily available p-position sulfonic acid groups can bound zirconyl ions, whereas the sulfonic acid groups in sterically unfavorable o-positions cannot do that.

IR spectroscopic characterization of X-DVB samples
The IR spectra of 2-DVB and 8-DVB resins showed the typical C-H and C-C/aromatic ring vibrational modes of the polystyrene-divinylbenzene copolymers [41,42] and the typical normal modes of aromatic sulfonic acids [42] as well. The evaluation of the spectra is a serious challenge due to strong overlaps of polystyrene network and sulfonic acid modes. The most important peaks around 667, 1007, 1036, 1126 and 1175 cm −1 agree well with the reference values [41][42][43][44][45][46][47][48][49] given for sulfonated polystyrene and polystyrene-DVB samples. The peaks at 1007 cm −1 are assigned to the in-plane bending vibration of the benzene ring but may contain the Zr = O vibrational modes as well. [45] The band that appears at 1126 cm −1 is a combined band, it contains the in-plane skeleton vibration of the benzene ring in the mono-sulfonated polystyrene rings (as a sign of di-substitution), but the mode of ν s (S-O) also appears in this band system confirmed by its shift to 1122 cm −1 when it forms a salt with zirconium(IV) ( Table E3, ESI Figure E1).
The peaks at 669 and 673 cm −1 represent aromatic CH out-of-plane bending (deformation) vibrations, whereas the absence of peaks at around 699 and 759 cm −1 confirms complete sulfonation of all the styrene rings at least once. There are no peaks at 1680 cm −1 associated with the non-polymerized vinyl groups of cross-linkers. Thus, there are no unreacted pendant vinyl groups. A strong peak at 833 cm −1 together with no peak at 775 cm −1 indicates that the substitution on the styrene ring is mainly in the para position [46].
Although the reactivity of the DVB ring toward sulfonation is higher than that of the styrene ring due to the inductive effect of the two alkyl groups, no tri-sulfonated derivatives can be formed because of the low reactivity of the di-sulfonated species toward further sulfonation under the conditions used in ion-exchanger preparation (sulfonation with 96% sulfuric acid).
Thus, the 2-DVB and 8-DVB resins contain only monoand di-sulfonated aromatic rings (STY and DVB). The first sulfonic acids are built into para-and ortho-positions of styrene and DVB rings, respectively. The favored ring positions for incorporation of the second sulfonic acid group are expected to be the ortho-positions of both rings due to ortho-induction of alkyl and m-induction of sulfonic acid groups of the mono-sulfonated derivatives. The sulfonation pathway is illustrated in Fig. 1. The average sulfonation degree of rings in 2-DVB and 8-DVB was found to be 1.63 and 1.69, respectively.

Vibrational spectroscopic characterization of Zr-loaded samples
A free sulfonate ion has a tetrahedral structure with pyramidal distortion (C3V local symmetry) with doubly degenerate antisymmetric stretching vibrations of S-O bonds in the sulfonate groups. The coordination of zirconyl ion to sulfonate groups may induce a change in local symmetry, depending on the coordination mode. A possible symmetric tridentate coordination mode (κ 1 O, κ 1 O j , κ 1 O jj -SO 3 -Zr(= O) does not change the original C 3 v symmetry, whereas upon monoor chelate-forming bidentate coordination of a zirconium ion, the C 3 v symmetry of the -SO 3 − anions is distorted and becomes Cs symmetry. It results in the appearance of two vibrational bands of ν as (SO 3 ) (Table E3).
In the experimental IR spectra of the X-DVB-Zr sample, we assigned the doubly degenerated ν as mode as a band and as a shoulder at room temperature. It excludes the presence of the symmetric (C 3v ) tridentate coordination mode of zirconium to the sulfonate groups and confirms the presence of a strong cation-anion interaction with a mono-or a bidentate coordination mode (Fig. 2) [42,46].
The IR spectra of zirconyl salts containing mono-or bidentate coordinated sulfonate groups should be considerably different, but other factors, e.g., the contribution of ν (C-S) vibrational modes to ν (S-O) modes can make their assignation unambiguous. Nonetheless, DFT calculations on model compounds of sulfonated STY polymers with monoand bidentate sulfonate groups showed that the bidentate sulfonate coordination mode always gave larger splitting of ν as (S-O) bands than monodentate coordination modes. The experimentally found ~ 65 cm −1 splittings of ν as (S-O) bands (decomposition of E degeneration) in the IR spectra of 2-DVB-Zr and 8-DVB-Zr suggest bidentate coordination modes of sulfonate groups to the zirconium in both samples. The coordination sphere of zirconium may contain hydroxyl ions (κ 1 O, κ 1 O j -O-SO 2 Zr(OH) 2 + or water (κ 1 O, κ 1 O j -OSO 2 Zr(O)(OH 2 ) + ). The surplus + 1 charge can be compensated by the neighboring uncoordinated sulfonate, hydroxyl or hydrogen carbonate (due to the carbon dioxide from the air) ions. The evolution of CO 2 was detected by MS (Fig. 3, m/z = 44) on heating our wet samples to around 100 °C; therefore, we tried to detect the presence of the HCO 3 − ion precursor. Unfortunately, the position of the only intense stretching mode (ν as (C-O)) and all the deformation modes in the IR spectra strongly coincide with bands of polymer network or water (ESI Table E1, Figures E2 and E3). The most intensive expected IR band of the hydrogen carbonate ion (ν as (C-O)) is covered by other bands of the X-DVB samples. The intensity of ν s (C-O) mode is very low in the IR spectra, but the highest in the Raman spectra of hydrogen carbonates. Thus, we tried to compare the Raman spectra of wet (we assumed the samples contained HCO 3 − ions) and dry (we assumed they contained no HCO 3 − due to thermal decomposition during drying at 105 °C for 2 h) X-DVB-Zr samples (ESI Figures E2 and E3) to detect differences between the spectra. Unfortunately, the excitation laser of Raman spectroscopy causes local heating and decomposition of the carbonate species so the Raman spectra of the wet and dry samples were identical (ESI Figures E3). Thus, we cannot assign hydrogen carbonate as a CO 2 source. However, in the IR spectra, there were shoulders in the wet samples around 1200 cm −1 , which disappeared in the dry samplethis may be assigned to chemisorbed CO 2 bound as bidentate carbonate ions [50,51]. We have no more information about the thermal stability of carbonates bidentately bound to zirconium, so the origin of CO 2 signals in the MS (Fig. 3) is left unexplained.

Thermal decomposition of X-DVB-Zr compounds
The analogs of styrene-DVB copolymers decompose on heating with the formation of solid carbonaceous materials and liquid/gaseous hydrocarbons [51][52][53]. The hardness and ratio of the solid carbon/liquid/gaseous decomposition materials depend on the DVB content. The increasing DVB content will increase the hardness and carbon content of the composite [24,27,[52][53][54].
As metals are loaded into sulfonic acid groups, decomposition becomes more controlled, more solid product is formed, and the sulfonate groups are converted into gaseous SO 2 or, in the case of sulfidofilic metals, incorporate into solid sulfides [24,25]. Zirconium does not belong to sulfidofilic elements. Therefore, the expected (and found) decomposition product (intermediate) is the oxide (ZrO 2 ), while carbide formation can be expected only at high temperatures [26] with reduction of Zr-oxide phases in the carbon-rich environment.
Thermal decomposition of 2-DVB-Zr and 8-DVB-Zr samples started analogously, namely with losing the swelling water from the pores of the hydrophilic polymer network phase around ~ 100 °C (ESI Figure E4a, g). No degradation of the organic part was observed in this temperature range. However, carbon dioxide formation was detected by TG-MS (m/z = 44, Fig. 3), which suggests that not only hydroxide ions but also hydrogen carbonate or carbonate ions may neutralize the charges of the zirconium cation.
The TG-MS curves for the main decomposition fragments (e.g., m/z = 64 and 48, SO 2 + and SO + , respectively) belonging to each decomposition step of 2-DVB-Zr and 8-DVB-Zr samples can be seen in Table E4 and ESI Figure E4 (c, i). The evolution of sulfur as SO 2 below the temperatures of ZrO 2 (or ZrC) formation show that the sulfur (sulfonic acid) content of the starting polymer does not directly affect ZrO 2 (ZrC) formation, and the sulfonate content plays a role only in adjusting the Zr/C ratio. Due to the appearance of organic fragments up to 800 °C (ESI Figure E5), the ZrO 2 @C samples were prepared above this temperature (1000, 1200, and 1400 °C) with isothermal heating under N 2 in a tube furnace.
The shape of m/z = 64 and 48 (SO 2 + and SO + , respectively) curves (ESI Figure E4 (c, i)) is different for 2-DVB-Zr and 8-DVB-Zr. Sulfur dioxide evolution appears at a much lower temperature (289 and 410 °C) in the decomposition of 8-DVB-Zr than in that of 2-DVB-Zr (~ 400 and 490 °C). The coinciding decomposition peak temperatures of m/z = 64 (SO 2 ) and m/z = 76, 77 (C 6 H 4 + and C 6 H 5 + ) fragments show that these are formed in the same decomposition process, namely in the cleavage of aromatic rings with the formation of SO 2 , water and organic fragments. No styrene but phenyl fragment formation was observed at 300 °C (ESI Fig. 5f, l) (there is no depolymerization at this temperature) in the case of 8-DVB-Zr (Fig. 5d,j). Styrene appears as a decomposition product, however, at 440 °C for both X-DVB-Zr samples (Fig. 5f,l). These results suggest that cross-linking has a great influence on the decomposition mechanism of the polymer network (Fig. 4). SO 3 was lost in trace amounts (ESI Figure E5), and although the main desulfurization reaction is the degradation of aromatic rings with a simple loss of SO 2 and the formation of aromatic hydrocarbons, the sulfone compounds also form [47] as a minor by-reaction and the presence of residual sulfur can be observed by XPS (Figs. 5, 6), and some Zr-oxysulfide may also be detected by XRD (ESI Figures E6, E7, E8).

Preparation of nano-ZrO 2 @C composites by high-temperature isothermal heat treatment of X-DVB-Zr samples
The 2-DVB-Zr and 8-DVB Zr samples were heated isothermally in a tube furnace between 1000 and 1400 °C. First, the swelling water content of X-DVB-8 samples released and surface melting forms ball-shaped drops, whose surface is carbonized fast due to the heat shock, even before complete liquefaction of the resin particles. The solid carbonaceous layer on the surface of the beads keeps their ball-like shape during carbonization/carbide formation [24][25][26].
Carbonization resulted in ZrO 2 formation due to the oxygen-rich coordination environment of zirconium. The formation of ZrO 2 (monoclinic or tetragonal), its relative amounts and the size of each ZrO 2 polymorph, and that of carbon (graphite or amorphous) prepared by carbonization of 2-and 8-DVB-Zr samples depend on temperature, DVB content and reaction time ( Table 2). Prolonged heating (8 h) at 1400 °C resulted in an unknown cubic phase (a = 0.76 nm) in a minor amount. The particle size of nano-ZrO 2 @C composites prepared from 2-DVB resin between 1000 and 1400 °C was smaller (~ 30-40 nm) than the size of particles prepared from 8-DVB samples (~ 60-70 nm) ( Table 2) under the same conditions.
In all experiments, two types of carbons-graphite and amorphous carbon-were produced. Their relative amounts strongly depend on the above-mentioned conditions. In general, the amount of graphite decreases, whereas the amorphous C content increases with increasing carbonization temperature and reaction time ( Table 2). Due to the low reactivity of ZrO 2 polymorphs toward carbon, no zirconium carbide formed up to 1400 °C.  The direct RF plasma treatment of 2-and 8-DVB-Zr samples was not successful due to their high water content and the short residence time of the injected grains in the plasma, which ultimately led to incomplete carbonization. However, RF thermal plasma treatment offers suitable conditions for the dehydrated intermediates for an in-flight reaction between ZrO 2 and carbon [22,55]. For those tests, the 8-DVB-ZrO 2 -1000-2 sample was selected and plasma processing was performed under an inert (He) and reducing atmosphere (H 2 ), respectively.

Surface characterization of the nano-ZrO 2 @C and nano-(ZrC,ZrO 2 )@C composites
The surface of the nano-ZrC-and nano-ZrO 2 -decorated carbon (graphite) balls was characterized by physisorption method, XPS, and TEM. The BET surface area of the tube furnace-treated samples was generally below 50 m 2 g −1 (ESI Table E5), except for the sample 8-DVB-1400-8, whose surface area was 136 m 2 g −1 . This increased value could probably be attributed to the presence of a cubic phase, which probably considerably contributed to the large surface area of the 8-DVB-Zr-1400-8 composite.
The surface characteristics of the nano-ZrO 2 @C (8-DVB-Zr-1400-2) and nano-(ZrC, ZrO 2 )@)C (H 2 ) composites were determined. Zirconium and carbon are present in various chemical states, as illustrated by the Zr 3d and C 1s spectra in Figs. 5 and 6. Zirconium is found in its carbide form (Zr3d5/2 179.4 eV) only in the 8-DVB-ZrC(H 2 ) sample (the appropriate carbon signal of ZrC is strongly covered by the signal of graphite). Zirconium is also present in oxide (Zr3d5/2 182.8 eV) form. The monoclinic, tetragonal and cubic ZrO 2 phases can be distinguished in the Zr3d spectrum. The line positions of the corresponding carbon, C1s 182.0 eV, and oxygen, O1s 531.2 eV, agree with the literature data [56]. The appearance of the ZrO 2 minor phase can be attributed to the high affinity between Zr and O [57,58]. Part of the oxygen is chemically bound to carbon (O1s 533.0 ± 0.2 eV, C1s 286.7 ± 0.2 eV).
As can be seen from Table 3, the zirconium carbide did not form in a tube furnace, even in 2 h at 1400 °C. Only the two (monoclinic and tetragonal) polymorphs of ZrO 2 and graphite were detected on the surface of the sample. The plasma-treated 8-DVB-ZrC(H 2 ) sample had nearly twice  (Table 3).

TEM analysis of plasma-treated samples
We investigated the microstructure of 8-DVB-Zr-1400-8 and 8-DVB-ZrC (H 2 ) composites by TEM. In BFTEM images, the carbon matrix can be recognized by its bright contrast, whereas the particles containing Zr are dark ( Fig. 7). In contrast, in HAADF images, the contrast is reversed (Fig. 8). As a result of heat and plasma treatments, the original amorphous carbon matrix became poorly crystalline graphite and both samples contained a high amount of unreacted ZrO 2 (blue arrows).
Only the plasma-treated sample contains ZrC (Fig. 7c,  d), which we associate with the effective decomposition of zirconia during the plasma treatment. Although, in theory, the precursor ZrO 2 particles passing through the plasma flame should evaporate, only a small portion dissociated due to the insufficient dwell time in the hottest zone of the plasma. The appropriate residence time varied depending on both the size and the actual thermal trajectory of the individual particles. The EDX composition maps (Fig. 8) suggest that the grains containing Zr are embedded in the carbon matrix. Surprisingly, the maps indicate sulfur content in all the samples. Sulfur dominantly occurs in environments containing oxygen. Thus, it is probably a TiS x O y -type compound (oxysulfide), whereas the sulfur content detected in the 8-DVB-ZrO 2 -1400-8 sample can be associated with Zr x S y .

Conclusions
We have developed a simple method to prepare nano(ZrC 0.93 , ZrO 2 polymorph)@carbon composites with adjustable Zr/C ratio in two-step (tube furnace and plasma-assisted heat treatment) pyrolysis of zirconium-loaded sulfonated styrene(STY)-divinylbenzene(DVB) copolymers. Prepyrolysis of zirconium-loaded sulfonated STY-DVB ion exchangers with 2 and 8 mass % DVB at temperatures between 1000 and 1400 °C (in a tube furnace) in 2 h led to nano-ZrO 2 @C with particle sizes of ~ 30-60 nm without ZrC formation. Plasma processing produced nano-(ZrC 0.93 , ZrO 2 )@C composites with 11% (under He) (C/Zr = 73) or 13% (under H 2 ) (C/Zr = 58) ZrC 0.93 content. Three polymorphs of the zirconium dioxide (tetragonal, monoclinic and cubic with 18-27 nm sizes) were found in the products. The average particle size of ZrC 0.93 prepared in this way was 21-23 nm. The BET surface area of the nano-(ZrC 0.93 , ZrO 2 )@C(graphite) composites prepared in He and H 2 was over 250 and 300 m 2 g −1 , respectively. We developed a reproducible and easy method to prepare nano-(ZrC, ZrO 2 )@C products by setting the DVB content, sulfonation degree, Zrloading and the thermal treatment conditions, which have an influence on the ZrC and graphite/amorphous carbon content of nano-ZrO 2 @C intermediates.