Introduction

N-containing heterocycles (NCH) are a group of compounds with wide occurrence in the environment and broad application in chemistry [1,2,3]. They have been detected at industrial and agricultural sites [4, 5], groundwater [6,7,8,9,10], and surface water [11, 12]. They were also found in coal liquids, shale oil, surface waters such as lakes, and marine sediments [13]. The features of chemicals such as pharmaceuticals, industrial solvents, ionic liquids, and pesticides are based on NCH moieties [14,15,16]. Important moieties of, e.g., pharmaceuticals, are piperidine, pyridine, and imidazole [17]. Aromatic N-heterocycles such as pyrazole, pyridine, and pyrazine are also by-products of chemical processes such as acrylonitrile manufacturing [18].

Many NCHs are mutagenic, carcinogenic, and thus, hazards to the environment and human health [19], e.g., pyridine is toxic for several bacterial species and a strong odor compound [20]. It can cause weakness, lung and liver damage, and gastrointestinal inflammation in humans by acute exposure [15]. NCHs such as pyridine and imidazole are categorized as water hazard category 2 (significantly hazardous to the aquatic environment), and pyrazole is categorized as water hazard category 3 (highly hazardous to the aquatic environment) by German Environment Agency (Umweltbundesamt – UBA) (Table 1) [22]. Hence, UBA has set a provisionary water standard of 3 μg/L for pyrazole [24].

Table 1 Physical and chemical properties of selected N-heterocycles

Due to their basic properties (Table 1), NCHs are readily water soluble as cations at typical pH values of surface and groundwaters and thus mobile in the aquatic environment. NCHs are not directly photo-oxidized due to low absorption maxima [8], and some of them, such as pyrazole, are poorly biodegradable [25]).

To abate NCHs, different methods can be utilized [12, 15, 19]. Ozonation is one of the methods that is extensively used to eliminate NCHs, whether in wastewater [12] or process studies [26, 27]. Monitoring these methods and processes requires specific and sensitive measurement of NCHs.

The high affinity of NCHs towards the aqueous medium results in a challenge for classical reversed-phase liquid chromatography (RPLC) separation techniques [28]. For polar and small compounds such as imidazole, piperidine, and pyridine, there is hardly any retention on C18 columns [14, 29, 30]. This causes co-elution with non-retarded matrix components present in the sample, such as salts, and thus interferences in the detection of these NCHs (e.g., spectral interference (UV-vis detection) or ionization (MS-detection)) [31]. Moreover, solid-phase extraction (SPE) hardly enriches polar compounds such as NCHs [32,33,34]. Using ion chromatography (IC) [35] or ion pair chromatography [36] retention of NCHs was accomplished. However, the use of nonvolatile buffers, necessity for strong ion-pairing agents, and high water percentage render hyphenation with MS challenging, which is the detection method of choice.

Recently the use of mixed-mode liquid chromatography (MMLC) was suggested for the separation of polar and charged compounds as an alternative to hydrophilic interaction liquid chromatography (HILIC) [14, 37]. MMLC provides different modes of interactions with the analytes, such as hydrophobic interaction and ion exchange [38]. This can be used for separation of analytes with very different properties. For example, MMLC columns with anion or cation exchangers embedded in their hydrophobic alkyl chains [37] were used for the simultaneous measurement of cationic, zwitterionic, and neutral compounds previously [39]. Moreover, obstacles such as long equilibration time and use of organic solvent in the sample do not exist in MMLC compared to HILIC [33, 34]. However, up to now, MMLC was not used for the determination of NCHs.

The present work aims to develop an MMLC-MS method for the analysis of NCHs in water samples. It furthermore characterizes the separation mechanisms and employs the method in measurement of selected NCHs in the river water matrix. Applicability of the method in process monitoring will be tested for ozonation of the same water matrix.

Material and methods

Chemicals

