Uncertainty Propagation in High-Resolution Mass Spectrometry Non-Target Screening for Food Contact Conformity Verification

Spike

Simulant Contact Protocols and Extraction Recoveries

Analytical measurements for food contact materials begin long before eluate reaches a chromatographic column. Non-target screening by liquid chromatography or gas chromatography coupled to high-resolution mass spectrometry evaluates complex chemical mixtures migrating from polymers, adhesives, coatings, and printing inks. The raw measurement signal reflects not only the concentration of an analyte in the liquid food simulant but also the thermodynamic partitioning and chemical stability of that compound during exposure.

Food simulants prescribed under European Union Regulation 10/2011, such as ten percent ethanol, three percent acetic acid, or vegetable oil, extract migrating species under defined time and temperature regimes. A ten-day contact period at sixty degrees Celsius alters labile migrants through hydrolysis, oxidation, or thermal rearrangement. When a polyolefin film containing ester-based antioxidant breakdown products undergoes migration testing in three percent acetic acid, acid-catalyzed hydrolysis degrades parent molecules into smaller fragment species, shifting the apparent molecular weight distribution observed during mass analysis.

Extraction efficiency across diverse chemical structures introduces the primary source of analytical uncertainty in non-target workflows. Unlike targeted quantitation where isotopically labeled internal standards correct for recovery losses of known target compounds, non-target screening covers hundreds of unknown non-intentionally added substances possessing uncharacterized polarities, partition coefficients, and pKa values. Solutes with high octanol-water partition coefficients adsorb onto glass collection vials, PTFE tubing, or syringe filter membranes during sample transfer steps.

Solid-phase extraction protocols used to concentrate migrants from aqueous simulants like ten percent ethanol display variable retention efficiencies, ranging from under ten percent for highly hydrophilic organic acids to over ninety percent for lipophilic oligomers.

Solid-phase extraction on polymeric sorbents returns recoveries below twenty percent for polar primary aromatic amines extracted from three percent acetic acid at forty degrees Celsius.

Sources of Systematic Extraction Variance

Quantifying recovery bias across unknown migrants necessitates evaluating physical losses occurring during sample pre-treatment. Liquid-liquid extractions using dichloromethane or hexane selectively enrich non-polar compounds while leaving ionic species in the aqueous phase. Evaporation of extraction solvents under nitrogen streams drives off volatile migrant species, including monomer residues such as vinyl chloride, butadiene derivatives, or short-chain aldehyde decomposition products.

The accumulated error from extraction recovery introduces an uncertainty factor of two to ten in final concentration estimates prior to mass spectrometric detection.

• Adsorptive Surface Depletion losses occur when non-polar migrants bind irreversibly to active silanol sites on glassware during solvent concentration steps.
• Thermal Degradation Artefacts generate synthetic signals when heat-sensitive additives decompose inside hot injector ports or during high-temperature simulant exposures.
• Volatilization Losses deplete low-boiling migrants during solvent evaporation, distorting the relative peak area ratios of light compounds against heavier homologues.
• Differential Sorption Efficiency leaves hydrophilic degradation products behind in aqueous food simulants during liquid-liquid extractions using non-polar organic solvents.

Standard additions using isotopically labeled analogues remain impossible when peak identities remain unknown during data acquisition. Surrogate standards spiked into food simulants before exposure provide a benchmark for recovery corrections, yet a surrogate only reflects compounds sharing its precise physical property space. Spiking deuterium-labeled bisphenol A measures recovery for structural analogues but fails to represent the recovery dynamics of cyclic polyester oligomers or slip agent degradation products present in the same migration solution.

Suppliers frequently defend missing or suppressed migrant peaks by asserting that analytical methods validated for target monomers cover all degradation chemistry originating from their resin matrices.

Drift

Mass Accuracy and Resolving Power Stability

High-resolution mass spectrometry relies on precise mass-to-charge ratios to assign chemical formulas to observed molecular ions. Instrument architectures such as Quadrupole Time-of-Flight and Orbitrap mass analyzers achieve resolving powers exceeding thirty thousand and one hundred thousand, respectively. Mass accuracy degrades over long analytical runs due to ambient laboratory temperature fluctuations, power supply instability, and space-charge accumulation within ion traps or flight tubes.

Drift in mass calibration by two parts per million expands the set of plausible elemental compositions for an unknown ion from three candidates to dozens, directly expanding identification ambiguity.

Continuous internal mass calibration using lock masses mitigates systemic drift across analytical batches. Lock mass compounds infused continuously or injected alongside chromatographic peaks correct mass scale shifts in real time. Instrument lock mass signals experience ion suppression when co-eluting high-concentration migrants dominate the ionization space.

