Determination of 60 Migrant Substances in Plastic Food Contact Materials by Vortex-Assisted Liquid-Liquid Extraction and GC-Q-Orbitrap HRMS

A GC-HRMS analytical method for the determination of 60 migrant substances, including aldehydes, ketones, phthalates and other plasticizers, phenol derivatives, acrylates, and methacrylates, in plastic food contact materials (FCM) has been developed and validated. The proposed method includes migration tests, according to Commission Regulation (EU) 10/2011, using four food simulants (A, B, C, and D1), followed by vortex-assisted liquid–liquid extraction (VA-LLE) and GC-Q-Orbitrap HRMS analysis in selected ion monitoring (SIM) mode, with a resolving power of 30,000 FWHM and a mass accuracy ≤5 ppm. The method was validated, showing satisfactory linearity (R2 ≥ 0.98 from 40 to 400 µg L−1), limits of quantification (40 µg L−1), precision (RSD, 0.6–12.6%), and relative recovery (81–120%). The proposed method was applied to the analysis of field samples, including an epoxy-coated tin food can, a drinking bottle made of Tritan copolyester, a disposable glass made of polycarbonate, and a baby feeding bottle made of polypropylene, showing that they were in compliance with the current European regulation regarding the studied substances.


Introduction
In Europe, the 'Commission Regulation (EU) No 10/2011 of 14 January 2011 on plastic materials and articles intended to come into contact with food' [1] (henceforth, Regulation 10/2011) contains the 'Union list of authorized substances' (IAS), which establishes the specific migration limits (SML) that they should comply with in order to ensure the safety of the consumers. Moreover, non-intentionally added substances (NIAS) could also be present in plastic food contact materials (FCM) as impurities of the authorized substances, reaction, or degradation products, formed during the manufacture, storage, or consumption processes [2]. In this sense, the development and validation of analytical methods for the quantification of migrant substances, both IAS and NIAS, in plastic FCM using food simulants are a matter of importance and a research field of growing interest [3][4][5][6].
While most of the published studies are focused on the analysis of particular families or groups of substances with a reduced number of target analytes, few multi-residue or multi-analyte methods are available in the analytical literature. The most significant are the determination of 48 contaminant residues in food contact plastic products by microwave-assisted extraction, using methanol as an alternative 'non-regulated' food simulant, followed by LC-HRMS, published by Zhang et al. [34]; the determination of 84 migrant substances using food simulants A and C, followed by liquid-liquid extraction using dichloromethane as extraction solvent and GC-MS, published by Tsochatzis et al. [35]; and the determination of 75 migrant substances using food simulants A, C, and D1, followed by salt-assisted liquid-liquid extraction using dichloromethane as extraction solvent and GC-MS/MS, also published by Tsochatzis et al. [36]. Moreover, GC presents the handicap of not being directly compatible with the hydro-ethanolic food simulants of Regulation 10/2011 [1]. For this reason, it is required to use alternative food simulants [17,18,23], such as organic solvents, or to perform a liquid-phase [15,35,36] or a solid-phase extraction [24,27] in order to enable the GC analysis. With regard to the detection technique, currently there are no published methods that use GC coupled to HRMS, which enhances the identification and confirmation confidence of the migrant substances by means of their exact masses.
To the best of our knowledge, this is the first work addressing a simple and reliable analytical procedure for the determination of 60 migrant substances in plastic FCM using four food simulants (A, B, C, and D1), according to Regulation 10/2011 [1], followed by vortex-assisted liquid-liquid extraction (VA-LLE), using only n-hexane as an extraction solvent, and GC-Q-Orbitrap HRMS analysis.
However, in the current market, there is a wide variety of FCM made of different plastics and polymer materials, such as epoxy-coatings, polypropylene (PP), polyethylene (PE), polyethylene terephthalate (PET), polycarbonate (PC), and polyester copolymers, among others. Furthermore, not all the studied substances, both IAS and NIAS, may be relevant for all the plastic FCM. In this sense, a fast and reliable multi-analyte method to quantify a large set of GC-suitable migrant substances, such as those proposed in the present manuscript, could be useful for food control authorities in order to evaluate compliance of plastic FCM with current European regulations as a part of a comprehensive analytical control plan [1,37].