All chemicals were commercially available and used as received (purity is presented in parenthesis): piperidine (ReagentPlus®, 99%) Sigma-Aldrich, pyridine (anhydrous, 99.8%) Sigma-Aldrich, imidazole (ACS reagent, ≥ 99%) Sigma-Aldrich, pyrazole (98%) Sigma-Aldrich, pyridazine (98 + %) Alfa Aesar. Water for chromatography (LC-MS Grade) LiChrosolv Merck Millipore, acetonitrile (ACN) HiPerSolv CHROMANORM® for LC-MS (≥ 99.9%,) VWR Chemicals, formic acid Suprapur (98–100%) Merck Millipore, nitrogen gas (> 99%) and oxygen gas (99.998%) from Air Liquid, Oberhausen.

Instrumentation

The following instruments were used: pH meter from Metrohm 827pH, TOC analyzer from Shimadzu TOC-L with ASI-L autosampler, and TOC Control-L Software. HPLC from Agilent Technologies (1100 Series) with Autosampler coupled with a Quadrupole LC/MS 6120 Mass spectrometer. LC-MS online software was used for data acquisition and processing. ESI was operated as follows: API-ES-positive mode, drying gas (nitrogen) temperature 250 °C, drying gas flow rate 11.0 L/min, nebulizer pressure 35 psi, capillary voltage 3.2 kV. The temperature of the column oven was kept constant at 30 °C in all measurements, and the flow rate was 300 μL/min with 5 μL injection volume. A Primesep 200 column (2.1 × 150 mm 5-μm particle size and 100-A pore size) and guard column (2.1 × 5 mm) from SIELC technology were used for separation. Single ion monitoring was used to measure real water samples. For further details on MS parameters, see Electronic Supplementary Material (ESM) Table S4.

Ozone was generated by continuous purging of ozone enriched O2 stream into ultrapure water, using a Philaqua 802 x ozone generator from BMT Messtechnik Berlin. The concentration of ozone was determined by the determination of UV-absorption (εozone at 258 nm = 2950 M−1 cm−1 [40]) using a Shimadzu UV-1800 spectrometer. The used absorption coefficient is different from the more recently published one of 3200 M−1 cm−1 [41]; however, it includes losses of ozone happening during the transfer of the ozone stock solution to the UV-vis spectrometer.

Chromatography

For optimization of the separation using MMLC, one has to determine the main separation modes. Based on reversed-phase (RP) interaction, an increase of the organic phase volume fraction (here ACN), and thus a decrease of aqueous volume fraction, results in a decrease of the retention factor (Eq. 1) [42].

$$\log \mathrm{k}^{\prime }=-\mathrm{S}\times {\updelta}_{\mathrm{A}\mathrm{CN}}+\log \mathrm{k}^{\prime}_{\mathrm{A}}$$
(1)

where k′ is the observed retention factor, δACN is the volume fraction of ACN, S is the slope the linear line, k′A is the retention factor when the volume fraction of ACN is zero, while for ion exchange (IE) chromatography, the retention factor is based on Eq. 2 [42].

$$\log \mathrm{k}^{\prime }=-\mathrm{Z}\times \log\ {C}_s+\log \mathrm{k}^{\prime}_z$$
(2)

where Cs is the concentration of counter ion (hydrogen ion), Z is the slope of the linear line, and k′z is a constant related to ion exchange properties of eluent and stationary phase [43].

To investigate the role of both IE and RP interactions for retention of different NCHs, a 1000 μg/L sample was analyzed using isocratic conditions with different δACN. In these experiments, the water was acidified with formic acid (FA), but ACN was not acidified. By mixing not acidified ACN with acidified water up to 60%, the final pH of the eluent will increase linearly according to Eq. 3 [44].

$${}_{\mathrm{w}}^{\mathrm{s}}\mathrm{pH}={}_{\mathrm{w}}^{\mathrm{w}}\mathrm{pH}+{m}_{\mathrm{pH}}\times {\updelta}_{\mathrm{ACN}}$$
(3)

where \({}_{\mathrm{w}}^{\mathrm{s}}\mathrm{pH}\) is the pH of the mixture, \({}_{\mathrm{w}}^\mathrm{w}\mathrm{pH}\) is the initial pH of the water before mixing with ACN, and mpH is the proportionality coefficient for the pH variation. mpH also depends on the concentration of acid and initial pH, which is discussed by Subirats et al. [44]. As a result of increasing pH, (decrease of formic acid concentration in the final mixture of eluent by increasing δACN and decrease of hydrogen ion as the counter ion), the IE interactions should become weaker by the increase of δACN in these experiments, i.e., the retention behavior will be contrary to RP interaction. The investigation of the separation performance was based on parameters shown in ESM Table S1.

