Framework to Estimate Total Particulate Mass and Nicotine Delivered to E-cig Users from Natural Environment Monitoring Data

A framework describing the joint effect of user topography behavior and product characteristics of one exemplar device on the total particulate mass (TPM) and aerosol constituent yield delivered to a user is presented and validated against seven user-specific ‘playback’ emissions observations. A pen-style e-cig was used to collect emissions across puff flow rates and durations spanning the range observed in the natural environment. Emissions were analyzed with GC-MS and used to construct empirical correlations for TPM concentration and nicotine mass ratio. TPM concentration was demonstrated to depend upon both puff flow rate and duration, while nicotine mass ratio was not observed to be flow-dependent under the conditions presented. The empirical model for TPM and nicotine yield demonstrated agreement with experimental observations, with Pearson correlation coefficients of r = 0.79 and r = 0.86 respectively. The mass of TPM and nicotine delivered to the mouth of an e-cig user are dependent upon the puffing behavior of the user. Product-specific empirical models of emissions may be used in conjunction with participant-specific topography observations to accurately quantify the mass of TPM and nicotine delivered to a user.


Approach
We define the cumulative yield, Y AC , of individual Aerosol Constituents (AC) of emissions delivered from an ENDS to the mouth of a user as the integral of the product of the time dependent mass ratio of the aerosol constituent, the Total Particulate Matter (TPM) concentration of the whole aerosol, and the user's volumetric flow rate: The constituents in the ENDS aerosol may be present in the un-puffed e-liquid, or generated as decomposition products. We normalize all aerosol constituents (including vapor phase constituents, compounds originating in the e-liquid, and thermal decomposition products) by the mass of TPM emissions to facilitate separation of variables between the f AC and C TPM terms. We posit the TPM concentration and mass ratio of constituents may be expressed as linearly independent functions of Product Characteristics (PC) and User Behavior Characteristics (UBC). Numerous ENDS PC may affect the C TPM including but not limited to the device operating power (reflected in coil wattage, amperage, or temperature), flow path geometry, coil design, and aspiration features. Additionally, the solvent composition of the e-liquid consumable (such as the PG/VG ratio which directly impacts the saturation temperature of the e-liquid) may impact C TPM . Additional consumable PC impacting f AC may include nicotine concentration, flavor additives, viscosity, and pH. Furthermore, a variety of UBC may affect the C TPM and/or f AC including puff duration, d, flow rate, q, volume, v, and interval, i: AC AC For the current study we limit variability in PC by selecting a single ENDS with no user-adjustable settings and a single e-liquid. We thus focus on the interaction between UBC and the flow path PC, reflected by the topography parameters q and d, and consider a single AC, nicotine, to illustrate the approach. Prior work 4 The experimental observations of C TPM can be computed as the ratio of the mass emissions captured on a filter pad per measured volume of aerosol passing through the pad. The coefficients in Eq. 4 may be estimated using ordinary least squares (OLS), weighted least squares (WLS), or other regression techniques. OLS regression is employed in the current work. Since the ENDS device chosen for the study does not have a user-selectable power setting, and the nicotine concentration of the e-liquid is held constant across all trials, we hypothesize a first order linear model, Eq. 5, for the model-predicted nicotine mass ratio, f NIC , as a function of puff flow rate, q.

