Next Article in Journal
Effects of Reducing Sugars on the Structural and Flavor Properties of the Maillard Reaction Products of Lycium barbarum Seed Meal
Next Article in Special Issue
Development of Pectin Particles as a Colon-Targeted Marjoram Phenolic Compound Delivery System
Previous Article in Journal
Nutritional and Technological Properties of Albino Peach Palm (Bactris gasipaes) from the Amazon: Influence of Cooking and Drying
Previous Article in Special Issue
Characterization, Thermal Stability and Antimicrobial Evaluation of the Inclusion Complex of Litsea cubeba Essential Oil in Large-Ring Cyclodextrins (CD9–CD22)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Microencapsulation of Juniper and Black Pepper Essential Oil Using the Coacervation Method and Its Properties after Freeze-Drying

by
Alicja Napiórkowska
1,
Arkadiusz Szpicer
1,
Iwona Wojtasik-Kalinowska
1,
Maria Dolores Torres Perez
2,
Herminia Dominguez González
2 and
Marcin Andrzej Kurek
1,*
1
Department of Technique and Food Development, Warsaw University of Life Sciences, 02-787 Warsaw, Poland
2
CINBIO, Departamento de Ingeniería Química, Campus Ourense, Universidade de Vigo, 32004 Ourense, Spain
*
Author to whom correspondence should be addressed.
Foods 2023, 12(23), 4345; https://doi.org/10.3390/foods12234345
Submission received: 16 October 2023 / Revised: 16 November 2023 / Accepted: 27 November 2023 / Published: 1 December 2023

Abstract

:
Essential oils are mixtures of chemical compounds that are very susceptible to the effects of the external environment. Hence, more attention has been drawn to their preservation methods. The aim of the study was to test the possibility of using the classical model of complex coacervation for the microencapsulation of essential oils. Black pepper (Piper nigrum) and juniper (Juniperus communis) essential oils were dissolved in grape seed (GSO) and soybean (SBO) oil to minimize their loss during the process, and formed the core material. Various mixing ratios of polymers (gelatin (G), gum Arabic (GA)) were tested: 1:1; 1:2, and 2:1. The oil content was 10%, and the essential oil content was 1%. The prepared coacervates were lyophilized and then screened to obtain a powder. The following analyses were determined: encapsulation efficiency (EE), Carr index (CI), Hausner ratio (HR), solubility, hygroscopicity, moisture content, and particle size. The highest encapsulation efficiency achieved was within the range of 64.09–59.89%. The mixing ratio G/GA = 2:1 allowed us to obtain powders that were characterized by the lowest solubility (6.55–11.20%). The smallest particle sizes, which did not exceed 6 μm, characterized the powders obtained by mixing G/GA = 1:1. All powder samples were characterized by high cohesiveness and thus poor or very poor flow (CI = 30.58–50.27, HR = 1.45–2.01).

1. Introduction

In response to consumer expectations, food producers have begun to focus on the production of high-quality and less processed products. These products also come with the longest shelf life. Moreover, a commitment to environmental stewardship has prompted manufacturers to prioritize the reduction or elimination of synthetic preservatives and stabilizers [1]. This has led to a growing emphasis on plant-based alternatives, such as essential oils, to meet these evolving expectations. This is because essential oils have a strong antimicrobial effect, and their use minimizes the risk of foodborne pathogens acquiring resistance to antibiotics. In addition, essential oils have several health-promoting properties [1,2,3].
The main components of juniper essential oil (Juniperus communis L.) are α-pinene, limonene, and myrcene [4,5], while the main ingredients of black pepper essential oil (Piper nigrum L.) are α-pinene, sabinene, β-pinene, δ-3-carene, limonene, and β-caryophyllene [6,7,8]. Both plants are known worldwide as spices but also as pharmaceutical raw materials [9,10]. They have strong antioxidant and antimicrobial properties [7,11], but also have anti-inflammatory and antispasmodic properties [5].
The growing interest in essential oils and the escalating demand for natural additives in the food industry have prompted research into their powdered forms as potential food ingredients [12]. Encapsulating essential oils has become imperative as these substances are sensitive to external environmental factors such as light, temperature, pH, and oxygen. However, their challenges extend beyond environmental sensitivity. Their lipophilic nature, resulting in low water solubility, poor bio-accessibility, and limited bioavailability, hinders their widespread application in food products [13,14,15]. The application of microencapsulation techniques has been proven to be essential in safeguarding the active ingredients of essential oils, enabling their incorporation into various products, such as fruit and vegetable drinks, yogurts, and other dairy products.
Crucially, the choice of the wall material composition for each system is paramount, as the encapsulating material influences the numerous physicochemical properties of the powdered product and dictates its storage behavior [2]. An intriguing method that facilitates microencapsulation is complex coacervation. This process holds a distinctive advantage as it can occur at room temperature, a particularly noteworthy feature when dealing with essential oils [16].
Complex coacervation is one of the most important microencapsulation methods used to protect sensitive substances such as aromas, omega-3 fatty acids, vegetable oils, antioxidants, and essential oils [9,10,15,17,18]. This phenomenon takes place in an aqueous solution when two polymers with opposite charges are drawn together by electrostatic forces. Subsequently, this interaction leads to phase separation, creating two distinct phases: a concentrated coacervated phase referred to as the “continuous phase,” and a less concentrated liquid phase known as the “equilibrium solution” [11,19,20]. The coacervated phase consists of wall material deposited in a thin layer of the core material. The mixtures of proteins and polysaccharides are commonly used, and the most widely studied coacervation system is gelatin (G): gum Arabic (GA) [11,21,22]. This technique allows for high encapsulation efficiency (up to 99%) [14] and enables the encapsulation of more core material per unit mass of wall material compared with other microencapsulation techniques [12,16,17,19,20]. This process has been used to microencapsulate essential oils [9,10,14,16,19]. However, so far, it has been used to do so without first dissolving EO in the oil. Complex coacervation allows for one to obtain a high efficiency of oil encapsulation [23]. However, essential oils, as highly volatile compounds, may suffer significant losses during the spray drying or freeze drying processes [9,10,12,24].
The main aim of this research was to determine whether the classical model of complex coacervation (G:GA) could be used for the microencapsulation of essential oils. The aim of this study was to examine how varying mixing ratios of wall materials (specifically, gelatin and gum Arabic) impact encapsulation efficiency. Additionally, the research sought to explore the physicochemical properties of the resulting powders. In order to minimize the loss of essential oil (due to the increased temperature of the process), the essential oils were dissolved in the oil before being microencapsulated. The literature review shows that, so far, such a solution has not been used. The amount of research into the application of the classical model of complex coacervation to microencapsulation of essential oils is also limited.
The obtained results indicate that the use of the classical model of complex coacervation for microencapsulation of essential oils is possible. This allows a large part of EO to be retained in the microcapsule structures and gives hope for their possible use in food.

2. Materials and Methods

2.1. Materials

Gelatin (Agnex, Białystok, Poland), and gum Arabic (Warchem, Warsaw, Poland) were used as wall materials. Juniper berry (Juniperus communis) essential oil and black pepper peppercorn (Piper nigrum) essential oil (Ancient Wisdom, Sheffield, Great Britain) were firstly dissolved in soybean oil (Dary Natury, Koryciny, Poland) or grape seed oil (Basso Fedele e figli s.r.l., Avellino, Italy) and used as the core material. Essential oils were dissolved in oil at a concentration of 1% v/v to reduce the risk of their evaporation during the freeze drying process.

2.2. Preparation of Coacervates

As a wall material, 5% gelatin solution with 5% gum Arabic solution was used. Gelatin (G) (Agnex, Białystok, Poland) and gum Arabic (GA) (Warchem, Warsaw, Poland) were dissolved at the temperature of 50 °C. The solutions were mixed in different mixing ratios: 1:1; 1:2, and 2:1. The mass of each system was 300 g. After the solutions were mixed together, they were subjected to high shear homogenization using Ultra turrax (IKA T18 basic, Staufen, Germany) for 5 min at 15,000 rpm/min at room temperature. In the homogenization process, every variant received an addition of 10% (30 g) of either soybean oil (SBO) or grape seed oil (GSO), within which 1% of the system mass (3 g) had juniper essential oil (JEO) or black pepper essential oil (BPEO) dissolved in SBO or GSO previously. After emulsification, the pH was adjusted to 4.0 (under the isoelectric point) using 1M HCl solution (Firma Chempur, Piekary Śląskie, Poland). Each emulsion underwent a storage process at 4 °C for 24 h, followed by transfer to −20 °C for an additional 24 h. Subsequently, the samples were further transferred to −60 °C for an additional 24 h. The frozen samples were then subjected to lyophilization for a duration of 72 h at −80 °C. After this period, the lyophilized products were sieved using a laboratory sieve with a mesh size of 710 μm. The screened material was vacuum packed and stored at 4 °C for subsequent analyses.
The samples were coded for easier identification (Table 1). The mixing ratio was marked as 1—1:1, 2—1:2, 3—2:1, juniper essential oil—J, black pepper oil—B, soybean oil—S, grape seed oil—G. For example, the sample in which G/GA = 1:1 and which contains juniper essential oil dissolved in soybean oil was coded as SJ1.

2.3. Complex Coacervation Yield, Solid Yield, and Encapsulation Efficiency

To determine the efficiency of the complex, coacervation yield (CY) was calculated according to the equation [25]:
CY = CM SM 100 %
  • CM—coacervate mass collected after the process,
  • SM—total mass of the freeze-dried powder.
To determine the losses during the freeze drying process, the solid yield (SY %) was calculated according to the equation [25]:
SY = PM LCM 100 %
where:
  • PM—powder mass collected after freeze drying,
  • LCM—liquid coacervate mass, before the freeze drying process.
Encapsulation efficiency (EE) was measured according to the method described by Hernandez-Nava [26] with slight modifications. Briefly, it was calculated by measuring the surface oil (SO) and total oil (TO) of the freeze-dried powders. All measurements were completed in triplicate.
SO was determined by weighing 1 g of sample and dispersion in 30 mL of n-hexane with constant stirring (60 rpm) for 15 min. The oil phase with n-hexane was filtered into the pre-weighted round-bottom flask and evaporated on a rotary evaporator (R-100 Büchi, Switzerland). Then, the sample was stored in the oven at 105 °C for 30 min to ensure all n-hexane had evaporated. After this, the sample was left in the desiccator until it had cooled (1 h). The round-bottom flask was weighed, and the SO was calculated as:
S O = OM 1 OM 2
where:
  • OM1—oil mass after extraction and evaporation of the solvent,
  • OM2—the theoretical weight of oil from the sample.
TO was determined by weighing 1.5 g of the sample and spreading it in 4 mL of KCl. In the next step, 4 mL of acetone and 8 mL of chloroform were added to the sample. The sample prepared was constantly stirred (60 rpm) for 15 min and then centrifuged (5000 rpm, 5 min). The sample was separated during centrifugation. The top layer was poured, and the bottom layer containing chloroform and the oil phase was passed through a filter paper with anhydrous sodium sulfate into a weighed round-bottom flask, so the test sample remained in the falcon test tube. After that, 4 mL of double-distilled water, 4 mL of acetone, and 8 mL of chloroform were added to the residue. The sample was centrifuged again, the top layer was poured, and the bottom layer was soaked through filter paper. Then the solvent was evaporated on a rotary evaporator (R-100 Büchi, Switzerland). After that, the sample was stored in the oven at 105 °C for 30 min to evaporate the excess chloroform. After this, the sample was left in the desiccator until it was cooled down (1 h). The round-bottom flask was weighed, and the TO was calculated as:
T O = OM 1 OM 2
where:
  • OM1—oil mass after extraction and evaporation of the solvent,
  • OM2—the theoretical weight of oil from the sample.
Having the results for SO and TO, EE (%) was calculated according to the formula:
EE = T O S O T O 100 %

2.4. Bulk and Tapped Density, Carr Index (CI), and Hausner Ratio (HR)

The bulk density (ρbulk) was determined by measuring the volume occupied by a known quantity of powder in a 10 mL measuring cylinder, measuring its weight before and after adding the powder [26]. Similarly, the tapped density (ρtap) was determined in a 10 mL measuring cylinder by repeatedly tapping manually by lifting and dropping the cylinder under its own weight at a vertical distance of 14 ± 2 mm high for one minute (one tap per second) [27]. Determinations were made in triplicate. Bulk and tapped densities were expressed in g/cm3.
The compressibility index (CI) and Hausner ratio (HR) are measures of the tendency of a powder to be compressed. Determinations were made in triplicate. To calculate the CI and HR, ρbulk and ρtap were used in these equations [28]:
CI = ρ tap ρ bulk ρ tap 100 %
HR = ρ tap ρ bulk

