DES/Extraction and Process Optimization of Callistemon citrinus (Curtis) Skeels Essential Oil

TL;DR

This study optimizes the extraction of essential oil from Callistemon citrinus using ultrasound and metal ions, improving yield and chemical composition.

Contribution

A novel combination of ultrasonic-assisted hydrodistillation and metal ion pretreatment is introduced to enhance essential oil extraction efficiency and composition.

Findings

Optimized extraction conditions yielded 0.8 mL/kg of essential oil with specific heating and ultrasonic parameters.

Metal ions like Sr²⁺ and Ca²⁺ significantly influenced the relative abundance of monoterpenoid compounds in the extracted oil.

PCA confirmed distinct chemical profiles of oils extracted under different ion pretreatment conditions.

Abstract

The essential oil of Callistemon citrinus (Curtis) Skeels is a valuable natural product with a complex volatile composition. In this study, an ultrasonic-assisted hydrodistillation process was systematically optimized using a combination of the steepest ascent method and response surface methodology (RSM) to refine key extraction parameters and identify the central point of the experimental design. Heating time, heating power, and ultrasonic treatment time were selected as critical variables influencing extraction efficiency. The optimized conditions were determined to be a heating time of 57 min, a heating power of 563 W, and an ultrasonic treatment time of 23 min, under which an essential oil yield of 0.8 mL/kg (dry weight basis) was achieved. Gas chromatography–mass spectrometry (GC–MS) analysis indicated that pretreatment with different metal ions resulted in pronounced differences in the chemical profiles of the extracted essential oils, particularly in major constituents such as cineole and α-terpineol. Among the tested ions, Sr2+ and Ca2+ pretreatments were associated with a higher relative abundance of monoterpenoid compounds, which suggests a selective influence of divalent metal ions on the extraction behavior of volatile components. Principal component analysis (PCA) further revealed clear discrimination among the essential oils obtained under different ion pretreatment conditions, confirming the compositional variability induced by ion-assisted extraction.

Figures

Full text

1. Introduction

Callistemon citrinus is an aromatic species widely distributed in subtropical and temperate regions and has attracted increasing attention due to its rich phytochemical profile and broad biological potential. Extracts and essential oils derived from C. citrinus are known to be abundant in terpenoids and phenolic compounds, which have been associated with multiple bioactivities, including insecticidal, antioxidant, and neuroactive effects. In addition, Callistemon citrinus (Curtis) Skeels (C. citrinus) is widely cultivated as a landscape plant because of its distinctive bottlebrush-like inflorescences [1,2] (Abdelhady and Aly, 2012; Salem et al., 2013a). Previous studies have reported that its leaves contain a wide range of functional compounds [3,4,5]; however, these phytochemical resources are often underutilized and insufficiently explored. Despite increasing reports on its chemical composition and general bioactivities, most existing studies remain largely descriptive, with limited insight into the molecular interactions between individual bioactive constituents and their corresponding biological targets. Consequently, the mechanistic understanding of how specific components from C. citrinus exert their biological effects remains fragmented. Addressing this knowledge gap is, therefore, essential for providing a clearer rationale for species selection and for advancing the targeted functional or biocontrol applications of C. citrinus.

Essential oils extracted from C. citrinus have been reported to have antifungal and insecticidal activity [1,6], and most have been used as antimicrobial and antifungal agents [2,7]. However, although the raw plant material contains a relatively high proportion of essential oils, their direct utilization remains constrained by the lack of efficient and environmentally friendly separation strategies [8], so it is crucial to find an efficient, green extraction method that maximizes the yield and quality of essential oils. Methods commonly used for plant essential oil extraction include steam distillation, microwave-assisted extraction, supercritical fluid extraction, and ultrasound-assisted extraction, et al. [9,10]. Therefore, it is important to develop an efficient and environmentally friendly pretreatment technology to destroy cell walls and improve the extraction efficiency of essential oils. The previous research results of Khursheed et al. showed that metal salt solutions can effectively dissolve cellulose components [11] (Khursheed et al., 2020), and metal ions with different valence states may have different interaction mechanisms with cell wall components due to their differences in charge density and ion radius [12], thus selectively affecting the release of active ingredients [13]; this phenomenon provides a new idea for the directional extraction of plant essential oils through ion regulation.

