# The Potential of High Voltage Discharges for Green Solvent Extraction of Bioactive Compounds and Aromas from Rosemary (Rosmarinus officinalis L.)-Computational Simulation and Experimental Methods

Marinela Nutrizio, Jasenka Gajdoš Kljusurić, Zvonimir Marijanović, Igor Dubrović, Marko Viskić, Elena Mikolaj, Farid Chemat, Anet Režek Jambrak

► To cite this version:

Marinela Nutrizio, Jasenka Gajdoš Kljusurić, Zvonimir Marijanović, Igor Dubrović, Marko Viskić, et al. The Potential of High Voltage Discharges for Green Solvent Extraction of Bioactive Compounds and Aromas from Rosemary (Rosmarinus officinalis L.)-Computational Simulation and Experimental Methods. Molecules, 2020, 25 (16), pp.3711. 10.3390/molecules25163711. hal-03611566

HAL Id: hal-03611566

https://hal.inrae.fr/hal-03611566

Submitted on 17 Mar 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Article

# The Potential of High Voltage Discharges for Green Solvent Extraction of Bioactive Compounds and Aromas from Rosemary (Rosmarinus officinalis L.)—Computational Simulation and Experimental Methods

Marinela Nutrizio<sup>1,*<sup>ID</sup></sup>, Jasenka Gajdoš Kljusurić<sup>1<sup>ID</sup></sup>, Zvonimir Marijanović<sup>2</sup>, Igor Dubrović<sup>3<sup>ID</sup></sup>, Marko Viskić<sup>4</sup>, Elena Mikolaj<sup>1</sup>, Farid Chemat<sup>5<sup>ID</sup></sup> and Anet Režek Jambrak<sup>1,*<sup>ID</sup></sup>

Academic Editor: Marcello Locatelli

Received: 3 July 2020; Accepted: 12 August 2020; Published: 14 August 2020

<figure><img src="image_2.png" /></figure>

check for
updates

**Abstract:** Rosemary (Rosmarinus officinalis L.) is a Mediterranean medicinal and aromatic plant widely used due to valuable bioactive compounds (BACs) and aromas. The aim of the study was to evaluate the extraction of intracellular compounds from rosemary combining experimental procedure by means of high voltage electrical discharge (HVED), with a theoretical approach using two computational simulation methods: conductor-like screening model for real solvents and Hansen solubility parameters. The optimal HVED parameters were as follows: frequency 100 Hz, pulse width 400 ns, gap between electrodes 15 mm, liquid to solid ratio 50 mL/g, voltage 15 and 20 kV for argon, and 20 and 25 kV for nitrogen gas. Green solvents were used, water and ethanol (25% and 50%). The comparison was done with modified conventional extraction (CE) extracted by magnetic stirring and physicochemical analyses of obtained extracts were done. Results showed that HVED extracts in average 2.13-times higher total phenol content compared to CE. Furthermore, nitrogen, longer treatment time and higher voltage enhanced higher yields in HVED extraction. HVED was confirmed to have a high potential for extraction of BACs from rosemary. The computational stimulation methods were confirmed by experimental study, ethanol had higher potential of solubility of BACs and aromas from rosemary compared to water.

**Keywords:** high voltage electrical discharge; rosemary; COSMO-RS; Hansen solubility parameters; bioactive compounds; food aromas; extractions

## 1. Introduction

Rosemary (Rosmarinus officinalis L.) is an autochthonous Mediterranean herb from Lamiaceae family. From ancient times, rosemary has been used as flavoring agent and for medicinal purposes due to its intense aromatic odor and health benefits [1]. The biological activity of rosemary is mostly related to the phenolic compounds, such as carnosol, carnosic acid, and rosmarinic acid [2], and volatile compounds from essential oil like α-pinene, camphor, eucalyptol, or 1,8-cineole [3]. Due to valuable bioactive

<footnote n="1">1 Faculty of Food Technology and Biotechnology, University of Zagreb, 10000 Zagreb, Croatia; jasenka.gajdos@pbf.hr (J.G.K.); elena.mikolaj@yahoo.com (E.M.)</footnote>

<footnote n="2">2 Faculty of Chemistry and Technology, University of Split, 21000 Split, Croatia; zvonimir.marijanovic@ktf-split.hr</footnote>

<footnote n="3">3 Teaching Institute of Public Health of the Primorsko-goranska County, 51000 Rijeka, Croatia; igor.dubrovic@zzjzpgz.hr</footnote>

<footnote n="4">4 Faculty of Agriculture, University of Zagreb, 10000 Zagreb, Croatia; mviskic@agr.hr</footnote>

<footnote n="5">5 Université d'Avignon et des Pays du Vaucluse, 84000 Avignon, France; farid.chemat@univ-avignon.fr</footnote>

<footnote>* Correspondence: marinela.nutrizio@pbf.hr (M.N.); anet.rezek.jambrak@pbf.unizg.hr (A.R.J.); Tel.: +38-51-460-5287 (M.N. & A.R.J.); Fax: +38-51-483-6072 (M.N. & A.R.J.)</footnote>
compounds (BACs), rosemary possess antioxidant [2], anticancer [4], diuretic [5], antimicrobial [6], antiproliferative [7], anti-inflammatory [8], and anti-hyperglycemic properties [9]. Rosemary leaves, extracts, and essential oil have received recognition as generally recognized as safe (GRAS) for their intended use, from Food and Drug Administration [10] and according to Commission Directive 2010/67/EU and Commission Directive 2010/69/EU. Therefore, rosemary products can be a useful functional ingredient for the production of new functional foods [11].

Mostly used conventional extraction (CE) techniques are often associated with long extraction time, use of organic solvents in huge amounts and possible thermal degradation of thermosensitive compounds such as antioxidants [12]. For that reason, various innovative methods have been developed that are within the principles of green extraction processes. Apart from solvent extraction and steam distillation techniques, rosemary extracts have already been prepared by several green extraction methods: ultrasound-assisted extraction (UAE), microwave-assisted extraction (MAE) [13,14], supercritical fluid extraction (SFE) [15–17], and pressurized liquid extraction (PLE) [18–20]. For example, UAE was used for extraction of bioactive compounds from dried rosemary. Increased recovery of carnosic acid was obtained in ethanol, while rosmarinic acid was better extracted using methanol as a solvent. [21,22]. In a similar report, UAE of rosemary produced a three-fold increase in concentration of rosmarinic and carnosic acid when compared to the solid–liquid extraction [23]. In a work by Jacotet-Navaro et al., (2015), several extraction techniques were compared. Carnosic and ursolic acid extraction from rosemary was enhanced by ultrasound, while microwave extraction was more suitable for rosmarinic acid extraction. Both procedures were performed with reduced energy consumption and carbon emission when compared to heat reflux extraction. UAE was performed at 40 °C, while MAE was performed at 70–150 °C. The extraction yield increased with temperature but might have caused increased degradation and loss of volatile components, which makes thermal methods less suitable for extraction of aromatic compounds.

A green nonthermal extraction method that has not been investigated for extractions from rosemary is high voltage electrical discharge (HVED)—cold plasma treatment [24]. HVED extraction is a novel, eco-friendly extraction technique, that has been efficiently used in the extraction of BACs from various plant sources. Furthermore, as a nonthermal technology, extraction with HVED is performed in mild temperatures (usually at room temperature) and the temperature elevation after the extraction is low, so the thermal degradation of BACs is prevented [25,26]. The extraction of intracellular compounds by HVED is enhanced by the phenomenon of electrical breakdown in liquid that provokes cell structure damage and formation of pores (electroporation). The electrical breakdown is managed by the liquid ionization that is presented from the application of a high voltage between two electrodes with a gas flow. This leads to the liquid turbulence and intense mixing, emission of high-intensity ultraviolet light, generation of active radicals, production of shock waves and also a bubble cavitation [27].

Most conventional organic solvents are volatile, flammable, explosive, toxic, and have a negative environmental impact. For that reason, another principle of green extraction is directed towards use of alternative green solvents. In order to define green solvents, Gu and Jérôme (2013) proposed 12 criteria for green solvents that should be fulfilled related to availability, price, recyclability, grade, synthesis, toxicity, biodegradability, performance, stability, flammability, storage, and renewability. The solvent is considered greener when compared to conventional solvents that should be replaced with alternative solvent that fulfils at least some of the mentioned criteria [28]. Among green solvents, the water and agro- or bio-solvents play an important role for the replacement of organic solvents. Such solvents are derived from a renewable resource produced from biomass such as wood, starch, vegetable oils, or fruits. They have a high solvent power, are biodegradable, non-toxic, and non-flammable [29]. In a work by Barbieri et al., (2020), UAE of polyphenols from rosemary was performed in ethanol and natural deep eutectic solvents (NADES). High viscosity of NADES was reduced by the addition of 10% of water. The extraction efficiency was comparable to the results of ethanol extraction, while the antioxidant capacity of choline-based extracts was significantly improved [30].
It is important to choose an optimal solvent for the desired extraction, and to reduce solvent usage during experimentation. Traditional experimental procedures for solvent selection use high amount of solvents to choose which one is suitable for the extraction of compounds. For that reason, various computational simulations for solvent selection have been developed, such as Hansen solubility parameters (HSPs) and conductor-like screening model for real solvents (COSMO-RS). These models use theoretical predictions for assessing the solubility of targeted compounds in each solvent [31].

Since HVED treatment involves many processes, different reactive compounds are also formed. Beside generation of free radical species, there is also a possibility of electrodes abrasion and release of metals to the sample. For that reason, it is important to monitor preferred substances like phenols, antioxidants and volatile compounds, but also metals and other undesired compounds. Near infrared (NIR) spectroscopy has been utilized widely in the food and agribusiness ventures in the course of the last 20–30 years to decide significant parts in numerous agrarian items and plant materials [32–34]. In regard to determination of qualitative and quantitative characteristics of agricultural and food products, NIR offers a number of advantages over traditional analytical methods: it is a fast, physical, non-destructive, and non-invasive method, requires minimal or no sample preparation, no reagents are required, and no hazardous wastes are produced [35]. The wavelength range of the NIR spectroscopy from 750 to 2500 nm [36] is related to the vibration of molecules, especially the bands that are due to hydrogen (C-H, O-H, and N-H) vibrations [37,38]. But although the spectra recording is simple, user friendly (no additional sample preparation and use of chemicals) the interpretation of NIR spectra is very complex and chemometric methods are required to extract relevant information and reduce those that are less informative. The most common used tool is principal component analysis (PCA), quantitative analysis using multivariate calibration methods and qualitative analysis using multivariate classification techniques [34,38].

The aim of this work was to understand the green extraction of BACs and volatile compounds from rosemary using computational programs (HSPs and COSMO-RS) and experimental extractions. The experimental solvent extraction was carried out using HVED as a nonthermal technology. The obtained extracts were analyzed for physical and chemical (analytical) parameters. Furthermore, analyses of pesticides and metals in dry rosemary leaves and in HVED extracts were also performed.

## 2. Results

The extraction of BACs and volatile compounds from dried rosemary leaves were assessed by computational simulation methods and experimental analysis using HVED. Computational simulation methods were performed by HSPs and COSMO-RS where various green solvents were assessed in comparison with conventional solvent n-hexane in order to theoretically predict the probability of solution of BACs from rosemary. Experimental extraction method from rosemary was done using water and ethanol solutions (25 and 50%) as green solvents by means of HVED, and it was compared with modified CE under the same extraction conditions. The obtained extracts by CE and HVED were analyzed for physical parameters (pH, conductivity, temperature, and power), non-volatile compounds (total phenolic content (TPC), antioxidant capacity by 2-diphenyl-2-picrylhydrazyl (DPPH) free radical assay and ferric reducing antioxidant power (FRAP) method, NIR and ultra-performance liquid chromatography-tandem mass spectrometry (UPLC-MS/MS)), volatile compounds by headspace solid-phase microextraction/gas chromatography-mass spectrometry (HS-SPME/GC-MS), and metal content. Additionally, analyzes of pesticides and metals in dry rosemary leaves were performed to assess the safety of raw material for further processing and human consumption. The flowchart of all performed analysis from rosemary is presented in Figure 1.

### 2.1. Computational Simulation Methods for Assessing Solubility of Rosemary Compounds

The solubility parameters of green solvents for extraction of BACs and aromas from rosemary leaves have been studied by means of the HSP and COSMO-RS theoretical predictions. For both models, solubility results have been presented for various green solvents (ethyl acetate, methylacetate,
ethyl oleate, ethanol, 1-butanol, isopropanol, methanol, limonene, α-pinene, cymene, β-myrcene, cyclopentyl methyl ether (CPME), dimethyl carbonate, methyltetrahydrofuran (MeTHF), and water) compared to conventional *n*-hexane (first column). The HSPs for rosemary compounds were assessed at room temperature (20 °C) for different solvents and relative energy difference (RED) values are summarized in Table 1. RED results have been used to quantify the solutes–solvents interaction.

COSMO-RS combines quantum chemical considerations (COSMO) and statistical thermodynamics (RS) to determine and predict thermodynamic properties without experimental data. The computational simulation results derived by COSMO-RS for rosemary compounds at room temperature is presented in Table 2.

## 2.2. *Experimental Method for Extraction of Bioactive Compounds and Aromas from Rosemary by High Voltage Electrical Discharges*

Experimental method for assessing the solubility of BACs from rosemary was performed by green extraction method using HVED and green solvents water and ethanol (25 and 50%). For comparison of results, a modified CE by magnetic stirring was performed at same conditions as HVED: at room temperature, extraction times 3 and 9 min, ratio plant:solvent 1 g:50 mL. For all results in experimental procedure, “R” stands for rosemary, “N” for nitrogen, and “A” for argon. For HVED, numbers 1–12 show the order of conducted treatment. For CE treatments, 3 and 9 are referred to extraction time while 0, 25, and 50 stands for concentration of an ethanol solvent (%).