Selection of the Target Analytes
The abovementioned Regulation 10/2011 [1], in its 'Union list of authorized substances', contains approximately one thousand IAS and their applicable restrictions, including SML. Moreover, NIAS should also be controlled to ensure that their migration does not surpass the established limits [1].
From this enormous number of substances to be controlled, the selection of the target analytes was performed on the basis of two criteria: their analytical standards should be commercially available, and they should be suitable for determination with the selected analytical technique, GC-HRMS (mainly hydrophobic organic compounds). In this sense, inorganic substances, such as oxides, hydroxides, silicates, and metals; macromolecules and oligomers, such as sugars, waxes, and resins; and highly polar compounds, such as acids, alcohols, and amines were not considered. Finally, a representative set of 60 substances, including aldehydes, ketones, phthalates and other plasticizers, phenol derivatives, acrylates, and methacrylates, among others, were selected as target analytes (see Table 1). However, this list of target analytes was not meant to be exhaustive, and it could be extended if necessary. . For authorized substances (IAS) with no SML, the applicable limit is 60 mg kg −1 ; and for substances whose SML is ND (not detected) and NIAS, the applicable limit is 0.01 mg kg −1 .

Selection of the Quantification and Confirmation Ions
In order to select the quantification and confirmation ions of the target analytes for their acquisition in SIM mode (see Section 3.4.4), individual standard solutions containing the target analytes at 400 µg L −1 in n-hexane were prepared and injected directly into the GC-HRMS system in full scan mode with a resolving power of 60,000 FWHM, in order to obtain their HRMS spectrum. For each target substance, the most abundant fragment ions were selected as quantification and confirmation ions, respectively (see Table S1). If possible, the molecular formulas of the selected fragment ions were obtained, and their theoretical exact masses were considered for further acquisition.
As an example, the chromatographic peak and the selected fragment ions for Octocrylene, with their exact masses, are shown in Figure 1.

Study of the Experimental Variables Involved in the VA-LLE Procedure
In order to make the aqueous (or hydro-ethanolic) food simulants obtained from the migration test of the samples (see Section 3.4.1) compatible with the selected analytical technique (GC-HRMS), VA-LLE with an organic solvent was performed.
The experimental variables involved in the extraction process were studied, in triplicate, and those which provided the highest analytical signal (peak area) for a broad number of target compounds were selected. In this sense, extraction solvent composition and volume, ionic strength, and extraction time were studied. To perform these studies, a standard solution containing the target analytes at 160 µg L −1 in food simulant A was used. The volume of the food simulant was fixed at 20 mL.

Extraction Solvent Composition
The organic solvent to perform the extraction had to be immiscible and less dense than water in order to facilitate the collection of the extract as a supernatant phase. In this sense, n-hexane, toluene, and a mixture of acetone/n-hexane (1:1) were tested as extraction solvents. Moreover, chlorinated solvents were not considered due to their higher density and toxicity. The conditions for the other variables were as follows: extraction solvent volume, 5 mL; ionic strength, 5% w/v of NaCl; extraction time, 60 s. The obtained results are shown in Figure 2. It should be noted that only a representative set of nine target analytes (benzophenone, bis(2-ethylhexyl) phthalate, butylated hydroxytoluene, cyclohexyl methacrylate, dibutyl phthalate, diphenyl sulphone, methyl benzoate, octocrylene, and tri-n-butyl acetyl citrate) have been included in the figures below in order to simplify the visualization of the results. These substances were selected for being the most common and/or representative of the different chemical families studied in the present manuscript. Similar trends and conclusions were observed for the other target compounds.

Study of the Experimental Variables Involved in the VA-LLE Procedure
In order to make the aqueous (or hydro-ethanolic) food simulants obtained from the migration test of the samples (see Section 3.4.1) compatible with the selected analytical technique (GC-HRMS), VA-LLE with an organic solvent was performed.
The experimental variables involved in the extraction process were studied, in triplicate, and those which provided the highest analytical signal (peak area) for a broad number of target compounds were selected. In this sense, extraction solvent composition and volume, ionic strength, and extraction time were studied. To perform these studies, a standard solution containing the target analytes at 160 µg L −1 in food simulant A was used. The volume of the food simulant was fixed at 20 mL.