Deep eutectic solvent (DES) is a new type of green solvent formed by a hydrogen bond acceptor and a hydrogen bond donor through a hydrogen bond interaction [14]. Compared to traditional ionic liquids, deep eutectic solvent has designable physical and chemical properties, low toxicity and good biocompatibility [15]. It is widely used in the extraction of various active ingredients and shows great potential to replace traditional solvents [16]. Based on the above background, an innovative pretreatment strategy was proposed to couple the solvent function of DES with its pH-regulating ability and the cell wall destruction function of metal salt solution. We choose choline chloride and dilute hydrochloric acid or sodium hydroxide to construct a deep eutectic solvent system with different pH values and directly configure different types of metal salt (NaCl, KCl, SrCl_2_, CaCl_2_) pretreatment solutions in this deep eutectic solvent. Based on this, we employed an ultrasound-assisted steam distillation method to extract essential oil from C. citrinus leaves.

The purpose of this study was to systematically optimize the extraction process through single-factor experiments and the response surface method. The chemical composition differences of essential oils obtained by different pretreatments were analyzed by gas chromatography–mass spectrometry (G–MS), and their chemical profile characteristics were revealed by principal component analysis. Finally, the potential inhibitory activity of acetylcholinesterase of essential oil was predicted by molecular docking technology. It provides a theoretical basis and technical path for the directional and efficient extraction of C. citrinus essential oil and its application in functional product development.

2. Results and Discussion

2.1. Single-Factor Analysis

As shown in Figure 1, the yields of essential oils obtained from the treatment with different metal salt solutions are presented. Different inorganic salt solutions were used to treat millephrine within a unit of time, and the relationship between the extraction rate and different inorganic salt solutions was discussed. Through comparison, it was found that the SrCl_2_ pretreatment group had the highest essential oil yield, reaching 0.55 mL/kgDW, followed by CaCl_2_ and KCl, while the NaCl group had the lowest yield. This difference mainly stems from the physicochemical properties of different metal ions and their mechanisms of interaction with plant cell walls: divalent metal ions (such as $[eqn]$ , $[eqn]$ ), due to their high charge density and large ionic radius, can more effectively undergo cross-linking or displacement reactions with components such as pectin and cellulose in the cell wall [12] (Dudev and Lim, 2014), disrupting the integrity of the cell structure and thereby enhancing the release efficiency of essential oils. However, the effect of monovalent ions ( $[eqn]$ , $[eqn]$ ) is relatively weak. The results of this study provide a theoretical basis for further optimizing the essential oil extraction process, indicating that choosing an appropriate metal salt pretreatment agent can significantly improve extraction efficiency, which has particular application potential in the green extraction technology of natural products [17] (Modi et al., 2021).

2.2. Plackett–Burman Design (PBD) Experimental Analysis

PBD is an effective method to screen out key factors from all independent variables. A total of 12 sets of experiments were conducted according to the PBD method. Table 1 shows the essential oil yield and analysis of variance (ANOVA) of the 12 groups of experiments. The experimental data were fitted by multiple regression with Design Expert 13.0 software, and the second-order polynomial modulus of essential oil yield ( $[eqn]$ ) and each coding factor was obtained as follows:

$$[eqn]$$

where $[eqn]$ is the essential oil yield, X1 is pH effect, X2 is ion type, X3 is ion mass, X4 is ultrasonic time, X5 is ultrasonic power, X6 is heating power, and X7 is heating time. From Table 1, it can be seen that the F-value of the model is 6.30 and the p-value is 0.0472, which indicate that the model is significant at the 5% level and that the probability that this F-value is caused by noise is only 4.72%. The coefficient of determination (R^2^) of the model was 0.9169, indicating that the model could explain 91.69% of the variance of the response values, and the fitting degree was good. Among the seven factors investigated, only heating time had a significant effect on the yield of essential oils (p = 0.0049); other factors were not significant (p > 0.05), which indicated that heating time was the key factor affecting the yield of essential oils. By comparing the p-values of the seven factors, we concluded that the importance of the seven factors for the yield of essential oils was heating time > heating power > ultrasonic time > ion type > ultrasonic power > ion type. After comprehensive consideration, the factors of heating time, heating power and ultrasonic time were selected for further optimization.