2.5. Color Measurement

The color of the microcapsules was examined using a Minolta CR-400 colorimeter (Konica Minolta Inc., Tokyo, Japan) (D65 illuminate, measuring surface—8 mm, standard 2° observers). The measurement results are expressed following the system of the International Commission on Lighting (Comission Internationale de L’Eclarige) in the CIELab color space. The following parameters were determined and tested: L* (L = 0 (black), L = 100 (white)), a* (−a = green, +a = red), and b* (−b = blue, +b = yellow) [29]. Determinations were made in triplicate immediately after production.

2.6. Solubility, Hygroscopicity, and Moisture Content

The solubility of samples was measured according to the method described by Shaddel et al. [30]. Briefly, 0.5 g of freeze-dried sample was weighed into a falcon test tube containing 50 mL of double-distilled H2O. The sample prepared in this way was homogenized for 30 min and then centrifuged (5000 rpm, 5 min). After centrifugation, 25 mL aliquot of the supernatant was transferred into a dried and pre-weighed Petri dish and immediately oven dried at 105 °C for 6 h. Then, the solubility (%) was calculated by weight difference. Determinations were made in triplicate.
The hygroscopicity of samples was measured according to the method described by Shaddel et al. [30]. Briefly, 0.2 g of each freeze-dried sample was placed in a Petri dish and stored in a container containing the saturated solution of Na2SO4 for one week. Hygroscopicity was expressed as g of water absorbed per 100 g sample (%). Determinations were made in triplicate.
The moisture content of microcapsules was determined gravimetrically by drying in a Binder FP 115 drying oven (Binder, Tuttlingen, Germany) at 70 °C for at least 24 h, then cooling to room temperature in desiccators until a constant weight was reached. Determinations were made in triplicate.

2.7. Particle Size Distribution

The measurement was carried out with the Morphologi® G3SE apparatus (Malvern Instruments Ltd., Malvern, UK) equipped with a dispersion unit for dry samples. The particle size distribution was calculated as the relative volume of molecules in the band, shown as size distribution curves (Malvern Microsoft ware v.5.40, Malvern Instruments Ltd.). The examined parameters of the size distribution contained the largest particle size (D (v, 0.9)), mean particle volume (D (v, 0.5)), and the smallest particle size (D (v, 0.1)) [31]. The particle size distribution (Span index (SI)) was estimated using the following formula [32]:
SI = D 90 D 10 D 50

2.8. Fourier Transform Infrared Spectroscopy (FT-IR)

FT-IR spectra were recorded on a Nicolet™ iS™ 5 FTIR Spectrometer (Thermo Scientific, Waltham, MA, USA), with horizontal device for attenuated reflectance and diamond crystal, on a spectral window ranging from 4000 to 400 cm−1, at a spectral resolution of 2 cm−1. Spectra were recorded without any sample preparation and were processed with the OMNIC program (Thermo Scientific, Waltham, MA, USA).

2.9. Differential Scanning Calorimetry (DSC)

The thermal properties of the samples were evaluated using differential scanning calorimetry (DSC 1) from Mettler Toledo (Schwerzenbach) under an argon atmosphere at a flow rate of 100 cm3/min, as per the method described by Zhao et al. [33] with some modifications. The instrument was calibrated with pure indium and zinc. Each sample (5.0 ± 0.1 mg) was placed in a standard 40 μL aluminum crucible (ME-51119870) and covered with a lid (ME-51119871) using the Mettler Toledo Crucible Sealing Press. DSC scans were recorded from 10 °C to 230 °C at a rate of 10 °C/min (β). The thermograms were analyzed using STARe Software (Version 9.30) to determine the start (Ton), maximum (Tmax), and end (Tend) temperatures, as well as the areas under the peaks (ΔH).

2.10. Electronic Nose Analysis

Volatile compounds within the microcapsules were extracted using the Heracles II electronic nose (Alpha M.O.S., Toulouse, France) following the methodology outlined in the works of Wojtasik-Kalinowska et al. [34] and Górska-Horczyczak et al. [35]. This electronic nose employs ultra-fast gas chromatography with headspace and features a detection system comprising two metal columns of varying polarities (nonpolar MXT-5 and slightly polar MXT1701, diameter = 180 µm, length = 10 m) and two flame ionization detectors (FID). Kovats indexes were established using alkane standards (n-butane to n-hexadecane) (Restek) measured under the same conditions as the samples. Identification of volatile compounds was accomplished using the AroChemBase (Alpha MOS Co., Toulouse, France), containing 44,000 compounds and a base of sensory descriptors for each compound.
For analysis, 10% solutions of each sample were placed in 20 mL headspace vials sealed with a Teflon-faced silicon rubber cap. These vials were then incubated at 35 °C for 900 s under an agitation speed of 8.33 Hz. Hydrogen was used as the carrying gas at a constant flow rate of 1 mL/min. The injector temperature was set at 200 °C, with an injected volume of 3500 µL and an injection speed of 125 mL s−1. The analytes were collected in the trap at 15 °C and subsequently divided and simultaneously transferred into the two columns. The carrying gas was maintained at a constant pressure of 80 kPa, with a split flow rate of 10 mL min−1 at the column heads. The temperature program in the oven was as follows: 60 °C for 2 s, a ramp of 3 °C s−1 to 270 °C, held for 20 s, and FID1/FID2 at 280 °C. Flavor data profiles were presented through principal component analysis (PCA) using AlphaSoft Version 8.0 software. All samples were analyzed in triplicate.

2.11. Scanning Electron Microscopy Analysis (SEM)

A scanning electron microscope (Jeol JSM6010LA, Tokyo, Japan) was employed to analyze the samples’ surface morphologies after the different treatments. Each sample was fixed on an adhesive carbon tape and coated (approximately 15 nm) with gold. Secondary electron images were taken at several magnifications and an accelerating voltage of 10 kV.

2.12. Statistical Analysis

For statistical analysis, the STATISTICA computer program was used. To check significant differences between the results, one-way and multivariate analyses of variance ANOVA and the Fisher LSD test (p-value < 0.05, α = 95%) were used.

3. Results

3.1. Evaluation of Coacervation Efficacy

The critical factor influencing the complex coacervation process is the protein-polysaccharide mass ratio. Different ratios impact the intensity of interaction and complexation due to the charge balance between protein and polysaccharide [26,36,37]. The statistical analysis demonstrated that the obtained results had been affected by the mixing ratio, essential oil, the type of oil used, and their interactions. However, the mixing ratio had the most significant impact on the resulting coacervation yield (CY) and solid yield (SY), closely followed by the essential oil (EO) (Table 2). The mixing ratio exhibited a strong positive effect on CY but had a negative impact on SY. Similarly, the type of oil had a small effect, with a slight decrease in CY and a slight increase in SY values. Notably, the content of essential oil caused an opposite effect—it led to a slight decrease in CY and a slight increase in SY values. In this study, both CY and SY values were relatively low, not exceeding 50%. Samples GJ3, GB3, SJ3, and SB3 exhibited the highest CY values, with significant differences (p < 0.05) among them. Despite having the highest CY values, these samples did not reflect the highest SY. The highest SY values were observed for samples SB2, GB2, GJ2, and SB1, characterized by the lowest CY values. The specific molar mixing ratio for the achievement of maximum coacervate yield occurs when protein and polysaccharide have the exact opposite charge density, leading to charge neutralization [14,38,39]. This likely occurred at a mixing ratio of 2:1, where the ratio of gelatin to gum Arabic was the highest, resulting in the highest CY. Gelatin, rich in positively charged NH3+ ions, and gum Arabic, with primarily negatively charged COO- ions, undergo electrostatic attraction in the coacervation process [40]. Despite gelatin’s high water-binding capacity, these polymers form weak hydrogen bonds with water, which disintegrate and release water during the freeze drying process [41].
Encapsulation efficiency (EE) represents the percentage of the core material enclosed within the powder particles, a critical parameter for essential oils given their volatility and susceptibility to loss during the drying process. Statistical analysis revealed that the results had been significantly influenced by the interaction between the mixing ratio of polymers and oil (Table 2). The type of oil exhibited a notable positive impact on EE, leading to an increase in its values. The highest EE was observed for sample SJ1 (64.09% ± 0.09). Additionally, notable encapsulation efficiency was achieved for samples SB3 and GJ2 (61.92 ± 0.04 and 59.89 ± 0.01, respectively), with significant differences (p < 0.05) in EE. Samples GB2, SB1, and SJ2 did not show significant differences (p < 0.05) in EE, with values falling within the middle of the obtained range (53.4–55.25%). Notably, similarities were identified in samples SB1 and GB2, which shared the same oil and essential oil, as well as in samples GB2 and SJ2, featuring the same mixing ratio (1:2). The influence of the interaction between mixing ratio and oil on the EE value is evident in these cases. Samples GJ1 and GJ3 were also within the same statistical group (49.30 ± 0.07 and 49.65 ± 0.02, respectively), containing the same essential oil. The lowest EE value was recorded for sample GB1 (Table 2).

3.2. Bulk Density, Tapped Density, Carr Index (CI), and Hausner Ratio (HR)

Bulk density (ρbulk) is defined as the ratio of the mass to the volume (including the inter-particle void volume) of an untapped (loose) powder sample. If the bulk density is high, the volume of the powder is lower; when the bulk density is low, the same powder mass takes up a larger volume. The tapped density (ρtap) is defined as an increased bulk density attained after mechanically tapping a cylinder containing the powder sample. This property influences the appropriate selection of the size of the containers for packing the powder (e.g., barrels, bags). Analysis of bulk and tapped density of powder gives a possibility to calculate the compressibility index (Carr index) and Hausner ratio of the powder [28].
The Carr index (CI) serves as an indicator of powder bridge strength and stability, while the Hausner ratio (HR) reflects interparticulate friction [42]. These metrics are instrumental in assessing the flow characteristics of powders. Powder bridging takes place when particles interlock or bond, forming an arch above the container outlet. The cohesive strength and internal friction of individual particles together play a crucial role in determining the stability and strength of this arch. In instances where a powder experiences bridging, the arch effectively holds the remaining contents within the container, impeding the discharge of the remaining powder. A powder’s lower CI or lower HR indicates better flow properties than higher ones. A CI of <10 or HR of <1.11 is considered ‘excellent’ flow, whereas CI > 38 or HR > 1.60 is considered ‘very very poor’ flow [43]. In CI and HR, the obtained values were influenced by the mixing ratio and its interaction with the oil. The increase in the mixing ratio resulted in a decrease in the cohesiveness of the obtained powders, and thus a decrease in the CI and HR values. The samples differed statistically significantly (p < 0.05). The GB2, GB3, SJ2, and SB2 samples were in the same statistical group (Table 3), presenting the CI in the range of 31.81–34.74 and the HR in the range of 1.47–1.53, thus demonstrating high cohesiveness and very poor flow [44]. These were also the lowest CI and HR values obtained in this study. The highest values were recorded for the GJ1 and GB1 samples, which were also together in the same statistical group (p < 0.05). CI values were 50.27 ± 2.53 and 48.84 ± 7.73, respectively, and HR values were 2.01 ± 0.11 and 2.98 ± 0.28, respectively. These powders had the highest cohesiveness and virtually no flow [45].
The powders obtained in this study were characterized by a very low bulk density in the range of 0.11–0.17 and a tapped density from 0.22 to 0.28 (Table 3). The tapped density exceeded the bulk density because smaller particles filled the voids between larger particles, achieving a more compact packing condition due to the tapping process [46]. The greatest influence on the obtained bulk density was the mixing ratio and the interaction between the mixing ratio, oil, and essential oil. All factors caused an increase in bulk density. The samples showed statistically significant differences among themselves (p < 0.05). Similarly, in the case of tapped density, the mixing ratio and interactions between it and the oil and essential oil had the greatest impact on the obtained results. The influence of the mixing ratio was so positive that it caused the tapped density value to increase as it increased. Oil also had a slight positive effect while EO had a slight negative effect. The samples differed statistically significantly (p < 0.05).

3.3. Solubility, Hygroscopicity, and Moisture Content