Extraction Solvent Composition
The organic solvent to perform the extraction had to be immiscible and less dense than water in order to facilitate the collection of the extract as a supernatant phase. In this sense, n-hexane, toluene, and a mixture of acetone/n-hexane (1:1) were tested as extraction solvents. Moreover, chlorinated solvents were not considered due to their higher density and toxicity. The conditions for the other variables were as follows: extraction solvent volume, 5 mL; ionic strength, 5% w/v of NaCl; extraction time, 60 s. The obtained results are shown in Figure 2. It should be noted that only a representative set of nine target analytes (benzophenone, bis(2-ethylhexyl) phthalate, butylated hydroxytoluene, cyclohexyl methacrylate, dibutyl phthalate, diphenyl sulphone, methyl benzoate, octocrylene, and tri-n-butyl acetyl citrate) have been included in the figures below in order to simplify the visualization of the results. These substances were selected for being the most common and/or representative of the different chemical families studied in the present manuscript. Similar trends and conclusions were observed for the other target compounds. As can be seen, the best results were obtained with n-hexane, so it was selected for further analysis. As can be seen, the best results were obtained with n-hexane, so it was selected for further analysis.

Extraction Solvent Volume
The mass transfer of the analytes from the food simulant to the extraction solvent depends on the ratio between the food simulant and the extraction solvent volumes. In this sense, the volume of the food simulant was fixed at 20 mL, and 5, 10, and 20 mL were tested as extraction solvent volumes. The conditions for the other variables were as follows: extraction solvent, n-hexane; ionic strength, 5% w/v of NaCl; extraction time, 60 s. The obtained results are shown in Figure 3. From the point of view of the mass transfer of the analytes, a higher volume of extraction solvent implies that a higher amount of the analytes could be extracted, but the obtained extract would be more diluted. As can be seen, the analytical signal decreased with higher volumes due to the dilution effect, so 5 mL of extraction solvent were selected for further analysis. Moreover, the obtained sensitivity using the selected volume of extraction solvent allowed us to reduce the organic solvent consumption of the proposed method and to avoid further experimental steps, such as repeated extractions or evaporation and reconstitution steps.

Ionic Strength
The influence of the ionic strength of the food simulant was studied in order to evaluate the so-called 'salting-out effect'. For this, additions of 0, 2.5, 5, and 10% w/v of NaCl were tested. The conditions for the other variables were as follows: extraction solvent, n-hexane; extraction solvent volume, 5 mL; extraction time, 60 s. These results are shown in Figure 4.
As can be seen, the best results were obtained with the addition of 0% w/v of NaCl, so no adjustment of the ionic strength was selected for further analysis. From the point of view of the mass transfer of the analytes, a higher volume of extraction solvent implies that a higher amount of the analytes could be extracted, but the obtained extract would be more diluted. As can be seen, the analytical signal decreased with higher volumes due to the dilution effect, so 5 mL of extraction solvent were selected for further analysis. Moreover, the obtained sensitivity using the selected volume of extraction solvent allowed us to reduce the organic solvent consumption of the proposed method and to avoid further experimental steps, such as repeated extractions or evaporation and reconstitution steps.

Ionic Strength
The influence of the ionic strength of the food simulant was studied in order to evaluate the so-called 'salting-out effect'. For this, additions of 0, 2.5, 5, and 10% w/v of NaCl were tested. The conditions for the other variables were as follows: extraction solvent, n-hexane; extraction solvent volume, 5 mL; extraction time, 60 s. These results are shown in Figure 4.
As can be seen, the best results were obtained with the addition of 0% w/v of NaCl, so no adjustment of the ionic strength was selected for further analysis.

Extraction Time
The influence of the extraction time, as the time of vortex mixing, was studied. For this, 30, 60, 90, and 120 s were tested. The conditions for the other variables were as follows: extraction solvent, n-hexane; extraction solvent volume, 5 mL; ionic strength, not adjusted. These results are shown in Figure 5.
As can be seen, similar results were observed between 30 and 60 s, so the shortest extraction time of 30 s was selected for further analysis. With 90 and 120 s, a decrease in the analytical response was observed for some of the target analytes. This may be due to kinetic or thermodynamic effects that alter the partition coefficient of the analytes between the food simulant and the extraction solvent with longer extraction times. evaporation and reconstitution steps.