Method validation

To investigate the linearity of the method, three series of 10 samples with concentrations 1, 2, 5, 10, 20, 50, 100, 200, 500, 1000 μg/L were chosen based on DIN 38402–51 [45] and analyzed. Mandel’s test was calculated based on [46]. In brief, an F test was performed to compare linear and second-order polynomial calibration according to the following equation:

$$F=\frac{\left(n-2\right)\times {S_{y/x, lin}}^2-\left(n-3\right)\times {S_{y/x, non}}^2}{{S_{y/x, non}}^2}$$
(4)

where n is the number of calibration points; Sy/x,lin and Sy/x,non are the standard error of the linear and nonlinear regression, respectively.

The method detection limit was quantified according to EPA [47] by analyzing nine different samples containing a mixture of all compounds, each with a concentration of 25 μg/L. Measurements were performed in 3 different batches in 3 different days and method detection limit (MDL) was calculated according to Eq. 5:

$$\mathrm{MDL}={t}_{\left(n-1,\kern0.75em 1-\alpha =0.99\right)}\times {S}_S$$
(5)

where t(n − 1, 0.99) is Student’s t value for a single-tailed 99th percentile (2.896 for n = 9), and Ss is the standard deviation of measured samples.

The performance of the method was investigated in river water. The river water was filtered by 0.45-μm cellulose acetate membranes and stored in the refrigerator at ≈ 5 °C until further use. Loss of the polar NCHs, e.g., due to sorption in this step, is negligible [33]. The sample had a TOC concentration of 3.53 mg/L and a pH of 7.87. Portions of 10 mL surface water samples were spiked with different volumes of the stock solution of the NCHs (mixture, dissolved in water, concentration: 10 mg/L) to prepare six different concentration levels. These samples were measured in three non-consecutive days, each in triplicate using freshly made eluent by the same operator. Recoveries were calculated as the average of all recoveries within the whole method validation studies. By utilizing analysis of variance (ANOVA), the variation of different performance parameters within a day (repeatability) and between days (intermediate precision) were calculated according to the Eurachem Guide [48].

The matrix effect was assessed by comparing the slope of spiked river water samples with ultrapure water samples according to Eq. 6.

$$\mathrm{Matrix}\ \mathrm{effect}=\left(\frac{\mathrm{Slope}\ \mathrm{in}\ \mathrm{spike}\ \mathrm{samples}}{\mathrm{Slope}\ \mathrm{in}\ \mathrm{ultrapure}\ \mathrm{water}}-1\right)\times 100$$
(6)

Results and discussion

Method development

Void time was determined by injecting NaBr and measurement of bromide at different ACN volume fractions (δACN) in negative ionization mode. The experiment was performed at the same chromatographic conditions used for determination of chromatographic behavior (section 2.2). As expected, no trend in void time was observed regarding δACN, ranging between 1.19 and 1.20 min.

Using the void time and retention times, retention factors were calculated at different δACN (Fig. 1). It was observed that an increase of δACN leads to a decrease in retention time and log(k′) for pyridazine and pyrazole. Therefore, these two compounds mainly interact via RP mode. It is also worth mentioning that at lower δACN pyrazole elutes later than pyridazine, while at higher fractions, the order of elution changes. On the other hand, the stronger bases imidazole, pyridine, and piperidine (c.f. pKa values, Table 1) display a first decrease of the retention factor until δACN reaches 0.25. At δACN > 0.25 retention factor increases. This indicates that at δACN ≤ 0.25, these three compounds interact via RP interactions and at higher δACN ≥ 0.25 via ion exchange interactions. This can be explained as follows. Pyridazine and pyrazole are mainly present as neutral species during the whole range of δACN, and thus, hydrophobic interactions may be most important. In the case of imidazole, pyridine, and piperidine, cationic species prevail in the whole range of applied δACN and can thus be separated by ion exchange interactions. The minimum of log(k´) δACN at 0.25 indicates that imidazole, pyridine, and piperidine also interact by the RP mode of the MM-column, which are, however, suppressed with increasing ACN fractions and cannot be observed at δACN > 0.25. The exemplary chromatogram in ESM Fig. S1 illustrates how the retention time of imidazole and pyrazole increases and decreases by the increase of δACN, respectively.