<figure><img src="image_3.png" /><figcaption>**Figure 1.** Flowchart of performed analysis in the paper.</figcaption></figure>
<table><caption>**Table 1.** Hansen Solubility Parameter (HSP) values of relative energy difference (RED) of bioactive compounds from rosemary for different solvents.</caption><thead><tr><th>Compounds</th><th>n-Hexane</th><th>Ethyl Acetate</th><th>Methyl Acetate</th><th>Ethyl Oleate</th><th>Ethanol</th><th>1-Butanol</th><th>Isopropanol</th><th>Methanol</th><th>Limonene</th><th>α-Pinene</th><th>Cymene</th><th>β-Myrcene</th><th>CPME</th><th>Dimethyl Carbonate</th><th>MeTHF</th><th>Water</th></tr></thead><tbody><tr><td colspan="17"><strong>Monoterpenes</strong></td></tr><tr><td>β-myrcene</td><td>1.04</td><td>1.45</td><td>2.77</td><td>0.1</td><td>4.5</td><td>3.39</td><td>3.56</td><td>5.57</td><td>0.67</td><td>0.32</td><td>0.67</td><td>0</td><td>0.5</td><td>2.36</td><td>0.91</td><td>10.49</td></tr><tr><td>p-cymen-7-ol</td><td>3.01</td><td>1.37</td><td>1.68</td><td>1.93</td><td>2.97</td><td>1.96</td><td>2.15</td><td>4.13</td><td>1.48</td><td>2.17</td><td>1.82</td><td>1.99</td><td>1.59</td><td>1.66</td><td>1.32</td><td>8.84</td></tr><tr><td>α-pinene</td><td>0.94</td><td>1.75</td><td>3.05</td><td>0.4</td><td>4.77</td><td>3.64</td><td>3.82</td><td>5.86</td><td>0.72</td><td>0</td><td>0.59</td><td>0.32</td><td>0.77</td><td>2.63</td><td>1.18</td><td>10.75</td></tr><tr><td>β-pinene</td><td>0.83</td><td>1.79</td><td>3.09</td><td>0.43</td><td>4.82</td><td>3.69</td><td>3.87</td><td>5.9</td><td>0.83</td><td>0.11</td><td>0.7</td><td>0.34</td><td>0.82</td><td>2.69</td><td>1.23</td><td>10.8</td></tr><tr><td>camphene</td><td>0.83</td><td>1.79</td><td>3.09</td><td>0.43</td><td>4.82</td><td>3.69</td><td>3.87</td><td>5.9</td><td>0.83</td><td>0.11</td><td>0.7</td><td>0.34</td><td>0.82</td><td>2.69</td><td>1.23</td><td>10.8</td></tr><tr><td>sabinene</td><td>1.03</td><td>1.74</td><td>3.08</td><td>0.42</td><td>4.78</td><td>3.67</td><td>3.85</td><td>5.87</td><td>0.7</td><td>0.15</td><td>0.47</td><td>0.37</td><td>0.74</td><td>2.6</td><td>1.14</td><td>10.77</td></tr><tr><td>α-phellandrene</td><td>1.19</td><td>1.48</td><td>2.75</td><td>0.25</td><td>4.46</td><td>3.34</td><td>3.52</td><td>5.56</td><td>0.46</td><td>0.32</td><td>0.51</td><td>0.22</td><td>0.5</td><td>2.38</td><td>0.9</td><td>10.44</td></tr><tr><td>α-terpinene</td><td>1.35</td><td>1.46</td><td>2.66</td><td>0.4</td><td>4.34</td><td>3.21</td><td>3.4</td><td>5.46</td><td>0.3</td><td>0.47</td><td>0.51</td><td>0.38</td><td>0.52</td><td>2.35</td><td>0.87</td><td>10.32</td></tr><tr><td>δ-terpinene</td><td>1.33</td><td>1.42</td><td>2.63</td><td>0.35</td><td>4.33</td><td>3.2</td><td>3.1</td><td>5.44</td><td>0.34</td><td>0.46</td><td>0.54</td><td>0.34</td><td>0.48</td><td>2.32</td><td>0.84</td><td>15.5</td></tr><tr><td colspan="17"><strong>Oxygenated monoterpenes</strong></td></tr><tr><td>camphor</td><td>1.88</td><td>1.42</td><td>2.93</td><td>0.99</td><td>4.49</td><td>3.5</td><td>3.67</td><td>5.52</td><td>1.02</td><td>1.14</td><td>0.73</td><td>1.06</td><td>0.83</td><td>1.99</td><td>0.87</td><td>10.45</td></tr><tr><td>borneol</td><td>2.28</td><td>0.91</td><td>1.72</td><td>1.22</td><td>3.33</td><td>2.23</td><td>2.42</td><td>4.46</td><td>0.9</td><td>1.5</td><td>1.3</td><td>1.28</td><td>0.9</td><td>1.63</td><td>0.71</td><td>9.3</td></tr><tr><td>α-terpineol</td><td>2.37</td><td>0.98</td><td>1.68</td><td>1.32</td><td>3.26</td><td>2.15</td><td>2.34</td><td>4.4</td><td>0.97</td><td>1.58</td><td>1.38</td><td>1.37</td><td>1.01</td><td>1.67</td><td>0.82</td><td>9.22</td></tr><tr><td>piperitone</td><td>1.82</td><td>1.18</td><td>2.62</td><td>0.8</td><td>4.22</td><td>3.19</td><td>3.37</td><td>5.28</td><td>0.74</td><td>1.02</td><td>0.63</td><td>0.88</td><td>0.55</td><td>1.86</td><td>0.57</td><td>10.19</td></tr><tr><td colspan="17"><strong>Sesquiterpenes</strong></td></tr><tr><td>β-caryophyllene</td><td>1.04</td><td>1.73</td><td>2.97</td><td>0.42</td><td>4.68</td><td>3.54</td><td>3.73</td><td>5.78</td><td>0.61</td><td>0.14</td><td>0.57</td><td>0.33</td><td>0.74</td><td>2.63</td><td>1.15</td><td>10.65</td></tr><tr><td colspan="17"><strong>Diterpenes</strong></td></tr><tr><td>carnosol</td><td>2.49</td><td>1.55</td><td>2.37</td><td>1.45</td><td>3.78</td><td>2.73</td><td>2.93</td><td>4.95</td><td>0.87</td><td>1.57</td><td>1.12</td><td>1.49</td><td>1.21</td><td>2.1</td><td>1.12</td><td>9.65</td></tr><tr><td>carnosic acid</td><td>2.99</td><td>1.69</td><td>2.05</td><td>1.94</td><td>3.29</td><td>2.26</td><td>2.46</td><td>4.48</td><td>1.38</td><td>2.1</td><td>1.72</td><td>1.98</td><td>1.66</td><td>2.09</td><td>1.48</td><td>9.09</td></tr><tr><td>rosmanol</td><td>2.95</td><td>1.42</td><td>1.78</td><td>1.88</td><td>3.07</td><td>2.05</td><td>2.24</td><td>4.24</td><td>1.4</td><td>2.09</td><td>1.74</td><td>1.93</td><td>1.56</td><td>1.78</td><td>1.32</td><td>8.93</td></tr><tr><td>epirosmanol</td><td>2.95</td><td>1.42</td><td>1.78</td><td>1.88</td><td>3.07</td><td>2.05</td><td>2.24</td><td>4.24</td><td>1.4</td><td>2.09</td><td>1.74</td><td>1.93</td><td>1.56</td><td>1.78</td><td>1.32</td><td>8.93</td></tr><tr><td>rosm radial</td><td>2.6</td><td>1.25</td><td>2.42</td><td>1.53</td><td>3.8</td><td>2.88</td><td>3.05</td><td>4.85</td><td>1.25</td><td>1.75</td><td>1.28</td><td>1.61</td><td>1.19</td><td>1.56</td><td>0.94</td><td>9.69</td></tr><tr><td colspan="17"><strong>Triterpenes</strong></td></tr><tr><td>betulinic acid</td><td>2.05</td><td>1.51</td><td>2.63</td><td>1.04</td><td>4.15</td><td>3.08</td><td>3.27</td><td>5.3</td><td>0.49</td><td>1.11</td><td>0.62</td><td>1.07</td><td>0.88</td><td>2.2</td><td>0.94</td><td>10.08</td></tr><tr><td>ursolic acid</td><td>2.12</td><td>1.54</td><td>2.6</td><td>1.11</td><td>4.11</td><td>3.03</td><td>3.23</td><td>5.26</td><td>0.53</td><td>1.18</td><td>0.71</td><td>1.14</td><td>0.94</td><td>2.22</td><td>0.99</td><td>10.02</td></tr><tr><td>rosmarinic acid</td><td>4.56</td><td>2.95</td><td>2.76</td><td>3.49</td><td>3.11</td><td>2.55</td><td>2.69</td><td>4.2</td><td>2.93</td><td>3.65</td><td>3.18</td><td>3.53</td><td>3.17</td><td>2.85</td><td>2.91</td><td>8.29</td></tr><tr><td colspan="17"><strong>Flavonoids</strong></td></tr><tr><td>apigenin</td><td>4.37</td><td>2.78</td><td>2.61</td><td>3.30</td><td>3.07</td><td>2.43</td><td>2.59</td><td>4.19</td><td>2.74</td><td>3.46</td><td>3.01</td><td>3.35</td><td>2.99</td><td>2.74</td><td>2.74</td><td>8.35</td></tr><tr><td>hispidulin</td><td>4.44</td><td>2.88</td><td>2.71</td><td>3.37</td><td>3.14</td><td>2.52</td><td>2.67</td><td>4.26</td><td>2.8</td><td>3.52</td><td>3.07</td><td>3.42</td><td>3.07</td><td>2.85</td><td>2.83</td><td>8.38</td></tr><tr><td>diosmetin</td><td>4.13</td><td>2.63</td><td>2.6</td><td>3.07</td><td>3.21</td><td>2.5</td><td>2.66</td><td>4.35</td><td>2.5</td><td>3.21</td><td>2.75</td><td>3.11</td><td>2.77</td><td>2.65</td><td>2.53</td><td>8.59</td></tr><tr><td>hesperidin</td><td>4.92</td><td>2.87</td><td>2.43</td><td>3.82</td><td>2.38</td><td>2.19</td><td>2.26</td><td>3.27</td><td>3.42</td><td>4.09</td><td>3.67</td><td>3.9</td><td>3.44</td><td>2.4</td><td>3.08</td><td>7.53</td></tr><tr><td>cirsimaritin</td><td>4.13</td><td>2.64</td><td>2.6</td><td>3.07</td><td>3.21</td><td>2.49</td><td>2.65</td><td>4.36</td><td>2.5</td><td>3.22</td><td>2.76</td><td>3.12</td><td>2.78</td><td>2.68</td><td>2.55</td><td>8.59</td></tr><tr><td>genkwanin</td><td>4.04</td><td>2.49</td><td>2.45</td><td>2.97</td><td>3.11</td><td>2.36</td><td>2.53</td><td>4.26</td><td>2.41</td><td>3.12</td><td>2.67</td><td>3.01</td><td>2.66</td><td>2.52</td><td>2.42</td><td>8.55</td></tr></tbody></table>

HSP: Relative energy difference (RED) very good solubility 0–1 (green color); medium solubility 1–3 (yellow color); poor solubility >3 (red color). CPME—cyclopentyl methyl ether, MeTHF—methyltetrahydrofuran.
<table><caption>**Table 2.** Conductor-like screening model for real solvents (COSMO-RS) probability of solubility (%) of bioactive compounds from rosemary for different solvents.</caption><thead><tr><th>Compounds</th><th>n-Hexane</th><th>Ethyl Acetate</th><th>Methyl Acetate</th><th>Ethyl Oleate</th><th>Ethanol</th><th>1-Butanol</th><th>Isopropanol</th><th>Methanol</th><th>Limonene</th><th>α-Pinene</th><th>Cymene</th><th>β-myrcene</th><th>CPME</th><th>Dimethyl-Carbonate</th><th>MeTHF</th><th>Water</th></tr></thead><tbody><tr><td colspan="18" style="text-align: left; font-weight: bold;">Monoterpenes</td></tr><tr><td>β-myrcene</td><td>69.18</td><td>81.28</td><td>60.26</td><td>100.00</td><td>11.22</td><td>20.89</td><td>16.98</td><td>4.17</td><td>95.50</td><td>81.28</td><td>100.00</td><td>100.00</td><td>100.00</td><td>39.81</td><td>100.00</td><td>0.00</td></tr><tr><td>α-pinene</td><td>99.08</td><td>42.66</td><td>26.92</td><td>100.00</td><td>10.47</td><td>22.91</td><td>17.38</td><td>3.63</td><td>95.50</td><td>100.00</td><td>83.18</td><td>87.10</td><td>89.13</td><td>16.60</td><td>85.11</td><td>0.00</td></tr><tr><td>β-pinene</td><td>95.50</td><td>57.54</td><td>34.67</td><td>100.00</td><td>12.59</td><td>25.70</td><td>19.95</td><td>4.57</td><td>97.95</td><td>99.08</td><td>89.13</td><td>89.13</td><td>97.72</td><td>22.39</td><td>95.50</td><td>0.00</td></tr><tr><td>camphene</td><td>97.72</td><td>51.29</td><td>34.67</td><td>100.00</td><td>12.59</td><td>25.70</td><td>19.95</td><td>4.68</td><td>97.95</td><td>99.31</td><td>89.13</td><td>91.20</td><td>97.72</td><td>22.39</td><td>93.33</td><td>0.00</td></tr><tr><td>sabinene</td><td>87.10</td><td>66.07</td><td>46.77</td><td>100.00</td><td>13.49</td><td>25.70</td><td>20.42</td><td>5.13</td><td>99.98</td><td>93.97</td><td>95.50</td><td>95.50</td><td>100.00</td><td>30.90</td><td>100.00</td><td>0.00</td></tr><tr><td>α-phellandrene</td><td>87.10</td><td>63.10</td><td>44.67</td><td>100.00</td><td>12.30</td><td>24.55</td><td>19.50</td><td>4.57</td><td>100.00</td><td>95.50</td><td>95.50</td><td>95.50</td><td>100.00</td><td>28.84</td><td>100.00</td><td>0.00</td></tr><tr><td>β-phellandrene</td><td>83.56</td><td>69.18</td><td>48.98</td><td>100.00</td><td>12.88</td><td>25.12</td><td>19.95</td><td>4.79</td><td>100.00</td><td>93.33</td><td>95.50</td><td>97.72</td><td>100.00</td><td>32.36</td><td>100.00</td><td>0.00</td></tr><tr><td colspan="18" style="text-align: left; font-weight: bold;">Oxygenated monoterpenes</td></tr><tr><td>camphor</td><td>48.54</td><td>86.98</td><td>70.07</td><td>99.25</td><td>44.06</td><td>70.94</td><td>58.92</td><td>21.45</td><td>85.15</td><td>63.09</td><td>91.76</td><td>89.39</td><td>95.89</td><td>53.75</td><td>100.00</td><td>0.10</td></tr><tr><td>borneol</td><td>11.22</td><td>100.00</td><td>89.13</td><td>81.28</td><td>85.11</td><td>100.00</td><td>100.00</td><td>41.69</td><td>17.78</td><td>13.49</td><td>16.98</td><td>16.60</td><td>100.00</td><td>38.02</td><td>100.00</td><td>0.02</td></tr><tr><td>α-terpineol</td><td>11.22</td><td>75.86</td><td>57.54</td><td>50.12</td><td>60.26</td><td>87.10</td><td>77.62</td><td>30.90</td><td>20.42</td><td>14.45</td><td>20.89</td><td>20.42</td><td>97.72</td><td>31.62</td><td>100.00</td><td>0.02</td></tr><tr><td>piperitone</td><td>32.36</td><td>100.00</td><td>89.13</td><td>85.11</td><td>91.20</td><td>100.00</td><td>100.00</td><td>51.29</td><td>72.44</td><td>46.77</td><td>87.10</td><td>83.18</td><td>87.10</td><td>70.79</td><td>100.00</td><td>0.07</td></tr><tr><td colspan="18" style="text-align: left; font-weight: bold;">Sesquiterpenes</td></tr><tr><td>β-caryophyllene</td><td>99.95</td><td>53.70</td><td>33.11</td><td>100.00</td><td>7.94</td><td>18.62</td><td>14.13</td><td>2.34</td><td>100.00</td><td>100.00</td><td>89.13</td><td>89.54</td><td>100.00</td><td>18.20</td><td>100.00</td><td>0.00</td></tr><tr><td colspan="18" style="text-align: left; font-weight: bold;">Diterpenes</td></tr><tr><td>carnosol</td><td>2.86</td><td>100.00</td><td>100.00</td><td>60.81</td><td>100.00</td><td>100.00</td><td>100.00</td><td>58.48</td><td>10.59</td><td>4.67</td><td>13.46</td><td>11.99</td><td>100.00</td><td>100.00</td><td>100.00</td><td>0.00</td></tr><tr><td>carnosic acid</td><td>0.86</td><td>100.00</td><td>100.00</td><td>100.00</td><td>100.00</td><td>100.00</td><td>100.00</td><td>100.00</td><td>3.37</td><td>1.41</td><td>3.89</td><td>3.49</td><td>100.00</td><td>100.00</td><td>100.00</td><td>0.00</td></tr><tr><td>rosmanol</td><td>0.69</td><td>100.00</td><td>100.00</td><td>41.70</td><td>100.00</td><td>100.00</td><td>100.00</td><td>100.00</td><td>3.22</td><td>1.23</td><td>4.33</td><td>3.76</td><td>100.00</td><td>100.00</td><td>100.00</td><td>0.00</td></tr><tr><td>epirosmanol</td><td>0.69</td><td>100.00</td><td>100.00</td><td>41.70</td><td>100.00</td><td>100.00</td><td>100.00</td><td>100.00</td><td>3.22</td><td>1.23</td><td>4.33</td><td>3.76</td><td>100.00</td><td>100.00</td><td>100.00</td><td>0.00</td></tr><tr><td>rosmadal</td><td>0.58</td><td>100.00</td><td>100.00</td><td>38.89</td><td>100.00</td><td>100.00</td><td>100.00</td><td>96.82</td><td>4.06</td><td>1.22</td><td>6.35</td><td>5.58</td><td>100.00</td><td>100.00</td><td>100.00</td><td>0.00</td></tr><tr><td colspan="18" style="text-align: left; font-weight: bold;">Triterpenes</td></tr><tr><td>betulinic acid</td><td>6.64</td><td>100.00</td><td>100.00</td><td>48.56</td><td>68.99</td><td>93.73</td><td>90.54</td><td>22.83</td><td>17.90</td><td>9.08</td><td>19.53</td><td>17.24</td><td>100.00</td><td>45.69</td><td>100.00</td><td>0.00</td></tr><tr><td>ursolic acid</td><td>1.07</td><td>100.00</td><td>100.00</td><td>83.73</td><td>100.00</td><td>100.00</td><td>100.00</td><td>100.00</td><td>2.78</td><td>1.39</td><td>2.59</td><td>2.32</td><td>100.00</td><td>45.15</td><td>100.00</td><td>0.00</td></tr><tr><td>rosmarinic acid</td><td>0.00</td><td>100.00</td><td>100.00</td><td>34.92</td><td>100.00</td><td>100.00</td><td>100.00</td><td>100.00</td><td>0.07</td><td>0.01</td><td>0.13</td><td>0.11</td><td>100.00</td><td>100.00</td><td>100.00</td><td>0.02</td></tr><tr><td colspan="18" style="text-align: left; font-weight: bold;">Flavonoids</td></tr><tr><td>apigenin</td><td>0.00</td><td>100.00</td><td>100.00</td><td>25.70</td><td>100.00</td><td>100.00</td><td>100.00</td><td>100.00</td><td>0.02</td><td>0.00</td><td>0.03</td><td>0.03</td><td>100.00</td><td>100.00</td><td>100.00</td><td>0.18</td></tr><tr><td>hispidulin</td><td>0.02</td><td>100.00</td><td>100.00</td><td>25.13</td><td>100.00</td><td>100.00</td><td>100.00</td><td>100.00</td><td>0.26</td><td>0.05</td><td>0.47</td><td>0.41</td><td>100.00</td><td>100.00</td><td>100.00</td><td>0.00</td></tr><tr><td>diosmetin</td><td>0.00</td><td>100.00</td><td>100.00</td><td>23.44</td><td>100.00</td><td>100.00</td><td>100.00</td><td>100.00</td><td>0.02</td><td>0.00</td><td>0.03</td><td>0.03</td><td>100.00</td><td>100.00</td><td>100.00</td><td>0.07</td></tr><tr><td>hesperidin</td><td>0.00</td><td>100.00</td><td>100.00</td><td>0.72</td><td>100.00</td><td>100.00</td><td>100.00</td><td>100.00</td><td>0.00</td><td>0.00</td><td>0.00</td><td>0.00</td><td>100.00</td><td>100.00</td><td>100.00</td><td>0.00</td></tr><tr><td>cirsimaritin</td><td>0.03</td><td>100.00</td><td>100.00</td><td>13.12</td><td>100.00</td><td>100.00</td><td>100.00</td><td>100.00</td><td>0.38</td><td>0.09</td><td>0.71</td><td>0.63</td><td>100.00</td><td>100.00</td><td>100.00</td><td>0.00</td></tr><tr><td>genkwanin</td><td>0.03</td><td>100.00</td><td>100.00</td><td>13.48</td><td>100.00</td><td>100.00</td><td>100.00</td><td>100.00</td><td>0.35</td><td>0.08</td><td>0.65</td><td>0.57</td><td>100.00</td><td>100.00</td><td>100.00</td><td>0.00</td></tr></tbody></table>