Ionic Strength
The influence of the ionic strength of the food simulant was studied in order to evaluate the so-called 'salting-out effect'. For this, additions of 0, 2.5, 5, and 10% w/v of NaCl were tested. The conditions for the other variables were as follows: extraction solvent, n-hexane; extraction solvent volume, 5 mL; extraction time, 60 s. These results are shown in Figure 4.
As can be seen, the best results were obtained with the addition of 0% w/v of NaCl, so no adjustment of the ionic strength was selected for further analysis.

Extraction Time
The influence of the extraction time, as the time of vortex mixing, was studied. For this, 30, 60, 90, and 120 s were tested. The conditions for the other variables were as follows: extraction solvent, n-hexane; extraction solvent volume, 5 mL; ionic strength, not adjusted. These results are shown in Figure 5.
As can be seen, similar results were observed between 30 and 60 s, so the shortest extraction time of 30 s was selected for further analysis. With 90 and 120 s, a decrease in the analytical response was observed for some of the target analytes. This may be due to kinetic or thermodynamic effects that alter the partition coefficient of the analytes between the food simulant and the extraction solvent with longer extraction times.

Study of the Extraction Efficiency
In order to evaluate the extraction efficiency of the proposed VA-LLE under the selected conditions, a standard solution containing the target analytes at 160 μg L -1 in n-hexane was directly injected into the GC-HRMS system, and the obtained analytical signal was compared with that obtained for a standard solution containing the target analytes at 160 μg L −1 in food simulant after performing the extraction process (see Section 3.4.3).
Considering the ratio between the food simulant (20 mL) and the extraction solvent (5 mL) volumes, the extraction yield was determined. The obtained absolute recovery values (extraction yield) were different depending on the target analyte, ranging from 8 to 110% (Table S2). In this sense, the most polar analytes (such as 4-methylphenol, butyl lactate, and phenol, among others) showed the lowest extraction efficiencies. Moreover, the composition of the studied food simulants (A, B, C, and D1) also affected the observed extraction yield, thus showing that procedural standard calibration (i.e., performing the extraction process in the preparation of the calibration curve) would be

Study of the Extraction Efficiency
In order to evaluate the extraction efficiency of the proposed VA-LLE under the selected conditions, a standard solution containing the target analytes at 160 µg L −1 in n-hexane was directly injected into the GC-HRMS system, and the obtained analytical signal was compared with that obtained for a standard solution containing the target analytes at 160 µg L −1 in food simulant after performing the extraction process (see Section 3.4.3).
Considering the ratio between the food simulant (20 mL) and the extraction solvent (5 mL) volumes, the extraction yield was determined. The obtained absolute recovery values (extraction yield) were different depending on the target analyte, ranging from 8 to 110% (Table S2). In this sense, the most polar analytes (such as 4-methylphenol, butyl lactate, and phenol, among others) showed the lowest extraction efficiencies. Moreover, the composition of the studied food simulants (A, B, C, and D1) also affected the observed extraction yield, thus showing that procedural standard calibration (i.e., performing the extraction process in the preparation of the calibration curve) would be necessary in order to compensate for the variation of the extraction yields to obtain accurate results.

Study of Stability of the Target Analytes
In cases where it was necessary to perform the proposed analytical method, including migration tests, VA-LLE, and GC-HRMS analysis, in different working sessions, the stability of the target analytes in the food simulants and in the n-hexane extracts was studied.
In triplicate, a standard solution containing the target analytes at 160 µg L −1 in the food simulant was analyzed following the proposed method immediately after its preparation, and after 1, 3, and 5 days of fridge storage (4 • C). The variation in the analytical signal was evaluated, and the obtained results shown that the target analytes were stable during the first 24 h of fridge storage (RSD, <15%). After 1 day of fridge storage, a decrease of the analytical signal was observed for some of the target analytes, such as 3-(trimethoxysilyl)propyl methacrylate, lauryl acrylate, triethoxyvinyl silane, triethyl phosphite, and vinyltrimethoxysilane. Similar results were obtained for all the tested food simulants (A, B, C, and D1).
Moreover, the stability of the target analytes was also tested in the n-hexane extract obtained after the VA-LLE analysis. In triplicate, a standard solution containing the target analytes at 160 µg L −1 in the food simulant was subjected to VA-LLE, and the obtained n-hexane extract was analyzed immediately, and after 1, 3, and 5 days of freezer storage (−24 • C). The variation in the analytical signal was evaluated, and the obtained results showed that the target analytes were stable for at least 3 days of freezer storage after the VA-LLE analysis (RSD, <20%). Similar results were obtained for all the tested food simulants (A, B, C, and D1).