Fig. 1
figure 1

Logarithm of retention factor (log(k′)) in isocratic elution of NCHs at different ACN volume fractions and the pH of the final eluent mixture. Imidazole (triangles), pyrazole (circles), pyridine (stars), pyridazine (crosses), and piperidine (squares). Eluent A: water +0.1% FA, Eluent B: ACN. Flowrate 300 μL/min, 5 μL injection volume. Separation column: Primesep 200

Efficiency and asymmetry are almost constant for imidazole, pyridine, and piperidine for all δACN, except for 0.05 in which pyridine and piperidine show somewhat lower efficiency (ESM Figs. S2 and S3). On the other hand, for pyridazine and pyrazole, efficiency decreases with increasing δACN. Higher peak widths are also observed for imidazole, pyridine, and piperidine congruent with higher retention time (ESM Fig. S4). Moreover, by comparing the selectivity (ESM Fig. S5) with the resolution (ESM Fig. S6) one can observe that they follow the same trend. It appears that the role of selectivity in the calculation of resolution (formula presented in ESM Table S1) is greater than the retention factor (Fig. 1) and efficiency (ESM Fig. S2).

Figure 1 also represents the relation between log(k) and pH. According to Eq. 3 [44] the increase in the pH of final the eluent mixture is linear up to δACN = 0.6. Therefore, the retention of the analytes with IE interaction should also increase, as H+ competes with the cationic analytes for ion exchange moieties (Eq. 2). Moreover, pyrazole and pyridazine are mostly present as neutral species in the pH between 2.6 to 2.9 and thus interact with RP. On the other hand, imidazole, pyridine, and piperidine are present as cations in the pH range of the eluent mixture and thus, the role of IE interaction is relevant for them. The U-shaped correlation of pH vs. log(k) for piperidine, pyridine, and imidazole indicates that two interaction modes work at the same time (i.e., RP and IE interactions) while in the range of linear increase at higher δACN, IE interactions prevail. This type of U-shaped retention behavior is known to happen in HILIC phases [49] and was also previously observed in MMLC using other mixed-mode stationary phases [14, 50]. Montes et al. [38] described this as showing RP retention at low δACN and HILIC-like retention at high δACN for MMLC.

In another set of experiments, both components of the eluent, i.e., water and ACN, were acidified with 0.1% and 0.2% FA, respectively in order to investigate the suppression of IE at higher δACN. Compared to the eluent with non-acidified ACN, the increase of retention factor at δACN > 0.45 is less pronounced (Fig. 2). This can be explained by a constant elution strength for the IE mechanisms which corroborates that IE is indeed an important separation mechanism in the case of imidazole, pyridine, and piperidine. Substantial decrease of retention time by the increase of FA concentration for amines with pKa values of the corresponding acids above 8 was also previously observed by other researchers [39]. On the other hand, the hydrophilic interactions that are relevant in HILIC seem to be not completely suppressed, leading to the observed increase of the retention at δACN > 0.45. Contrary to RP and IE interactions (Eqs. 1 and 2), the relation between log(k′) and δACN in HILIC interactions is not first order. HILIC is modeled as having RP interaction in low δACN and normal phase (NP) interaction at high δACN according to [49].

$$\log \mathrm{k}^{\prime }={m}_1\times \left(1-{\updelta}_{\mathrm{ACN}}\right)-{m}_2\times \log\ \left(1-{\updelta}_{\mathrm{ACN}}\right)+\mathrm{constant}$$
(7)

where m1 and m2 are the empirical constants related to RP and NP interactions, respectively. Equation 4 can well describe the chromatographic behavior of all the compounds under investigation. Calculated parameters are presented in Table S2 and Fig. S7 (see ESM) represents the modeled data. All in all, it can be confirmed that NCHs that have basic properties follow a separation based on three interactions (RP/IE/HILIC). These multiple interactions provide many options to improve a separation from compounds which have fewer interactions with the stationary phase.

Fig. 2
figure 2