COSMO-RS: Low probability of solubility 0–20% (red color); medium probability of solubility 20–60% (yellow color); high probability of solubility 60–100% (green color). CPME—cyclopentyl methyl ether, MeTHF—methyltetrahydrofur.
Since the aim of the extraction with HVED was to achieve electrical discharge that is responsible for plant cell disruption and consequently the extraction of BACs from intracellular area. During the experiments, it was difficult to achieve discharge using nitrogen under 20 kV, therefore, for nitrogen treatments voltage of 20 and 25 kV were chosen, while lower voltages of 15 and 20 kV were obtained with argon treatments. All results were measured in duplicates and are presented as average ± standard deviation (SD).

The experiment design performed in STATGRAPHICS Centurion software is presented in Table 3.

<table><caption>**Table 3.** Denotation of samples, experimental design, and process parameters.</caption><thead><tr><th>Sample</th><th>Treatment Time (min)</th><th>Voltage (kV)</th><th>Ethanol Content (%)</th><th>Stirring Time (min)</th><th>Extraction Type</th></tr></thead><tbody><tr><td>3 R0</td><td>0</td><td>0</td><td>0</td><td>3</td><td rowspan="6">CE</td></tr><tr><td>9 R0</td><td>0</td><td>0</td><td>0</td><td>9</td></tr><tr><td>3 R25</td><td>0</td><td>0</td><td>25</td><td>3</td></tr><tr><td>9 R25</td><td>0</td><td>0</td><td>25</td><td>9</td></tr><tr><td>3 R50</td><td>0</td><td>0</td><td>50</td><td>3</td></tr><tr><td>9 R50</td><td>0</td><td>0</td><td>50</td><td>9</td></tr><tr><td>RN1</td><td>3</td><td>20</td><td>50</td><td>/</td><td rowspan="12">HVED</td></tr><tr><td>RN2</td><td>9</td><td>20</td><td>0</td><td>/</td></tr><tr><td>RN3</td><td>3</td><td>20</td><td>0</td><td>/</td></tr><tr><td>RN4</td><td>3</td><td>25</td><td>0</td><td>/</td></tr><tr><td>RN5</td><td>9</td><td>25</td><td>25</td><td>/</td></tr><tr><td>RN6</td><td>9</td><td>20</td><td>25</td><td>/</td></tr><tr><td>RN7</td><td>9</td><td>20</td><td>50</td><td>/</td></tr><tr><td>RN8</td><td>9</td><td>25</td><td>50</td><td>/</td></tr><tr><td>RN9</td><td>3</td><td>25</td><td>25</td><td>/</td></tr><tr><td>RN10</td><td>9</td><td>25</td><td>0</td><td>/</td></tr><tr><td>RN11</td><td>3</td><td>25</td><td>50</td><td>/</td></tr><tr><td>RN12</td><td>3</td><td>20</td><td>25</td><td>/</td></tr><tr><td>RA1</td><td>3</td><td>15</td><td>50</td><td>/</td></tr><tr><td>RA2</td><td>9</td><td>15</td><td>0</td><td>/</td></tr><tr><td>RA3</td><td>3</td><td>15</td><td>0</td><td>/</td></tr><tr><td>RA4</td><td>3</td><td>20</td><td>0</td><td>/</td></tr><tr><td>RA5</td><td>9</td><td>20</td><td>25</td><td>/</td></tr><tr><td>RA6</td><td>9</td><td>15</td><td>25</td><td>/</td></tr><tr><td>RA7</td><td>9</td><td>15</td><td>50</td><td>/</td></tr><tr><td>RA8</td><td>9</td><td>20</td><td>50</td><td>/</td></tr><tr><td>RA9</td><td>3</td><td>20</td><td>25</td><td>/</td></tr><tr><td>RA10</td><td>9</td><td>20</td><td>0</td><td>/</td></tr><tr><td>RA11</td><td>3</td><td>20</td><td>50</td><td>/</td></tr><tr><td>RA12</td><td>3</td><td>15</td><td>25</td><td>/</td></tr></tbody></table>

### 2.2.1. Determination of Physical Parameters of Rosemary Extracts

Results of physical parameters of CE and HVED treatment are given Figure 2, including pH, conductivity, power and temperature before and after the HVED treatment.
<figure><img src="image_4.png" /><figcaption>**Figure 2.** Values of pH, conductivity (mS/cm) for CE and HVED treated samples, starting temperature (°C), final temperature (°C) and power (kW) after HVED treatments.</figcaption></figure>

### 2.2.2. Determination of Phenols and Antioxidant Activity of Rosemary Extracts

The difference between CE and HVED extraction according to results of TPC, antioxidant parameters and extraction yields are presented in Figure 3. TPC results ranged from 7.21 to 31.64 mg GAE/g of sample. Antioxidant capacity was measured with two different methods—DPPH and FRAP. DPPH ranged from 25.85 to 32.92 µmol TE/g of sample, while FRAP ranged from 44.07 to 562.64 µmol FE/g of sample. Yield of extraction was calculated as g GAE/g of sample × 100 (%).

### 2.2.3. Near Infrared Spectroscopy and Qualitative Modeling

The recorded NIR spectra in wavelength range from 904 to 1699 nm were used in the qualitative modelling using principal component analysis (PCA) to identify potential grouping (Figure 4) where the sample 9 R50 seemed to be an outlier. Additionally, the Grubbs test was conducted and this sample was confirmed as an outlier.

Furthermore, modeling with different wavelength ranges was performed (Table S1) using the partial least squares regression (PLSR). Four different wavelength ranges were used for four models: Model 1: $\lambda$ = 904–1699 nm; Model 2: $\lambda$ = 1349–1699 nm; Model 3: $\lambda$ = 904–932 and 1349–1699 nm; and Model 4: $\lambda$ = 904–932 nm. Model efficacy was evaluated using the coefficient of determination ($R^2$) and the regression point displacement that is the ratio of the standard error of performance (RPD) and the ratio of the range of reference chemistry values to standard error of prediction (RER). RPD is the ratio of standard deviation of the validation data set (SDv) and the standard error of prediction (SEP). How efficient this quantitative prediction of TPC, FRAP, and DPPH is on the rest of 40 % of the samples is presented in Figure 5.

### 2.2.4. Determination of Color of Rosemary Extracts

The colorimetric analyses of rosemary extracts were measured under International Commission on Illumination (CIE)—L*a*b* color system. Results are shown in Table 4. Results of $\Delta C$, $\Delta E$, and $\Delta H$ are
presented for HVED extracts as a result of total color difference, difference in tone color and difference in saturation, respectively, compared to same conditions with CE.

<figure><img src="image_5.png" /><figcaption>**Figure 3.** Determination of bioactive compounds—total phenolic compounds (TPC) values, antioxidant activity (2-diphenyl-2-picrylhydrazyl (DPPH) and ferric reducing ability of plasma (FRAP) and yield of extraction—measurements for CE and HVED treated samples. TPC = total phenolic content, DPPH = 2,2-diphenyl-2-picrylhydrazyl free radical assay, FRAP = ferric reducing ability of plasma; Treatment type: CE—conventional extraction, N₂—HVED treatment with nitrogen, Ar—HVED treatment with argon.</figcaption></figure>

<figure><img src="image_6.png" /><figcaption>**Figure 4.** The principal component analysis (PCA) of rosemary extracts: (a) The data is denoting untreated rosemary samples, and rosemary samples treated with HVED using nitrogen; (b) The data is denoting untreated rosemary samples, and rosemary samples treated with HVED using argon.</figcaption></figure>
<figure><img src="image_7.png" /><figcaption>**Figure 5.** Control chart for predicted and experimental data of analyzed samples, using principal component regression (PCR) models.</figcaption></figure>

<table><caption>**Table 4.** International Commission on Illumination (CIE)—L*a*b* color parameters of CE and HVED treated rosemary extracts.</caption><thead><tr><th>Sample</th><th>L*</th><th>a*</th><th>b*</th><th>C</th><th>h</th><th>ΔC</th><th>ΔE</th><th>ΔH</th><th>Extraction Type</th></tr></thead><tbody><tr><td>3 R0</td><td>86.01 ± 2.01</td><td>3.87 ± 0.52</td><td>35.77 ± 1.79</td><td>35.98 ± 1.03</td><td>1.46 ± 0.02</td><td>/</td><td>/</td><td>/</td><td rowspan="5">CE</td></tr><tr><td>9 R0</td><td>82.09 ± 3.64</td><td>7.09 ± 1.07</td><td>44.15 ± 0.61</td><td>44.72 ± 1.56</td><td>1.41 ± 0.43</td><td>/</td><td>/</td><td>/</td></tr><tr><td>3 R25</td><td>95.13 ± 1.97</td><td>-0.77 ± 0.03</td><td>15.26 ± 3.72</td><td>15.28 ± 0.82</td><td>-1.52 ± 0.06</td><td>/</td><td>/</td><td>/</td></tr><tr><td>9 R25</td><td>93.55 ± 2.27</td><td>-0.41 ± 0.04</td><td>20.56 ± 4.08</td><td>20.56 ± 1.35</td><td>-1.55 ± 0.15</td><td>/</td><td>/</td><td>/</td></tr><tr><td>3 R50</td><td>95.32 ± 4.31</td><td>-1.16 ± 0.37</td><td>13.79 ± 1.06</td><td>13.84 ± 1.10</td><td>-1.49 ± 0.09</td><td>/</td><td>/</td><td>/</td></tr><tr><td>9 R50</td><td>94.11 ± 2.55</td><td>-0.77 ± 0.15</td><td>15.79 ± 2.35</td><td>15.81 ± 2.74</td><td>-1.52 ± 0.07</td><td>/</td><td>/</td><td>/</td></tr></tbody></table>
<table><caption>Table 4. *Cont.*</caption><thead><tr><th>Sample</th><th>L*</th><th>a*</th><th>b*</th><th>C</th><th>h</th><th>ΔC</th><th>ΔE</th><th>ΔH</th><th>Extraction Type</th></tr></thead><tbody><tr><td>RN1</td><td>88.56 ± 1.96</td><td>-1.81 ± 1.06</td><td>43.33 ± 1.71</td><td>43.37 ± 2.67</td><td>-1.53 ± 0.16</td><td>29.53</td><td>30.31</td><td>1.03</td><td rowspan="14">HVED</td></tr><tr><td>RN2</td><td>80.89 ± 1.74</td><td>6.58 ± 1.21</td><td>49.65 ± 2.09</td><td>50.08 ± 1.65</td><td>1.44 ± 0.27</td><td>5.37</td><td>5.65</td><td>1.30</td></tr><tr><td>RN3</td><td>83.86 ± 0.72</td><td>4.45 ± 0.47</td><td>43.24 ± 2.52</td><td>43.47 ± 1.79</td><td>1.47 ± 0.20</td><td>7.49</td><td>7.79</td><td>0.21</td></tr><tr><td>RN4</td><td>84.41 ± 1.82</td><td>4.38 ± 0.69</td><td>43.40 ± 1.79</td><td>43.62 ± 3.40</td><td>1.47 ± 0.13</td><td>7.64</td><td>7.81</td><td>0.28</td></tr><tr><td>RN5</td><td>92.19 ± 4.18</td><td>-1.70 ± 0.82</td><td>28.03 ± 0.67</td><td>28.08 ± 1.06</td><td>-1.51 ± 0.06</td><td>7.52</td><td>7.70</td><td>0.98</td></tr><tr><td>RN6</td><td>92.22 ± 2.33</td><td>-2.38 ± 0.04</td><td>28.18 ± 1.38</td><td>28.28 ± 1.75</td><td>-1.49 ± 0.00</td><td>7.72</td><td>7.98</td><td>1.55</td></tr><tr><td>RN7</td><td>92.42 ± 1.79</td><td>-5.22 ± 0.57</td><td>36.22 ± 1.64</td><td>36.59 ± 0.89</td><td>-1.43 ± 0.14</td><td>20.79</td><td>20.98</td><td>2.27</td></tr><tr><td>RN8</td><td>92.52 ± 0.64</td><td>-4.07 ± 0.03</td><td>32.48 ± 0.82</td><td>32.73 ± 2.07</td><td>-1.45 ± 0.01</td><td>16.93</td><td>17.09</td><td>1.73</td></tr><tr><td>RN9</td><td>92.40 ± 2.87</td><td>-5.28 ± 0.50</td><td>37.19 ± 2.07</td><td>37.56 ± 1.46</td><td>-1.43 ± 0.03</td><td>22.28</td><td>22.55</td><td>2.17</td></tr><tr><td>RN10</td><td>81.96 ± 2.91</td><td>6.35 ± 0.73</td><td>47.41 ± 4.23</td><td>47.83 ± 2.37</td><td>1.44 ± 0.12</td><td>3.12</td><td>3.35</td><td>1.21</td></tr><tr><td>RN11</td><td>91.92 ± 5.01</td><td>-4.42 ± 0.06</td><td>36.02 ± 3.16</td><td>36.29 ± 1.68</td><td>-1.45 ± 0.07</td><td>22.45</td><td>22.72</td><td>0.86</td></tr><tr><td>RN12</td><td>95.30 ± 2.69</td><td>-1.86 ± 0.07</td><td>18.18 ± 1.06</td><td>18.27 ± 0.69</td><td>-1.47 ± 0.06</td><td>3.00</td><td>3.12</td><td>0.86</td></tr><tr><td>RA1</td><td>92.23 ± 3.15</td><td>-2.04 ± 0.00</td><td>28.69 ± 2.74</td><td>28.76 ± 1.41</td><td>-1.50 ± 0.04</td><td>14.92</td><td>15.24</td><td>0.26</td></tr><tr><td>RA2</td><td>82.44 ± 1.82</td><td>5.79 ± 1.14</td><td>43.08 ± 1.95</td><td>43.47 ± 2.38</td><td>1.44 ± 0.19</td><td>-1.25</td><td>1.72</td><td>1.13</td></tr><tr><td>RA3</td><td>86.00 ± 2.47</td><td>3.43 ± 0.97</td><td>36.88 ± 0.56</td><td>37.04 ± 1.06</td><td>1.48 ± 0.02</td><td>1.06</td><td>1.19</td><td>0.55</td></tr><tr><td>RA4</td><td>84.95 ± 0.53</td><td>3.72 ± 0.38</td><td>38.07 ± 2.03</td><td>38.25 ± 0.76</td><td>1.47 ± 0.14</td><td>2.27</td><td>2.54</td><td>0.38</td></tr><tr><td>RA5</td><td>90.25 ± 2.74</td><td>-0.73 ± 0.05</td><td>30.17 ± 1.73</td><td>30.18 ± 0.19</td><td>-1.55 ± 0.12</td><td>9.61</td><td>10.17</td><td>0.11</td></tr><tr><td>RA6</td><td>90.89 ± 3.17</td><td>-0.76 ± 0.14</td><td>28.81 ± 2.49</td><td>28.82 ± 1.03</td><td>-1.54 ± 0.10</td><td>8.26</td><td>8.68</td><td>0.16</td></tr><tr><td>RA7</td><td>92.36 ± 1.56</td><td>-2.91 ± 0.09</td><td>29.96 ± 0.86</td><td>30.10 ± 2.07</td><td>-1.47 ± 0.06</td><td>14.29</td><td>14.44</td><td>1.05</td></tr><tr><td>RA8</td><td>91.24 ± 2.40</td><td>-2.84 ± 0.16</td><td>33.36 ± 1.78</td><td>33.48 ± 1.56</td><td>-1.49 ± 0.15</td><td>17.67</td><td>17.92</td><td>0.83</td></tr><tr><td>RA9</td><td>92.84 ± 3.77</td><td>-1.43 ± 0.62</td><td>24.87 ± 0.93</td><td>24.91 ± 0.59</td><td>-1.51 ± 0.08</td><td>9.63</td><td>9.90</td><td>0.14</td></tr><tr><td>RA10</td><td>81.04 ± 1.82</td><td>5.72 ± 0.11</td><td>43.50 ± 2.00</td><td>43.87 ± 2.24</td><td>1.44 ± 0.20</td><td>-0.84</td><td>1.84</td><td>1.26</td></tr><tr><td>RA11</td><td>94.56 ± 3.09</td><td>-2.81 ± 0.23</td><td>23.50 ± 1.37</td><td>23.67 ± 1.38</td><td>-1.45 ± 0.04</td><td>9.83</td><td>9.88</td><td>0.63</td></tr><tr><td>RA12</td><td>93.93 ± 2.66</td><td>-1.57 ± 0.07</td><td>21.44 ± 1.08</td><td>21.50 ± 1.59</td><td>-1.50 ± 0.13</td><td>6.22</td><td>6.35</td><td>0.41</td></tr></tbody></table>

L*—lightness from black to white; a* from green to red, and b* from blue to yellow; C—chroma; h—hue angle; ΔE—total color difference, ΔC—difference in chroma; ΔH—difference in hue.
### 2.2.5. Principal Component Analysis (PCA) of Rosemary Extracts

Data matrix used to identify similarities and/or differences in the data set presenting the samples as well as the physical and chemical properties. The PCAs ability to reduce the dimensionality and increasing the interpretability with minimal information lost was used on the data matrix of experimental data. The matrix included as active variables the physical parameters, total phenolic content and antioxidant capacity, as well as the parameters of color, while the supplementary data set were the experimental conditions (extraction type, HVED treatment time, stirring time, treatment time, voltage, and the ethanol content). The Extraction type (CE or HVED) were included in the analysis as qualitative variables (Figure 6). To present the differences in the samples which are extracted by CE or HVED, in the Supplementary Information are added boxplots for the TPC and AOX by use of DPPH and FRAP method and the yield, as well as the PCA biplot for different extractions and the physical-chemical properties Figures S3 and S4).