In the food industry, the solubility of powders in water is a crucial factor to expand their range of applications. However, complex coacervation aims to produce water-insoluble microcapsules, allowing for controlled release of the core material. In this study, all samples exhibited low solubility in water (>26%) (Table 3). The solubility of the powders was primarily influenced by the mixing ratio of polymers, as well as the interaction between the mixing ratio and the oil, and the interaction between the mixing ratio and essential oil (p < 0.05). These interactions led to a decrease in the solubility of the powders. The least soluble samples were GJ3, GB3, and SB3, forming a homogeneous statistical group with solubility below 8%. The impact of the mixing ratio and the oil used to dissolve the essential oil is evident here, with sample SB2 showing the highest solubility (25.86% ± 6.81).
Hygroscopicity is a critical factor affecting the stability of food products, storage conditions, and packaging materials. Powders with hygroscopic values exceeding 15.1% are considered highly hygroscopic [47]. However, the hygroscopicity of the powders in this study ranged from 1.49% to 7.29% (Table 3), indicating low hygroscopicity. The examined factors did not have a statistically significant influence on the hygroscopicity of the obtained powders (p < 0.05), with no statistical differences found between samples.
Moisture content is another important property related to storage stability and shelf life. The moisture content in the obtained powders was low, ranging from 0.05% to 0.56% (Table 3). Similar to hygroscopicity, the examined factors did not have a statistically significant influence on the moisture content of the powders (p < 0.05), and no statistical differences were observed between samples.

3.4. Color Measurement

The color of the resulting powders was primarily influenced by the oil used to dissolve the essential oil, with grape seed oil imparting a slight greenish tint and soybean oil exhibiting a vivid yellow hue. However, the powder color was also affected by the mixing ratio, essential oil, and the interactions among these factors.
Statistical analysis revealed significant differences (p < 0.05) in each of the color components among the samples (Table 4). Regarding the brightness of the samples (L*), the variations were not substantial, with L* parameter values ranging from 78.72 ± 0.34 to 88.52 ± 0.25 for samples containing soybean oil and grape seed oil. Statistical homogeneity was determined for GJ3, GB3, and SJ1 samples, as well as for SJ2 and SB2. The a* parameter, indicative of the green color, exhibited negative values for all samples. For samples containing grape seed oil, the a* values ranged from −1.19 ± 0.02 to −3.83 ± 0.15. The mixing ratio played a significant role in influencing the a* value, with an increase corresponding to a higher proportion of green color. Samples containing soybean oil showed a* values ranging from −7.06 ± 0.27 to −7.66 ± 0.05, and a homogeneous group was identified for samples SB1 and SB2. In terms of the b* parameter, representing the yellow color, positive values were observed for all samples. For samples containing grape seed oil, b* ranged from 5.61 ± 0.04 to 12.3 ± 0.47, while for those with soybean oil, the range was 41.87 ± 1.05 to 51.71 ± 0.47.

3.5. Particle Size Distribution

The size of microcapsules should be in the range of 0.1 μm to 100 μm [12]. The particle sizes in all tested trials met these assumptions. The obtained results depend on all examined factors individually and their interactions. However, the mixing ratio of polymers had the greatest influence (p < 0.05). The highest values of the SI characterized the GJ3, GB3, SB3, and SJ3 samples. This means they had the highest diversity in terms of particle size (Table 5) [48]. In sample GJ3, the SI was the closest to 1.0 (SI = 0.96), so this powder contained particles with the most uniform size distribution in the tested ranges [48]. The smallest range of particle sizes and the smallest SI value were found for samples SB2 and GB2 (D10 = 59.04 ± 1.18 μm, D90 = 96.72 ± 0.5 μm, SI = 0.43 ± 0.01 and D10 = 60.13 ± 0.74 μm, D90 = 96.86 ± 0.41 μm, SI = 0.44 ± 0.01, respectively). The smallest powder particles were obtained for a mixing ratio of 1:1 (GB1, SJ1, SB1). However, no clear relationship was observed between the increase in the content of one polymer and the size of the obtained powder particles.

3.6. FT-IR

The obtained FT-IR spectra for JEO showed characteristic peaks at 2917.56 cm−1, 2878.29 cm−1, 1446.07 cm−1, 887.14 cm−1, 786.52 cm−1, and 418.96 cm−1. For BPO, characteristic peaks occurred at 2954.86 cm−1, 2922.95 cm−1, 2867.06 cm−1, 1446.24 cm−1, 885.83 cm−1, 875.20 cm−1, 786.32 cm−1, 543.65 cm−1, 421.49 cm−1, and 442.11 cm−1. Based on the obtained results, the compositional similarity of both essential oils can be seen, which was confirmed by the analysis of volatile compounds using the e-nose (3.8).
In the FT-IR spectra for wall materials, characteristic peaks were observed at wavenumber 3300.05 cm−1, 2893.28 cm−1, 1596.88 cm−1, 1413.05 cm−1, 1017.17 cm−1, 448.52 cm−1, 437.30 cm−1 and 416.39 cm−1 for gum Arabic and at 1627.16 cm−1, 1526.30 cm−1, 1444.40 cm−1, 1332.90 cm−1, 1235.97 cm−1, 1078.27 cm−1, 597.00 cm−1, 507.22 cm−1, 490.10 cm−1, 469.20 cm−1, 448.82 cm−1, 416.66 cm−1, 409.46 cm−1 and 402.23 cm−1 for gelatin. Regular soy oil contains approximately 54% linoleic acid (18:2), 23% oleic acid (18:1), 11% palmitic acid (16:0), 8% linolenic acid (18:3), and 4% steric acid (18:0) [49]. Grape seed oil contains the same fatty acids in similar amounts, hence the similarity of the obtained FT-IR spectra [50]. Characteristic peaks were observed at wavenumber 3008.09 cm−1, 2922.44 cm−1, 2852.97 cm−1, 1742.88 cm−1, 1457.23 cm−1, 1377.09 cm−1, 1159.43 cm−1, 1097.51 cm−1 and 721.59 cm−1 for grapeseed oil. For soybean oil characteristic peaks were found at 3008.58 cm−1, 2922.40 cm−1, 2852.84 cm−1, 1742.92 cm−1, 1456.95 cm−1, 1377.02 cm−1, 1159.01 cm−1, 1097.76 cm−1, and 720.99 cm−1.
For gum Arabic a broad absorption band was observed at wavenumber 3300.05 cm−1 which corresponds to hydrogen bond. Because this was followed by the presence of spectra at the 1600–1300 cm−1 (1596.88 cm−1, 413.05 cm−1) frequencies, we can confirm the existence of a hydroxyl (-OH) group. The next narrow band was found at 2893.28 cm−1 followed by peaks between 1470–720 cm−1, which can correspond with absorption band for long-chain linear aliphatic compounds. A characteristic peak at 1017.17 cm−1 corresponds to the fingerprint of carbohydrates, which, along with peaks at 448.52 cm−1 and 437.30 cm−1, is characteristic for CCO, COC, symmetrical, and asymmetrical ring breathing vibration. For gelatin, characteristic peaks were found at 1627.16 cm−1, which corresponds to COO asymmetric stretching, and at 1526.30 cm−1 and 1444.40 cm−1 which can be associated with COO symmetric stretching. Peaks at 1332.90 cm−1, 1235.97 cm−1, and 1078.27 cm−1 can be attributed to C=O stretching. The remaining characteristic peaks at 597.00 cm−1, 507.22 cm−1, 490.10 cm−1, 469.20 cm−1, 448.82 cm−1, 416.66 cm−1, 409.46 cm−1 and 402.23 cm−1 are CCO, COC, symmetrical and asymmetrical ring breathing vibrations [51].
Figure 1 and Figure 2 show the FT-IR spectra for individual microcapsules. The spectra of the microcapsules were similar to those of the wall materials. Most of the peaks corresponding to the essential oils disappeared in these spectra. This phenomenon can be related to the overlapping of the peaks of the matrix, oils, and essential oils, which is due to the low concentration of essential oils in the total weight of the microcapsules [24]. Nevertheless, the obtained results confirm that both oils were successfully encapsulated. This is evident from the characteristic peaks at 1455.95 cm−1, 1454.10 cm−1, 1454.84 cm−1 and 1453.80 cm−1 seen in samples GJ3, GB3, SJ1 and SB3, respectively. It is also worth noting that samples SJ1 and SB3 have the highest encapsulation efficiency among the remaining samples.
In all samples, a characteristic broad absorption band between 3305.44–3308.15 cm−1 can be seen. This is characteristic for gum Arabic and hydroxyl groups [52]. For samples SB2, SJ3 and GJ3 this peak is followed by characteristic narrow peaks between wavenumber 3008.19–3008.43 cm−1 which can be linked with unsaturated compounds or aromatic rings characteristic for essential oils. Those peaks are not visible in other sample spectra. For all samples, peaks between 2922.73–2924.19 cm−1 and 2853.23–2854.13 cm−1 were found. Since it was found in the single bond area those peaks can contribute to long-chain aliphatic compounds [53]. There were no characteristic peaks in the triple-bond region (2000–2500 cm−1). In the double-bond region (1500–2000 cm−1) for all samples, peaks were found between 1743.15–1744.20 cm−1. Those peaks can describe carbonyl compounds such as aldehydes, ketones or esters characteristic for essential oils, at the same time being a proof for successful EOs encapsulation. These peaks for samples SB2, SJ3 and GJ3 were much more intense compared with the rest of the graphs, which may indicate a higher content of essential oils in these samples (according to the EE% results). In the same region we found peaks at 1634.11–1646.11 cm−1, which might be because of the presence of unsaturated bonds, probably C=C. The rest of the peaks are characteristic for the fingerprint region (600–1500 cm−1)—those peaks can contribute to the methylene C-H bond (1485–1445 cm−1), methyne C-H bond (1350–1330 cm−1), C-O stretch (~1200 cm−1, ~1150 cm−1), and CN stretch (1210–1150 cm−1). Again, only for samples SB2, SJ3 I GJ3 were characteristic peaks at 720.68–721.18 cm−1 found and these can be linked with the aromatic ring coming from essential oils [53].

3.7. DSC

DSC analysis was performed to verify the thermal behavior after encapsulation of the essential oils. In the curves showing the thermal behavior of the wall materials (gelatin, gum Arabic), an endothermic peak can be seen in each case—60.58 ± 0.002 °C (−1430.91 ± 0.001 mJ) for gelatin and 137.93 ± 0.001 °C (−696.52 ± 0.001 mJ) for gum Arabic. Those peaks probably refer to a glass transition. DSC curves for juniper and black pepper essential oils show an endothermic event at 24–25 ± 0.001 °C (−80.00 ± 0.001 mJ and −87.04 ± 0.002 mJ, respectively) related to the residual water and at 150 ± 0.001–158 ± 0.001 °C (−350.87 ± 0.001 mJ and −367.14 ± 0.001 mJ, respectively) related to the decomposition of JEO and BPO [24]. Table 6 shows DSC onset, peaks and endset temperatures with corresponding enthalpies for all designed microcapsules. Onset and endset temperatures for GJ1, GJ2, and GJ3 were found to be 63.67 ± 0.002–87.18 ± 0.002 °C, 56.93 ± 0.002–161.87 ± 0.002 °C and 68.75 ± 0.001–87.94 ± 0.002 °C, respectively, with different enthalpy values (−36.31 ± 0.001 mJ, −14.00 ± 0.001 mJ and −20.52 ± 0.002 mJ, respectively). Onset and endset temperatures for GB1, GB2 and GB3 were found to be 46.21 ± 0.001–113.29 ± 0.001 °C, 83.46 ± 0.001–99.21 ± 0.001 °C and 94.72 ± 0.001–151.31 ± 0.002 °C, respectively, with the following respective enthalpy values: −222.81 ± 0.001 mJ, −14.00 ± 0.001 mJ and −327.77 ± 0.001 mJ. For samples SJ1, SJ2 and SJ3 onset and endset temperatures were 55.23 ± 0.001–175.95 ± 0.001 °C, 82.32 ± 0.001–108.89 ± 0.001 °C and 68.86 ± 0.001–89.36 ± 0.001 °C, respectively. Each sample exhibited unique enthalpy values, indicative of the heat absorbed during thermal transitions. The respective enthalpy values were: −186.11 ± 0.001 mJ, −51.25 ± 0.001 mJ and −0.45 ± 0.002 mJ. Similarly, for samples SB1, SB2 and SB3, onset and endset temperatures were 59.94 ± 0.001–178.75 ± 0.001 °C, 70.13 ± 0.001–84.80 ± 0.001 °C and 71.79 ± 0.001–161.17 ± 0.001 °C, respectively, with enthalpy values of −285.82 ± 0.002 mJ, −14.50 ± 0.001 mJ and −210.46 ± 0.001 mJ, respectively. Due to the lack of peaks characteristic of JEO and BPO, it can be assumed that these essential oils have been successfully encapsulated in the wall material [24]. In terms of thermal stability, samples containing soybean oil were found to be the most stable. Additionally, the samples with the G/GA ratio of 1:1 were characterized by the highest onset temperature. From this it can be concluded that those samples would be the best variant of EO microcapsules, allowing their use in processes requiring the use of elevated temperature (not higher than 68 °C).