Method Validation
In order to validate the proposed methodology, the analytical figures of merit (i.e., linearity, limits of quantification (LOQ), precision, and relative recovery) were evaluated in the studied food simulants (A, B, C, and D1) following the guidelines of the 'Analytical quality control and method validation procedures for pesticide residues analysis in food and feed' (Document No. SANTE/12682/2019) [38]. The obtained analytical figures of merit using food simulant A are shown in Table 2, while those obtained for food simulants B, C, and D1 are shown in Tables S3-S5, respectively.
Linearity comprised at least from 40 to 400 µg L −1 , obtaining coefficients of determination (R 2 ) ≥ 0.98 in all the cases.
The LOQ, defined for practical reasons as the lowest concentration of the calibration curve, was fixed at 40 µg L −1 in the food simulant. Although many of the target analytes could be determined in concentrations below the reported LOQ, lower concentration limits were not pursued because the achieved method LOQ enables the evaluation of compliance with their SML [1]. The method LOQ, expressed in mg kg −1 of food, depends on the contact surface of the material exposed to the food simulant, considering the conventional ratio of 6 dm 2 of contact surface per kg of food (see Section 3.4.5). As an example, a migration test performed with 1 dm 2 of the plastic FCM exposed to 40 mL of food simulant would provide a method LOQ of 0.01 mg kg −1 , applying the proposed methodology. According to Regulation 10/2011 [1], the SML for IAS range from 0.05 to 60 mg kg −1 of food, except for those whose SML states ND (not detected) and for NIAS, which should not be found in concentrations higher than 0.01 mg kg −1 [1].
Precision was evaluated in terms of the relative standard deviation (RSD, %) obtained in the analysis of five replicates of standard solutions containing the target analytes at three concentration levels (40, 160, and 400 µg L −1 ). The obtained results ranged from 0.6 to 12.6%, thus showing that satisfactory precision was achieved. In order to validate the method accuracy, the matrix effects were studied by comparing the analytical response obtained in the analysis of solvent standards and matrix-matched standards (spiked samples). In this sense, the relative recovery values (%) obtained in the analysis of spiked and non-spiked field samples, applying the proposed methodology with the four food simulants (A, B, C, and D1), were calculated.
In triplicate, the sample extracts obtained in the migration tests of the field samples (see Section 3.4.1) were spiked with the target analytes at three concentration levels (40, 160, and 400 µg L −1 ), and they were quantified using the proposed procedural standard calibration. The relative recovery values (%) were evaluated by comparing the determined concentration of the target analytes in the spiked and non-spiked samples with the known spiked concentration (solvent standard). The obtained relative recovery values ranged from 81 to 120%, thus showing that there were no matrix effects and, therefore, the proposed procedural standard calibration provided accurate results.

Analysis of Field Samples
The proposed method was applied in the analysis of four field samples made of different plastic FCM (see Section 3.2). The migration tests were performed according to Regulation 10/2011 [1] (see Section 3.4.1).
In the analysis of the field samples, no target analytes were found in concentrations above the method LOQ, thus showing that the tested samples were in accordance with the abovementioned regulation with regard to the studied substances [1]. Moreover, the study of matrix effects discussed above (see Section 2.6.1) demonstrated that there were no matrix effects that could affect the accuracy of the obtained results.
However, it could be interesting to apply the proposed methodology to a larger number of field samples of plastic FCM to obtain data regarding the appearance rate of the studied substances as market surveillance.

