## 2 Materials and Methods
### Measurements  
The structure of the obtained polymers was confirmed by spectroscopic methods, using the ATR FTIR technique, using the BRUKER FTIR spectrophotometer TENSOR 27 with a resolution of 4 cm\({}^{-1}\). The measurements were carried out in the frequency range of 4000-600 cm\({}^{-1}\). Elemental analyses (CHN) were performed using a Perkin Elmer CHN 2400 analyzer (Palo Alto, CA, USA). \begin{table}
\begin{tabular}{c c c c c c} \hline \hline
**Polymer** & **DVB [mL]** & **GMA [mL]** & **t [mL]** & **d [mL]** & **b [mL]** \\ \hline pDVB t & 30 & - & 45 & - & - \\ DVB-GMA 2:1 t & 20 & 10 & 45 & - & - \\ DVB-GMA 1:1 t & 15 & 15 & 45 & - & - \\ DVB-GMA1:2 t & 10 & 20 & 45 & - & - \\ \hline pDVB db & 30 & - & - & 10 & 35 \\ DVB-GMA 2:1 db & 20 & 10 & - & 10 & 35 \\ DVB-GMA 1:1 db & 15 & 15 & - & 10 & 35 \\ DVB-GMA 1:2 db & 10 & 20 & - & 10 & 35 \\ \hline \hline \end{tabular} DVB–divinylbenzene, GMA–glycidyl methacrylate, t—toluene, d—decanol, b—benzyl alcohol.  
\end{table} 
Table 1: Volumes of the reagents used to the syntheses. Figure 1: Polymerization of DVB and GMA 1:1 [18]. TG/DTG curves were determined using a STA 449 F1 Jupiter thermal analyzer ([PERSON], Selb, Germany). Measurements were made in alumina crucibles, under a helium or synthetic air atmosphere, with a gas flow of 30 mL min\({}^{-1}\). Dynamic scans were carried out in the temperature range of 20-1000 \({}^{\circ}\)C, and the heating rate was 10 \({}^{\circ}\)C min\({}^{-1}\). The sample weight was about 7 mg and the measurements were made in the presence of an empty crucible as a reference. Images of polymeric microspheres were taken with optical microscope (Malvern Morphologi G3 with, Morphologi G3 software). The magnification was set on 10\(\times\). The porous structure of the studied polymers was investigated using the ASAP 2405 adsorption analyzer (Micromerities Inc., Norcross, GA, USA). Nitrogen adsorption-desorption measurements were performed at \(-\)190 \({}^{\circ}\)C, after degassing the samples at 70 \({}^{\circ}\)C. The specific surface area was calculated by the standard BET method (S\({}_{\text{BET}}\)), while the total pore volume (V\({}_{\text{tot}}\)) was evaluated as the volume of the liquid adsorbed at 0.99 relative pressure. The pore size distribution (PSD) was calculated from the adsorption branch of the isotherm using the 2D-NLDFT model, available as free version from the SAIEUS 3.0 software (Micromerities) on www.nldft.com. The efficiency of removing methylene blue from an aqueous solution was performed using a dynamic sorption method. For this purpose, SPE columns with a diameter of 1 cm and a height of 5 cm were used. The columns were filled with a 100 mg of the tested polymer placed between two Teflon disks. The columns were conditioned with 2.5 cm\({}^{3}\) of methanol and then 5 cm\({}^{3}\) of distilled water. Sorption studies were carried out at room temperature. An aqueous solutions of methylene blue (MB) of stock concentration equal to 1 mg/L were passed through the bed. The eluate was collected in separate fractions with a volume of 10 cm\({}^{3}\). The MB content in each fraction was determined photometrically using a UV-1800 spectrophotometer from Schimadzu (Canby, OR, USA). Changes in the concentration of methylene blue were determined on the basis of a calibration curve (\(\lambda\) = 666 nm). An increment of the MB concentration in each fraction expressed as a dependency C/C\({}_{0}\) = f(V) was used to plot the breakthrough curves. The sorption capacity \(q_{a}\) [mg/g] was calculated on the basis of Equation (1), where \(C_{0}\) is the concentration of stock solution, \(C_{i}\) is the concentration of the solution leaving the column, \(V_{i}\) is the volume [L] of a single fraction of methylene blue solution. \[q_{a}=\frac{\sum(C_{0}-C_{i})}{m}, \tag{1}\]  
Total mass adsorbed was calculated as a percentage value of mass adsorbed until the total breakthrough of the column in relation to the mass introduced into the column, as is shown in Equation (2).