2.3. Optimization of Essential Oil Yields with Box–Behnken Design

Table 2 shows 17 groups of experiments designed to study the interaction of heating time, heating power and ultrasonic time on essential oil yield by Design Expert 13.0 software. Multivariate regression fitting was performed on the experimental data to obtain a second-order polynomial model between essential oil yield ( $[eqn]$ ) and each coding factor as follows:

$$[eqn]$$

where $[eqn]$ is the yield and $[eqn]$ , $[eqn]$ , $[eqn]$ are ultrasonic time, heating power, and heating time, respectively.

The results of the ANOVA analysis of C. citrinus essential oil yield based on response surface design are shown in Table 2. The F-value (10.19) and low p-value (0.0029) of the model showed that the model was statistically significant. The probability that the F-value of the model is caused by error is only 0.29%. In this study, the lack of fit was not significant (p > 0.05), indicating that the model fit was good and no obvious lack-of-fit phenomenon was found. The coefficient of determination ( $[eqn]$ ) of the model was 0.9291, indicating that 92.91% of the essential oil yield values matched the predicted values of the model. The adjusted R^2^ (Adj. $[eqn]$ ) was 0.8379, indicating that the regression equation was in good agreement with the actual experiments.

Table 2 shows that the linear factors ultrasonic time ( $[eqn]$ ), heating power ( $[eqn]$ ), heating time ( $[eqn]$ ) and quadratic factors were significant (p < 0.05) for essential oil yield. Among them, the linear effect of heating time ( $[eqn]$ ) was the most significant (p < 0.001), indicating that heating time was the most important factor affecting the yield of essential oil. The heating power (X6) and ultrasonic time (X4) also showed significant linear effects. The interaction factors were not significant (p > 0.05), indicating that the interaction between various factors had little effect on the yield of essential oil. Among the quadratic factors, only $[eqn]$ was significant, indicating a clear nonlinear relationship between ultrasound time and yield. The factors affecting the volatile oil extraction rate were sorted according to the F-value—heating time> ultrasonic time> heating power—indicating that the heating time had a greater influence on the volatile oil extraction rate and that the influence of ultrasonic time and heating power on the yield of volatile oil was not much different.

According to the regression equation, the response surface and contour plot between heating time, heating power and ultrasonic time can be plotted (Figure 2). Figure 2a,d show the interaction between heating time and heating power. The results showed that the yield of essential oil increased from 0.68% to 0.82% with the increase in heating time from 30 min to 50 min but decreased to 0.75% at 60 min after the heating time exceeded 50 min, indicating that there was an optimal heating time. At the same time, when the heating power was increased from 350 W to 450 W, the yield decreased from 0.78% to 0.74%, which indicates that too high power may have an inhibitory effect. Figure 2b,e reflect the interaction between heating time and ultrasonic time, and the results show that when the heating time is extended from 30 min to 50 min and the ultrasound time is increased from 15 min to 25 min, the yield continues to increase and reaches 0.86%. However, the growth slowed down after more than 25 min and only increased to 0.87% at 30 min, indicating that prolonging ultrasound could not further significantly improve the extraction efficiency. Figure 2c,f show the interaction between heating power and ultrasonic time. The results show that the yield increased significantly with ultrasonic time from 15–25 min, from 0.70% to 0.84%, while the heating power dropped slightly to 450 W, from 0.80% to 0.76%. Finally, ultrasound time is the main driving factor, while excessive heating power will inhibit the yield of essential oil. In conclusion, the optimal extraction conditions determined by model optimization were 57 min heating time, 563 W heating power and 23 min ultrasonic time, and the yield of essential oil reached 0.80 mL/kg DW.

2.4. GC–MS Analysis of Essential Oil Components