Reagents
Analytical standards of the target analytes were purchased from Merck KGaA (Darmstadt, Germany). The migrant substances studied in this paper, and their SML, are summarized in Table 1.

Samples
An epoxy-coated tin food can (Sample 1), a drinking bottle made of Tritan copolyester (Sample 2), a disposable glass made of polycarbonate (Sample 3), and a baby feeding bottle made of polypropylene (Sample 4) were analyzed. The samples were obtained directly from the manufacturer, and they had not previously been in contact with food.

Preparation of Samples. Migration Tests
According to Regulation 10/2011 [1], migration tests have to be performed by exposing the plastic FCM in the appropriate food simulant under standardized test conditions, according to its worst foreseeable use.
The food simulants that were considered and can be applied in the proposed method are as follows: food simulant A (ethanol, 10% v/v), food simulant B (acetic acid, 3% w/v), food simulant C (ethanol, 20% v/v), and food simulant D1 (ethanol, 50% v/v).
The intended use of the analyzed samples and the selected conditions for the migration tests are summarized in Table 3. Briefly, Sample 1, intended to contain preserved vegetables or animal products in an oily medium at room temperature for a long period of time (above 30 days), was tested using food simulant A for 10 days at 60 • C; Sample 2, intended to contain clear drinks, such as water, fruit juice, or energy drinks, was tested using food simulant B for 24 h at 40 • C; Sample 3, intended to contain alcoholic or non-alcoholic beverages, such as wine or beer, was tested using food simulant C for 2 h at 40 • C; and Sample 4, intended to contain milk or milk-based drinks, including reconstituted milkbased on infant formula, was tested using food simulant D1 for 2 hours at 70 • C, followed by 24 h at 40 • C.
Three individual units of each sample were filled with the selected food simulant, sealed with aluminum foil to prevent solvent evaporation, and placed into a laboratory oven for its incubation. After the migration tests, the food simulants were collected and submitted to the subsequent VA-LLE procedure.
It should be noted that the specific conditions for the migration test (food simulant composition, contact time, and contact temperature) depend on the worst foreseeable use of the plastic FCM. The methodology proposed in this paper can be applied to the analysis of any plastic FCM that uses any of the described food simulants (A, B, C, and D1) for its migration test.

Preparation of Standards
Stock individual standard solutions of the target analytes and the internal standards were prepared at 1600 mg L −1 in methanol. From these solutions, stock multicomponent solutions containing the target analytes and the internal standards, respectively, at 16 mg L −1 were prepared in acetone. These stock standard solutions were stored in a freezer (−24 • C) when necessary until use.
The standard solutions used to prepare the calibration curve, containing the target analytes from 40 to 400 µg L −1 , were prepared using the four food simulants. Within the working range, a 6-point calibration curve was constructed using the following concentration levels: 40, 80, 160, 240, 320, and 400 µg L −1 . For this, aliquots of the stock multicomponent standard solution, ranging from 50 to 500 µL, were placed into a series of 20 mL volumetric flasks, and they were filled to the line using food simulant as a solvent. These standard solutions were freshly prepared in each analytical sequence.

Vortex-Assisted Liquid-Liquid Extraction (VA-LLE)
After the migration test, 20 mL of the food simulant were collected and placed into a 50 mL glass centrifuge tube, and 5 mL of n-hexane were added. The mixture was stirred vigorously for 30 s in a vortex mixer (ca. 2000 rpm) in order to enhance the extraction effectiveness, and it was then centrifuged at 3000 rpm for 5 min for phase separation. Finally, 1 mL of the n-hexane extract (supernatant phase) was collected, placed into an injection vial, and spiked with 20 µL of the internal standards solution for subsequent GC-HRMS analysis. Additionally, sample blanks consisting of blank food simulant subjected to VA-LLE were also prepared in triplicate following the same procedure.
The VA-LLE was also applied to the calibration curve (see Section 3.4.2) following the described extraction procedure under the same conditions prior to GC-HRMS analysis (procedural standard calibration).