A heavy slip additive like erucamide eluting at high concentration saturates detector channels, depressing the lock mass intensity and causing localized mass assignment errors exceeding five parts per million precisely during the elution window of critical unknown migrants.

Where Does Mass Drift Impact Formula Assignment Accuracy?

Determining whether an observed peak originates from a regulated packaging component or an uncharacterized contaminant hinges on isotopic pattern matching and narrow mass tolerances. At a mass tolerance of one part per million, a monoisotopic mass of 285.1122 yields a unique formula of C17H16O4. Expanding the tolerance window to five parts per million due to thermal instrument drift introduces alternate formulas incorporating nitrogen and sulfur atoms, such as C14H17N2O4S.

This shift alters the predicted toxicity profile and regulatory classification of the substance, transforming a benign ester compound into a potential structural alert requiring targeted toxicological evaluation.

Chromatographic retention time drift compounds mass accuracy errors during feature alignment across multi-sample studies. Gradient elution high-performance liquid chromatography experiences slight retention time shifts due to mobile phase evaporation, column temperature variation, or matrix-induced stationary phase modification. When aligning chromatographic features from ten consecutive migration extracts, a retention time drift window exceeding three seconds risks misaligning identical chemical entities or merging distinct isomers into a single composite feature, skewing peak abundance calculations.

Signal intensity drift across analytical runs further complicates semi-quantitative screening. Ion source contamination caused by non-volatile matrix components eluting from food simulants reduces ionization efficiency over time, causing identical migrant concentrations to yield lower peak areas at the end of an analytical sequence compared to the start. Regular injection of quality control pools throughout sequence runs allows software algorithms to model and correct intensity decay across run sequences, provided signal decay remains linear.

A calibrated instrument operated within a controlled environment limits mass error to predictable boundaries over routine sequence durations.

Deconvolution

Algorithm Architecture and Feature Extraction Errors

Raw files generated by mass spectrometers contain millions of individual data points spanning retention time, mass-to-charge ratio, and ion intensity. Deconvolution software converts continuous signals into discrete chemical features representing discrete molecular species. Algorithms must distinguish real analyte signals from background noise, electronic spikes, chemical contaminants in solvents, and mobile phase adducts.

Setting peak detection threshold values too high omits low-concentration migrants present near the regulatory screening threshold of ten parts per billion, whereas setting thresholds too low introduces thousands of false-positive features originating from baseline noise.

Adduct formation, in-source fragmentation, and isotopic distributions convolute feature lists by generating multiple peaks for a single chemical entity. Electrospray ionization yields protonated molecules alongside sodium, potassium, and ammonium adducts. In recycled polyolefin extracts, a single synthetic antioxidant like Irganox 1010 produces protonated ions, sodiated ions, dimer species, and multiple in-source fragment ions.

If deconvolution software fails to group these related signals into a single compound feature, the software reports five distinct chemical species, undercounting the total integrated mass of the parent antioxidant while creating phantom non-intentionally added substances in conformity dossiers.

Standardized evaluation under ISO/IEC 17025 requires documented signal-to-noise ratios exceeding ten to one prior to reporting unassigned chromatographic features in compliance dossiers.

Isotopic Pattern Matching and Formula Exclusion

Isotopic distributions provide a secondary filter for validating formula assignments derived from high-resolution mass measurements. Chlorine, bromine, sulfur, and silicon contain distinctive polyisotopic ratios that constrain chemical composition choices. Carbon-13 natural abundance patterns confirm the number of carbon atoms present in an unknown molecule.

High intensity uncertainty in small chromatographic peaks distorts observed isotopic ratios, leading deconvolution algorithms to reject correct chemical formulas due to theoretical pattern mismatch.

Co-eluting compounds possessing overlapping mass spectra distort background subtraction routines. When a plasticizer like diethylhexyl phthalate co-elutes with an unknown adhesive degradation product, shared fragment ions mask the true spectrum of the lower-abundance migrant. Automatic spectral deconvolution algorithms attempt to unmix composite spectra using ion chromatographic profiles, but highly correlated elution profiles result in incomplete spectral extraction and corrupted library searches against databases like NIST or Wiley.

1. Signal Threshold Definition establishes the baseline noise floor to exclude electronic artifacts while retaining trace migrant peaks above ten parts per billion.
2. Chromatographic Peak Detection fits mathematical models to total ion chromatograms to isolate eluting chemical species from mobile phase noise.
3. Adduct and Isotope Grouping clusters sodium, potassium, and protonated species alongside natural isotope patterns into unified molecular features.
4. Background Subtraction Routine eliminates chemical noise contributions originating from migration cell components and extraction solvents.
5. Library Mass Spectral Matching compares deconvoluted fragmentation spectra against certified reference databases to calculate match confidence scores.