Methods
Aerosol generation, capture and analysis. The ENDS used for the current study was an NJOY vape pen e-cigarette with a top-fill (top-coil) e-liquid tank and no user adjustable settings such as power or coil temperature. The same model of ENDS was used for the emissions testing and natural environment monitoring of all results presented here. Similarly, the e-liquid used was AVAIL brand Tobacco Row (TR) having a manufacturer's labeled nicotine concentration of 1.8% (measured using GC-MS to be 0.0140 ± 0.00014 [mg/mg] in the un-puffed e-liquid) and a measured solvent mixture ratio of 50:50 propylene glycol to glycerin. The nominal coil resistance of the vape pens used in the study ranged from 2.3 ± 0.3 (95% CI) ohms across nine coils prior-to and following the test series. The battery was recharged at the beginning of each test series, so device power was comparable between experimental conditions and repeated trials.
The PES-1 Programmable Emissions System was used to generate and collect aerosol emissions. The machine can be configured to collect particulate, vapor or liquid phase emissions from ENDS. A vacuum box can accommodate pressures as low as −25 [kPA], and hold sampling bags up to 5.0 [L] for mixed gas/liquid phase collection. The current study used a particulate phase collection mode with Cambridge style filter pads. The machine is driven by a vacuum pump (Model DOA-F704-AA, GAST) and flow is controlled by a proportional valve with response time of 10 milliseconds and range of 0-20 SLPM (KPIH-VP-20-156-25, Kelly Pneumatic Inc.). The system flow rate is measured using a precision gas flow meter (M-50SLPM-D-30PSIA/5M, Alicat Scientific, Inc.) third-party calibrated over its rated range of 20 SLPM. The emissions system is operated under direct digital control, enabling operation across a variety of flow conditions including the instantaneous puff profile, puff volume, flow rate, duration, and interval so that the machine can reproduce homogeneous (repeated) puff profiles or mimic heterogeneous (non-repeating) human puffing topographies for a range of behaviors.
The Total Particulate Matter (TPM) collected on filter pads following each emissions trial was determined gravimetrically using a Mettler AE240 Analytical Balance, with a protected weighing cell with a measurement accuracy of 0.0001 [g] and a repeatability of 0.1 [mg].
Nicotine mass ratio was determined using methods similar to those previously used in our lab 13 . In brief, pads were spiked with quinolone as an internal standard and submerged in methanol and broken up using orbital and wrist shaking followed by filtering through a 0.045 μm cellulose filter prior to GCMS analysis. In addition, calibration standards solutions of nicotine, that span the concentration range observed, and quinolone internal standard were prepared. The standard solutions were subjected to the same procedure as all samples, including introduction of a filter pad in each. Concentrations of nicotine were determined using a Shimadzu 2020 GCMS equipped with an AOC 20i autosampler. Triplicate runs of all samples and standards were run, 5 μL samples and standards were sequentially injected into the GC and passed through a Restek Rxi-5 ms fused silica column (30 m × 0.25 mm I.D. × 0.25 μm). Helium carrier gas was run through a split injector (50:1 split) at a temperature of 230C. The GC oven temperature was increased from 60C to 200C at 20C/min and then held for 3 minutes. The MS source and transfer line were kept at 180C and 280C, respectively, and the MS was run in the single-ion-monitoring mode (SIM) at m/z values of 133 (quantitation) and 162 (confirmation) for nicotine and 102 (quantitation) and 129 (confirmation) for quinolone. Peaks were integrated and the ratio of nicotine's integration to its corresponding surrogate's integration was used to establish calibration curves and determine nicotine concentrations in the samples. The mass of nicotine measured in each sample was divided by the mass of TPM measured on the corresponding filter pad to compute the measured nicotine mass ratio, f NIC .
Method for quantifying the model coefficients. The empirical model coefficients for Ĉ TPM (Eq. 4) and f NIC (Eq. 5) were quantified by characterizing the ENDS and e-liquid over the range of puffing behaviors as observed in the natural environment [14][15][16] . Figure 2 illustrates the design of experiments used for the emissions testing experiment. The emissions testing began with a 'screening' experiment 4 of 10 flow conditions wherein the puff duration and puff flow rate were varied inversely to maintain a nominally constant puff volume and cumulative session volume. The previous study 4 presented screening results for cig-alike, vape-pen, box-mod and pod-style electronic cigarettes under varying conditions of coil location, operating power, and e-liquid flavors. The current work further investigates the joint effects of puff flow rate and puff duration on emissions for a single device, and thus controlling coil location, operating power and e-liquid. Using knowledge from the screening experiment, 24 additional 'full' flow conditions were tested wherein puff flow rate and duration were varied. Each flow condition was repeated for six trials, resulting in a total of 204 trials across 34 flow conditions. Each trial consisted of a fixed number of repeated puffs such that the average cumulative aerosol volume generated per trial was approximately 740 [mL], ranging from 522 [mL] to 1000 [mL]. All results are presented in terms of observed actual flow conditions as opposed to the nominal programmed command conditions. The TPM mass was measured and divided by the cumulative observed aerosol volume associated with each trial to determine the value of C TPM for each trial at each flow condition. The resulting observations of C TPM as a function of q and d are used in conjunction with an ordinary linear least squares algorithm to estimate the coefficients … b b 1 6 appearing in Eq. 4. The filter pad from selected 'screening' trials were analyzed using GC-MS to determine the mass ratio of nicotine to TPM in conjunction with a linear least squares algorithm to estimate the coefficients β 1 and β 2 appearing in Eq. 5.
Method for validating the model. www.nature.com/scientificreports www.nature.com/scientificreports/ One natural environment puffing session 'playback' data set was selected for each of the seven participants in a previously reported two week switching study 14 who elected to use e-liquid with a nicotine concentration of 1.8%. Each session listed in Table 1 was selected as representative of the topography behavior exhibited by the participant during six full days of ENDS use with TR e-liquid. Figure 3 illustrates exemplar puffing session data collected for seven participants in the natural environment during a recent two-week switching study. These observations were used to generate 'playback' command profiles to generate emissions and estimate the participant-specific yields of TPM and NIC.
Validation of the model consisted of the following sequence of operations.
1. Seven sessions representative of natural environment topography data (Table 1 and Fig. 3) were used to generate 'playback' puffing profiles as input to the PES-1 system. 2. The emissions system was operated with the subject ENDS device and e-liquid to mimic the behavior of each study participant. Each 'playback' puffing profile was repeated for six trials, with the device battery regularly recharged and the e-liquid tank confirmed to be greater than 1/2 full at the beginning of each trial. 3. The filter pads were measured after each of the 42 trials to determine the experimentally observed TPM yield, Y TPM . 4. Each filter pad was analyzed using GC-MS to determine the experimentally observed nicotine yield, Y NIC . 5. Apply the empirical model for Ĉ TPM (Eq. 4) to each 'playback' puff, based on the actual q and d measured in the PES-1, and sum the individual puff results to estimate the TPM yield for each 'playback' session, wherein N puff is the number of discrete puffs during each session:

Results
Empirical coefficients in the framework. Figure 4 shows the results of the emissions study conducted using NJOY e-cig operated across a wide range of puff flow rates and puff durations, using Avail brand Tobacco Row flavored e-liquid with a manufacturer's labeled nicotine concentration of 1.8%. Each trial is represented by a sphere, and the shaded surface illustrates the empirical model used to describe the C TPM as a function of puff flow rate, q, and duration, d.  with a root mean square error of 0.012 [mg/mL]. Figure 5 shows the results of the GC-MS analysis for nicotine mass ratio present in the TPM emissions samples. Also shown is the empirical curve fit for the nicotine mass ratio as a function of puff flow rate. The linear least squares regression coefficients in the empirical model for f NIC (Eq. 5) are shown in Table 3. The coefficient of determination for the model is R 2 = 0.0002, reflecting the significant scatter in the observations, with a root mean square error of 0.0025 [mg/mg].
Validation results. Figure 6 illustrates the correlation of the empirical model estimates of Ŷ TPM and Ŷ NIC versus machine puffed playback measurements for six repeated trials of seven exemplar playback sessions. The slope for Ŷ TPM is m = 1.08, R 2 = 0.58, and Pearson's correlation coefficient of r = 0.79. The slope for Ŷ NIC is m = 0.904, R 2 = 0.71 and a Pearson's correlation coefficient of r = 0.86. The scale of the two validation figures varies by a factor of 100 due to the circumstance that f NIC = 0.01 and is independent of flow rate for the e-cigarette and e-liquid tested here. The difference in the r value reflects experimental differences in the observations of Y TPM and Y NIC .