To clarify the chemical composition of C. citrinus essential oil (EO), this study employed gas chromatography–mass spectrometry (GC–MS) for the separation and identification of oil components. The relative area percentage (RA) in GC–MS was used to represent the relative content of the essential oil compounds. The results are shown in Table 3. The relative peak area (RA) of each component was calculated using the area normalization method. During the detection process, four different ions (Sr^2+^, Na^+^, Ca^2+^, K^+^, corresponding to EO_a,_ EO_b_, EO_c_, and EO_d_ in the data table) were used in the pretreatment experiments, along with one blank group (EO_e_). The GC–MS chromatogram clearly shows that the peaks of key substances in different samples are concentrated at similar retention times, indicating that they have similar compositions.

The results reveal that 79 compounds were separated and identified from C. citrinus essential oil, encompassing various types such as alcohols, ketones, esters, alkenes, and phenols. Specifically, these can be categorized into monoterpenes (16 species), sesquiterpenes (17 species), and oxygenated compounds (46 species). The primary components of the essential oil include eucalyptol, 1R-α-Pinene, Leptospermone, and L-α-Terpineol. The compounds with relatively high concentrations in the essential oil are mainly focused among monoterpenes, sesquiterpenes, and oxygenated monoterpenes. The content proportion of terpenoids is relatively stable with a small fluctuation range. In terms of relative content, the total relative content of terpenoids (including olefinic monoterpenes, sesquiterpenes, alcoholic terpenoid derivatives, etc.) in the five samples are EO_a_ 71.81%, EO_b_ 59.78%, EO_c_ 79.16%, EO_d_ 79.03%, and EO_e_ 81.62%. This indicates that terpenoids are the core active components of the essential oil and that the addition of salts has a slight promoting effect on the increase in the total content of terpenoids.

Monoterpenoid compounds are the primary contributors to the aroma of C. citrinus leaf essential oil, while also exhibiting physiological activities such as antibacterial and anti-inflammatory properties. Variations in their content directly affect the sensory quality and application value of the essential oil. Limonene mainly imparts a fresh citrus fragrance to the oil. The cyclic cluster heatmap of the compounds in the essential oils is shown in Figure 3. In the figure, the EO_a_ ring exhibits a darker color in the Limonene region, indicating a relatively high RA value. It is inferred that Sr^2+^ may promote the release of Limonene within the cells by regulating the permeability of the cell membranes of leaves, thereby reducing volatilization losses during the extraction process.

By contrast, the corresponding area in the EO_e_ ring appears lighter, with a lower RA value, as pure water cannot alter the cell membrane structure, which results in a relatively lower content of Limonene compared to the strontium chloride treatment group. 1R-α-Pinene imparts a fresh, pine-needle-like woody aroma to essential oils and exhibits strong antioxidant activity. In the EO_c_ zone, the 1R-α-Pinene region shows a darker color and higher RA values, which may be attributed to Ca^2+^ acting as a key ion in plant cell signal transduction. Ca^2+^ can activate enzymes related to 1R-α-Pinene biosynthesis, such as terpene synthase, within the leaves and also enhance cell rupture efficiency, which leads to a more complete release of 1R-α-Pinene. By contrast, in the EO_b_ zone, the RA values of 1R-α-Pinene are lower because Na^+^ has a weaker regulatory effect on cell membrane permeability and cannot effectively activate terpene synthase, which results in no significant promotion of 1R-α-Pinene synthesis or release [18] (Schilcher, 1993).

Based on the specific chemical composition of each sample, the biological activities of different essential oils also vary. Regarding antimicrobial activity, EO_a_ may exhibit the strongest antibacterial properties due to its higher content of monoterpene alcohol compounds, which are generally known for their excellent antibacterial efficacy. In terms of antioxidant activity, EO_c_ and EO_d_ may hold greater advantages, as they are rich in ketone compounds, which generally demonstrate potent antioxidant capabilities. Regarding anti-inflammatory activity, EO_c_ may demonstrate more pronounced effects, closely linked to its abundant sesquiterpene content. Sesquiterpene compounds have been extensively documented to exhibit significant anti-inflammatory effects [19] (Cho et al., 2000).