Calculating the true confidence level of a deconvoluted feature requires propagating mass accuracy errors, isotopic fit scores, and fragmentation match metrics into a unified identification probability score. How can laboratories mathematically bound the false-discovery rate of unknown migrants when spectral libraries contain fewer than fifteen percent of all commercially used plastic additives and transformation products?

Surrogate

Semi-Quantitative Response Factor Variance

Determining the concentration of an unidentified non-intentionally added substance requires converting raw chromatographic peak areas into mass fractions. In targeted analysis, pure chemical reference standards generate compound-specific calibration curves that correlate ion current directly with concentration. Non-target screening lacks reference standards for newly discovered migrants.

Laboratories adopt semi-quantitative approaches, applying the response factor of an added internal standard or surrogate compound to calculate concentrations for all unknown features present in the extract.

Ionization efficiency in electrospray mass spectrometry varies by up to three orders of magnitude across different chemical structures. A highly ionizable compound like a tertiary amine produces a signal intensity one thousand times higher than a weakly ionizable compound like a saturated fatty acid ester at identical molar concentrations. Assuming a uniform response factor across unknown migrants introduces substantial quantitative error, causing a laboratory to underestimate the concentration of a poorly ionizable migrant by a factor of one hundred, or overestimate a highly ionizable contaminant by a similar margin.

Uncertainty Propagation Mathematics in Concentration Estimation

Quantifying the combined uncertainty of semi-quantitative non-target screening requires integrating multiple variance terms. The total relative standard uncertainty of a calculated migrant concentration depends on the uncertainty of the extraction recovery, the repeatability of the mass spectrometer peak integration, the uncertainty in the internal standard concentration, and the probability distribution of response factors across the chemical domain under evaluation.

Models utilizing response factor distributions constructed from hundreds of target compounds estimate the probability that an unknown feature’s true concentration lies below a specific regulatory limit. By fitting log-normal distribution curves to response factor databases, analysts determine a conversion factor corresponding to the eightieth or ninety-fifth percentile of ionization efficiency. Applying an upper-percentile response factor yields a conservative concentration estimate, minimizing the risk of false-negative compliance decisions while increasing the rate of unnecessary follow-up investigations for benign migrants.

Ionization response factors across unknown migrants span three orders of magnitude in electrospray ionization, making single-standard semi-quantification a screening estimate rather than an absolute measurement.

The mathematical propagation of uncertainty through a log-normal response factor distribution proceeds through defined sequential steps during data processing.

1. Calculates raw peak area ratios between the unknown non-target feature and the nearest eluting deuterated surrogate standard.
2. Applies the known concentration of the surrogate standard to establish a baseline nominal concentration value.
3. Multiplies nominal concentration by the extraction recovery correction factor derived from matrix validation spikes.
4. Incorporates the log-normal variance parameter of the compound-class response factor model to compute upper confidence bounds.
5. Compares the calculated upper confidence limit concentration directly against toxicological screening thresholds.

Deploying uncalibrated semi-quantitative screening data directly to declare food contact compliance creates immediate legal exposure when border control authorities perform targeted testing with authentic reference standards and discover restricted migrants exceeding legal limits.

Benchmark

Toxicological Threshold Mapping and Schymanski Levels

Interpreting non-target screening data for regulatory compliance requires combining chemical identification confidence with toxicological risk criteria. The Schymanski communication scale categorizes mass spectrometric identification confidence into five distinct tiers. Level 1 represents confirmed structures verified by authentic reference standards with matching mass spectra, retention times, and fragmentation profiles.

Level 2 provides probable structures supported by library spectrum matching or diagnostic fragment interpretation. Level 3 indicates tentative candidate structures sharing specific molecular formulas or sub-structural fragments. Level 4 establishes molecular formulas without structural resolution, while Level 5 defines exact mass features lacking formula assignment.

Regulatory evaluation of unknown migrants relies on the Threshold of Toxicological Concern framework when empirical toxicological data remains unavailable. Under European Food Safety Authority guidelines, uncharacterized substances lacking structural alerts for genotoxicity operate under a default exposure threshold of 0.015 micrograms per kilogram body weight per day, corresponding to a concentration limit of 0.01 milligrams per kilogram in food or food simulants. If an unknown migrant contains structural alerts for genotoxicity, such as aromatic amines, epoxides, or alkyl hydrazines, the toxicological threshold drops to 0.0025 micrograms per kilogram body weight per day, requiring analytical detection thresholds down to 0.15 parts per billion in food simulants.