Discussion
Natural environment observations of user topography behavior have been used to create 'playback' command profiles and demonstrate estimates of user-specific exposure to harmful and potentially harmful constituents. The mass of particulate matter and nicotine delivered to the mouth of an electronic cigarette user are dependent upon the puffing topography behavior of the user. The total particulate matter concentration and nicotine mass ratio may be determined as a function of product characteristics and user behavior characteristics and used to create a predictive emissions model specific to each tobacco product. The predictive emissions model may be used in conjunction with user-specific topography observations to estimate the mass of particulate matter and aerosol constituents delivered to the mouth of an individual user. The proposed framework may be used to assess the effect of e-liquid and e-cigarette product characteristics on emissions. The proposed framework may be used to  www.nature.com/scientificreports www.nature.com/scientificreports/ examine the relationship between product characteristics, user behavior characteristics, and user-specific exposure to whole aerosol mass and individual constituents there of.
A framework was introduced and validated. Total Particulate Matter (TPM) emissions were collected from a machine puffing system across the range of puff duration and puff flow rate typical of observed natural environment behavior. The TPM concentration was demonstrated to be a function of both puff flow rate, q, and puff    www.nature.com/scientificreports www.nature.com/scientificreports/ duration, d. The nicotine mass ratio was determined across the range of puff duration and flow rate and found to be relatively independent of flow conditions for the products studied. The emissions data was used to create empirical models for Ĉ TPM and f NIC . The model predictions of Ŷ TPM and Ŷ NIC yield were validated against experimental observations of Y TPM and Y NIC representative of natural environment behavior of ENDS users. The framework was demonstrated to be predict the yield of machine 'playback' emissions profiles with Pearson's correlation coefficients of r = 0.79 (TPM) and r = 0.86 (NIC).
The framework provides a context to enhance the understanding of the functional relationship between topography and consumption behavior and product characteristics of ENDS and e-liquids to estimate the total particulate and HPHC constituent yield delivered to a user. The framework defines the constituent mass ratio and TPM concentration as a means of separating the influence of ENDS flow path and power characteristics from the composition of e-liquids, offering the potential for more robust data sharing between research groups and enabling data re-use. The ability to quantify the yield of TPM and HPHC delivered to an ENDS user enables the investigation of a causal relationship between the time-history product use and observed bio-markers of exposure (indicators of uptake). The approach may lead to the ability to predict a-priori the effects of changes in product characteristics on resulting emissions, yield, and uptake. It is recommended that TPM emissions test results be reported as the mass of TPM per volume of whole aerosol. It is recommended that HPHC emissions test results be reported as the ratio of mass of the HPHC per mass of TPM. The actual flow conditions, device parameters and e-liquid composition should be reported with all emissions and yield data.
The data and analysis presented in this manuscript are in terms of the "puffed aerosol" collected under various topography (puff flow rate and duration) conditions. The study employed e-liquids directly from manufacturer's retail bottles, drawn from the same production lots as the e-liquids provided to participants in the preceding human subject study 14 . The emissions test results (Figs 4 and 5) and the experimental validation results (Fig. 6) were all conducted by refilling the vape pen tank from manufacturer's packaging, reflecting the actual usage of the previous study participants.
While the approach is generally applicable to a wide variety of ENDS and e-liquids, the numerical results are not. The coefficients computed for the empirical functions describing Ĉ TPM (q, d) and f NIC (q) are specific to emissions collected with an NJOY e-cigarette filled with AVAIL brand Tobacco Row e-liquid having a labeled nicotine concentration of 1.8% and a solvent mixture ratio of 50:50 propylene glycol to glycerin. The functional form of Ĉ TPM (q, d) and f NIC (q) are anticipated to be device and consumable dependent, particularly in the case of ENDS with adjustable power and/or flow path settings. The single model of vape pen presented in this study has no user-selectable power adjustment options, so any variation in device power which may have been induced by variations in coil resistance between test articles and over time are reflected in the repeated trial variability present in the reported experimental data. A full investigation of the electronic principles of operation and power management of electronic cigarettes is beyond scope of the current work and is deferred to a future study.