In summary, different ionic treatments may significantly influence the RA of compounds such as monoterpenes, sesquiterpenes, and alcohols in essential oils by regulating the permeability of C. citrinus leaf cell membranes, activating relevant synthases, or inhibiting oxidative degradation processes. Changes in the content of these components not only directly determine the aroma characteristics of essential oils (such as citrus, woody, and sweet floral scents) but also relate to their biological activity (such as antibacterial and anti-inflammatory properties). Further comparative analysis of different metal salts revealed that, compared to the control group, strontium chloride (EO_a_) and sodium chloride (EO_b_) favored the accumulation of monoterpenes, while potassium chloride (EO_d_) and calcium chloride (EO_c_) promoted the release of sesquiterpenes and oxygen-containing compounds. This indicates that ion treatment methods can be selected based on the requirements of target bioactive compounds to selectively enhance functional properties such as antioxidant and antibacterial activity in essential oils. Furthermore, for easily oxidized components like terpenes, their relative content and color changes under different ion treatments suggest potential stability differences. Because terpenes readily oxidize, leading to essential oil degradation and off-flavor formation, these results also indicate that ion treatment methods significantly impact the storage stability of essential oils.

2.5. Principal Component Analysis

To systematically evaluate the overall impact of different metal salt pretreatments on the chemical composition of C. citrinus leaf essential oil, principal component analysis (PCA) was performed on the GC–MS relative content data of five essential oil samples (EO_a–e_). The analysis results are presented in Figure 4, which is a PCA biplot showing both sample distribution and compound contributions. The PCA results indicated that the cumulative variance contribution rate of the first two principal components reached 97.8%. This demonstrates that these two principal components can fully capture and explain most of the variation information in the original GC–MS data, confirming high model reliability.

In Figure 4, both EO_a_ and EO_b_ are located on the positive semi-axis of PC1 but separated along the PC2 axis, which indicates that the effects of the Sr^2+^ and Na^+^ pretreatments on the chemical composition of the essential oil have both similarities and specificity. The three samples (EO_c_, EO_d_, and EO_e_) cluster on the negative semi-axis of PC1, which suggests that their chemical compositions are relatively similar but significantly different from those of EO_a_ and EO_b_, especially EO_e_, which has the largest projection in the negative direction of PC2. Key compounds in the right half of the figure, such as eucalyptol and α-pinene, have loading vectors pointing to the positive direction of PC1, which is fully consistent with the GC–MS data showing that EO_a_ and EO_b_ have a higher eucalyptol content and that EO_a_ is rich in α-pinene. Compounds in the upper and left halves of the figure make a significant negative contribution to PC1 and PC2, which indicates that the Ca^2+^ and K^+^ pretreatments, as well as the control treatment, are more likely to promote the dissolution of oxygenated compounds. Therefore, Sr^2+^ and Na^+^ pretreatments tend to enhance the release of monoterpene olefins, while Ca^2+^ and K^+^ pretreatments are closely related to the content of specific oxygenated terpenoids and ketone compounds, with the chemical profile of the blank control group falling between the two.

2.6. Microscopic Change Analysis

2.6.1. Electron Scanning Microscopy (SEM) Analysis

Figure 5 shows the microscopic morphology of raw material leaf powder and extracted filter residue visualized with scanning electron microscopy. Figure 5a,c,e show micrographs extracted without ultrasound at different scales, and Figure 5b,d,f correspond to SEM images of C. citrinus leaf powder that has undergone ultrasonic extraction at different scales. Under SEM, the surface of the sample after ultrasonic extraction showed significant wrinkles, cracks, and holes, compared to the relatively intact and smooth surface of the sample without sonication. This was due to the powerful ultrasonic cavitation effect that produces instantaneous high pressure and shock waves [20] (Shen et al., 2023), which physically break the solid cell wall and effectively destroy the microstructure of the C. citrinus leaves, which increases the mass transfer area and channels [21] (Shi et al., 2025), making the essential oil more easily soluble and extracted by solvents.

2.6.2. Fourier Transform Infrared Spectroscopy Analysis