Modeling that has followed included the steps of calibration, validation and prediction. In the modelling were included the NIR spectra what is in detail explained in the Section 4.5. Calibration model was developed by principal component regression (PCR). From the data matrix 2/4 of it was used for the calibration and 1/4 was used for the validation while the remaining ¼ was used for the prediction. Validation was done by K-fold cross-validation.

### 2.2.6. Ultra-Performance Liquid Chromatography–Tandem Mass Spectrometry (UPLC–MS/MS) Analysis of Phenolic Compounds from Rosemary Extracts

UPLC-MS/MS analysis was performed to quantify individual phenolic compounds (apigenin, carnosol, diosmetin, hydroxytyrosol, luteolin, oleanolic acid, quercetin, rosmarinic acid, p-cymene, camphor, thymol, and carvacrol) from rosemary extracts (Table 5).

### 2.2.7. Determination of Volatile Compounds from Rosemary Extracts

Determination of main volatile compounds from rosemary extracts was performed by HS-SPME/GC-MS and analysis included eucalyptol, camphor, borneol and linalool. Results of performed analysis are presented in Table 6.

### 2.2.8. Determination of Pesticides and Metals

Results of pesticides and metals measured from dried rosemary leaves are shown in Table 7. Also, table includes results of heavy metals measured in selected rosemary extracts. Maximum residue levels (MRLs) are given according to European Commission (EC) Regulations EC No. 396/2005 for pesticides for a rosemary plant and No. 1881/200 for metals for food supplements with rosemary, since rosemary is not listed in this Regulations as a plant.
<figure><img src="image_8.png" /><figcaption>**Figure 6.** Biplot of the principal component analysis applied on all samples including the active variables (physical parameters, phenols and antioxidants, and parameters of color) and parameters of the experiment design as supplementary variables.</figcaption></figure>
<table><caption>**Table 5.** Ultra-Performance Liquid Chromatography–Tandem Mass Spectrometry (UPLC–MS/MS) analysis of extractive compounds from rosemary (measurements for CE and HVED treated samples) (ng/mL).</caption><thead><tr><th>Sample</th><th>Apigenin</th><th>Carnosol</th><th>Diosmetin</th><th>Hydroxytyrosol</th><th>Luteolin</th><th>Oleanolic Acid</th><th>Quercetin</th><th>Rosmarinic Acid</th><th>p-Cymene</th><th>Camphor</th><th>Thymol</th><th>Carvacrol</th><th>Extraction Type</th></tr></thead><tbody><tr><td>3 R0</td><td>44.460</td><td>0.940</td><td>115.897</td><td>0.394</td><td>180.406</td><td>/</td><td>/</td><td>13.030</td><td>0.009</td><td>0.602</td><td>0.002</td><td>0.013</td><td rowspan="6">CE</td></tr><tr><td>9 R0</td><td>32.818</td><td>0.869</td><td>111.415</td><td>0.104</td><td>152.254</td><td>/</td><td>/</td><td>0.408</td><td>0.059</td><td>0.038</td><td>0.002</td><td>0.001</td></tr><tr><td>3 R25</td><td>29.996</td><td>1.849</td><td>80.044</td><td>3.241</td><td>154.296</td><td>/</td><td>0.035</td><td>0.767</td><td>/</td><td>0.007</td><td>0.032</td><td>0.002</td></tr><tr><td>9 R25</td><td>27.342</td><td>2.410</td><td>82.583</td><td>2.983</td><td>147.022</td><td>/</td><td>/</td><td>0.761</td><td>0.070</td><td>0.003</td><td>0.002</td><td>0.002</td></tr><tr><td>3 R50</td><td>66.946</td><td>69.323</td><td>140.454</td><td>15.749</td><td>107.979</td><td>307.057</td><td>0.481</td><td>4756.226</td><td>0.033</td><td>0.949</td><td>0.031</td><td>0.029</td></tr><tr><td>9 R50</td><td>90.244</td><td>34.363</td><td>179.350</td><td>59.951</td><td>236.985</td><td>390.762</td><td>0.711</td><td>5100.455</td><td>0.729</td><td>0.251</td><td>0.020</td><td>0.300</td></tr><tr><td>RN1</td><td>159.160</td><td>207.346</td><td>310.578</td><td>68.141</td><td>305.866</td><td>288.807</td><td>11.271</td><td>5797.821</td><td>0.001</td><td>0.066</td><td>0.001</td><td>0.035</td><td rowspan="12">HVED</td></tr><tr><td>RN2</td><td>80.659</td><td>2.548</td><td>255.637</td><td>0.362</td><td>291.207</td><td>/</td><td>/</td><td>23.421</td><td>/</td><td>0.003</td><td>1.576</td><td>/</td></tr><tr><td>RN3</td><td>50.146</td><td>1.143</td><td>166.246</td><td>0.583</td><td>239.840</td><td>/</td><td>/</td><td>5.173</td><td>0.012</td><td>0.214</td><td>0.062</td><td>0.224</td></tr><tr><td>RN4</td><td>60.764</td><td>1.317</td><td>218.043</td><td>0.469</td><td>399.846</td><td>/</td><td>/</td><td>2.544</td><td>0.001</td><td>0.069</td><td>0.001</td><td>0.000</td></tr><tr><td>RN5</td><td>137.663</td><td>27.816</td><td>376.440</td><td>96.537</td><td>415.194</td><td>/</td><td>1.608</td><td>4228.058</td><td>0.010</td><td>0.238</td><td>/</td><td>/</td></tr><tr><td>RN6</td><td>107.933</td><td>8.988</td><td>314.863</td><td>61.995</td><td>326.021</td><td>/</td><td>0.269</td><td>3591.086</td><td>0.025</td><td>0.005</td><td>0.094</td><td>0.004</td></tr><tr><td>RN7</td><td>119.723</td><td>349.797</td><td>177.469</td><td>39.265</td><td>126.156</td><td>2091.128</td><td>1.510</td><td>5950.966</td><td>0.001</td><td>0.003</td><td>0.002</td><td>0.001</td></tr><tr><td>RN8</td><td>164.683</td><td>117.627</td><td>335.963</td><td>72.962</td><td>304.784</td><td>325.866</td><td>7.824</td><td>5745.552</td><td>0.002</td><td>0.004</td><td>0.003</td><td>0.005</td></tr><tr><td>RN9</td><td>123.606</td><td>303.095</td><td>191.720</td><td>37.747</td><td>122.569</td><td>2053.066</td><td>1.380</td><td>6002.350</td><td>/</td><td>0.006</td><td>0.205</td><td>0.000</td></tr><tr><td>RN10</td><td>95.125</td><td>3.532</td><td>362.800</td><td>3.845</td><td>600.262</td><td>4.907</td><td>/</td><td>29.489</td><td>0.001</td><td>0.153</td><td>0.001</td><td>0.000</td></tr><tr><td>RN11</td><td>116.828</td><td>195.651</td><td>194.323</td><td>38.968</td><td>124.223</td><td>1464.630</td><td>1.115</td><td>5700.140</td><td>0.001</td><td>0.015</td><td>0.001</td><td>0.000</td></tr><tr><td>RN12</td><td>50.827</td><td>6.240</td><td>125.214</td><td>11.977</td><td>167.060</td><td>14.068</td><td>/</td><td>68.983</td><td>0.000</td><td>0.007</td><td>0.000</td><td>0.000</td></tr><tr><td>RA1</td><td>112.850</td><td>286.709</td><td>206.963</td><td>37.852</td><td>167.721</td><td>920.212</td><td>3.707</td><td>5648.074</td><td>0.002</td><td>0.054</td><td>0.005</td><td>0.069</td></tr><tr><td>RA2</td><td>71.536</td><td>2.555</td><td>179.178</td><td>0.377</td><td>246.834</td><td>/</td><td>/</td><td>32.988</td><td>0.032</td><td>1.776</td><td>7.304</td><td>0.011</td></tr><tr><td>RA3</td><td>42.925</td><td>1.156</td><td>141.806</td><td>0.298</td><td>151.249</td><td>/</td><td>/</td><td>16.758</td><td>0.009</td><td>0.239</td><td>0.000</td><td>0.000</td></tr><tr><td>RA4</td><td>42.412</td><td>0.521</td><td>140.039</td><td>0.268</td><td>181.444</td><td>/</td><td>/</td><td>1.106</td><td>0.002</td><td>0.003</td><td>0.000</td><td>0.687</td></tr><tr><td>RA5</td><td>95.701</td><td>11.270</td><td>266.810</td><td>58.260</td><td>285.150</td><td>/</td><td>0.124</td><td>236.826</td><td>0.001</td><td>0.030</td><td>0.013</td><td>0.059</td></tr></tbody></table>
<table><caption>**Table 5.** *Cont.*</caption><thead><tr><th>Sample</th><th>Apigenin</th><th>Carnosol</th><th>Diosmetin</th><th>Hydroxytyrosol</th><th>Luteolin</th><th>Oleanolic Acid</th><th>Quercetin</th><th>Rosmarinic Acid</th><th>p-Cymene</th><th>Camphor</th><th>Thymol</th><th>Carvacrol</th><th>Extraction Type</th></tr></thead><tbody><tr><td>RA6</td><td>73.333</td><td>7.219</td><td>207.054</td><td>18.538</td><td>197.301</td><td>/</td><td>/</td><td>9.584</td><td>0.001</td><td>0.007</td><td>0.000</td><td>0.000</td><td rowspan="6">HVED</td></tr><tr><td>RA7</td><td>90.353</td><td>236.740</td><td>164.949</td><td>40.098</td><td>127.597</td><td>954.465</td><td>2.300</td><td>5829.363</td><td>0.001</td><td>/</td><td>0.000</td><td>0.000</td></tr><tr><td>RA8</td><td>111.501</td><td>291.279</td><td>193.902</td><td>39.888</td><td>159.881</td><td>1001.253</td><td>4.224</td><td>5872.906</td><td>0.018</td><td>0.019</td><td>0.000</td><td>0.000</td></tr><tr><td>RA9</td><td>56.649</td><td>10.003</td><td>173.879</td><td>11.916</td><td>203.449</td><td>/</td><td>/</td><td>30.652</td><td>0.019</td><td>0.137</td><td>0.000</td><td>0.000</td></tr><tr><td>RA10</td><td>34.136</td><td>1.326</td><td>111.703</td><td>0.276</td><td>239.407</td><td>/</td><td>/</td><td>5.206</td><td>0.000</td><td>0.006</td><td>0.000</td><td>0.001</td></tr><tr><td>RA11</td><td>76.580</td><td>157.254</td><td>131.357</td><td>26.408</td><td>98.402</td><td>757.572</td><td>1.479</td><td>5531.217</td><td>0.002</td><td>/</td><td>0.002</td><td>0.001</td></tr><tr><td>RA12</td><td>38.301</td><td>7.088</td><td>129.120</td><td>5.393</td><td>148.453</td><td>/</td><td>0.022</td><td>21.810</td><td>0.000</td><td>0.015</td><td>0.001</td><td>/</td></tr></tbody></table>

/—not detected.