Comparison between logarithm of retention factor (log(k′)) using two different eluent compositions. (1) Eluent A: water +0.1% FA, Eluent B: ACN, imidazole (triangles) pyridine (stars), and piperidine (squares). (2) Eluent A: water + 0.1% FA, Eluent B: ACN + 0.2% FA imidazole (filled triangles), pyridine (crosses), and piperidine (filled squares). Flowrate 300 μL/min, 5 μL injection volume. Separation column: Primesep 200

Based on the above data, a method was developed for the detection of these NCHs in a surface water matrix to achieve proper retention and adequate separation from other constituents of the sample to avoid interferences during ESI [51]. The chromatographic method optimization was performed to achieve the best peak shape in a short runtime (k′, 2 to 5). As can be seen from Fig. 2, k′ of 5 (log(k′) ≈ 0.7) can be achieved for stronger bases only with the second set of eluents and at the middle δACN. Using different compositions of the mobile phase and gradient elution, the final method was developed according to ESM Table S3. MS measurements were done in two modes: in single ion monitoring (SIM) to improve the sensitivity and in scan mode (ESM Table S4).

Using the abovementioned conditions, a separation with parameters shown in Table 2 was achieved. Under these conditions, the efficiency of pyridine and piperidine separation largely increased, while the other compounds did not change much compared to the previous measurements. The method performance is very good regarding peak width (13 to 18 s) and asymmetry factors of maximum 1.65, particularly taking into account that the best peak shapes were obtained for the compounds with no retention in RP (i.e., imidazole, pyridine, piperidine).

Table 2 Method performance parameters for measurement of NCHs. Flowrate 300 μL/min, 5 μL injection volume. Separation column: Primesep 200, Gradient according to ESM Table S3. MS conditions according to ESM Table S4

Method validation

Table 3 shows Mandel’s test results for linearity and method detection limit (MDL). For all the compounds, calibrating by consideration of the highest investigated concentration (1000 μg/L) leads to a lack of linearity according to Mandel’s test. Therefore, 500 μg/L is the highest linear point in the logarithmic range.

Table 3 Method performance parameters for measurement of NCHs. Flowrate 300 μL/min, 5 μL injection volume. Separation column: Primesep 200, Gradient according to ESM Table S3. MS condition according to ESM Table S4

All the measured MDLs are between 3 and 6 μg/L, according to Table 3. It should be noted that no peaks were observed in blank samples; hence, no correction of MDL was required. For pyridine the detection limit of a previously developed method was reported to be 4.6 μg/L achieved by extraction of 1 L water into 1 mL methylene chloride followed by GC-MS [52]. Using electrical field–stimulated liquid-phase microextraction followed by HPLC-UV [53], 0.01 μg/L MDL was also achieved for pyridine. Pyrazole has been measured using vacuum-assisted headspace solid-phase microextraction followed by GC-MS (LOQ = 0.04 μg/L) [54], or by injecting sample in volumes as high as 90 μL in LC followed by MS in multiple reaction monitoring (0.05 μg/L) [55]. Utilizing conductivity detector and RP columns provided detection limits for ionic liquids of piperidinium, pyridinium, and imidazolium in the milligram per liter range [56, 57]. It can be stated that the MMLC method is similarly sensitive, considering that it could be improved by a factor of 1000 using an appropriate extraction method. The method standard deviation (residual standard deviation divided by slope) is similar for most of the investigated compounds. Pyridazine, however, showed an approximately three times higher method standard deviation than the average.

Table 4 shows method performance in river samples. The recoveries for imidazole, pyridine, and piperidine are more than 90%, while pyridazine showed a recovery of ca. 30%. In a similar manner, a relatively small matrix effect is observed in the measurement of NCHs except for pyridazine (ESM Fig. S8). One reason for the small matrix effect can be the low flow rate in this method, which can reduce matrix effects in ESI [51]. This can be further improved by post-column splitting [58] or using smaller dimensions of the column with proportionally lower flow rate. Moreover, using other ionization sources such as APCI [59] and APPI [60] might also help to reduce the matrix effect as they are less prone to interferences. The high matrix effect for pyridazine affected the overall recovery of the compound. It can be stated that the method is not performing well for the detection of pyridazine in the tested river water matrix.