Fourier transform infrared spectroscopy (FTIR) analysis confirmed the successful extraction of essential oils by ultrasound-assisted steam distillation. As shown in Figure 6, the structural analysis of raw leaf powder (a) and extracted filter residue (b) of C. citrinus by Fourier transform infrared spectroscopy shows that the new absorption peaks observed at 2918.8 $[eqn]$ and 2847.0 $[eqn]$ are attributed to C-H telescopic vibrations, indicating the presence of fat chain structures unique to terpenoids. The strong and wide absorption bands at 3305.8 $[eqn]$ and 3329.4 $[eqn]$ are attributed to the O-H telescopic vibration, indicating the presence of phenolic or alcoholic components in the essential oil. The characteristic peak at 1686.9 $[eqn]$ corresponds to the C=O telescopic vibration, indicating the presence of carbonyl-containing compounds such as aldehydes, ketones, or esters. The displacement of the C-O telescopic vibration peak to 1031.4 $[eqn]$ further demonstrates the successful extraction of oxygen-containing terpenoids. These spectral results indicate that essential oils contain characteristic functional groups of terpenes and their oxygen-containing derivatives, confirming the effectiveness of the extraction process.

3. Materials and Methods

3.1. Materials and Reagents

Fresh plant materials were collected in May 2025 from areas surrounding the campus of Jiaying University, Guangdong, China. The following reagents were used in this study: choline chloride (C_5_H_14_ClNO, purity ≥ 98%; Yuanye Bio-Technology Co., Ltd. (Shanghai, China), catalog No. S30169-500 g), sodium hydroxide (NaOH, purity ≥ 97%; Shanghai Zhanyun Chemical Co., Ltd. (Shanghai, China)), and hydrochloric acid (HCl, purity 36–38%; Xilong Scientific Co., Ltd. (Guangzhou, China)). These reagents were employed for the preparation of deep eutectic solvents (DESs).

The inorganic salts used for pretreatment included sodium chloride (NaCl, purity ≥ 99.5%; Sinopharm Chemical Reagent Co., Ltd. (Shanghai, China)), strontium chloride (SrCl_2_, purity ≥ 99.5%; Tianjin Zhonglian Chemical Reagent Co., Ltd. (Tianjin, China)), calcium chloride (CaCl_2_, purity ≥ 96%; Tianjin Beichen Fangzheng Reagent Co., Ltd. (Tianjin, China)), and potassium chloride (KCl, purity ≥ 99.5%; Tianjin Baishi Chemical Co., Ltd. (Tianjin, China)). The deep eutectic solvents (DESs) used in this study were prepared by combining choline chloride, which acted as the hydrogen bond acceptor, with dilute hydrochloric acid or sodium hydroxide as hydrogen bond donors to construct eutectic systems with different pH values. The components were mixed at predetermined ratios and heated at 80 °C under continuous magnetic stirring until a transparent and homogeneous liquid was obtained, which indicated the formation of a stable eutectic phase.

After cooling to room temperature, inorganic salts (NaCl, KCl, SrCl_2_, or CaCl_2_) were directly dissolved into the prepared DES to obtain ion-containing pretreatment solutions with controlled pH and ionic composition.

3.2. Instruments and Equipment

The main instruments and equipment used in this research include: electronic analytical balance (JM-B5003, Zhuji Chaoze Weighing Instrument Equipment Co., Ltd., Shaoxing, China), grinding crusher (2500C, Yongkang Hongtaiyang Electromechanical Co., Ltd., Yongkang, China), temperature regulating heating mantle (DZTW 3000 mL, Shanghai Lizhan Bangxi Instrument Technology Co., Ltd., Shanghai, China), numerical control ultrasonic cleaning instrument (Shenzhen Peninsula Medical Co., Ltd., Shenzhen, China), distillation apparatus (1000 mL, Sichuan Shubo Co., Ltd., Sichuan, China), scanning electron microscope (Zeiss, Carl Zeiss AG, Shanghai, China), gas chromatography–mass spectrometry (Agilent 6890GC, Agilent Technologies, Inc., Beijing, China) and Fourier transform infrared spectrometer (Nicolet Magna-IR 560, Thermo Fisher Scientific, Inc., Shanghai, China).

3.3. Optimization of the Experimental Process