3.8. Electronic Nose Analysis

Figure 3 presents the classification of scent profiles in relation to their experimental groups. Samples are represented in a two-dimensional plane with reference to selected components: principal component 1 and principal component 2. When considering essential oils separately, it can be concluded that, in the case of grape seed with juniper as an essential oil, the total contribution variances of PC1 and PC2 for direct electronic nose measurements were 60.6% and 14.27%, respectively. The combination of grape seed and black pepper as essential oil resulted in values explaining 65.23% of the data variance, where 10.81% was intercepted by the horizontal axis, explaining differences among samples along that axis. For soybean with juniper, the total contribution variances of PC1 and PC2 for direct electronic nose measurements were 63.21% and 14.10%, respectively. Soybean with black pepper caused a 67.56% data variance, with 10.51% intercepted by the horizontal axis. Additionally, there was no overlap between the mixing ratios of the single samples for all essential oils.
Table 7 lists the Kovats indexes along with the volatiles identified in each sample and their corresponding sensory descriptors. According to literature data, juniper (Juniperus communis) essential oil contains mainly α-pinene (25.3%), β-pinene (2.98%), sabinene (5.99%), β-myrcene (12.6%), α-terpinene (0.21%), γ-terpinene (0.80%), and limonene (2.99%) [7]. All of the abovementioned compounds except for sabinene and α-terpinene were identified in all of the samples containing juniper essential oil. These two compounds are highly volatile, and their absence in the samples can be explained by their probable evaporation during the freeze drying process [54]. Samples containing EO from juniper dissolved in soybean oil were also characterized by the content of 3-methyl-octane, 2-butylfuran and dimethyl disulfide. In turn, those with grape seed oil were the only ones that contained 3-methyl−2-butene-1-thiol, geranial and thymol. All JEO capsules contained nonane and propyl nonanoate. From this it can be seen that the oil in which the EO is dissolved has an effect on the composition of the EO volatiles after freeze drying. Similarly, in the samples containing black pepper essential oil, the compounds characteristic of this oil were identified as α-pinene, β-pinene, myrcene, and limonene. Sabinene and β-caryophyllene were not identified [55]. In addition, 5-pentanolide and isoamylacetate were identified in samples containing BPO.
In each of the samples, ethanol and 2-propanol were identified, which are probably residues from the extraction process [56].

3.9. SEM

Figure 4 shows SEM images of lyophilized microcapsules. As described earlier, the encapsulation efficiency was not very high (42.7–64.09%). This fact can be explained by the morphology of the microcapsules. The obtained powder was characterized by an irregular, very porous structure with a large surface area, suggesting that the core material was not completely covered. For this reason, oxygen availability can be high, causing oxidation of both the oil and the essential oil [51]. In addition, the highly porous surface of the microcapsules may make it easier for the essential oil to evaporate during storage. An increase in the concentration of polymers can make the wall of the microcapsule thicker, which can have a direct effect on the size of the resulting capsules [54]—in figures 12–15 it is clear that the microcapsules containing G:GA = 1:1 were smaller than those where the concentration of G:GA was 1:2 or 2:1.

4. Discussion

4.1. Evaluation of Coacervation Efficacy

Karagozlu [57], using a complex coacervation between gelatin and acacia, obtained a SY of about 26.66%. Hernandez-Nava [26] used G and GA as wall material (mixing ratio 2:1) and obtained an SY in the 72–88% range. These differences may be caused primarily by the difference in pH during the coacervation process and the difference in the selection of the drying method.
The most common method of encapsulating essential oils is spray drying with the use of various types of carriers—whey protein (WP), maltodextrin (MD), inulin (IN), gelatin, gum Arabic, various types of modified starches (MS), etc. The use of this method allows for encapsulation efficiencies at levels of 47.13% for WPI, 61.64% for WPI:MD, 48.14% for WPI:IN [12], 87.5–94.5% for G:GA (depending on the air temperature used at the inlet) [26] or even 98.87% for MD:MS [42]. The encapsulation efficiency for spray drying depends largely on the wall materials used, as is the case with complex coacervation. Here, however, the mixing ratio of the polymers is also important. It is well known that encapsulation efficiency increases with increasing polymer concentration, however, this is not the only factor affecting EE. The interactions between the polymers used and the amount of core material also play an important role, as shown by the results of our research. Manaf et al. [44] used simple and complex coacervation to encapsulate citronella EO. For simple coacervation they used GA (12,5% w/v), for complex coacervation they used G:GA (mixing ratio 1:1, 12,5% v/w), and in both cases they used 1% of core material. The encapsulation efficiencies of the system prepared in this way were 93.91% and 94.42%, respectively. The team did not use any method of drying the obtained coacervates, liquid forms were used for the research. Various concentrations of wall and core material were tested by Samakradhamrongthai et al. [45]. A complex coacervation with spray drying was performed between G and GA at pH = 4.0 with a mixing ratio of 1:1; 1:2, and 2:1. The results show that the lowest EE was obtained with a mixing ratio of 1:1 and a core material content of 10% (38.56%). The highest was obtained with a mixing ratio of 2:1 and a content of core material of 5% (95.15%). Much less research is available where lyophilization has been used to dry the coacervates. This method gives good results. Freeze drying of complex coacervates made between WPI and carbomethylcellulose (CMC), WPI and sodium alginate (SA) as well as WPI and chitosan (CH), has resulted in EE at the levels of 83.94%, 79.28% and 82.88%, respectively [19]. A complex coacervation with freeze drying has also been performed between G and GA [57]. The encapsulation efficiency reported was 64.31%. In our research, the encapsulation efficiency was similar to that obtained by Karagozlu et al. [57] and slightly lower compared with results obtained by Rojas-Moreno et al. [19] (max. 64.09%, min. 42.7%). To increase the efficiency of encapsulation, in our case, the amount of core material should probably be reduced (less than 10%) or the concentration of wall materials should be increased (more than 5%).

4.2. Bulk density, Tapped Density, Carr Index (CI), and Hausner Ratio (HR)

Bulk and tapped density depend not only on composition of the wall material, but also on the drying method. The use of WPI:IN (3:1, 1:1, 1:3) as a wall material for rosemary EO encapsulation by spray drying (170 °C) resulted in powders with bulk densities of 0.20 g/cm3, 0.24 g/cm3, and 0.26 g/cm3, respectively. Tapped densities were slightly higher—0.30 g/cm3, 0.34 g/cm3, and 0.41 g/cm3, respectively [12]. On the other hand, Yekdane et al. [58], using spray drying (170 °C), microencapsulated pomegranate seed oil and reported that the bulk density, depending on the content of pomegranate juice in the wall material, ranged from 0.4 to 0.5 g/cm3. Hernandez-Nava et al. [26], used the complex coacervation with spray drying between G:GA and gelatin and chia mucilage (CM) to microencapsulate oregano EO. Depending on the spray drying parameters (160 °C or 180 °C), the bulk density in their study was in the range of G:GA 0.157–0.202 g/cm3 and 0.234–0.282 g/cm3 for G:CM. In contrast, the tapped density for G:GA was 0.314–0.404 g/cm3, and for G:CM 0.335–0.403 g/cm3. However, a low bulk density may indicate a high air content in the microcapsules, which in turn may favor the negative effect of lipid oxidation. The authors of [57], investigating the bulk density of freeze-dried coacervates (G:GA) containing oregano essential oil, indicated that the bulk density of the powders thus obtained was 0.527 g/cm3. In our case, the bulk density was lower, probably due to different solid content.
The use of MD:WPI as wall materials for spray drying and thus encapsulation of olive oil allowed us to obtain a CI in the range of 25.61–54.70 [59]. Such a large dispersion of values was caused by different mixing ratios of polymers. Hernandez-Nava et al. [26], have reported that the powders obtained by spray drying G:GA coacervates were characterized by a Carr index of approximately 50, while the HR was approximately 2. Similarly, for the G:CM system, CI was approximately 30 and HR approximately 1.43. Another team used complex coacervation between soy protein and GA to microencapsulate chia seed oil. The obtained coacervates were spray dried (130 °C) and were characterized by CI = 35 and HR = 1.25 [23]. In turn, freeze drying of MD-based emulsions (10, 15 and 20% v/w) allowed for the obtaining of powders characterized by a CI in the range of 11.62 to 17.64 and HR in the range of 1.13–1.21 [60]. These values were much lower than in our case, where we freeze dried coacervates (CI min. 30.58, max. 50.27, HR min. 1.45, max. 2.01). Therefore, it can be concluded that the powders obtained as a result of complex coacervation are characterized by very high cohesiveness and virtually no flow.

4.3. Solubility, Hygroscopicity, and Moisture Content

Fernandez et al. [12] investigated the solubility of spray-dried powders containing ginger EO. Depending on the wall material used (WPI, WPI:MD, WPI:IN), the solubility of the powders was in the range of 76.94–81.58%. Another team [19], using MD:MS as a wall material for microencapsulation of orange EO (by spray drying), obtained powders with a solubility of 57.10%. Caparino et al. [27] investigated the water solubility of mango powders obtained from the emulsion with MD by spray and freeze drying. Powders obtained by spray drying had a solubility of 95.31%, while those obtained by freeze drying had a lower solubility—89.70%. Thus, the influence of the choice of method on the solubility of the obtained powder is visible—lyophilization makes the powders less soluble. In our case, the solubility was low and did not exceed 26%. These results correlate with those obtained by Shaddel et al. [30]. The powders obtained by freeze drying coacervates (5% v/v G/GA) were characterized by low water solubility, not exceeding 30%. Similarly, in the case of freeze-dried coacervates obtained by mixing G:GA in the ratio 1:1 with the addition of shrimp lipid extract [61]—the solubility of the obtained powder was low—approximately 9.6%. The above indicates that not only the method of drying, but also the type of dried material (emulsion, coacervates), has a significant impact on the solubility of powders. Powders obtained by the freeze drying of complex coacervates show much lower solubility in water than those obtained from the classical emulsion.
Zotarelli et al. [47] found that mango powders obtained by the spray drying of mango pulp with maltodextrin (MD) displayed high hygroscopicity in the range of 23.9–26.9%. In a study by Caparino et al. [27], the hygroscopicity of mango powders obtained by spray and freeze drying of mango puree with MD was reported to be around 16.5% for spray-dried samples and approximately 18% for freeze-dried samples. Similarly, in the case of the spray drying (at 170 °C) of an emulsion containing WPI:IN and rosemary essential oil, the hygroscopicity of the resulting powders varied from 15.7% to 17.1%, depending on the mixing ratio of the wall materials [32].
In contrast, powders obtained by the freeze drying of coacervates showed significantly lower hygroscopicity. Tavares et al. [62] have reported that coacervates with ginger essential oil exhibited lower hygroscopicity than coacervates without, ranging from 6.64% to 9.65%, which is below the threshold of 15.1%. This aligns more closely with the values obtained in our study. Gomez-Estaca [59] demonstrated that coacervates (G:GA, freeze dried) containing shrimp lipid extract had hygroscopicity at the level of 2.33 ± 0.12 g/100 g. Consequently, it can be inferred that both the addition of essential oil and the utilization of complex coacervation as a microencapsulation method have a positive effect on reducing the hygroscopicity of the obtained powders. Based on the results obtained and their comparison with literature data, it can be concluded that the obtained microcapsules could be stored for a longer time without compromising their properties.
The moisture content of powders obtained through various encapsulation methods and wall materials can vary significantly. Fernandes et al. [12] demonstrated that ginger essential oil (EO) powders obtained by spray drying had moisture content ranging from 1.05% to 1.98%, depending on the wall material used (WPI, WPI:MD, WPI:IN). Similar results were reported by Zotarelli et al. [47], where the moisture content of spray-dried mango powders with maltodextrin (MD) was 1.5%. Slightly higher hygroscopicity was observed in spray-dried coacervates, with moisture content ranging from 3.86% to 4.55% (G:GA) and 3.49% to 4.58% (G:CM) [26]. Rojas-Moreno et al. [19] found that the moisture content of spray-dried microcapsules containing orange essential oil varied from 1% to 4.5%, while freeze-dried samples exhibited slightly higher moisture content (2% to 5.5%). This difference in moisture content between spray-dried and freeze-dried samples could be attributed to the presence of oil droplets, which act as a vapor transport barrier, reducing the evaporation of water and increasing the hydrophobicity of the microcapsules. The freeze drying of coacervates containing ginger essential oil, with different wall materials (GA:CH, WPI:GA), resulted in powders with moisture contents of approximately 3.04% and 3.10%, respectively. The addition of essential oil significantly reduced these values, as coacervates without essential oil had moisture contents of approximately 6.01% and 8.02%, respectively [63]. The freeze drying of coacervates (G:GA) containing oregano essential oil produced powders with a water content of 3.39% [57].
The selection and concentration of wall materials play a crucial role in determining the moisture content of the obtained powders. The results obtained in our study (moisture content = 0.05% to 0.56%) are consistent with those obtained by Samakradhamrongthai et al. [45], who, by mixing gelatin and gum Arabic in a 1:1 ratio, obtained powders with moisture content ranging from 0.32% to 0.72% for different concentrations of core material. Based on the above, it can be concluded that the use of gelatin with gum Arabic in the coacervation process allows one to obtain powders with low moisture content, and that are thus suitable for long-term storage.