<table><caption>**Table 6.** Headspace solid-phase microextraction/gas chromatography-mass spectrometry (HS-SPME/GC-MS) analysis of volatile compounds from rosemary (measurements for CE and HVED treated samples) (%).</caption><thead><tr><th rowspan="2">Sample</th><th colspan="4">Area (%)</th><th rowspan="2">Extraction Type</th></tr><tr><th>Eucalyptol (RI = 1038)</th><th>Camphor (RI = 1150)</th><th>Borneol (RI = 1172)</th><th>Linalool (RI = 1103)</th></tr></thead><tbody><tr><td>3 R0</td><td>40.33</td><td>26.70</td><td>13.46</td><td>/</td><td rowspan="5">CE</td></tr><tr><td>9 R0</td><td>34.89</td><td>24.81</td><td>15.78</td><td>/</td></tr><tr><td>3 R25</td><td>/</td><td>/</td><td>/</td><td>/</td></tr><tr><td>9 R25</td><td>/</td><td>/</td><td>/</td><td>/</td></tr><tr><td>3 R50</td><td>/</td><td>/</td><td>/</td><td>/</td></tr><tr><td>9 R50</td><td>/</td><td>/</td><td>/</td><td>/</td><td rowspan="5">HVED</td></tr><tr><td>RN1</td><td>/</td><td>/</td><td>/</td><td>/</td></tr><tr><td>RN2</td><td>31.04</td><td>22.8</td><td>14.73</td><td>3.46</td></tr><tr><td>RN3</td><td>39.44</td><td>24.63</td><td>15.58</td><td>2.96</td></tr><tr><td>RN4</td><td>32.56</td><td>22.19</td><td>13.42</td><td>2.79</td></tr><tr><td>RN5</td><td>2.52</td><td>1.12</td><td>0.24</td><td>/</td><td></td></tr></tbody></table>
<table><caption>Table 6. *Cont.*</caption><thead><tr><th rowspan="2">Sample</th><th colspan="4">Area (%)</th><th rowspan="2">Extraction Type</th></tr><tr><th>Eucalyptol (RI = 1038)</th><th>Camphor (RI = 1150)</th><th>Borneol (RI = 1172)</th><th>Linalool (RI = 1103)</th></tr></thead><tbody><tr><td>RN6</td><td>2.28</td><td>0.99</td><td>0.40</td><td>/</td><td rowspan="12">HVED</td></tr><tr><td>RN7</td><td>/</td><td>/</td><td>/</td><td>/</td></tr><tr><td>RN8</td><td>/</td><td>/</td><td>/</td><td>/</td></tr><tr><td>RN9</td><td>/</td><td>/</td><td>/</td><td>/</td></tr><tr><td>RN10</td><td>25.66</td><td>20.88</td><td>15.99</td><td>/</td></tr><tr><td>RN11</td><td>/</td><td>/</td><td>/</td><td>/</td></tr><tr><td>RN12</td><td>/</td><td>/</td><td>/</td><td>/</td></tr><tr><td>RA1</td><td>/</td><td>/</td><td>/</td><td>/</td></tr><tr><td>RA2</td><td>28.36</td><td>25.01</td><td>15.91</td><td>1.56</td></tr><tr><td>RA3</td><td>34.58</td><td>24.14</td><td>11.60</td><td>/</td></tr><tr><td>RA4</td><td>30.27</td><td>24.64</td><td>9.51</td><td>/</td></tr><tr><td>RA5</td><td>3.92</td><td>2.07</td><td>0.46</td><td>/</td></tr><tr><td>RA6</td><td>3.20</td><td>1.68</td><td>0.61</td><td>/</td></tr><tr><td>RA7</td><td>/</td><td>/</td><td>/</td><td>/</td></tr><tr><td>RA8</td><td>/</td><td>/</td><td>/</td><td>/</td></tr><tr><td>RA9</td><td>/</td><td>/</td><td>/</td><td>/</td></tr><tr><td>RA10</td><td>30.65</td><td>37.83</td><td>6.22</td><td>1.19</td></tr><tr><td>RA11</td><td>/</td><td>/</td><td>/</td><td>/</td></tr><tr><td>RA12</td><td>/</td><td>/</td><td>/</td><td>/</td></tr></tbody></table>

/—not detected.
<table><caption>**Table 7.** Residue levels and maximum residue levels (MRL) of pesticides (mg/kg) and metals (mg/kg) in rosemary samples.</caption><thead><tr><th rowspan="2">Component</th><th rowspan="2">MRL (mg/kg)</th><th rowspan="2">Content (mg/kg)</th><th colspan="4">HVED Extracts</th></tr><tr><th>RA8</th><th>RN7</th><th>RN9</th><th>RN11</th></tr></thead><tbody><tr><td>Alachlor</td><td>0.02</td><td>&lt;0.005</td><td>/</td><td>/</td><td>/</td><td>/</td></tr><tr><td>Aldrin and Dieldrin (Aldrin and dieldrin combined expressed as dieldrin)</td><td>0.01</td><td>&lt;0.002</td><td>/</td><td>/</td><td>/</td><td>/</td></tr><tr><td>Captan (Sum of captan and THPI, expressed as captan)</td><td>0.06</td><td>&lt;0.020</td><td>/</td><td>/</td><td>/</td><td>/</td></tr><tr><td>DDT (sum of p,p'-DDT, o,p'-DDT, p-p'-DDE and p,p'-TDE (DDD) expressed as DDT)</td><td>0.05</td><td>&lt;0.004</td><td>/</td><td>/</td><td>/</td><td>/</td></tr><tr><td>Endosulfan (sum of alpha- and beta-isomers and endosulfan-sulphate expressed as endosulfan)</td><td>0.05</td><td>&lt;0.002</td><td>/</td><td>/</td><td>/</td><td>/</td></tr><tr><td>Endrin</td><td>0.01</td><td>&lt;0.004</td><td>/</td><td>/</td><td>/</td><td>/</td></tr><tr><td>Heptachlor (sum of heptachlor and heptachlor epoxide expressed as heptachlor)</td><td>0.01</td><td>&lt;0.002</td><td>/</td><td>/</td><td>/</td><td>/</td></tr><tr><td>Hexachlorobenzene</td><td>0.01</td><td>&lt;0.002</td><td>/</td><td>/</td><td>/</td><td>/</td></tr><tr><td>Hexachlorocyclohexane (HCH), alpha-isomer</td><td>0.01</td><td>&lt;0.002</td><td>/</td><td>/</td><td>/</td><td>/</td></tr><tr><td>Hexachlorocyclohexane (HCH), beta-isomer</td><td>0.01</td><td>&lt;0.002</td><td>/</td><td>/</td><td>/</td><td>/</td></tr><tr><td>Iprodione</td><td>0.02</td><td>&lt;0.010</td><td>/</td><td>/</td><td>/</td><td>/</td></tr><tr><td>Lindane (Gamma-isomer of hexachlorocyclohexane (HCH))</td><td>0.01</td><td>&lt;0.002</td><td>/</td><td>/</td><td>/</td><td>/</td></tr><tr><td>Methoxychlor</td><td>0.01</td><td>&lt;0.010</td><td>/</td><td>/</td><td>/</td><td>/</td></tr><tr><td>Tolylfluanid (Sum of tolylfluanid and dimethylaminosulfotoluidide expressed as tolylfluanid)</td><td>0.05</td><td>&lt;0.020</td><td>/</td><td>/</td><td>/</td><td>/</td></tr><tr><td>Vinclozolin</td><td>0.02</td><td>&lt;0.004</td><td>/</td><td>/</td><td>/</td><td>/</td></tr><tr><td>Lead (Pb)</td><td>3.00</td><td>&lt;0.050</td><td>/</td><td>/</td><td>/</td><td>/</td></tr><tr><td>Cadmium (Cd)</td><td>1.00</td><td>&lt;0.006</td><td>/</td><td>/</td><td>/</td><td>/</td></tr><tr><td>Mercury (Hg)</td><td>0.10</td><td>0.026</td><td>/</td><td>/</td><td>/</td><td>/</td></tr><tr><td>Chromium (Cr)</td><td>/</td><td>0.240</td><td>55.3</td><td>66.1</td><td>71.0</td><td>60.5</td></tr><tr><td>Nickel (Ni)</td><td>/</td><td>0.322</td><td>2.10</td><td>1.10</td><td>1.20</td><td>0.950</td></tr><tr><td>Manganese (Mn)</td><td>/</td><td>21.00</td><td>7.10</td><td>5.10</td><td>5.45</td><td>6.20</td></tr><tr><td>Iron (Fe)</td><td>/</td><td>163</td><td>23.6</td><td>17.6</td><td>17.0</td><td>19.8</td></tr><tr><td>Copper (Cu)</td><td>/</td><td>6.40</td><td>3.00</td><td>3.75</td><td>3.95</td><td>6.90</td></tr><tr><td>Zinc (Zn)</td><td>/</td><td>26.0</td><td>6.65</td><td>9.10</td><td>10.7</td><td>20.5</td></tr></tbody></table>

/—no data available.
# 3. Discussion

The aim of the present study was to provide the potential of green solvents to extract BACs from rosemary leaves comparing theoretical and experimental methods.

## 3.1. Computational Simulation Methods for Assessing Solubility of Rosemary Compounds

For extractions of BACs, *n*-hexane is one of the most used solvents due to its low polarity, optimal boiling point, easy removal from the product by evaporation and stability. On the other side, *n*-hexane is a solvent of petrochemical origin, which are nowadays strictly regulated by European Directives. For that reason, industries are forced to replace such solvents with more sustainable alternative solvents [39]. Therefore, various green solvents have been chosen for assessing the potential for extraction of BACs from rosemary (Tables 1 and 2). Every solvent showed different theoretical solubility for rosemary BACs, and that can be explained by the differences in the solvent polarities. Generally, the optimum HSPs for very good solubility of solutes in solvents are presented with green color (Table 1). It can be concluded that many alternative (green) solvents are capable for extraction of BACs from rosemary, some even with higher affinity for extraction, compared to conventionally used *n*-hexane. According to RED results, by evaluating compounds with most green color (very good solubility), followed by yellow (medium solubility) and red color (poor solubility), the potential for extractions of BACs from rosemary was in the following order: CPME > limonene > cymene > MeTHF > ethyl oleate > β-myrcene > α-pinene > ethyl acetate > *n*-hexane > methylacetate > dimethyl carbonate > isopropanol > 1-butanol > methanol > ethanol > water. Results showed that α-pinene, cymene, β-myrcene, CPME, and ethyl oleate had high potential for extraction of volatile compounds such as monoterpenes and sesquiterpenes. For that reason, these solvents should be used for extraction of volatile compounds of essential oil from rosemary. Moreover, water and ethanol showed low potential for extraction of most evaluated compounds from rosemary.

Results of COSMO-RS solubility assessment (Table 2) showed similar trend for solvents like HSPs, although some differences have been noticed. COSMO-RS results showed following order according to most results with high probability of solubility (green color): MeTHF > CPME > ethyl acetate > methylacetate > 1-butanol > isopropanol > ethanol > ethyl-oleate > dimethylcarbonate > methanol > limonene > cymene > β-myrcene > α-pinene > *n*-hexane > water. It is clear that all green solvents have higher potential for extraction of rosemary compounds compared to conventional *n*-hexane, except water. According to results, ethyl acetate, methylacetate, ethanol, 1-butanol, isopropanol, methanol, CPME, dimethylcarbonate and MeTHF, showed high probability of solubility for diterpenes, triterpenes and flavonoids, which was not showed with HSPs assessment.

## 3.2. Experimental Analysis of Extraction of BACs and Aromas from Rosemary Using HVED

### 3.2.1. Physical Parameters of Rosemary Extracts

Results of pH, conductivity, power and temperature during extraction of rosemary BACs using HVED and CE are given in Figure 2. Temperature was measured before and after the HVED treatment and the maximum elevation of 5.9 °C was noted with sample RA5. In average, extracts treated for 9 min had 1 °C higher final temperature compared to extracts treated for 3 min. Also, maximum temperature was 29.5 °C after HVED treatment and it can be concluded that no significant elevation in temperature was noted during HVED treatment. Since all temperatures were under 30 °C, HVED was confirmed as a non-thermal extraction method. Results of pH for all extracts varied between 5.57 and 6.54 and no significant changes in pH were observed during the HVED treatment. pH and conductivity significantly depended (p ≤ 0.05) only on ethanol content. With higher content of ethanol in the solution, pH increased and conductivity decreased for both extraction types that is expected according to literature data [40].

The power used for the treatment with HVED changed from 7.0 to 27.0 kW. The highest power was noted for sample RN10 (27.0 ± 3.0 kW) which is expected since it was the sample treated for longer
time 9 min) and the highest voltage (25 kV). Accordingly, power significantly increased with higher voltage used for the treatment and decreased with higher ethanol content in the sample.

### 3.2.2. Total Phenolic Content and Antioxidant Activity of Rosemary Extracts

The goal of the extraction with HVED was to extract highest yield of phenolic compounds and antioxidants and compare results with same conditions by modified CE (at room temperature) (Figure 3). HVED treated samples showed 0.76–3.39-times higher yield of phenolic content for same extraction parameters than CE. On average, HVED had 2.13-times better yield of extraction of phenols than CE. The highest TPC value was noted for sample RN8 that was treated with HVED for 9 min at 25 kV with 50% of ethanol as a solvent. Nitrogen, longer treatment time and higher voltage yielded higher results of phenolic compounds, but only the treatment time had a statistically significant influence to TPC score ($p \le 0.05$). Ethanol content showed different trend for each treatment, with CE, the highest scores were obtained with water, HVED treatment with nitrogen was highest with 50% of ethanol, while treatment with argon had highest yield of TPC with 25% of ethanol. Bellumori et al., (2016) performed extractions from dried rosemary leaves by ultrasound-assisted extraction and microwave-assisted extraction for 10 min. Their results also showed higher amounts of TPC with ethanol as a solvent, compared to water, but their maximum obtained results were slightly higher compared to HVED, 35.0 mg/g for ultrasound and 36.6 mg/g DL for microwave-assisted extraction with ethanol (ACS grade, ≥99%), but different method for calculation was used [23].

Results of DPPH did not vary notably between extracts. However, the significant correlation was observed for DPPH with treatment time and ethanol content, it was higher with shorter treatment time and higher percentage of ethanol in the solution. Some differences in DPPH results could be due to generation of free radicals that is characteristic for HVED treatment that could influence to the antioxidant activity of extracts and have a possibility to interact with DPPH radical. FRAP results showed similarities with TPC values for HVED extraction, it was higher with longer treatment time, higher voltage used and with treatment using nitrogen, compared to argon. In average, HVED extraction showed 2.39-times higher antioxidant capacity, according to FRAP results, when compared with CE.

In total, higher results for phenolic content and antioxidants was noted with ethanol than with water. These results are in line with theoretical predictions with HSPs and COSMO-RS, although pure ethanol was used for calculations, while in experimental procedure 25 and 50% ethanol was used. Additionally, other green solvents were used for experimental assessment for extraction of BACs using HVED including limonene, α-pinene, glycerol, ethyl acetate and dimethyl carbonate. However, no electrical discharge was achieved during the extraction with HVED when mentioned solvents were used. This could be explained by high viscosity and density of these solvents [41,42]. Also, terpenes, such as limonene and pinene are oil solvents that are not able to mix with water and are therefore suitable for extraction of volatile compounds and not water-soluble compounds [41].

### 3.2.3. Near Infrared Spectroscopy of Rosemary Extracts and Modeling

The recorded NIR spectra showed grouping in the range of $\lambda = 904–1699$ nm which is the result of different vibrations of molecules as C–H, O–H, and N–H bonds. NIR spectra was specific in two regions: 904–925 nm and 1350–1699 nm. Overlapping of NIR scans is visible with specific differences in the range of 904–925 nm, indicating absorption detecting differences of the third overtone region and detecting different vibration of C–H bonds as well as HOH region at 1400 nm (water $\lambda = 1400–1460$ nm) and continues to the end of the recorded spectra indicating differences in the vibrations of the first overtone of C–H3, Ar–CH, C–H, and C–H2 bonds and second overtone of O–H, N–H, Ar–CH, and R–OH [43,44].

From those findings it was clear that the content of TPC and the antioxidative activity could be predicted. A part of the input spectral data was used for the model training (60%). After the training followed the PLSR model evaluation and its testing on unknown samples. Our aim was to predict those parameters not only qualitative but also quantitatively. In order to gain more accurate models,
four different wavelength ranges were used (Model 1: $\lambda$ = 904–1699 nm; Model 2: $\lambda$ = 1349–1699 nm; Model 3: $\lambda$ = 904–932 and 1349–1699 nm; and Model 4: $\lambda$ = 904–932 nm). The first model (Table S1) which included the total range of NIR spectra resulted with the best model efficiency parameters (RPDs > 3 and RER > 10) what is an indication of very good quantitative model prediction [45]. The efficacy of model for prediction of TPC, DPPH, and FRAP using NIR is presented in Figure 5.

### 3.2.4. Colorimetric Analysis of Rosemary Extracts

The color of the product is an important aspect that specifies the commercial quality of the product and has an effect to the consumer's final purchase decision [46]. The color measurements of the different extracts prepared by conventional and HVED extraction are given in terms of L*, a*, b*, C, and h values under CIE—L*a*b* color system. According to the results, HVED treated samples had lower values of L* (darker), a* (more green), and higher values of parameter b* (more yellow) when compared to CE in same conditions (Table 4). The HVED extraction also caused increasing in parameter C resulting with increase and discernible difference in color intensity and slight increase in hue of extracts. Differences with untreated extracts are expressed as ΔE—total color difference, ΔC—difference in tone color; ΔH—saturation. Based on these data with more intense coloring of HVED extracts (mass and the solvent used were equally as with CE), it can be concluded that HVED treatment had damaged the cell structure and that the pigments contained within the plant cells exited the cell surface causing changes in extract color [47]. Since the cavitation is caused, the process of extraction of BACs is also facilitated [48]. In our study, ethanol content had significant effect on color parameters (L*, a* and b*) in rosemary HVED treated samples with both argon and nitrogen. All colorimetric parameters, except ΔH, had statistically significant dependence (p ≤ 0.05) on ethanol content—with higher ethanol content, parameters L, ΔC, and ΔE increased, while a*, b*, C, and h decreased. Saturation significantly depended on treatment time and gas used, more saturation was noted with longer treatment time and use of nitrogen during the treatment.