Table 4 Recovery, repeatability, and intermediate precision for measurement NCHs in spiked surface water. Flowrate 300 μL/min, 5 μL injection volume. Separation column: Primesep 200, Gradient according to ESM Table S3. MS condition according to ESM Table S4

Moreover, higher intermediate precision in Table 4 can be attributed to a higher deviation from actual values in different days due to further variations of solvents, (room) temperatures, performance of the instrument (e.g., condition of seals, capillaries), etc. in comparison to triplicate injection within a day. While the intermediate precision for NCHs except pyridazine was ca. 5.2%, repeatability had a clear better performance for compounds with small or no tailing, i.e., imidazole, pyridine, and piperidine.

Pyridazine also showed worse performance in terms of retention time, peak width, and peak asymmetry with over 20% RSD (Table 5). Moreover, a difference between repeatability and intermediate precision of retention times was observed for all compounds due to the decrease in retention time. However, both within-a-day and between-days variation are very small in absolute terms, except for pyridazine. The highest intermediate precision (0.4% for imidazole) originates from the highest decrease of 4 s in retention time in three non-consecutive days. The peak width and asymmetry had no significant difference in intraday and interday measurements for NCHs other than pyridazine.

Table 5 Repeatability and intermediate precision of retention time, peak width, and peak asymmetry measured for NCHs in spiked surface water Flowrate 300 μL/min, 5 μL injection volume. Separation column: Primesep 200, Gradient according to ESM Table S3. MS condition according to ESM Table S4

Application in ozonation of NCH-containing matrix

As mentioned before, NCHs have to be abated during wastewater treatment and in drinking water since NCHs may contaminate drinking water resources and ozonation is one of the methods [11, 18]. A set of experiments was designed to follow the reduction of NCHs using ozone, in which river water was spiked to 100 μg/L NCHs each. Afterwards, the solution was divided into 10 mL samples in which ozone was dosed. The samples were analyzed with the developed LC-MS method at minimum 24 h after ozone was dosed for assuring full ozone depletion (note that typical lifetimes of ozone in ozonation are below 1 h). Figure 3 shows that ozonation has different efficiencies in the reduction of NCHs. Imidazole is the most efficiently eliminated compound, followed by pyrazole. Both compounds may be degraded at typical conditions of drinking water ozonation. The other compounds seem to be much more ozone refractory. Pyrazole, piperidine, and pyridine require elevated ozone dosages, which could result in undesired by-product formation. Pyridazine was not degraded in the present experiment, which covers a wide range of ozone dosages applied in drinking water. The results show that the developed method can be readily applied to analyze NCHs in real water samples and monitor their degradation in oxidation processes.

Fig. 3
figure 3

Relative degradation of NCHs with ozone in a surface water matrix spiked with a mixture of compounds with a concentration of 100 μg/L of each measured by LC-MS method. Imidazole (triangles, dotted line), pyrazole (circles, dashed line), pyridine (stars, long dashed line), pyridazine (crosses, long dash-dotted line), and piperidine (squares, long dash double dotted line). Error bars are the standard deviation of triplicate samples. Flowrate 300 μL/min, 5 μL injection volume. Separation column: Primesep 200, Gradient according to ESM Table S3. MS condition according to ESM Table S4. Note that recovery of pyridazine was 30%

Conclusion

The present work has shown that the MMLC is a complementary alternative for HILIC in the separation of polar compounds such as NCHs. Decent knowledge on the separation mechanisms largely helps to develop and optimize methods for their separation. The successful measurement of the selected NCHs should also facilitate the use of MMLC for the quantification of other NCHs and polar compounds that are difficult to retain in RPLC, such as 1,2,4-Triazole. Moreover, short run time, no sample preparation, and internal standard make this method, time-efficient, robust, and cheap. However, to further increase sensitivity for environmental monitoring, higher injection volumes or enrichment methods that were recently reported such as multi-layer solid-phase extraction [32], vacuum-assisted evaporative [61], or freeze-drying [62] might be necessary. It was demonstrated that the present method could be used to determine NCHs in trace concentrations in real water samples and to monitor their abatement during water treatment processes such as oxidation processes.