4.4. Color Measurement

Syed et al. [63] conducted color measurements of, inter alia, unrefined soybean oil. The results obtained in the study correlate with those obtained in this study in terms of the b* parameter—the researchers indicate the value of b* as 40.0 ± 0.60. There is some discrepancy in the case of the a* parameter. The team indicated positive values (4.0 ± 0.65), while the results we obtained clearly indicate negative a* values. Such a change in the a* parameter should be attributed to the wall materials used for encapsulation and to the presence of EO. Bruhl and Unbehend [64] carried out measurements of the color of grape seed oil. The L* parameter obtained by the researchers was similar to the results obtained in this study L* = 99.44. In turn, parameters a* and b* had higher values, although with the same sign: a* = −4.68, b* = 12.49. Again, the discrepancy in the results was largely caused by wall materials and essential oil.

4.5. Particle Size Distribution

Several studies have investigated particle size and similar parameters for powders produced by various methods, such as spray drying and complex coacervation. Rojas-Moreno et al. [19] examined powders obtained by spray drying coacervates (WPI:CMC, WPI:SA, WPI:CH), reporting particle sizes of 6.31 μm, 10.59 μm, and 9.44 μm, respectively, along with surface index (SI) values of 0.95 μm, 1.03 μm, and 10.6 μm, respectively.
Hernandez-Nava et al. [26] conducted similar studies involving complex coacervation between gelatin (G) and gum Arabic (GA). Depending on the parameters, they observed particle sizes (D10, D50, and D90) ranging from 6.08 μm to 30.31 μm. It is worth noting that spray drying generally allows for smaller particle sizes [25]. For instance, the particle size of spray-dried smoke powder food flavoring has been reported to be 6.3–6.9 μm, significantly smaller than freeze-dried preparations (134.7–580.4 μm) [65]. Freeze-dried G:GA coacervates had a particle size of 41.26 µm. Tavares et al. [63] conducted complex coacervation between WPI:CH and GA:CH, followed by freeze drying. Their results include D10 of approximately 25.06 μm, D50 of approximately 103.39 μm, D90 of approximately 293.51 μm, and an SI of approximately 2.60 for WPI:CH. For GA:CH, the results were D10 ≈ 19.18 μm, D50 ≈ 74.23 μm, D90 ≈ 189.21 μm, and SI ≈ 2.29. In our study, the screening of the obtained lyophilizates on a laboratory sieve played a crucial role in achieving size uniformity and obtaining relatively small particles, consistent with the standard for lyophilizates. This method allowed for the control and standardization of particle sizes.
Particle size stands as a crucial quality parameter in the determination of microcapsules’ application areas. It significantly influences delivery properties and the flowability of the powder, playing a pivotal role in the selection process. Complex coacervation between CH and GA [66] in the ratio 1:1, 1:2 and 2:1 has shown that the ratio of the wall material has a significant effect on the particle size of the microcapsules. The particle size for the samples was 16.39 μm, 28.98 μm, and 49.53 µm, respectively. Particle size is significantly affected by the proportion of polymer used in the wall material, and as the particle size decreases, particle–particle interactions become easier due to the increase in surface area [66].

4.6. SEM

Marfil et al. [67], using the same mixing ratios G:GA (1:1, 1:2, 2:1), obtained microcapsules with a completely different surface—smooth, without pores—and individual microcapsules were connected with each other via solid bridges, characteristic of freeze-dried products. This type of surface suggests that the entire core material was covered with wall material and the encapsulation efficiency should be high. Meanwhile, in the case of capsules that had the same concentration of wall materials and core material as in the present study, the encapsulation efficiency was quite similar, ranging from 50.08 to 83.5%. Differences in the obtained structure may be due to the difference in lyophilization temperature, which in the case of studies conducted by Marfil et al. [68] was higher (−40 °C) than in this study (−80 °C). As described by Krokida et al. [69]—where the porosity of the obtained product depends on the lyophilization conditions, including the temperature—the lower the temperature, the higher the porosity. The same was later confirmed by Barresi et al. [70]. Perhaps increasing the lyophilization temperature would allow for higher encapsulation efficiency in the case of this study.

5. Conclusions

Complex coacervation between gelatin and gum Arabic can be used to microencapsulate essential oils. However, not all of the combinations of mixing ratios were equally effective. The highest encapsulation efficiency characterized the SJ1 sample.
A relatively low EE was obtained, influenced by the interaction between the mixing ratio and the oil. A low efficiency was probably caused by a too-high concentration of core material in relation to the wall material. The highest EE value was obtained for sample SJ1 (64.09% ± 0.09). Similarly, high encapsulation efficiency was obtained for samples SB3 and GJ2 (61.92 ± 0.04 and 59.89 ± 0.01, respectively). Thus, no clear influence of the mixing ratio between gelatin and gum Arabic was demonstrated on the obtained efficiency of the encapsulation process. The powders obtained because of lyophilization of the coacervates with a 1:1 mixing ratio were characterized by the smallest particle sizes. However, this did not reflect an increase in bulk and tapped density, and did not reflect an increase in their hygroscopicity. However, the powders containing G/GA = 1:1 were the most soluble among the powders tested. The least dissolving microcapsules were obtained with a mixing ratio of G/GA = 2:1.
The aim of this research was to determine whether the classical model of complex coacervation (G: GA) could be used for microencapsulation of essential oils. These findings suggest that such a combination of wall materials can be used for this purpose. Further research should be carried out to find the optimal mixing ratio between gelatin and gum Arabic that would allow for higher encapsulation efficiency. Future investigations may also focus on optimizing all process parameters for the microencapsulation of essential oils to obtain powders with better properties.

Author Contributions