In the process of optimizing the extraction process, the Plackett–Burman design (PBD) was used to screen the key influencing factors based on the results of the univariate test. Seven factors—pH (X1), ion type (X2), ion mass (X3), ultrasonic time (X4), ultrasonic power (X5), heating power (X6) and heating time (X7)—were selected as the investigation variables. All the factors were set into two levels: high (+1) or low (−1). The factor codes and levels are shown in Table S1. Subsequently, the steepest climbing experiment was carried out, starting from the center point of PBD and moving along a path that gradually increased in heating time, heating power and ultrasonic time. The experimental design and results are shown in Table S2. The optimal parameters for the extraction process were a heating time of 50 min, heating power of 566 W, and ultrasonic time of 20 min, while the extraction efficiency approaches the maximum response range. On this basis, the three-factor three-level Box–Behnken design (BBD) was further adopted, and the heating time, heating power and ultrasonic time were selected as independent variables, and the volatile oil yield was used as the response value, and the response surface was analyzed and optimized using Design-Expert 13.0 software. The factor coding and level are shown in Table S3.

3.4. The Characterization of the Samples

3.4.1. GC–MS Characterization of Essential Oil Components

The essential oil components were characterized by gas chromatography–mass spectrometry (GC–MS) under different ionic treatments using a DB-5 capillary column (30 m × 0.25 mm, 0.25 μm). The carrier gas (He) was maintained at a flow rate of 1.0 mL/min, and the column temperature for injecting 1 μL of sample in split mode was set to an initial temperature of 50 °C, warmed up to 180 °C at a rate of 3 °C/min for 10 min, and then warmed to 280 °C at a rate of 10 °C/min. The injection temperature was set at 280 °C. Compound identification was achieved by matching mass spectrometry to NIST libraries.

3.4.2. Micromorphology of Extracts

In order to explore the failure mechanism of ultrasound-assisted pretreatment of different inorganic salts on plant cell structure, a morphology analysis of dried leaf powder before extraction and filter residue after steam distillation was performed. After drying and sieving, the sample was fixed on the conductive glue and sprayed with gold, and the structure of the sample was observed by scanning electron microscopy (SEM) with an acceleration voltage of 3 kV and a working distance of 7.0 mm or 7.2 mm.

3.4.3. Fourier Transform Infrared Spectroscopy Analyzes the Chemical Structure of Extracted Samples

Liquid nitrogen was used to homogenize the freeze-dried samples (leaf powder raw material and post-extraction residue) into fine powders, and approximately 1 mg of each powder sample was thoroughly mixed with 100 mg of spectral-grade potassium bromide (KBr) and pressed into transparent pellets under hydraulic pressure. Fourier transform infrared spectroscopy (FTIR) measurements were performed using a Nicolet Magna-IR 560 Fourier transform infrared spectrometer. The test conditions were as follows: the wavenumber range was set to 4000–400 cm^−1^, the resolution was set to 4 cm^−1^, and each sample was scanned 64 times to ensure that the signal-to-noise ratio met the analysis requirements. The characteristic functional groups were identified by comparing the obtained absorption peaks with the standard spectra reported in the existing literature. All assays were repeated three times independently to ensure the reproducibility and reliability of the experimental data.

4. Conclusions

It should be noted that the novelty of the present work lies primarily in the process–composition relationship rather than in the extraction technique itself. While the general extraction framework has been previously described, this study emphasizes how ion-assisted pretreatment, coupled with multivariate optimization, influences the distribution of key volatile components in C. citrinus essential oil. Such an approach extends conventional yield-focused optimization toward a more composition-driven perspective. In this study, an ultrasound-assisted metal salt pretreatment strategy was optimized to enhance the extraction efficiency of C. citrinus essential oil and to elucidate the regulatory effects of different metal ions on its chemical composition. Response surface methodology indicated that heating time was the dominant factor influencing essential oil yield. Furthermore, GC–MS combined with PCA analysis demonstrated that metal ions with different valence states selectively modulated the chemical profile of the extracted essential oils. Specifically, Sr^2+^ and Ca^2+^ favored the enrichment of monoterpenoids, whereas Na^+^ and K^+^ promoted the release of certain sesquiterpenes and oxygenated compounds. Overall, this study establishes a green and efficient approach for both the extraction and compositional regulation of C. citrinus essential oil and demonstrates that its chemical profile can be tailored through appropriate salt pretreatment. These findings provide an experimental basis for the targeted production of essential oils with specific aroma characteristics or potential functional properties.