Integration of QSAR Tools and Cramer Classification

Assigning a Threshold of Toxicological Concern to a Level 2 or Level 3 non-target feature requires running molecular structures through Quantitative Structure-Activity Relationship models like ToXsuite or Derek Nexus. These computational tools screen candidate structures against databases of chemical carcinogens, mutagens, and endocrine disruptors. If an analytical feature maps to five candidate structures at Level 3, and one candidate contains a structural alert for mutagenicity, the analyst must apply the most conservative toxicological threshold across all plausible candidates.

Uncertainty in identification confidence propagates directly into toxicological evaluation. Misinterpreting an isotope peak as a protonated molecular ion leads to an incorrect molecular weight calculation, which subsequently generates false structural inputs for toxicological alert screening software. A false negative during structural alert screening assigns a high toxicological threshold of 0.05 milligrams per kilogram to a potent carcinogen, allowing a hazardous migrant to pass compliance evaluation undetected.

Mandatory documentation requirements govern the submission of non-target screening data for compliance dossiers supporting declarations of conformity.

• Mass Accuracy Traceability Records document daily instrument calibration routines and internal lock mass stability across analytical batches.
• Schymanski Tier Classifications assign explicit confidence levels to every reported chemical feature exceeding screening thresholds.
• Ionization Mode Complementarity Data demonstrates acquisition using both positive and negative electrospray modes alongside gas chromatography mass spectrometry for non-polar fractions.
• Extraction Blank Profiles identify matrix spikes, solvent impurities, and laboratory system contaminants to prevent false non-intentionally added substance assignments.

Article 16 of Regulation EC 1935/2004 demands that supporting documentation demonstrate compliance through test reports backed by verified analytical methodologies.

Discrepancy

Inter-Laboratory Variability and Legal Defensibility

When two independent accredited laboratories analyze identical samples of migrating packaging film using high-resolution non-target screening, reported chemical feature lists frequently diverge. Discrepancies arise from differences in sample preparation protocols, instrument front-end geometries, chromatographic stationary phase selectivity, deconvolution algorithms, and spectral library coverage. Laboratory A running a liquid chromatography time-of-flight mass spectrometer with a methanol-water gradient identifies forty distinct migrant features, whereas Laboratory B operating an Orbitrap system with an acetonitrile-water gradient reports sixty-five features from the same simulant extract.

Variability between analytical facilities creates significant commercial friction during cross-border trade inspections. A regulatory authority performing official controls may flag an unlisted non-intentionally added substance exceeding ten parts per billion based on an in-house non-target screening protocol, while the manufacturer’s supporting dossier contains a test report from a contract laboratory showing no signals above background noise. Resolving these conflicting measurements requires establishing standardized workflows for uncertainty estimation that account for instrument-specific variance, ionization suppression, and algorithm settings.

Quantifying measurement uncertainty in non-target screening transforms raw mass spectrometry counts into defensible legal evidence. Combining individual uncertainty components, extraction recovery, mass measurement drift, peak deconvolution variance, and response factor spread, yields an expanded measurement uncertainty that routinely reaches fifty to two hundred percent of the estimated concentration. A reported non-target concentration of eight parts per billion possessing a two hundred percent expanded uncertainty spans a true value range from two point six to twenty-four parts per billion, straddling the regulatory action limit of ten parts per billion.

Regulatory authorities faced with data straddling compliance limits evaluate the statistical probability of non-compliance based on confidence intervals. Importers placing packaging materials on the market bear the responsibility of demonstrating that migrant concentrations remain safely below legal limits even at the upper bound of analytical uncertainty. Downstream brand owners increasingly demand that converted films and multi-layer laminates come backed by non-target screening reports that explicit report lower and upper quantitative confidence intervals rather than single-point semi-quantitative values.

Standardizing non-target screening protocols across enforcement laboratories requires establishing certified reference materials containing known mixtures of non-intentionally added substances at trace concentrations. Synthetic mixtures containing cyclic oligomers, degradation products of secondary antioxidants, and ink photoinitiators allow testing facilities to benchmark feature detection thresholds, formula assignment accuracy, and semi-quantitative response models against standardized ground truth data. Method validation frameworks developed under European Committee for Standardization protocols aim to harmonize reporting formats, signal-to-noise calculations, and identification confidence assignments, ensuring that non-target screening results stand up to judicial review during market enforcement actions.