A.N.—conceptualization, methodology, formal analysis, writing—original draft, investigation, A.S.—investigation, I.W.-K.—investigation, M.D.T.P.—investigation, H.D.G.—investigation, M.A.K.—supervision, writing—review and editing. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Data is contained within the article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kunicka-Styczyńska, A. Olejki eteryczne jako alternatywa dla syntetycznych konserwantów żywności—Praca przeglądowa. In Innowacyjne Rozwiązania w Technologii Żywności i Żywieniu Człowieka; Tarko, T., Drożdż, I., Najgebauer-Lejko, D., Duda-Chodak, A., Eds.; Oddział Małopolski Polskiego Towarzystwa Technologów Żywności: Kraków, Poland, 2016; pp. 175–184. [Google Scholar]
  2. Tajkarimi, M.M.; Ibrahim, S.A.; Cliver, D.O. Antimicrobial herb and spice compounds in food. Food Control 2010, 21, 1199–1218. [Google Scholar] [CrossRef]
  3. Maruyama, S.; Streletskaya, N.A.; Lim, J. Clean label: Why this ingredient but not that one? Food Qual. Prefer. 2020, 87, 1–9. [Google Scholar] [CrossRef]
  4. Hojjati, F.; Sereshti, H.; Hojjati, M. Leaf essential oils and their application in systematics of Juniperus excelsa complex in Iran. Biochem. Syst. Ecol. 2019, 84, 29–34. [Google Scholar] [CrossRef]
  5. Ghorbanzadeh, A.; Ghasemnezhad, A.; Sarmast, M.K.; Ebrahimi, S.N. An analysis of variations in morphological characteristics, essential oil content, and genetic sequencing among and within major Iranian Juniper (Juniperus spp.) populations. Phytochem 2021, 186, 1–10. [Google Scholar] [CrossRef] [PubMed]
  6. Nikolić, M.; Stojković, D.; Glamočlija, J.; Ćirić, A.; Marković, T.; Smiljković, M.; Soković, M. Could essential oils of green and black pepper be used as food preservatives? J. Food Sci. Technol. 2015, 52, 1–9. [Google Scholar] [CrossRef]
  7. Zheljazkov, D.V.; Kacaniova, M.; Dincheva, I.; Radoukova, T.; Semerdjieva, I.B.; Astatkie, T.; Schlegel, V. Essential oil composition, antioxidant and antimicrobial activity of the galbuli of six juniper species. Ind. Crops Prod. 2018, 124, 449–458. [Google Scholar] [CrossRef]
  8. Dosoky, N.S.; Satyal, P.; Barata, L.M.; da Silva, J.K.R.; Setzer, W.N. Volatiles of Black Pepper Fruits (Piper nigrum L.). Molecules 2019, 24, 4244. [Google Scholar] [CrossRef]
  9. Bastos, L.P.H.; de Sá Costa, B.; Siqueira, R.P.; Garcia-Rojas, E.E. Complex coacervates of β-lactoglobulin/sodium alginate for the microencapsulation of black pepper (Piper nigrum L.) essential oil: Simulated gastrointestinal conditions and modeling release kinetics. Int. J. Biol. Macromol. 2020, 160, 861–870. [Google Scholar] [CrossRef]
  10. Bastos, L.P.H.; Corrêa dos Santos, C.H.; de Carvalho, M.G.; Garcia-Rojas, E.E. Encapsulation of the black pepper (Piper nigrum L.) essential oil by lactoferrin-sodium alginate complex coacervates: Structural characterization and simulated gastrointestinal conditions. Food Chem. 2020, 316, 861–870. [Google Scholar] [CrossRef]
  11. Amalraj, A.; Haponiuk, J.T.; Thomas, S.; Gopi, S. Preparation, characterization and antimicrobial activity of polyvinyl alcohol/gum arabic/chitosan composite films incorporated with black pepper essential oil and ginger essential oil. Int. J. Biol. Macromol. 2020, 151, 366–375. [Google Scholar] [CrossRef]
  12. Fernandes, R.V.d.B.; Silva, E.K.; Borges, S.V.; de Oliveira, C.R.; Yoshida, M.I.; da Silva, Y.F.; do Carmo, E.L.; Azevedo, V.M.; Botrel, D.A. Proposing Novel Encapsulating Matrices for Spray-Dried Ginger Essential Oil from the Whey Protein Isolate-Inulin/Maltodextrin Blends. Food Bioproc. Technol. 2017, 10, 115–130. [Google Scholar] [CrossRef]
  13. Napiórkowska, A.; Kurek, M. Coacervation as a Novel Method of Microencapsulation of Essential Oils—A Review. Molecules 2022, 27, 5142. [Google Scholar] [CrossRef] [PubMed]
  14. Bakry, A.M.; Shabbar, A.; Barkat, A.; Majeed, H.; Abouelwafa, M.Y.; Mousa, A.; Li, L. Microencapsulation of Oils: A Comprehensive Review of Benefits, Techniques, and Applications. CRFSFS 2016, 15, 143–182. [Google Scholar] [CrossRef] [PubMed]
  15. Delshadi, R.; Bahrami, A.; Tafti, A.G.; Barba, F.J.; Williams, L.L. Micro and nano-encapsulation of vegetable and essential oils to develop functional food products with improved nutritional profiles. Trends Food. Sci. Technol. 2020, 104, 72–83. [Google Scholar] [CrossRef]
  16. Timilsena, Y.P.; Taiwo, O.A.; Nauman, K.; Benu, A.; Colin, J.B. Complex coacervation: Principles, mechanisms and applications in microencapsulation. Int. J. Biol. Macromol. 2019, 121, 1276–1286. [Google Scholar] [CrossRef]
  17. Kralovec, J.; Zhang, S.; Zhang, W.; Barrow, C. A review of the progress in enzymatic concentration and microencapsulation of omega-3 rich oil from fish and microbial sources. Food Chem. 2020, 131, 639–644. [Google Scholar] [CrossRef]
  18. Ogilvie-Battersby, J.D.; Nagarajan, R.; Mosurkal, R.; Orbey, N. Microencapsulation and controlled release of insect repellent geraniol in gelatin/gum arabic microcapsules. Colloids Surf. A Physicochem. Eng. 2022, 640, 1–11. [Google Scholar] [CrossRef]
  19. Rojas-Moreno, S.; Osorio-Revilla, G.; Gallardo-Velázquez, T.; Cárdenas-Bailón, F.; Meza-Márquez, G. Effect of the cross-linking agent and drying method on encapsulation efficiency of orange essential oil by complex coacervation using whey protein isolate with different polysaccharides. J. Microencapsul. 2018, 35, 165–180. [Google Scholar] [CrossRef]
  20. Pakzad, H.; Alemzadeh, I.; Kazemi, A. Encapsulation of Peppermint Oil with Arabic Gum-gelatin by Complex Coacervation Method. Int. J. Eng. 2013, 26, 807–814. [Google Scholar] [CrossRef]
  21. Devi, N.; Sarmah, M.; Khatun, B.; Tarun, K.M. Encapsulation of active ingredients in polysaccharide-protein complex coacervates. Adv. Colloid Interface Sci. 2017, 239, 136–145. [Google Scholar] [CrossRef]
  22. Warnakulasuriya, S.N.; Nickerson, M.T. Review on plant protein-polysaccharide complex coacervation, and the functionality and applicability of formed complexes. J. Sci. Food Agric. 2018, 98, 5559–5571. [Google Scholar] [CrossRef] [PubMed]
  23. Bordón, M.G.; Paredes, A.J.; Camacho, N.M.; Penci, M.C.; González, A.; Palma, S.D.; Ribotta, P.D.; Martinez, M.L. Formulation, spray-drying and physicochemical characterization of functional powders loaded with chia seed oil and prepared by complex coacervation. J. Powder Technol. 2021, 391, 479–493. [Google Scholar] [CrossRef]
  24. Mehran, M.; Masoum, S.; Memarzadeh, M. Microencapsulation of Mentha spicata essential oil by spray drying: Optimization, characterization, release kinetics of essential oil from microcapsules in food models. Ind. Crops Prod. 2020, 154, 1–8. [Google Scholar] [CrossRef]
  25. Xin, X.; Essien, S.; Dell, K.; Woo, M.W.; Baroutian, S. Effects of Spray-Drying and Freeze-Drying on Bioactive and Volatile Compounds of Smoke Powder Food Flavouring. Food Bioproc. Technol. 2022, 15, 785–794. [Google Scholar] [CrossRef]
  26. Hernandez-Nava, R.; Lopez-Malo, A.; Palou, E.; Ramírez-Corona, N.; Jimenez-Munguía, M.T. Encapsulation of oregano essential oil (Origanum vulgare) by complex coacervation between gelatin and chia mucilage and its properties after spray drying. Food Hydrocoll. 2020, 109, 1–36. [Google Scholar] [CrossRef]
  27. Caparino, O.A.; Tang, J.; Nindo, C.I.; Sablani, S.S.; Powers, J.R.; Fellman, J.K. Effect of drying methods on the physical properties and microstructures of mango (Philippine ‘Carabao’ var.) powder. J. Food Eng. 2012, 111, 135–148. [Google Scholar] [CrossRef]
  28. Akseli, I.; Hilden, J.; Katz, J.M.; Kelly, R.C.; Kramer, T.T.; Mao, C.; Osei-Yeboah, F.; Strong, J.C. Reproducibility of the Measurement of Bulk/Tapped Density of Pharmaceutical Powders Between Pharmaceutical Laboratories. J. Pharm. Sci. 2019, 108, 1081–1084. [Google Scholar] [CrossRef]
  29. Kurek, M.A.; Moczkowska-Wyrwisz, M.; Wyrwisz, J.; Karp, S. Development of Gluten-Free Muffins with β-Glucan and Pomegranate Powder Using Response Surface Methodology. Foods 2021, 10, 2551. [Google Scholar] [CrossRef]
  30. Shaddel, R.; Hesari, J.; Azadmard-Damirchi, S.; Hamishehkar, H.; Fathi-Achachlouei, B.; Huang, Q. Use of gelatin and gum Arabic for encapsulation of black raspberry anthocyanins by complex coacervation. Int. J. Biol. Macromol. 2018, 107, 1800–1810. [Google Scholar] [CrossRef]
  31. Pieczykolan, E.; Kurek, M.A. Use of guar gum, gum arabic, pectin, beta-glucan and inulin for microencapsulation of anthocyanins from chokeberry. Int. J. Biol. Macromol. 2019, 129, 665–671. [Google Scholar] [CrossRef]
  32. Fernandes, R.V.d.B.; Borges, S.V.; Botrel, D.A.; Oliveira, C. Physical and chemical properties of encapsulated rosemary essential oil by spray drying using whey protein–inulin blends as carriers. Int. J. Food Sci. 2014, 49, 1–8. [Google Scholar] [CrossRef]
  33. Zhao, Y.; Khalid, N.; Shu, G.; Neves, M.A.; Kobayashi, I.; Nakajima, M. Complex Coacervates from Gelatin and Octenyl Succinic Anhydride Modified Kudzu Starch: Insights of Formulation and Characterization. Food Hydrocoll. 2019, 86, 70–77. [Google Scholar] [CrossRef]
  34. Wojtasik-Kalinowska, I.; Guzek, D.; Górska-Horczyczak, E.; Brodowska, M.; Sun, D.W.; Wierzbicka, A. Diet with linseed oil and organic selenium yields low n-6/n-3 ratio pork Semimembranosus meat with unchanged volatile compound profiles. Int. J. Food Sci. 2018, 53, 1838–1846. [Google Scholar] [CrossRef]
  35. Górska-Horczyczak, E.; Wojtasik-Kalinowska, I.; Guzek, D.; Sun, D.W.; Wierzbicka, A. Differentiation of chill-stored and frozen pork necks using electronic nose with ultra-fast gas chromatography. J. Food Process Eng. 2017, 40, 1–9. [Google Scholar] [CrossRef]
  36. Yuan, Y.; Li, M.F.; Chen, W.S.; Zeng, Q.Z.; Su, D.X.; Tian, B.; He, S. Microencapsulation of shiitake (Lentinula edodes) essential oil by complex coacervation: Formation, rheological property, oxidative stability and odour attenuation effect. Int. J. Food Sci. 2018, 53, 1681–1688. [Google Scholar] [CrossRef]
  37. Klemmer, K.J.; Waldner, L.; Stone, A.; Low, N.H.; Nickerson, M.T. Complex coacervation of pea protein isolate and alginate polysaccharides. Food Chem. 2012, 130, 710–715. [Google Scholar] [CrossRef]
  38. Rohman, A.; Windarsih, A.; Erwanto, Y.; Zakaria, Z. Review on analytical methods for analysis of porcine gelatine in food and pharmaceutical products for halal authentication. Trends Food Sci. Technol. 2020, 101, 122–132. [Google Scholar] [CrossRef]
  39. Atgié, M.; Garrigues, J.C.; Chennevière, A.; Masbernat, O.; Roger, K. Gum Arabic in solution: Composition and multi-scale structures. Food Hydrocoll. 2019, 91, 319–330. [Google Scholar] [CrossRef]
  40. Muhoza, B.; Xia, S.; Wang, X.; Zhang, X.; Li, Y.; Zhang, S. Microencapsulation of essential oils by complex coacervation method: Preparation, thermal stability, release properties and applications. Crit. Rev. Food Sci. Nutr. 2020, 62, 1363–1382. [Google Scholar] [CrossRef]
  41. Kavoosi, G.; Rahmatollahi, A.; Dadfar, S.M.M.; Purfard, A.M. Effects of essential oil on the water binding capacity, physico-mechanical properties, antioxidant and antibacterial activity of gelatin films. LWT 2014, 57, 556–561. [Google Scholar] [CrossRef]
  42. de Melo Ramos, F.; Silveira Júnior, V.; Prata, A.S. Assessing the Vacuum Spray Drying Effects on the Properties of Orange Essential Oil Microparticles. Food Bioproc. Technol. 2019, 12, 1917–1927. [Google Scholar] [CrossRef]
  43. Shah, R.B.; Tawakkul, M.A.; Khan, M.A. Comparative Evaluation of Flow for Pharmaceutical Powders and Granules. AAPS Pharmscitech 2008, 9, 250–258. [Google Scholar] [CrossRef] [PubMed]
  44. Manaf, M.A.; Subuki, I.; Jai, J.; Raslan, R.; Mustapa, A.N. Encapsulation of Volatile Citronella Essential Oil by Coacervation: Efficiency and Release Study. In IOP Conference Series: Materials Science and Engineering, 358, Proceedings of the 3rd International Conference on Global Sustainability and Chemical Engineering (ICGSCE) Putrajaya, Malaysia, 15–16 February 2017; IOP Publishing: Bristol, UK, 2018. [Google Scholar]
  45. Samakradhamrongthai, R.S.; Angeli, P.T.; Kopermsub, P.; Utama-ang, N. Optimization of gelatin and gum arabic capsule infused with pandan flavor for multi-core flavor powder encapsulation. Carbohydr. Polym. 2019, 226, 2584–2591. [Google Scholar] [CrossRef] [PubMed]
  46. Mitra, H.; Pushpadass, H.A.; Franklin, M.E.E.; Ambrose, R.K.; Ghoroi, C.; Battula, S. Influence of moisture content on the flow properties of basundi mix. J. Powder Technol. 2017, 312, 133–143. [Google Scholar] [CrossRef]
  47. Zotarelli, M.F.; da Silva, V.M.; Durigon, A.; Hubinger, M.D.; Laurindo, J.B. Production of mango powder by spray drying and cast-tape drying. Powder Technol. 2017, 305, 447–454. [Google Scholar] [CrossRef]
  48. Drozińska, E.; Kanclerz, A.; Kurek, M.A. Microencapsulation of sea buckthorn oil with β-glucan from barley as coating material. Int. J. Biol. Macromol. 2019, 131, 1014–1020. [Google Scholar] [CrossRef] [PubMed]
  49. Walsh, M.K. Immobilized enzyme technology for food applications. In Novel Enzyme Technology for Food Applications; Rastall, R., Ed.; Woodhead Publishing: Sawston, UK, 2007; pp. 60–84. [Google Scholar] [CrossRef]
  50. Kapcsándi, V.; Lakatos, E.H.; Sik, B.; Linka, L.A.; Székelyhidi, R. Characterization of fatty acid, antioxidant, and polyphenol content of grape seed oil from different Vitis vinifera L. varieties. OCL 2021, 28, 1–6. [Google Scholar] [CrossRef]
  51. Cui, S.W.; Phillips, G.O.; Blackwell, B.; Nikiforuk, J. Characterisation and properties of Acacia senegal (L.) Willd. var. senegal with enhanced properties (Acacia (sen) SUPERGUM™): Part 4. Spectroscopic characterisation of Acacia senegal var. senegal and Acacia (sen) SUPERGUM™ Arabic. Food Hydrocoll. 2007, 21, 347–352. [Google Scholar] [CrossRef]
  52. Nandiyanto, A.B.D.; Oktiani, R.; Ragadhita, R. How to read and interpret FTIR spectroscope of organic material. Indones. J. Sci. Technol. 2019, 4, 97–118. [Google Scholar] [CrossRef]
  53. Musa, H.H.; Ahmed, A.A.; Musa, T.H. Chemistry, Biological, and Pharmacological Properties of Gum Arabic. In Bioactive Molecules in Food. Reference Series in Phytochemistry; Mérillon, J.M., Ramawat, K.G., Eds.; Springer: Cham, Switzerland, 2019; pp. 797–814. [Google Scholar] [CrossRef]
  54. Marković, M.S.; Radosavljević, D.B.; Pavićević, V.P.; Ristić, M.S.; Milojević, S.Z.; Bošković-Vragolović, N.M.; Veljković, V.B. Influence of common juniper berries pretreatment on the essential oil yield, chemical composition and extraction kinetics of classical and microwave-assisted hydrodistillation. Ind. Crops Prod. 2018, 122, 402–413. [Google Scholar] [CrossRef]
  55. Bastos, L.P.H.; Vicente, J.; Corrêa dos Santos, C.H.; Geraldo de Carvalho, M.; Garcia-Rojas, E.E. Encapsulation of black pepper (Piper nigrum L.) essential oil with gelatin and sodium alginate by complex coacervation. Food Hydrocoll. 2019, 102, 1–8. [Google Scholar] [CrossRef]
  56. Roohinejad, S.; Koubaa, M.; Barba, F.J.; Leong, S.Y.; Khalefa, A.; Greiner, R.; Chemat, F. Extraction methods of essential oils from herbs and spices. In Essential Oils in Food Processing: Chemistry, Safety and Applications; Hashemi, S.M.B., Khaneghah, A.M., de Souza Sant’Ana, A., Eds.; John Wiley & Sons Ltd.: Hoboken, NJ, USA, 2017; pp. 21–55. [Google Scholar] [CrossRef]
  57. Karagozlu, M.; Ocak, B.; Özdestan-Ocak, Ö. Effect of Tannic Acid Concentration on the Physicochemical, Thermal, and Antioxidant Properties of Gelatin/Gum Arabic–Walled Microcapsules Containing Origanum onites L. Essential Oil. Food Bioproc. Technol. 2021, 14, 1231–1243. [Google Scholar] [CrossRef]
  58. Yekdane, N.; Goli, S.A.H. Effect of Pomegranate Juice on Characteristics and Oxidative Stability of Microencapsulated Pomegranate Seed Oil Using Spray Drying. Food Bioproc. Technol. 2019, 12, 1614–1625. [Google Scholar] [CrossRef]
  59. Koç, M.; Güngör, Ö.; Zungur, A.; Yalçın, B.; Selek, I.; Ertekin, F.K.; Ötles, S. Microencapsulation of Extra Virgin Olive Oil by Spray Drying: Effect of Wall Materials Composition, Process Conditions, and Emulsification Method. Food Bioproc. Technol. 2015, 8, 301–318. [Google Scholar] [CrossRef]
  60. Rafiq, S.; Sofi, S.A.; Kumar, H.; Kaul, R.K.; Mehra, R.; Awuchi, C.G.; Okpala, C.O.R.; Korzeniowska, M. Physicochemical, antioxidant, and polyphenolic attributes of microencapsulated freeze-dried kinnow peel extract powder using maltodextrin as wall material. J. Food Process. Preserv. 2022, 46, 1–10. [Google Scholar] [CrossRef]
  61. Gomez-Estaca, J.; Comunian, T.A.; Montero, P.; Favaro-Trindade, C.S. Physico-Chemical Properties, Stability, and Potential Food Applications of Shrimp Lipid Extract Encapsulated by Complex Coacervation. Food Bioproc. Technol. 2018, 11, 1596–1604. [Google Scholar] [CrossRef]
  62. Tavares, L.; Norena, C.P.Z. Encapsulation of Ginger Essential Oil Using Complex Coacervation Method: Coacervate Formation, Rheological Property, and Physicochemical Characterization. Food Bioproc. Technol. 2020, 13, 1405–1420. [Google Scholar] [CrossRef]
  63. Syed, N.; Mahesar, S.A.; Sherazi, S.T.H.; Soylak, M. Quality assessment and safety measurement of different industrial processing stages of soybean oil. TURJFAS 2020, 1, 28–33. [Google Scholar] [CrossRef]
  64. Bruhl, L.; Unbehemd, G. Precise Color Communication by Determination of the Color of Vegetable Oils and Fats in the CIELAB 1976 (L*a*b*) Color Space. Eur. J. Lipid Sci. Technol. 2021, 123, 1–9. [Google Scholar] [CrossRef]
  65. Sanchez, C.; Nige, M.; Meji, T.; Doco, T.; Williams, P.; Amine, C.; Renard, D. Acacia gum: History of the future. Food Hydrocoll. 2018, 78, 140–160. [Google Scholar] [CrossRef]
  66. Ocak, B. Gum arabic and collagen hydrolysate extracted from hide fleshing wastes as novel wall materials for microencapsulation of Origanum onites L. essential oil through complex coacervation. Environ. Sci. Pollut. Res. 2020, 27, 42727–42737. [Google Scholar] [CrossRef] [PubMed]
  67. Marfil, P.H.M.; Paulo, B.B.; Alvim, I.D.; Nicoletti, R. Production and characterization of palm oil microcapsules obtained by complex coacervation in gelatin/gum Arabic. J. Food Process Eng. 2018, 41, 1–11. [Google Scholar] [CrossRef]
  68. Tavares, L.; Barros, H.L.; Barbosa, J.C.P.; Vaghetti, C.P.Z. Noreña, Microencapsulation of Garlic Extract by Complex Coacervation Using Whey Protein Isolate/Chitosan and Gum Arabic/Chitosan as Wall Materials: Influence of Anionic Biopolymers on the Physicochemical and Structural Properties of Microparticles. Food Bioproc. Technol. 2019, 12, 2093–2106. [Google Scholar] [CrossRef]
  69. Krokida, M.K.; Karathanos, V.T.; Maroulis, Z.B. Effect of freeze-drying conditions on shrinkage and porosity of dehydrated agricultural products. J. Food Eng. 1998, 35, 369–380. [Google Scholar] [CrossRef]
  70. Barresi, A.A.; Ghio, S.; Fissore, D.; Pisano, R. Freeze Drying of Pharmaceutical Excipients Close to Collapse Temperature: Influence of the Process Conditions on Process Time and Product Quality. Dry Technol. 2009, 27, 805–816. [Google Scholar] [CrossRef]