GC-HRMS Analysis
One microliter of the standard or sample extract was injected into the GC-HRMS system (splitless injection). The inlet temperature was set at 280 • C. The GC operated in constant flow mode at 1.2 mL min −1 with helium as a carrier gas, using the following oven temperature programs: 40 • C, held for 5 min; 5 • C min −1 up to 315 • C, held for 10 min. The MS transfer line was set at 300 • C. The EI ion source operated at 70 eV, and its temperature was 250 • C. The acquisition was performed in selected ion monitoring (SIM) mode with a resolving power of 30,000 FWHM and a mass accuracy ≤ 5 ppm. The selected ions for signal acquisition and the retention times of the target analytes are shown in Table S1.

Quantification and Confirmation
The quantification of the target analytes was performed by procedural standard calibration with internal standards in the working range, from 40 to 400 µg L −1 (6-point calibration curve). The calibration curves were obtained by simple linear regression, with the analytical signal (i.e., ratio between the target analyte peak area and the internal standard peak area) being the dependent variable and the analyte concentration (µg L −1 ) being the independent variable. The internal standard bis(2-ethylhexyl) phthalate-D 4 was used to quantify benzyl butyl phthalate, bis(2-ethylhexyl) adipate, bis(2-ethylhexyl) phthalate, diallyl phthalate, dibutyl adipate, dibutyl maleate, dibutyl phthalate, diethyl phthalate, diisobutyl phthalate, dimethyl isophthalate, dimethyl terephthalate, and diphenyl phthalate, and the internal standard phenol-13 C 6 was used to quantify the other target analytes.
The results, expressed in µg L −1 of food simulant, were obtained by interpolation in the calibration curves. In order to evaluate compliance with the SML of Regulation 10/2011 [1], the results needed to be converted to mg kg −1 of food. For this, the conventional ratio of 6 dm 2 of food contact surface per kg of food was considered using the following expression: M = C × (6 V/S) × 10 −3 ; 'M' being the migration of the target analyte in mg kg −1 of food; 'C' the obtained concentration of the target analyte in µg L −1 of food simulant; 'V' the volume of food simulant in L or dm 3 ; and 'S' the contact surface of the FCM in dm 2 .
The identification and confirmation of the target analytes was performed by comparing the chromatographic peaks obtained for the standards and the samples, considering the following criteria: retention time should not differ by more than 0.1 min, the two selected ions should be detected (quantification ion and confirmation ion) (see Table S1) with mass accuracy ≤2 ppm, and the analyte peaks in the extracted ion chromatogram of the two selected ions must fully overlap.

Conclusions
In this work, a useful analytical method for the determination of 60 migrant substances, including aldehydes, ketones, phthalates and other plasticizers, phenol derivatives, acrylates, and methacrylates, from plastic FCM using food simulants has been developed and validated. The proposed methodology included migration tests of plastic FCM using food simulants A, B, C, and D1, followed by VA-LLE using only n-hexane as an extraction solvent and GC-Q-Orbitrap HRMS analysis. Migration tests were performed following standardized conditions, according to current Regulation 10/2011, and the experimental variables involved in the vortex-assisted liquid-liquid extraction were studied. The proposed method was then validated using procedural standard calibration, showing satisfactory linearity, method limits of quantification, precision, and relative recovery; it was applied in the analysis of four field samples made of different plastic FCM, showing that the target analytes were in compliance with their SML according to Regulation 10/2011. In this sense, the analytical method addressed in this paper constitutes a fast, useful, and reliable procedure for the analysis of plastic FCM using food simulants for the control of compliance with current European regulations. It should be noted that, in the current market, there is a wide variety of FCM made of different plastics and polymer materials and not all the studied substances, both IAS and NIAS, may be relevant for all the plastic FCM. In this sense, a multi-analyte method such as the one proposed in the present manuscript, could be useful for food control authorities in order to evaluate compliance of plastic FCM with current European regulations as a part of a comprehensive analytical control plan.
Supplementary Materials: The following are available online. Table S1: Selected ions (m/z) for HRMS acquisition and retention time (RT) of the target analytes, Table S2: Extraction efficiencies of the target analytes from the studied food simulants (A, B, C, and D1), Table S3: Analytical figures of merit of the proposed method using food simulant B, Table S4: Analytical figures of merit of the proposed method using food simulant C, Table S5: Analytical figures of merit of the proposed method using food simulant D1.