In general, results indicated that HVED extracts had dark greenish-brown color with darker color when compared to CE which is associated also with higher phenolic compounds and antioxidant content. Plasma reactive species induce release of some BACs that are covalently bonded to the plant matrix which accordingly results with increase in TPC and greater antioxidant capacity and consequently with changes in extract colors. The reason is that phenolic compounds result with different color in free and bound forms [49,50].

### 3.2.5. Principal Component Analysis (PCA) of Rosemary Extracts

Based on the PCA biplot showed in Figure 6, the extraction type divided samples on the left and right part of chart quadrants. Extraction with CE positioned the samples in the second and third quadrant while the HVED extraction has spread the samples in all four quadrants with the main sample concentration in the fourth quadrant. The physical composition, content of TPC and antioxidant capacity as well the color parameters of the extracts. The experiment conditions are dominantly positioned near to the first principal component, PC1 (time (of HVED treatment, stirring and total treatment), voltage) with the exception—ethanol content.

The subjected table to the PCA chart-squared cosines of the variables asserts the HVED treatment time, stirring time, total treatment time and voltage used as experiment condition as values for which the squared cosine is the largest.

Biplot in the form where the active variables (physical properties, TPC, antioxidant capacity (DPPH and FRAP), and color parameters) are related to the supplementary which are the experimental design parameters, show the correlations between them and for which samples are they dominant. So accordingly, in the fourth quadrant the almost overlap of parameters as the DPPH content and the ethanol content indicating that the antioxidant capacity by DPPH method was higher in those samples where the ethanol content was 25% or 50% (Samples RA1; RA7; RA8; and all other positioned in the fourth quadrant, with the higher DPPH values of 31.28% (RA1 and RA7) and 31.31% (RA8),
respectively. Applying the same rule of data relation analysis in the third quadrant, samples RN12; RA11; and RA12 as well as 3 R25; 3 R50; 9 R25; and 9 R50 will have the highest values of the parameter L*, what confirmed the results presented in Table 4.

### 3.2.6. Analysis of Individual Bioactive Compounds from Rosemary Extracts

Data of UPLC-MS/MS analysis for individual phenolic compounds (Table 5) showed that main phenolic compounds in rosemary are: apigenin, diosmetin and rosmarinic acid. Comparison of phenolic compound in different type of extraction (CE and HVED) and conditions in extraction type showed that: apigenin, carnosol, diosmetin, hydroxytyrosol, luteolin, oleanolic acid, oleuropein, quercetin, and rosmarinic acid were higher in HVED extracts compared with CE extracts. The example of chromatograms for selected extracts is presented in Supplementary Materials (Figure S1). For this purpose, the extract with highest content of phenolic compounds detected with UPLC-MS/MS was chosen (RN9, Figure S1b) and compared with extract extracted with CE (3 R25, Figure S1a) under same conditions (3 min, 25 % of ethanol).

Statistical analysis showed that most of the measured compounds significantly depended only on ethanol content (apigenin, carnosol, diosmetin, hydroxytyrosol, luteolin, oleanolic acid, oleuropein, quercetin, and rosmarinic acid) and additionally, apigenin, carnosol, diosmetin, hydroxytyrosol, and luteolin depended on treatment time.

HVED is considered to be energy- and cost- saving method for successful extraction of phenolic compounds from rosemary. However, further analysis of energy and environmental impact should be performed. Hirondart et al., (2020) have obtained rosmarinic acid, carnosic acid, and carnosol by PLE in hydroalcoholic solution and conventional Soxhlet extraction. Extract yields of bioactive compounds were similar with both methods, energy consumption was lower for PLE extraction because less solvent had to be heated, and the cost was reduced with a smaller amount of waste generated [51].

When compared with theoretical results, it is clear that a similar trend for solution in water and ethanol was noticed. With higher ethanol content, solubility of most of extracted compounds increased, as well as a sum of all phenolic compounds. According to HSP results, all compounds that were analyzed with UPLC-MS/MS showed poor solubility (red color in Table 1) in both water and ethanol, but better results (higher solubility) was presented with ethanol as a solvent. However, COSMO-RS results gave better solubility results for extraction with ethanol. Results showed that carnosol, rosmarinic acid, apigenin and diosmetin have 100% of solubility in ethanol (green color in Table 2), while camphor has 44.06% solubility (yellow color). Experimental results confirmed these results since apigenin and carnosol showed the highest results with 50% of ethanol as a solvent, 164.68 and 349.80 ng/mL respectively, and diosmetin and rosmarinic acid showed highest measured results 25% of ethanol as an extraction solvent, 376.44 and 6002.35 ng/mL, respectively. Camphor was found in small amounts in all extracts (<1 ng/mL) except in water extract RA2 (1.78 ng/mL). Results are in line with one previously reported for oregano [52].

### 3.2.7. Analysis of Volatile Compounds from Rosemary Extracts

HS-SPME is a rapid, simple, inexpensive, solvent-free and highly sensitive technique [53]. Volatile organic compounds composition is strongly dependent on the extraction method. Results of the chemical composition of HS is presented in Table 6. The predominant HS compound was the cyclic monoterpene ether eucalyptol (40.33–2.28%). A second compound was cyclic monoterpene ketone camphor (26.70–0.99%), followed by bicyclic monoterpene borneol (15.99–0.24%). Terpene alcohol linalool (3.46–1.19%) was found in a smaller percentage. The concentration of these terpenes depends on the treatment of the plant with ethanol content or gas used (nitrogen or argon). From Table 6, it is notable that when water is used as an extraction solvent, the percentage of all three monoterpenes is high except for linalool. It was difficult to characterize volatile compounds in extracts with ethanol since overlapping profile of peaks happened with ethanol peak. Therefore, no results for most samples with ethanol were presented. The traceability is similar in the RN2-RN6 sample as in the RA2-RA6
sample, except that the concentration of linalool monoterpene alcohol in these samples deviated. Linalool was not found in the RN10 sample, while it was found in RA10. Chromatograms for two extracts with highest concentrations obtained by CE (3 R0, Figure S2a) and HVED (RA10, Figure S2b) are given in Supplementary Materials.

It was not possible to compare data with theoretical predictions for ethanol and water solubility of volatile compounds, since incomplete results are provided for ethanol extracts. However, the comparison between camphor and borneol can be provided to compare theoretical and experimental data. HSPs showed poor solubility in both ethanol and water for both compounds, but borneol had slightly higher chances for solubility (RED = 3.33 in ethanol and 9.3 in water), compared to camphor (RED = 4.49 for ethanol and 10.45 for water). From COSMO-RS results, it was also predicted that camphor and borneol have low chances for solubility in water, 0.1% and 0.02%, respectively, but different results were given for ethanol: camphor has medium probability of solubility (44.06%), while borneol has high probability (85.11%). Experimental results were opposite and higher concentrations of camphor were found in all extracts, compared to borneol. Moreover, results were closer to COSMO-RS results than HSPs since both compounds were extracted in higher concentrations: 37.83% of camphor was extracted in sample RA10 and 15.99% of borneol in sample RN10. Mena et al., (2016) showed similar results in conventionally performed acetone-based extraction from rosemary, they have extracted more camphor (41.52 ± 6.00 µg/g) than borneol (11.92 ± 2.01 µg/g) in their extracts [54].

SFE is a procedure recently used for the extraction of bioactive compounds and purification of rosemary essential oil. Mouahid et al., (2017) have compared the efficiency of SFE and hydrodistillation of rosemary leaves. The essential oil obtained by SFE had an increase in yield of monoterpenoids for 37%, while yields of individual monoterpenes varied [16]. Pereira et al., (2007) have performed a cost analysis of the extraction of rosemary essential oil by SFE and steam distillation. The manufacturing cost with SFE was lowered, while the lower profitability of steam distillation was a consequence of higher energy consumption and lower content of essential oil in the extracts [55]. SFE was used to extract the volatiles from rosemary and combined with SWE for the recovery of polyphenols from the produced extract. A combination of these processes resulted in a 28% reduction of operating cost when compared to the separate use of these techniques [56]. Since HVED was presented as a method that is more efficient for extraction of non-volatile than volatile compounds, it could be considered to be used in a combination with some other techniques for better extraction efficiency as well.

### 3.2.8. Analysis of Pesticides and Metals

Although rosemary could be considered as a nutritional supplement, there are still no categories in European Commission (EC) Regulation, herbs or plant tea, therefore the high levels of pesticides and heavy metals that can be found in its dried leaves could possess serious toxicological effects on human health. For that reason, the analysis of pesticides and heavy metals in the rosemary samples were measured and are presented in Table 7. This analysis is important for preparing healthy extracts from dried rosemary leaves that could be further used for new functional food. Trace analysis of pesticides residues were analyzed in dried leaves and for heavy metals residues were analyzed in dried leaves and HVED extracts. Residues levels of all pesticides were lower than limit of quantitation of method which is quite below maximum residue level (MRL) according to EC Regulations No. 396/2005 for rosemary as a plant. In EC regulation MRL of some pesticides with similar structures (such as DDT, endosulfan, aldrin, heptachlor, etc.) are grouped.

Residue levels of Lead (Pb) and Cadmium (Cd) in dried rosemary leaves were also below the limit of quantitation of the method, and only the level of Mercury (Hg) was slightly higher than limit of quantitation. However, all these values were quite below MRL according to EC No. 1881/2006. Furthermore, other metals (Chromium (Cr), Nickel (Ni), Manganese (Mn), Iron (Fe), Copper (Cu), and Zinc (Zn)) were measured in selected HVED extracts with high phenolic and antioxidant content (RA8, RN7, RN9, and RN11) as well and results are presented in Table 7. These metals are not included in EC Regulations so no MRL data were provided. Although the data should not be compared since
data before the HVED extraction are given per g of dried herb and data after HVED extraction are given per gram of extract, it is clear that content of Cr and Ni increased in extracts, while level of other metals decreased after HVED extraction. From this data it is notable that during the HVED treatment, levels of some toxic metals are increasing that could be the result of abrasion of electrodes during the treatment.

Rosemary plant was compliant regarding content of contaminants, pesticide residues and toxic metals. With respect to this, obtained results showed that rosemary samples are safe for use in human dietary. On the other hand, some changes in levels of metals could happen during the treatment and further detailed analyses should be done to assess this issue.

# 4. Materials and Methods

The concept of this work is presented in Figure 1, where all analysis performed for rosemary are presented as a flowchart.

## 4.1. Plant Materials

Dried rosemary leaves (*Rosmarinus officinalis* L.) were provided by local drugstore (Suban d.o.o., Samobor). Herbs were stored in polyethylene bags in a dark and dry place until extractions. Dried rosemary leaves were grinded to plant particle size distribution of d(0.1) ≤ 3 9.683 µm; d(0.5) ≤ 224.816 µm; d(0.9) ≤ 425.819 µm measured by the laser particle size analyzer Mastersizer 2000 (Malvern Instruments GmbH, Herrenberg, Germany). For the extraction, 1 g of herb material was weighted into the beaker of 100 mL and mixed with 50 mL of extracting solvent at room temperature (22 °C). Extraction was carried out using distilled water, 25% and 50% aqueous ethanol (v/v) as extraction solvents.

## 4.2. Computational Simulation Methods

### 4.2.1. Hansen Solubility Parameters (HSPs)

HSP provide a convenient and efficient way for characterization of solute-solvent interactions according to the classical “like dissolves like” rule. A detailed concept of HSPs is described in Aissou et al., (2017) [57]. For solvent optimization, a simple composite affinity parameter, the RED number, has been calculated to determine the solubility between solvents and solutes.

$$ \text{RED} = \frac{\text{Ra}}{\text{Ro}} \qquad (1) $$

where R₀ is the radius of a Hansen solubility sphere and Rₐ is the distance of a solvent from the center of the Hansen solubility sphere.

A potentially good solvent has RED number smaller than 1 (the compound has similar properties and will dissolve), while medium and poor solvents have RED values of from one to three and more than 3, respectively. The chemical structures of the solvents and solutes discussed in this article could be mutually transformed by JChemPaint version 3.3 (GitHub Pages, San Francisco, CA, USA) to their simplified molecular input line entry syntax (SMILES) notations, which were subsequently used to calculate the solubility parameters of the solvents and compounds (HSPiP Version 4.0, Hansen Solubility, Hørsholm, Denmark).

### 4.2.2. COSMO-RS Software

The COSMO-RS was developed by Klamt and co-workers as a statistical thermodynamic method for molecular description and solvent screening based on a quantum-chemical approach [58]. COSMO-RS prediction is a two-step procedure—microscopic and macroscopic. The procedure was explained in details by Aissou et al., (2017) [57]. The COSMOthermX program (version C30 release 13.01, COSMOlogic, Leverkusen, Germany) was used to calculate the relative solubility between the
solid compound and the liquid solvent in terms of the logarithm of the solubility in mole fractions ($\log_{10}(x_{solub})$). The logarithm of the best solubility was set to 0 and all other solvents were given relative to the best solvent. Also, the logarithm was transformed into probability of solubility (%). The calculation was performed at room temperature (20 °C) and at boiling temperature for each solvent.

### 4.3. High Voltage Electrical Discharge (HVED) and Conventional Extraction (CE)

HVED was performed with “IMP-SSPG-1200” generator (Impel group d.o.o., Zagreb, Croatia) that generated rectangular pulses using direct current and achieving high voltage. Maximum adjustable current was 30 mA and voltage up to 25 kV. Based on conducted preliminary experiments with different HVED parameters (frequency, voltage, pulse length, distance between electrodes, and ratio mass to solvent), fixed HVED parameters were chosen as follows: frequency of 100 Hz, pulse width 0.4 microseconds, voltage 15 and 20 kV for argon gas, and 20 and 25 kV for nitrogen gas, the gap between electrodes of 15 mm, treatment duration 3 and 9 min, and ratio mass to solvent 1 g:50 mL (according to pharmacopoeia). Mixture of herb material and solvent was transferred to beaker shaped reactor of 100 mL. This reactor that is opened on both sides was fitted with silicone tops with 1 cm in diameter. Silicone tops were used due to easier mounting of the electrode from the top and needle form the bottom. Gases (argon or nitrogen) were flowed in through the needle with the flow 0.5–1 L/min. Set-up of generator and reactor are shown in Figure 7. For measuring the output voltage (data not shown), oscilloscope Hantek DS05202BM (Tektronix, Inc., Beaverton, OR, USA) connected to the high voltage probe Tektronix P6015A (Hantek Electronic Co., Ltd., Qingdao, China) was used.

<figure><img src="image_9.png" /><figcaption>**Figure 7.** Set-up of generator and reactor for HVED treatments: (a) HVED and plasma generator, (b) Beaker shaped reactor: (1)—ground electrode; (2)—high voltage electrode (needle with empty interior for argon and nitrogen flow) during treatments; and (3)—discharge (plasma).</figcaption></figure>
For comparison, modified CE (untreated samples) was performed at room temperature as well, by dissolving the dried rosemary material in the solvent with light magnetic stirring during 3 or 9 min. Both extractions, HVED and CE, were performed in duplicates.

## 4.4. Analytical Methods

### 4.4.1. Determination of Total Phenolic Content (TPC)

For determination of TPC of rosemary extracts, a Folin–Ciocalteu method was used [59] with slight modifications. A volume of 0.1 mL of extract (appropriately diluted) was mixed with 0.2 mL of Folin–Ciocalteu reagent. After 3 min 1 mL of 20% Na₂CO₃ (m/v) was added. After thorough mixing by vortex, the reaction mixtures were incubated at 50 °C for 25 min, followed by absorbance reading at 765 nm against blank (instead of an extract, extraction solvent was used). The calibration curve was prepared using 50 to 500 mg/L of gallic acid in ethanol as a standard. The concentration of TPC was expressed in mg of gallic acid equivalents per gram of sample (mg GAE/g of sample).

### 4.4.2. 2-Diphenyl-2-Picrylhydrazyl (DPPH) Free Radical Assay

The antioxidant activity of rosemary extracts determined by DPPH method was determined as reported by Shortle et al., (2014) [59] with slight modifications. An aliquot (0.75 mL) of rosemary extracts or methanol solution of Trolox (25–200 mM) was mixed with 1.5 mL of 0.5 mM DPPH methanolic solution. After mixing, the solutions were stored in the dark for 20min at room temperature and the absorbance was measured at 517 nm against 100% methanol as a blank. The results were calculated using calibration curve for Trolox and expressed as µmol of Trolox equivalents per gram of sample (µmol TE/g of sample).

### 4.4.3. Ferric Reducing Antioxidant Power (FRAP) Assay

The FRAP assay was conducted according to literature [59] with modifications. The FRAP reagent was prepared by mixing 0.3 M acetate buffer (pH 3.6) with 10 mM TPTZ solution and 20 mM FeCl₃ solution in ratio 10:1:1. An aliquot (80 µL) of rosemary extract was mixed with 240 µL of water and 2080 µL of FRAP reagent. Following incubation at 37 °C for 5 min, the absorbance was measured at 595 nm. FRAP values were calculated according to the calibration curve for FeSO₄·7H₂O and expressed as µmol of Fe²⁺ equivalents (FE) per g of sample (µmol FE/g of sample).

### 4.4.4. Near Infrared Spectroscopy (NIR)

NIR spectroscopy was conducted using the NIR-128-1.7-USB/6.25/50µm (Control Development Inc., South Bend, IN, USA) to record sample spectra using the SPEC 32 Control Development software. NIR spectra was recorded in the wavelength range from 904 to 1699 nm. Each sample was recorded in triplicate and the average spectrum was calculated which was used for further processing.

### 4.4.5. Colorimetric Evaluation of Rosemary Extracts

Color parameters for all trials was measured by Konica Minolta colorimeter CM 3500d (Konica Minolta, Tokyo, Japan) at CIE Standard Illuminant D65 by 8 mm thick plate. All measurements were conducted in the Specular Component Included (SCI) mode as previously reported [60]. The color measurements of the different extracts prepared by CE and HVED extraction are given in terms of L*, a*, b*, C, and h values under CIE—L*a*b* color system (L*—lightness from black to white; a*—from green to red, and b*—from blue to yellow; C—chroma; and h—hue angle). Differences compared to CE were expressed as ΔE—total color difference, ΔC—difference in tone color; and ΔH—saturation.
### 4.4.6. Ultra-Performance Liquid Chromatography-Tandem Mass Spectrometry Characterization of Phenolic Compounds (UPLC-MS/MS)

UPLC-MS/MS Eksigent Expert Ultra LC 110, SCIEX 4500 QTRAP (SCIEX, USA) method for reference conditions [61] was conducted using Luna Omega 3 µm Polar C18 100 Å, 100 × 4.6 mm (column), thermostat column temperature 40 °C, automatic sampling temperature 4 °C, and injection volume of 10 µL. Mobile phases consisted of: A 100% H₂O with 0.1% HCOOH (v/v) and B 100% acetonitrile with 0.1 % HCOOH (v/v) with mobile phase flow 0.40 mL/min. Gradient was set as follows: 1 min 10% B, 2 min 10% B, 15 min 90% B, 25 min 90% B, 27 min 10% B, 30 min 10% B. Determination conditions for MS/MS detector were: ionization -negative ionization mode atmospheric pressure (API)—negative ionization at atmospheric pressure; ionization temperature: 500 °C, i.e., gas temperature combining the mobile phase at the exit from the capillary before ionization. Voltage on the electrode after capillary and next to ionization (Ion Spray Voltage) was -4500 V.

### 4.4.7. Headspace Solid-Phase Microextraction (HS-SPME) Followed by Gas Chromatography and Mass Spectrometry Analysis (GC-MS)

HS-SPME was performed with a manual SPME holder using three fiber covered with DVB/CAR/PDMS obtained from Supelco Co. (Bellefonte, PA, USA). For HS-SPME, the finely samples 2 mL were placed separately in 10 mL glass vials and hermetically sealed. The vials were maintained at 60 °C during equilibration (15 min) and extraction (45 min). Thereafter, the SPME fiber was withdrawn and inserted into GC-MS injector (250 °C) for 6 min for thermal desorption. The procedure was similar as in previous paper [62].

Gas chromatography and mass spectrometry (GC-MS) analyses were done on an Agilent Technologies (Palo Alto, CA, USA) gas chromatograph model 7890A equipped with a mass spectrometer (MSD) model 5977E (Palo Alto, CA, USA) and HP-5MS capillary column (5% phenyl-methylpolysiloxane, Agilent J & W). The GC conditions were the same as reported previously [62]. In brief, the oven temperature was set at 70 °C for 2 min, then increased from 70 to 200 °C (3 °C/min) and held at 200 °C for 18 min; the carrier gas was helium (1.0 mL/min). The compounds identification was based on the comparison of their retention indices (RI), determined relatively to the retention times of *n*-alkanes (C₉–C₂₅), with those reported in the literature [63] and those from Wiley 9 (Wiley, New York, NY, USA) and NIST 14 (National Institute of Standards and Technology; Gaithersburg, MD, USA) mass spectral database. The percentage composition of the samples was computed from the GC peak areas using the normalization method (without correction factors).

### 4.4.8. Determination of Pesticides and Metals in Rosemary Samples

The contents of the pesticides were performed by modified procedures with following national regulations HRN EN ISO 12393-1, 12393-2, and 12393-3: 2013, i.e., extraction with petroleum ether/dichloromethane and determination using the GC-ECD Varian CP-3800 instrument (Varian, Inc., Walnut Creek, CA, USA). Metal trace content was determined according to the HRN EN ISO 14084: 2005 procedure, or by wet sample digestion by HNO₃ (microwave digestion) with microwave reaction system Multiwave 3000 (Anton Paar GmbH, Graz, Austria). Determination of metals were conducted on the Perkin Elmer AAS Analyst 800 and ICP-MS Perkin Elmer NexION 300× (PerkinElmer, Inc., Waltham, MA, USA), while Hg traces were determined by the Leco AMA254 Hg analyzer (LECO Inc., St. Joseph, MI, USA).

## 4.5. *Experimental Design and Statistical Analysis*

The experiment was designed in STATGRAPHICS Centurion (StatPoint Technologies Inc., Warrenton, VA, USA) software. Multifactorial design consisting of 12 experimental trials using per gas (argon and nitrogen). The three chosen independent variables for HVED assisted extraction were: treatment time (3 and 9 min), voltage applied and gas type (15 or 20 kV for argon, and 20 or 25 kV for nitrogen) and concentration of ethanol (0%, 25%, or 50%). For CE, the independent
variables included: concentration of ethanol (0%, 25%, or 50%) and treatment time (3 and 9 min). The experimental design is presented in Table 3. A total of 30 extracts were prepared in duplicates.

In order to provide information about experimental results, a descriptive statistic was used. Analysis of covariance (ANCOVA) was used for assessment of correlation and all parameters according to dependent variables (gas, treatment time, voltage, ethanol content). The *p*-values present the statistical significance of each of the factor and it was significant at *p* ≤ 0.05. Statistics was performed using XLStat (MS Excel 2010) (data not shown).

Information investigation of NIR spectra includes preprocessing and calibration modeling [64] where preprocessing will minimize commotions and undesirable components in spectra, which are subjected to construction of calibration models. Data matrix which included NIR spectra and physical-chemical properties of the samples consisted of 128 rows and 797 columns. This matrix was used for the identification of qualitative differences, by use of principal component analysis (PCA) as well as in the modelling that followed. The NIR spectra were pre-treated to enhance the prediction accuracy. Several spectra pre-treatment methods were arranged to the original absorbance spectra such as multiplicative scatter correction (MSC), standard normal variate (SNV), Smoothing (Moving Average, Gaussian and Median filter, and Savitzky-Golay), first and second derivative absorbance (d1a and d2a), Savitzky-Golay first and second derivative absorbance (S-G d1a and S-G d2a), and the combination of the MSC and SNV + d1a or d2a, but as the most effective was the Savitzky-Golay first derivation (S-G d1a). Calibration model was developed by principal component regression (PCR). From the data matrix 2/4 of it was used for the calibration and 1/4 was used for the validation while the remaining 1/4 was used for the prediction. Validation was done by K-fold cross-validation. [65]. Then followed the application of the multivariate tool mostly used in model prediction, the PLSR and PCR [45]. As in the case of pretreatments, one or more multivariate tools can be used in the calibration, which implies quantitative or qualitative analysis. Model efficacy is evaluated using R² and the regression point displacement that is the ratio of the standard RPD and the ratio of the range of reference chemistry values to RER. RPD is the ratio of SDv and SEP, while RER range, min and max, and minimum and maximum values of the validation set. All data analyses were conducted in MS Excel and its additional statistical tool pack: XLStat.

## 5. Conclusions

In this study, the potential of high voltage discharges for green solvent extraction of BACs and aromas from rosemary leaves was assessed by computational simulation and experimental method by means of HVED. The experimental results were compared with untreated samples (modified CE) and HVED was presented to yield 2.13-times higher TPC and 2.39-times higher antioxidant capacity. Nitrogen, longer treatment time, and higher voltage yielded higher results of phenolic compounds and antioxidants. Also, NIR spectra and modelling with analytical data were shown as an extremely useful tools that can help in assessing whether there is a "cost-effectiveness" of extracting phenols or antioxidants from specific samples, in a quick and easy way. The results presented that NIR spectroscopy combined with chemometrics approach gave accurate TPC, FRAP, and DPPH content prediction, showing that indicates the potential of the method in estimating the quantitative expected antioxidant potential as well as the content of total phenols. Generally, HVED extracts had a dark greenish-brown color with darker color when compared to CE which is associated also with higher phenolic compounds and antioxidant content. An UPLC-MS/MS showed that main phenolic compounds in rosemary were apigenin, diosmetin, and rosmarinic acid, while the predominant volatile compounds in rosemary extracts was eucalyptol. Altogether, results showed that HVED confirmed high potential for extraction of BACs and food aromas from rosemary with increased yield of individual compounds and total phenolic and antioxidant properties, compared to untreated samples. Furthermore, rosemary was presented as safe raw material for further processing in human nutrition in terms of pesticides and metals. The computational stimulation methods were confirmed by experimental study, ethanol had higher potential of solubility of BACs and aromas from rosemary compared to water. Therefore,
these theoretical prediction methods present a new approach for assessment of solubility of individual compounds in selected solvents that could impact to lower solvent usage during experimentation and lower environmental impact.

**Supplementary Materials:** The following are available online. Table S1: Model statistics for prediction of compositional parameters of rosemary extracts based on the NIR spectra, Figure S1: UPLC-MS/MS chromatograms of representative extracts: (a) CE (sample 3 R25), and (b) HVED (sample RN9), Figure S2: UPLC-MS/MS chromatograms of representative extracts: (a) CE (sample 3 R0), and (b) HVED (sample RA10), and (c) chemical structure of main detected compounds, Figure S3: Box plots for the (A) content of total phenols; antioxidant activity of the samples conducted by the (B) DPPH and (C) FRAP method and the (D) yield for samples treated by CE (R) and HVED (RN & RA), Figure S4: PCA biplot for different extraction types (CE: R; HVED; RN & RA).

**Author Contributions:** Conceptualization, A.R.J. and M.N.; Methodology, M.N and A.R.J.; Software, M.V. and F.C.; Validation, A.R.J.; Formal Analysis. M.N., J.G.K., Z.M., I.D., M.V. and E.M.; Investigation, M.N. and E.M.; Resources, A.R.J.; Data Curation, M.N. and A.R.J.; Writing—Original Draft Preparation, M.N.; Writing—Review & Editing, J.G.K., Z.M., I.D. and A.R.J.; Visualization, M.N.; Supervision, A.R.J.; Project Administration, A.R.J.; Funding Acquisition, A.R.J. and M.N. All authors have read and agreed to the published version of the manuscript.

**Funding:** This work was supported by grant from the Croatian Science Foundation: "High voltage discharges for green solvent extraction of bioactive compounds from Mediterranean herbs (IP-2016-06-1913)".

**Acknowledgments:** The work of doctoral student Marinela Nutrizio has been fully supported by the "Young researchers' career development project—training of doctoral students" of the Croatian Science Foundation funded by the European Union from the European Social Fund. Authors would like to thank Renata Fujfer and Paula Žurga from Teaching Institute for Public health of Primorje-Gorski Kotar County for help with analyses of pesticides and metals.

**Conflicts of Interest:** All authors declare that they have no conflict of interest. This article does not contain any studies with human or animal subjects.

## References

1. Ali, A.; Chua, B.L.; Chow, Y.H. An insight into the extraction and fractionation technologies of the essential oils and bioactive compounds in *Rosmarinus officinalis* L.: Past, present and future. *TrAC-Trends Anal. Chem.* **2019**, *118*, 338–351. [CrossRef]

2. Babovic, N.; Djilas, S.; Jadranin, M.; Vajs, V.; Ivanovic, J.; Petrovic, S.; Zizovic, I. Supercritical carbon dioxide extraction of antioxidant fractions from selected Lamiaceae herbs and their antioxidant capacity. *Innov. Food Sci. Emerg. Technol.* **2010**, *11*, 98–107. [CrossRef]

3. Viuda-Martos, M.; Ruiz Navajas, Y.; Sánchez Zapata, E.; Fernández-López, J.; Pérez-Álvarez, J.A. Antioxidant activity of essential oils of five spice plants widely used in a Mediterranean diet. *Flavour Fragr. J.* **2010**, *25*, 13–19. [CrossRef]

4. Moore, J.; Yousef, M.; Tsiani, E. Anticancer effects of rosemary (*Rosmarinus officinalis* L.) extract and rosemary extract polyphenols. *Nutrients* **2016**, *8*, 731. [CrossRef]

5. Haloui, M.; Louedec, L.; Michel, J.B.; Lyoussi, B. Experimental diuretic effects of *Rosmarinus officinalis* and *Centaurium erythraea*. *J. Ethnopharmacol.* **2000**, *71*, 465–472. [CrossRef]

6. Rožman, T.; Jeršek, B. Antimicrobial activity of rosemary extracts (*Rosmarinus officinalis* L.) against different species of Listeria. *Acta Agric. Slov.* **2009**, *93*, 51–58. [CrossRef]

7. Sánchez-Camargo, A.P.; Mendiola, J.A.; Valdés, A.; Castro-Puyana, M.; García-Cañas, V.; Cifuentes, A.; Herrero, M.; Ibáñez, E. Supercritical antisolvent fractionation of rosemary extracts obtained by pressurized liquid extraction to enhance their antiproliferative activity. *J. Supercrit. Fluids* **2016**, *107*, 581–589. [CrossRef]

8. Takaki, I.; Bersani-Amado, L.E.; Vendruscolo, A.; Sartoretto, S.M.; Diniz, S.P.; Bersani-Amado, C.A.; Cuman, R.K.N. Anti-inflammatory and antinociceptive effects of *Rosmarinus officinalis* L. essential oil in experimental animal models. *J. Med. Food* **2008**, *11*, 741–746. [CrossRef]

9. Naimi, M.; Vlavcheski, F.; Shamshoum, H.; Tsiani, E. Rosemary extract as a potential anti-hyperglycemic agent: Current evidence and future perspectives. *Nutrients* **2017**, *9*, 968. [CrossRef]

10. FDA Code of Federal Regulations (CFR). Title 21: Food and Drugs. Chapter I–Food and Drug Administration, Department of Health and Human Services (Continued), Subchapter B–Food for Human Consumption, Part 182–Substances Generally Recognized as Safe (GRAS). Available online: https://www.accessdata.fda.gov/scripts/cdrh/cfdocs/cfcfr/CFRSearch.cfm?fr=182.10 (accessed on 9 June 2020).
11. del Pilar Sánchez-Camargo, A.; Herrero, M. Rosemary (*Rosmarinus officinalis*) as a functional ingredient: Recent scientific evidence. *Curr. Opin. Food Sci.* **2017**, *14*, 13–19. [CrossRef]

12. Azmir, J.; Zaidul, I.S.M.; Rahman, M.M.; Sharif, K.M.; Mohamed, A.; Sahena, F.; Jahurul, M.H.A.; Ghafoor, K.; Norulaini, N.A.N.; Omar, A.K.M. Techniques for extraction of bioactive compounds from plant materials: A review. *J. Food Eng.* **2013**, *117*, 426–436. [CrossRef]

13. Jacotet-Navarro, M.; Rombaut, N.; Fabiano-Tixier, A.S.; Danguien, M.; Bily, A.; Chemat, F. Ultrasound versus microwave as green processes for extraction of rosmarinic, carnosic and ursolic acids from rosemary. *Ultrason. Sonochem.* **2015**, *27*, 102–109. [CrossRef] [PubMed]

14. Švarc-Gajić, J.; Stojanović, Z.; Segura Carretero, A.; Arráez Román, D.; Borrás, I.; Vasiljević, I. Development of a microwave-assisted extraction for the analysis of phenolic compounds from Rosmarinus officinalis. *J. Food Eng.* **2013**, *119*, 525–532. [CrossRef]

15. Zermane, A.; Larkeche, O.; Meniai, A.H.; Crampon, C.; Badens, E. Optimization of Algerian rosemary essential oil extraction yield by supercritical CO2 using response surface methodology. In *Proceedings of the Comptes Rendus Chimie*; Elsevier: Amsterdam, The Netherlands, 2016; Volume 19, pp. 538–543.

16. Mouahid, A.; Dufour, C.; Badens, E. Supercritical CO2 extraction from endemic Corsican plants; Comparison of oil composition and extraction yield with hydrodistillation method. *J. CO2 Util.* **2017**, *20*, 263–273. [CrossRef]

17. Tena, M.T.; Valcárcel, M.; Hidalgo, P.J.; Ubera, J.L. Supercritical Fluid Extraction of Natural Antioxidants from Rosemary: Comparison with Liquid Solvent Sonication. *Anal. Chem.* **1997**, *69*, 521–526. [CrossRef] [PubMed]

18. Nieto, A.; Borrull, F.; Pocurull, E.; Marcé, R.M. Pressurized liquid extraction: A useful technique to extract pharmaceuticals and personal-care products from sewage sludge. *TrAC-Trends Anal. Chem.* **2010**, *29*, 752–764. [CrossRef]

19. Smith, R.M. Extractions with superheated water. *J. Chromatogr.* **A** **2002**, *975*, 31–46. [CrossRef]

20. Hossain, M.B.; Barry-Ryan, C.; Martin-Diana, A.B.; Brunton, N.P. Optimisation of accelerated solvent extraction of antioxidant compounds from rosemary (*Rosmarinus officinalis* L.), marjoram (*Origanum majorana* L.) and oregano (*Origanum vulgare* L.) using response surface methodology. *Food Chem.* **2011**, *126*, 339–346. [CrossRef]

21. Albu, S.; Joyce, E.; Paniwnyk, L.; Lorimer, J.P.; Mason, T.J. Potential for the use of ultrasound in the extraction of antioxidants from *Rosmarinus officinalis* for the food and pharmaceutical industry. In *Proceedings of the Ultrasonics Sonochemistry*; Elsevier: Amsterdam, The Netherlands, 2004; Volume 11, pp. 261–265.

22. Paniwnyk, L.; Cai, H.; Albu, S.; Mason, T.J.; Cole, R. The enhancement and scale up of the extraction of anti-oxidants from *Rosmarinus officinalis* using ultrasound. *Ultrason. Sonochem.* **2009**, *16*, 287–292. [CrossRef]

23. Bellumori, M.; Innocenti, M.; Binello, A.; Boffa, L.; Mulinacci, N.; Cravotto, G. Selective recovery of rosmarinic and carnosic acids from rosemary leaves under ultrasound-and microwave-assisted extraction procedures. *Comptes Rendus Chim.* **2016**, *19*, 699–706. [CrossRef]

24. Chemat, F.; Abert Vian, M.; Fabiano-Tixier, A.-S.; Nutrizio, M.; Režek Jambrak, A.; Munekata, P.E.S.; Lorenzo, J.M.; Barba, F.J.; Binello, A.; Cravotto, G. A review of sustainable and intensified techniques for extraction of food and natural products. *Green Chem.* **2020**, *22*, 2325–2353. [CrossRef]

25. Li, Z.; Fan, Y.; Xi, J. Recent advances in high voltage electric discharge extraction of bioactive ingredients from plant materials. *Food Chem.* **2019**, *277*, 246–260. [CrossRef] [PubMed]

26. Almohammed, F.; Koubaa, M.; Khelfa, A.; Nakaya, M.; Mhemdi, H.; Vorobiev, E. Pectin recovery from sugar beet pulp enhanced by high-voltage electrical discharges. *Food Bioprod. Process.* **2017**, *103*, 95–103. [CrossRef]

27. Boussetta, N.; Vorobiev, E.; Reess, T.; De Ferron, A.; Pecastaing, L.; Ruscassié, R.; Lanoisellé, J.L. Scale-up of high voltage electrical discharges for polyphenols extraction from grape pomace: Effect of the dynamic shock waves. *Innov. Food Sci. Emerg. Technol.* **2012**, *16*, 129–136. [CrossRef]

28. Gu, Y.; Jérôme, F. Bio-based solvents: An emerging generation of fluids for the design of eco-efficient processes in catalysis and organic chemistry. *Chem. Soc. Rev.* **2013**, *42*, 9550–9570. [CrossRef]

29. Chemat, F.; Vian, M.A.; Cravotto, G. Green extraction of natural products: Concept and principles. *Int. J. Mol. Sci.* **2012**, *13*, 8615–8627. [CrossRef]

30. Barbieri, J.B.; Goltz, C.; Batistão Cavalheiro, F.; Theodoro Toci, A.; Igarashi-Mafra, L.; Mafra, M.R. Deep eutectic solvents applied in the extraction and stabilization of rosemary (*Rosmarinus officinalis* L.) phenolic compounds. *Ind. Crops Prod.* **2020**, *144*, 112049. [CrossRef]
31. Sicaire, A.-G.; Vian, M.A.; Fine, F.; Carré, P.; Tostain, S.; Chemat, F. Experimental approach versus COSMO-RS assisted solvent screening for predicting the solubility of rapeseed oil. *OCL* **2015**, *22*, D404. [CrossRef]

32. Belščak-Cvitanović, A.; Valinger, D.; Benković, M.; Tušek, A.J.; Jurina, T.; Komes, D.; Gajdoš Kljusurić, J. Integrated approach for bioactive quality evaluation of medicinal plant extracts using HPLC-DAD, spectrophotometric, near infrared spectroscopy and chemometric techniques. *Int. J. Food Prop.* **2017**, *20*, S2463-S2480. [CrossRef]

33. Valinger, D.; Kušen, M.; Jurinjak Tušek, A.; Panić, M.; Jurina, T.; Benković, M.; Radojčić Redovniković, I.; Gajdoš Kljusurić, J. Development of near infrared spectroscopy models for quantitative prediction of the content of bioactive compounds in olive leaves. *Chem. Biochem. Eng. Q.* **2018**, *32*, 535-543. [CrossRef]

34. Kljusurić, J.G.; Jurina, T.; Valinger, D.; Benković, M.; Tušek, A.J. NIR spectroscopy and management of bioactive components, antioxidant activity, and macronutrients in fruits. In *Fruit Crops*; Elsevier: Amsterdam, The Netherlands, 2020; pp. 95-109.

35. Wang, W.; Paliwal, J. Near-infrared spectroscopy and imaging in food quality and safety. *Sens. Instrum. Food Qual. Saf.* **2007**, *1*, 193-207. [CrossRef]

36. Oliveira, L.S.; Franca, A.S. Applications of near infrared spectroscopy (NIRS) in food quality evaluation. In *Food Quality: Control, Analysis and Consumer Concerns*; Medina, D.A., Laine, A.M., Eds.; Nova Science Publishers: Hauppauge, NY, USA, 2011.

37. Dyrby, M.; Engelsen, S.B.; Nørgaard, L.; Bruhn, M.; Lundsberg-Nielsen, L. Chemometric Quantitation of the Active Substance (Containing C≡N) in a Pharmaceutical Tablet Using Near-Infrared (NIR) Transmittance and NIR FT-Raman Spectra. *Appl. Spectrosc.* **2002**, *56*, 579-585. [CrossRef]

38. Patel, H. Near Infrared Spectroscopy: Basic principles and use in tablet evaluation. *Int. J. Chem. Life Sci.* **2017**, *6*, 2006-2015. [CrossRef]

39. Yara-Varón, E.; Fabiano-Tixier, A.S.; Balcells, M.; Canela-Garayoa, R.; Bily, A.; Chemat, F. Is it possible to substitute hexane with green solvents for extraction of carotenoids? A theoretical versus experimental solubility study. *RSC Adv.* **2016**, *6*, 27750-27759. [CrossRef]

40. Maeda, K.; Maeno, K.; Fukui, K.; Moritoki, M.; Kuramochi, H. Electrical conductivity of aqueous ethanol solutions containing ammonium salts under high pressure at 298 K. *J. Chem. Eng. Data* **2013**, *58*, 264-270. [CrossRef]

41. Tanzi, C.D.; Vian, M.A.; Ginies, C.; Elmaataoui, M.; Chemat, F. Terpenes as green solvents for extraction of oil from microalgae. *Molecules* **2012**, *17*, 8196-8205. [CrossRef]

42. Wolfson, A.; Dlugy, C.; Shotland, Y. Glycerol as a green solvent for high product yields and selectivities. *Environ. Chem. Lett.* **2007**, *5*, 67-71. [CrossRef]

43. Neveu, V.; Perez-Jiménez, J.; Vos, F.; Crespy, V.; du Chaffaut, L.; Mennen, L.; Knox, C.; Eisner, R.; Cruz, J.; Wishart, D.; et al. Phenol-Explorer: An Online Comprehensive Database on Polyphenol Contents in Foods-PubMed. *Database* **2010**, *2010*, bap024. [CrossRef]

44. Rothwell, J.A.; Perez-Jimenez, J.; Neveu, V.; Medina-Remón, A.; M'hiri, N.; García-Lobato, P.; Manach, C.; Knox, C.; Eisner, R.; Wishart, D.S.; et al. Phenol-Explorer 3.0: A major update of the Phenol-Explorer database to incorporate data on the effects of food processing on polyphenol content | Database | Oxford Academic. *Database* **2013**, bat070. [CrossRef]

45. Kujundžić, D.; Režek Jambrak, A.; Vukušić, T.; Stulić, S.; Gajdoš Kljusurić, J.; Banović, M.; Herceg, Z. Near-infrared spectroscopic characterization of steviol glycosides extracted from Stevia rebaudiana Bertoni using high-power ultrasound and gas-phase plasma. *J. Food Nutr. Res.* **2017**, *56*, 109-120.

46. Hertwig, C.; Reineke, K.; Ehlbeck, J.; Erdoğdu, B.; Rauh, C.; Schlüter, O. Impact of remote plasma treatment on natural microbial load and quality parameters of selected herbs and spices. *J. Food Eng.* **2015**, *167*, 12-17. [CrossRef]

47. Mahnič-Kalamiza, S.; Vorobiev, E.; Miklavčič, D. Electroporation in Food Processing and Biorefinery. *J. Membr. Biol.* **2014**, *247*, 1279-1304. [CrossRef]

48. Fincan, M.; DeVito, F.; Dejmek, P. Pulsed electric field treatment for solid-liquid extraction of red beetroot pigment. *J. Food Eng.* **2004**, *64*, 381-388. [CrossRef]

49. Rodríguez, Ó.; Gomes, W.F.; Rodrigues, S.; Fernandes, F.A.N. Effect of indirect cold plasma treatment on cashew apple juice (Anacardium occidentale L.). *LWT* **2017**, *84*, 457-463. [CrossRef]

50. Massini, L.; Rico, D.; Martin-Diana, A.B. Quality Attributes of Apple Juice: Role and Effect of Phenolic Compounds. *Fruit Juices* **2018**, *45-57*.
51. Hirondart, M.; Rombaut, N.; Fabiano-Tixier, A.S.; Bily, A.; Chemat, F. Comparison between pressurized liquid extraction and conventional Soxhlet extraction for rosemary antioxidants, yield, composition, and environmental footprint. *Foods* **2020**, *9*, 584. [CrossRef]

52. Nutrizio, M.; Maltar-Strmečki, N.; Chemat, F.; Duić, B.; Jambrak, A.R. High-Voltage Electrical Discharges in Green Extractions of Bioactives from Oregano Leaves (*Origanum vulgare* L.) Using Water and Ethanol as Green Solvents Assessed by Theoretical and Experimental Procedures. *Food Eng. Rev.* **2020**, *1–14*. [CrossRef]

53. Carrillo, J.D.; Tena, M.T. Determination of volatile compounds in antioxidant rosemary extracts by multiple headspace solid-phase microextraction and gas chromatography. *Flavour Fragr. J.* **2006**, *21*, 626–633. [CrossRef]

54. Mena, P.; Cirlini, M.; Tassotti, M.; Herrlinger, K.; Dall'Asta, C.; Del Rio, D. Phytochemical Profiling of Flavonoids, Phenolic Acids, Terpenoids, and Volatile Fraction of a Rosemary (*Rosmarinus officinalis* L.) Extract. *Molecules* **2016**, *21*, 1576. [CrossRef]

55. Pereira, C.G.; Meireles, M.A.A. Economic analysis of rosemary, fennel and anise essential oils obtained by supercritical fluid extraction. *Flavour Fragr. J.* **2007**, *22*, 407–413. [CrossRef]

56. Zabot, G.L.; Moraes, M.N.; Carvalho, P.I.N.; Meireles, M.A.A. New proposal for extracting rosemary compounds: Process intensification and economic evaluation. *Ind. Crops Prod.* **2015**, *77*, 758–771. [CrossRef]

57. Aissou, M.; Chemat-Djenni, Z.; Yara-Varón, E.; Fabiano-Tixier, A.-S.; Chemat, F. Limonene as an agro-chemical building block for the synthesis and extraction of bioactive compounds. *Comptes Rendus Chim.* **2017**, *20*, 346–358. [CrossRef]

58. Klamt, A. Prediction of the mutual solubilities of hydrocarbons and water with COSMO-RS. *Fluid Phase Equilib.* **2003**, *206*, 223–235. [CrossRef]

59. Shortle, E.; O'Grady, M.N.; Gilroy, D.; Furey, A.; Quinn, N.; Kerry, J.P. Influence of extraction technique on the anti-oxidative potential of hawthorn (Crataegus monogyna) extracts in bovine muscle homogenates. *Meat Sci.* **2014**, *98*, 828–834. [CrossRef] [PubMed]

60. Žuntar, I.; Putnik, P.; Bursać Kovačević, D.; Nutrizio, M.; Šupljika, F.; Poljanec, A.; Dubrović, I.; Barba, F.J.; Režek Jambrak, A. Phenolic and Antioxidant Analysis of Olive Leaves Extracts (*Olea europaea* L.) Obtained by High Voltage Electrical Discharges (HVED). *Foods* **2019**, *8*, 248.

61. Arranz, E.; Santoyo, S.; Jaime, L.; Fornari, T.; Reglero, G.; Guri, A.; Corredig, M. Improved Bioavailability of Supercritical Rosemary Extract Through Encapsulation in Different Delivery Systems After In Vitro Digestion. *Food Dig.* **2015**, *6*, 30–37.

62. Jerković, I.; Marijanović, Z.; Kuš, P.M.; Tuberoso, C.I.G. Comprehensive Study of Mediterranean (Croatian) Propolis Peculiarity: Headspace, Volatiles, Anti-Varroa-Treatment Residue, Phenolics, and Antioxidant Properties. *Chem. Biodivers.* **2016**, *13*, 210–218. [CrossRef]

63. El-Sayed, A.M. The Pherobase: Database of Insect Pheromones and Semiochemicals. 2014. Available online: http://www.pherobase.com/ (accessed on 24 June 2020).

64. Cortés, V.; Blasco, J.; Aleixos, N.; Cubero, S.; Talens, P. Monitoring strategies for quality control of agricultural products using visible and near-infrared spectroscopy: A review. *Trends Food Sci. Technol.* **2019**, *85*, 138–148. [CrossRef]

65. Wang, H.; Peng, J.; Xie, C.; Bao, Y.; He, Y. Fruit quality evaluation using spectroscopy technology: A review. *Sensors* **2015**, *15*, 11889–11927. [CrossRef]

**Sample Availability:** Samples of the compounds for UPLC-MS/MS and GC-MS are available from the authors.

<figure><img src="image_10.png" /></figure>

© 2020 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).