Figure 1. FT-IR spectra for samples GJ1, GJ2, GJ3, GB1, GB2 and GB3.
Figure 1. FT-IR spectra for samples GJ1, GJ2, GJ3, GB1, GB2 and GB3.
Foods 12 04345 g001
Figure 2. FT-IR spectra for samples SJ1, SJ2, SJ3, SB1, SB2 and SB3.
Figure 2. FT-IR spectra for samples SJ1, SJ2, SJ3, SB1, SB2 and SB3.
Foods 12 04345 g002
Figure 3. Classification of scent profiles of microcapsules containing juniper and black pepper essential oils.
Figure 3. Classification of scent profiles of microcapsules containing juniper and black pepper essential oils.
Foods 12 04345 g003
Figure 4. SEM photos.
Figure 4. SEM photos.
Foods 12 04345 g004
Table 1. Coding of samples.
Table 1. Coding of samples.
OilEssential OilMixing Ratio G/GACode
Grape seedJuniper1:1GJ1
1:2GJ2
2:1GJ3
Black pepper1:1GB1
1:2GB2
2:1GB3
SoybeanJuniper1:1SJ1
1:2SJ2
2:1SJ3
Black pepper1:1SB1
1:2SB2
2:1SB3
Table 2. Complex coacervation yield (CY), solid yield (SY), and encapsulation efficiency (EE) values [%].
Table 2. Complex coacervation yield (CY), solid yield (SY), and encapsulation efficiency (EE) values [%].
SampleCY%SY%EE%
GJ130.93 ± 0.44 b29.79 ± 1.21 c49.3 ± 0.07 abe
GJ228.01 ± 0.16 a37.11 ± 0.49 e59.89 ± 0.01 cdf
GJ339.22 ± 0.16 g21.32 ± 0.25 a49.65 ± 0.02 abe
GB130.49 ± 0.56 b29.25 ± 1.29 c42.7 ± 0.11 e
GB229.45 ± 0.08 d36.23 ± 0.08 de55.25 ± 0.01 abcd
GB341.64 ± 0.15 i23.01 ± 0.36 b47.21 ± 0.04 ae
SJ133.78 ± 0.46 f26.58 ± 0.36 g64.09 ± 0.09 f
SJ231.91 ± 0.12 e33.8 ± 0.02 h53.4 ± 0.04 abcd
SJ346.58 ± 0.19 j22.02 ± 0.18 ab56.8 ± 0.01 bcdf
SB123.95 ± 0.16 c35.44 ± 0.47 d54.14 ± 0.05 abcd
SB228.34 ± 0.12 a38.45 ± 0.34 i52.33 ± 0.01 abc
SB341.04 ± 0.51 h20.21 ± 0.61 f61.92 ± 0.04 df
S.E.M0.3690.09626.524
EffectOil**NS**
Essential oil****NS
Mixing ratio****NS
Oil × essential oil****NS
Essential oil × mixing ratio******
Oil × mixing ratio****NS
Oil × essential oil × mixing ratio****NS
Results in this table are expressed as mean ± standard deviation. Mean values with the same superscript letters within a row are not significantly different at p ≤ 0.05. S.E.M.—standard error of the mean. * * = p 0.001 ; NS non - significant   effect = p > 0.05 .
Table 3. Bulk density, tapped density, Carr index, Hausner ratio [g/cm3], and solubility, hygroscopicity and moisture content (g/100 g).
Table 3. Bulk density, tapped density, Carr index, Hausner ratio [g/cm3], and solubility, hygroscopicity and moisture content (g/100 g).
SampleBulk Density g/cm3Tapped Density g/cm3Carr Index
g/cm3
Hausner Ratio g/cm3Solubility
g/100 g
Hygroscopicity g/100 gMoisture Content g/100 g
GJ10.12 ± 0.01 c0.24 ± 0.00 abc50.27 ± 2.53 d2.01 ± 0.11 d19.4 ± 2.92 be2.53 ± 0.32 a0.27 ± 0.04 a
GJ20.16 ± 0.01 a0.23 ± 0.01 ab30.58 ± 4.93 a1.45 ± 0.10 a12.54 ± 3.76 acd7.29 ± 8.02 a0.18 ± 0.12 a
GJ30.16 ± 0.01 ab0.24 ± 0.01 abc31.81 ± 4.45 a1.47 ± 0.09 ab7.85 ± 3.54 a2.87 ± 2.15 a0.56 ± 0.62 a
GB10.11 ± 0.01 c0.22 ± 0.02 ab48.84 ± 7.73 d1.98 ± 0.28 d18.22 ± 3.35 bd3.77 ± 2.51 a0.29 ± 0.42 a
GB20.17 ± 0.01 b0.26 ± 0.00 de33.56 ± 3.46 ab1.51 ± 0.08 ab20.7 ± 4.15 be6.98 ± 3.85 a0.05 ± 0.05 a
GB30.16 ± 0.01 a0.24 ± 0.01 abc34.64 ± 2.13 ab1.53 ± 0.05 ab6.58 ± 2.69 a6.00 ± 4.84 a0.36 ± 0.29 a
SJ10.12 ± 0.00 c0.22 ± 0.03 a44.35 ± 6.21 cd1.81 ± 0.21 cd16.39 ± 4.46 bcd4.92 ± 4.18 a0.30 ± 0.21 a
SJ20.16 ± 0.00 ab0.25 ± 0.00 cd34.74 ± 0.64 ab1.53 ± 0.02 ab18.72 ± 5.65 bd4.32 ± 3.06 a0.44 ± 0.66 a
SJ30.17 ± 0.01 ab0.28 ± 0.00 e40.28 ± 4.6 bc1.68 ± 0.13 bc11.14 ± 3.01 ac2.71 ± 1.97 a0.36 ± 0.30 a
SB10.14 ± 0.02 d0.23 ± 0.01 ab38.05 ± 8.82 abc1.64 ± 0.22 abc10.74 ± 2.41 ac2.63 ± 0.17 a0.24 ± 0.25 a
SB20.16 ± 0.00 ab0.25 ± 0.01 bcd34.33 ± 2.05 ab1.52 ± 0.05 ab25.86 ± 6.81 e2.86 ± 0.82 a0.28 ± 0.24 a
SB30.15 ± 0.01 ad0.24 ± 0.01 abcd37.75 ± 2.98 abc1.61 ± 0.08 abc7.47 ± 1.34 a1.48 ± 0.48 a0.49 ± 0.08 a
S.E.M0.0000.00023.1520.02015.46612.1050.113
EffectOilNSNSNSNSNSNSNS
Essential oilNSNSNSNSNSNSNS
Mixing ratio*********NSNS
Oil × Essential oilNS*NSNSNSNSNS
Essential oil × mixing ratio*****NSNS
Oil × mixing ratioNS*NSNS*NSNS
Oil × essential oil × mixing ratio**NSNSNSNSNS
Results in this table are expressed as mean ± standard deviation. Mean values with the same superscript letters within a row are not significantly different at p ≤ 0.05. S.E.M.—standard error of the mean. * = p 0.05 ;   * * = p 0.001 ; NS non - significant   effect = p > 0.05 .
Table 4. Color parameters of obtained powders.
Table 4. Color parameters of obtained powders.
SampleL*a*b*
GJ186.99 ± 0.32 cd−1.74 ± 0.14 h6.3 ± 0.19 bc
GJ288.52 ± 0.25 g−2.09 ± 0.02 g6.56 ± 0.13 c
GJ384.46 ± 0.97 a−3.05 ± 0.08 e9.47 ± 0.37 d
GB187.08 ± 0.23 c−1.19 ± 0.02 i5.61 ± 0.04 b
GB287.00 ± 0.11 c−2.91 ± 0.01 e9.24 ± 0.02 d
GB384.10 ± 0.87 a−3.83 ± 0.15 f12.30 ± 0.47 f
SJ184.52 ± 0.05 a−7.66 ± 0.05 c46.52 ± 0.13 ae
SJ285.64 ± 0.33 b−7.38 ± 0.14 ab45.82 ± 0.31 e
SJ378.72 ± 0.34 d−7.06 ± 0.27 d51.71 ± 0.47 h
SB186.07 ± 1.03 bc−7.37 ± 0.21 ab41.87 ± 1.05 g
SB285.70 ± 0.30 b−7.49 ± 0.11 ac46.86 ± 0.38 a
SB382.29 ± 0.56 f−7.2 ± 0.07 bd47.06 ± 0.64 a
S.E.M0.3020.0170.201
EffectOil******
Essential oil*****
Mixing ratio******
Oil × essential oil******
Essential oil × mixing ratio****NS
Oil × mixing ratio******
Oil × essential oil × mixing ratio*****
Results in this table are expressed as mean  ±  standard deviation. Mean values with the same superscript letters within a row are not significantly different at p  ≤  0.05. S.E.M.—standard error of the mean. * = p 0.05 ;   * * = p 0.001 ;   NS non - significant   effect = p > 0.05 .
Table 5. Particle size and particle size distribution index—span (SI) [μm].
Table 5. Particle size and particle size distribution index—span (SI) [μm].
SampleParticle Size μmD10 μmD50 μmD90 μmSI μm
GJ120.48 ± 0.16 h56.79 ± 0.62 b83.42 ± 0.36 b96.84 ± 0.35 a0.48 ± 0.01 b
GJ29.18 ± 0.14 f53.66 ± 0.54 h81.98 ± 0.64 a95.55 ± 0.41 b0.47 ± 0.01 b
GJ34.51 ± 0.18 d51.9 ± 0.4 g81.7 ± 0.98 a95.37 ± 0.36 b0.96 ± 0.01 g
GB15.55 ± 0.29 a45.86 ± 0.33 f82.79 ± 0.75 ab95.7 ± 0.34 b0.6 ± 0.00 e
GB210.5 ± 0.25 c60.13 ± 0.74 c85.82 ± 1.46 d96.86 ± 0.41 a0.44 ± 0.01 a
GB312.74 ± 0.34 g43.14 ± 1.01 e77.77 ± 0.18 i95.38 ± 0.26 b1.17 ± 0.01 h
SJ15.57 ± 0.35 a22.2 ± 0.25 d72.98 ± 0.79 h92 ± 0.18 d0.53 ± 0.01 d
SJ26.67 ± 0.29 b57.59 ± 0.88 b84.11 ± 0.29 bc96.78 ± 0.55 a0.51 ± 0.01 c
SJ36.48 ± 0.46 b17.42 ± 0.5 a51.6 ± 1.16 e90.01 ± 0.14 b1.41 ± 0.02 j
SB15.71 ± 0.38 a17.37 ± 0.46 a68.25 ± 0.58 g97.24 ± 0.1 a0.67 ± 0.01 f
SB210.61 ± 0.35 c59.04 ± 1.18 c84.82 ± 0.22 cd96.72 ± 0.5 a0.43 ± 0.01 a
SB37.59 ± 0.34 e16.69 ± 0.25 a60.77 ± 0.94 f92.67 ± 0.14 e1.25 ± 0.02 i
S.E.M0.0940.4370.6300.1180.000
EffectOil**********
Essential oilNS*******
Mixing ratio**********
Oil × essential oil**********
Essential oil × mixing ratio**********
Oil × mixing ratio**********
Oil × essential oil × mixing ratio**********
Results in this table are expressed as mean  ±  standard deviation. Mean values with the same superscript letters within a row are not significantly different at p  ≤  0.05. S.E.M.—standard error of the mean. * = p 0.05 ;   * * = p 0.001 ;   NS non - significant   effect = p > 0.05 .
Table 6. Temperatures of the onsets, peaks and endsets of endothermic reactions with enthalpy.
Table 6. Temperatures of the onsets, peaks and endsets of endothermic reactions with enthalpy.
SampleOnset (°C)Peak (°C)Endset (°C)Enthalpy (mJ)
GJ163.67 ± 0.00276.61 ± 0.00187.18 ± 0.002−36.31 ± 0.001
GJ256.93 ± 0.00298.08 ± 0.001161.87 ± 0.002−353.88 ± 0.001
GJ368.75 ± 0.00179.31 ± 0.00287.94 ± 0.002−20.52 ± 0.002
GB146.21 ± 0.00198.60 ± 0.001113.29 ± 0.001−222.81 ± 0.001
GB283.46 ± 0.00192.76 ± 0.00299.21 ± 0.001−14.00 ± 0.001
GB394.72 ± 0.001121.90 ± 0.001151.31 ± 0.002−327.77 ± 0.001
SJ155.23 ± 0.001107.75 ± 0.001175.95 ± 0.001−186.11 ± 0.001
SJ282.32 ± 0.00189.61 ± 0.001108.89 ± 0.001−51.25 ± 0.001
SJ368.86 ± 0.00180.62 ± 0.00189.36 ± 0.001−0.45 ± 0.002
SB159.94 ± 0.001109.09 ± 0.002178.75 ± 0.001−285.82 ± 0.002
SB270.13 ± 0.00178.48 ± 0.00284.80 ± 0.001−14.50 ± 0.001
SB371.79 ± 0.001121.62 ± 0.002161.17 ± 0.001−210.46 ± 0.001
Table 7. Kovats indexes, identified volatile compounds and sensory descriptors assigned to them.
Table 7. Kovats indexes, identified volatile compounds and sensory descriptors assigned to them.
Kovats IndexIdentified Volatile CompoundSensory
Descriptors
GJ1GJ2GJ3SJ1SJ2SJ3GB1GB2GB3SB1SB2SB3
437ethanol
2-propanol
ethanethiol
alcoholic++++++++++++
508alcoholic++++++++++++
519earthy, fruity, garlic +
5412-methylpropanal
butanal
butan-2-one
aldehydic++++++++++++
567chocolate, green, malty++++++++++++
594acetone, butter++++++++++++
6082-methylfuran
methyl propanoate
acetic acid
acetone, burnt + + +
610apple, etheral, fresh +
617acetic, acidic + +
6181-propanethiol
1-butanamine
1-propanethiol
alliaceous, cabbage, onion, sween +
638ammoniacal, fish +++
643allioceous, cabbage, onion, sweet +
437ethanolalcoholic++++++++++++
5082-propanolalcoholic++++++++++++
519ethanethiolearthy, fruity, garlic +
5412-methylpropanalaldehydic++++++++++++
567butanalchocolate, green, malty++++++++++++
594butan-2-oneacetone, butter++++++++++++
6082-methylfuranacetone, burnt + + +
610methyl propanoateapple, etheral, fresh +
617acetic acidacetic, acidic + +
6181-propanethiolalliaceous, cabbage, onion, sween +
6381-butanamineammoniacal, fish +++
6431-propanethiolallioceous, cabbage, onion, sweet +
658n-butanolalcoholic, cheese, fermented + ++
667trichloroethanechloroform, etheral++++ +++++++
6731-methoxy-2-propanolmild +++
684pentan-2-oneacetone, banana, etheral +++ +++
700pentan-2-olalcoholic, etheral, fermented+++ +++
700heptanealkane, fruity +
735dimethyl disulfidecabbage, cheese, garlic +++
735butanethiolcoffee, garlic+++
742(e)-2-pentanalapple, fruity, green +
745propanoic acidcheese, fruity+++++++ +++
759pentanolalcoholic, anise, balsamic +
7673-methylbut-2-en-1-olherbaceous, lavender ++
780toluenecaramelized, etheral, fruity+++ ++++++
7672-methylpentane- +++
792hexanalacorn, aldehydic, fatty ++++++++
801ethyl butyrateacetone, banana++++ +
812octanealkane, fruity, fusel ++++ ++
810(e)-2-octene- ++
8474-ethylheptane- +++
8473-methyl-2-butene-1-thiolamine, leek, onion+++
852methyl pentanoateapple, etheral, fruity, green +
866ethyl isovalerateanise, apple, blackcurrant+++ +
8733-methyl-octane- +++
884nonanealkane, fusel++++++
8962-butylfuranfruity, mild +++
898isoamyl acetateapple, banana, ester, fresh++ +++++
907nonanealkane, fusel +
9221s-()-a-pinenefresh, herbaceous++++++++++++
9405-pentanolide- +++
962alpha-pinenecamphore, citrus++++++++++++
974beta-pinenedry, green, hay++++++++++++
966myrcenebalsamic, etheral, fruity++++++++++++
1033beta-phellandrenefruity, herbaceous++++++++++++
1049limonenecitrus, fruity++++++++++++
1076gamma-terpinenecitrus. etheral, fruity++++++++++++
1106methylacetophenone- +++
1107ethyl heptanoatefruity+++ ++++++
1135n-nonanalaldehydic, chlorine, citrus++++++++++++
1150ethyl cyclohexanecarboxylate- +
11532,3-diethyl-5-methylpyrazienfragrant, hazelnut +++
1153e-2-nonen-1-olgreen, melon+++ ++
1164benzyl acetateburnt, floral, fresh, fruity++++++++++
1187p-methylacetophenonealmond, bitter almond, cherry+++++ +++++
1205decanalaldehydic, burnt, citrus+++ +++++
12076-decenal- +++
1217dodecanealkane, fusel ++
1248ethyl phenylacetateanise, cinnamon, cocoa, flaral, rose +++
12652-butenoic acid, hexyl esterfruity, green, oily, walnut
1277tridecanealkane, citrus, fruity
1300ethyl nonanoatefruity, rose, rum +
1313nonyl acetatefruity, leafy, sweet +
1272tridecane- +
1278geranial-++
1293thymolaromatic, earthy, herbaceous++
1301tricdecanealkane, citrus, fruity +
1313butyl heptanoatefresh, fruity, grasssy, green +
1316anetholeanise, herbaceous,++
1330nonyl acetatefruity, leafy, sweet +
13653-ethyl dodecane- +++
1367eugenolbalsamic, camphore, floral ++
1371tetradecanealkane, fusel, herbaceous+++ +++
1386n-hexyl-hexanoateapple, fresh, fruity +++
1369butyl octanoatebutter, floral, fruity, green, oily+++ +++
1421propyl nonanoatefermented, melon++++++ +++
1468n-octylbenzene- +++
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Napiórkowska, A.; Szpicer, A.; Wojtasik-Kalinowska, I.; Perez, M.D.T.; González, H.D.; Kurek, M.A. Microencapsulation of Juniper and Black Pepper Essential Oil Using the Coacervation Method and Its Properties after Freeze-Drying. Foods 2023, 12, 4345. https://doi.org/10.3390/foods12234345

AMA Style

Napiórkowska A, Szpicer A, Wojtasik-Kalinowska I, Perez MDT, González HD, Kurek MA. Microencapsulation of Juniper and Black Pepper Essential Oil Using the Coacervation Method and Its Properties after Freeze-Drying. Foods. 2023; 12(23):4345. https://doi.org/10.3390/foods12234345

Chicago/Turabian Style

Napiórkowska, Alicja, Arkadiusz Szpicer, Iwona Wojtasik-Kalinowska, Maria Dolores Torres Perez, Herminia Dominguez González, and Marcin Andrzej Kurek. 2023. "Microencapsulation of Juniper and Black Pepper Essential Oil Using the Coacervation Method and Its Properties after Freeze-Drying" Foods 12, no. 23: 4345. https://doi.org/10.3390/foods12234345

APA Style

Napiórkowska, A., Szpicer, A., Wojtasik-Kalinowska, I., Perez, M. D. T., González, H. D., & Kurek, M. A. (2023). Microencapsulation of Juniper and Black Pepper Essential Oil Using the Coacervation Method and Its Properties after Freeze-Drying. Foods, 12(23), 4345. https://doi.org/10.3390/foods12234345

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop