Have a personal or library account? Click to login

Comparison of the Particle Size Distribution and Vapor Phase of Electronic Nicotine Delivery Systems Using Two Impactors

Open Access
|Sep 2023

Full Article

1.
INTRODUCTION

Partitioning of Electronic Nicotine Delivery Systems (ENDS) aerosol constituents between the liquid particle and vapor phase affects all aspects of their evaluation and measurement. The dynamic evolution of ENDS aerosols was recently highlighted and it was noted that aerosol dilution and residence time were critical parameters for measurement of their particle size distribution (1). Aerosol dilution and residence time primarily affect Brownian diffusion that is partially responsible, in addition to hydroscopic growth (condensation of water), for particle size growth of dense ENDS aerosols. This is consistent with the theoretical and experimental work demonstrating differences in the particle size distribution measurements among various aerosol instruments (cascade impactors, low pressure impactors, and optical particle counters) due to evaporation and condensation that occur within these aerosol instruments (2).

Particle size distribution measurement of concentrated, chemically complex, and temporally dynamic aerosols (e.g., ENDS products or mainstream/sidestream tobacco smoke) is difficult and requires compromises regardless of the method used (1, 3). Optical particle sizing methods use a light source (infrared, laser, etc.) to scan a small volume of aerosol that absorbs/scatters the light which is collected by one or more detectors. The pattern of light absorption/scattering is used to calculate the particle size distribution (4). As such, optical particle sizing methods do not directly interact with the aerosol and there are several challenges in using them reliably. Concentrated aerosols can present a unique challenge since coincidence counting is not only a function of particle concentration (number of particles per cm3) but also of the wavelength of the light source (in general shorter wavelength light has better resolution for the same volume). Coincidence counting occurs when two aerosol particles are too close together and are counted as one particle. Previous reports on ENDS aerosol particle concentrations have ranged from 1×109/cm3 to 4.39 × 109/cm3 (5, 6). This range of particle concentrations results in a coincidence counting concern that can only be resolved by diluting the aerosol sample. Without diluting the aerosol, this range of particle concentrations also can result in cloud behavior (7,8,9,10) within the small sensing volume of optical equipment that further complicates particle size distribution measurements. Another challenge in using optical particle sizing methods to measure particle size distribution arises from using the aerosol refractive index in the calculation. ENDS aerosols are temporally dynamic, which means depending on the relative humidity (RH), the humectants (propylene glycol and glycerin) are absorbing water (hygroscopicity) while simultaneously semi-volatile constituents are evaporating from the aerosol particles. Both processes, hygroscopicity and evaporation can change the refractive index of the aerosol resulting in a change in the particle size distribution. When comparing different ENDS products, the different chemical constituents will yield differences in the aerosol refractive index, and this will result in slight differences in the calculated particle size distribution among the ENDS products.

Although inertial particle sizing methods (e.g., cascade impactors, etc.) do directly interact with the aerosol by using the aerosol particle's inertia for size fractionation, there are also challenges to their use for measurement of the particle size distribution of ENDS aerosols. Cascade impactors are not generally affected by particle concentration unless they are overloaded and are not affected by refractive index, however they do require a constant airflow (4). One challenge is to match/interface the constant airflow of the cascade impactor with the ENDS puffing regimens. Unless numerous ENDS devices are sampled simultaneously some aerosol dilution is required to match the cascade impactor flow. Since cascade impactors operate under vacuum with a constant airflow, as the particle size cutoff of each impactor stage decreases, there is also a corresponding higher negative pressure, especially for low-pressure impactors (2). Since ENDS aerosols contain various semi-volatile constituents that have different vapor pressures (e.g., nicotine, propylene glycol, glycerin, and flavors) losses due to volatilization from the various impactor stages could occur during and after collection. Most impactor stages must be analyzed manually after aerosol collection either by gravimetric or chemical techniques, which could result in additional losses of very volatile constituents. One cascade impactor, Electrical Low Pressure Impactor (ELPI®+), electrically charges the aerosol particles and with electrometers on each stage counts the number of particles in real time, removing the need for post gravimetric or chemical analysis, however density of the aerosol is required to convert the measured count median diameter to a mass median diameter that is provided directly by other cascade impactors.

Recent work has attempted to identify (11), quantify (12, 13) and determine the evaporation time course (14) for some of the volatile constituents of ENDS aerosol. Berenguer et al. (11) identified 34 volatile compounds in the ENDS aerosol they tested using standard analytical chemistry techniques (Headspace solid phase extraction + gas chromatography-mass spectrometry). Lalo et al. (12) used mass recovery efficiency from cascade impactor measurements to estimate the amount of volatilized nicotine which increased when the power was switched from low to high on the specific ENDS product they used. Luo et al. (13) used an impactor stage to focus particle deposition on a transparent adhesive surface where digital pictures were recorded using a holographic microscope (2 frames/s) and reported rapid evaporation of propylene glycol that changed with the addition of nicotine or several flavors. Finally, laser-trapping (counter-propagating optical tweezers) and single particle Raman scattering was used to investigate the time course of nicotine and propylene glycol volatilization as well as water condensation on ENDS aerosol as a function of e-liquid pH (6). This work (14) reported that propylene glycol and nicotine (from basic solutions) evaporation, as well as water absorption started immediately and was completed in approximately 20–30 s when exposed to 65% relative humidity and 20° C air (1.95 s after generation).

As summarized by recent work (1), past studies measuring the particle size distribution of ENDS aerosols have used different ENDS products, measured the particle size distribution at different times after generation due to the use of different aerosol sizing equipment that has required various dilutions of ENDS aerosols (1, 3, 5, 6, 15,16,17,18,19,20,21,22,23). The different aerosol sizing equipment also used different measurement principles (inertial cascade impactors; light scattering - optical particle counters; electrical - mobility diameters), so when combined with different ENDS aerosol dilutions, different ENDS products and time post aerosol generation, it is not surprising that different particle size distributions have been reported. Based upon this past work, it seems clear that the temporal and chemical dynamics of ENDS aerosol requires consideration of the volatile constituents when measuring the particle size distribution.

Another recent publication (24) noted that the impact of the volatility of propylene glycol and glycerin in ENDS aerosol on the performance of the cascade impactor systems has not been reported in the literature. In a few studies measuring the particle size distribution of ENDS aerosols, an indirect measure of volatile constituents has been reported (1, 3, 5, 12). For example, when a low flow cascade impactor was used, the reported collection efficiency ranged from 74% to 91% indicating that, depending on the specific ENDS product tested, between 9% and 26% of the generated aerosol mass was either deposited in the sampling apparatus or was in the vapor phase and not collected in the cascade impactor (3). In another study, a collection efficiency of approximately 90% was obtained using a similar low flow impactor with an e-liquid formulation containing only a mixture (40/60) of propylene glycol and glycerin (1).

The determination of ENDS particle size is necessary to understand product performance and to comply with the regulatory requirements in some locations. The most widely used method for determining particle size of ENDS aerosol is with a cascade impactor and both the Mini MOUDI™ and ELPI®+ have been used in recent publications (1, 24). Based upon the above considerations, ENDS aerosol particle size distribution was determined using two different cascade impactors (hereafter referred to as impactors) each requiring a different dilution. For each impactor setup the average residence time post aerosol generation was determined, and the volatile component of ENDS aerosol was collected as part of the particle size distribution measurement with volatile chemicals (some quantitatively and others qualitatively) subsequently measured off-line.

2.
MATERIALS AND METHODS
2.1
Products tested

The two products tested (device + pods) were JUUL ENDS Mango flavor 5.0% nicotine by weight (w/w) and JUUL ENDS Mint flavor 5.0% nicotine by weight. This study started in 2020 shortly after Juul Labs, Inc. suspended commercial sales of all Mango (October 2019) and Mint (November 2019) flavored pods in the U.S. market. All pods containing the Mango- and Mint-flavored e-liquid were packaged in U.S. commercial packaging and were 11.5 to 12.5 months post manufacturing when testing was completed using the MSP-135-8 Mini MOUDI™ and ELPI®+, respectively. The primary constituents of both products are United States Pharmacopeia grade glycerin, propylene glycol, nicotine and benzoic acid combined with proprietary flavors. All devices were fully charged according to the manufacturer's instructions prior to testing.

2.2
Particle size distribution

The mass median aerodynamic diameter (MMAD) and geometric standard deviation (GSD) of the two products were determined using the MSP-135-8 Mini MOUDI™ impactor (MSP Corporation, Shoreview, MN, USA; currently a division of TSI® Shoreview, MN, USA) and the Electrical Low Pressure Impactor (ELPI®+; Dekati Ltd., Kangasala, Finland). The MSP-135-8 Mini MOUDI™ impactor has eight impaction cutoff sizes (10.0, 5.6, 3.2, 1.8, 1.0, 0.56, 0.32, and 0.18 μm), while the ELPI®+ has 14 impaction cutoff sizes (5.37, 3.66, 1.63, 0.911, 0.591, 0.315, 0.17, 0.097, 0.0496, 0.033, 0.0192, 0.0141, and 0.006 μm) and both have a filter to collect particles smaller than their respective minimum impaction cutoff size. The MSP-135-8 Mini MOUDI™ impactor uses replaceable aluminum substrates and a glass fiber filter. The ELPI®+ uses sintered metal plates with a small amount of vacuum oil on each plate and a pleated HEPA filter. The ISO puffing regimen for ENDS products (25) (55 mL over 3 s with 30 s interpuff interval) was used for both products and impactors.

To minimize potential dilution of the ENDs aerosol, the MSP-135-8 Mini MOUDI™ impactor was used as previously described (3). This method introduces the ENDS puff into a bias flow of ultra-zero air (2 ppm water; < 1% RH; Airgas) controlled by a mass flow meter (HFC-202, Teledyne Hastings Instruments, Hampton, VA, USA) for a total impactor flow of 2 L/min (3). In this setup, during the 3 s-puff, the 18.3 cm3/s of ENDS aerosol is mixed with an additional 15 cm3/s of ultra-zero (dilution) air resulting in 1 to 0.82 nominal dilution. The internal volume of the connection from the ENDS product to the impactor (3/8 inch Swagelok® “T”) was 1.1 cm3 resulting in a residence time of 0.06 s. In order to not exceed the maximum stage loading, the first five puffs using the ISO ENDS product puffing regimen (25) from three replicates of each product were used to determine the particle size distribution data. Previous studies with the same and similar JUUL products had demonstrated that the particle size distribution was independent of puff regimen (non-intense or intense puffing regimen), consistent throughout the pod life and over a twelve month controlled storage period (See Supplemental Tables 1–6). After completion of the fifth puff, flow through the impactor was maintained for 3–5 s so the aerosol from the fifth puff was pulled through the impactor. The total sampling time was 123–125 s. The size-fractionated particles deposited on pre-weighed 35mm aluminum foils (MSP Corporation) at each stage and filter were weighed post-sampling (Nimbus® analytical balance NBL-124i, Adam Equipment, Inc., Oxford, CT, USA) with the MMAD and GSD determined using the ISO standard (26). The collection efficiency (total mass on impactor stages and filter divided by pod mass loss) of each particle size distribution measurement was also determined. The mean of the three replicates of each product was used to determine the particle size distribution data.

Since the ELPI®+ has aerosol concentration limits due to stage mass loading limitations and requires approximately 10 L/min flowrate (ELPI®+ calibrated by Dekati Ltd. at 9.930 L/min), a custom dilution system was constructed that minimized the residence time by minimizing the total volume of the dilution system. The custom dilution device was assembled by connecting two 3/8 inch diameter Swagelock® stainless steel “Y”s (Figure 1).

Figure 1.

ELPI®+ custom dilution system consisting of two 3/8″ diameter Swagelock® stainless steel “Y”s with airflow connections.

The internal volume of the custom dilution system was measured by quantifying the amount of water required to fill it, resulting in a volume of 27.9 cm3. One end of the first “Y” was connected to the puffing machine (Burkhart Programmable Dual Syringe Pump, PDSP TE, Cerulean, MPRD Ltd., Milton Keynes, UK) and the other end was open to the laboratory environment. One end of the second “Y” supplied ultra-zero air (2 ppm water; < 1% RH; Airgas) as dilution controlled with a mass flow meter and the other end was connected to the ELPI®+ impactor. Although the ISO ENDS product puffing regimen was used (55 mL over 3 s with a 30 s interpuff interval) by the puffing machine, exhaust from the puffing machine into the dilution system occurred over 8 s. For each replicate run, the ELPI®+ and ultra-zero air dilution were set with mass flow meters and the air laboratory air entering the system was measured with flowmeter (Series 4000 flowmeter; TSI, Shoreview, MN, USA) to determine the amount of the puff sampled by ELPI®+. To make the amount of puff sampled calculation more accurate, aerosol deposition within the puffing machine was quantified by measuring the aerosol mass pre and post-puffing machine resulting in 29.045% deposition within the puffing machine that was used to adjust the amount of puff sampled. At the end of the 8 s puff, the exhaust/room air leg of the dilution system, as well as, the connection between the two “Y”s would be full of ENDS aerosol. The percent of puff sampled was calculated as a range including and excluding the volume of the exhaust/room air leg of the dilution system, which was measured at 6 cm3. Between puffs (22 s) the measured amount of room air entered the dilution system to provide the ELPI®+ with the required (approximately 10 L/min) flowrate. Based upon the required ELPI®+ dilution, ENDS puff sampling rate and limits of detection of the analytical methods for vapor phase constituents, a total of 100 puffs for each ENDS product replicate measurement were performed. The total sampling time was 49.5 min. Because the ELPI®+ counts the number of charged particles that land on each impactor stage, to determine the MMAD, GSD and collection efficiency in a similar manner as was done for the Mini MOUDI™ impactor, required using the ELPI®+ software to provide aerosol concentration in mg/m3 and then converting the mg/m3 output of the ELPI®+ to mg using the total amount of air sampled. An aerosol particle density of 1.15 g/cm3 was used by the ELPI®+ software to provide aerosol concentration in mg/m3 since the two most abundant primary ingredients were glycerin (1.26 g/cm3) and propylene glycol (1.04 g/cm3). After conversion to mg, the MMAD and GSD were determined using the ISO standard (26). Recovery efficiency of each particle size distribution measurement was calculated using the mass sampled by the ELPI®+ divided by the percent of puffs sampled, which then was divided by the pod mass loss. The ENDS devices were oriented at a 45° angle to gravity for puffing. A total of six replicate measurements were performed using the ELPI®+ for each ENDS product.

2.3
Vapor phase measurement

During particle size distribution measurements using the MSP-135-8 Mini MOUDI™ impactor, the 2 L/min total vacuum flow was split using mass flow controllers. A flow of 100 mL/min using a mass flow controller (MKS 179 MassFlo®, Andover, MA, USA) was pulled through a custom in-house constructed (Enthalpy Analytical, Richmond, VA, USA) thermal desorption tube (Camsco, Houston, TX, USA) using Tenax GR, Carbotrap B and Carboex 569 sorbents (Sigma-Aldrich, St. Louis, MO, USA) and subsequently analyzed using EPA and ISO Gas Chromatograph/Mass Spectroscopy (GC/MS) methods (27, 28) and the rest of the flow was exhausted to waste. Total amount of volatile compounds collected was based upon a ratio to an internal standard. Since the two ENDS products were tested on separate days, one room air blank sample was collected each day to correct samples for background levels of volatile compounds.

During particle size distribution measurements using the ELPI®+, the approximately 10 L/min vacuum flow was collected using a trap cooled with liquid nitrogen placed in line between the rough pump and ELPI®+. The trap was allowed to equilibrate for 30 min prior to volatile compound sampling. After particle size distribution samples were collected, the vacuum was shut off, one side of the trap was closed down completely, and the other side was disconnected from the ELPI®+, so a custom thermal desorption tube could be inserted. The trap was allowed to come to room temperature, allowing the contents to be trapped onto the custom thermal desorption tube at 200 mL/min flowrate controlled by a mass flow controller (MKS 179 MassFlo®). To ensure complete sampling, a small vacuum pump was attached to the end of the thermal desorption tube and vacuum was drawn for approximately 1 min before capping the thermal desorption tube. The thermal desorption tubes were subsequently analyzed using EPA and ISO Methods (27, 28).

The EPA and ISO methods (27, 28) provide quantitative determination of 86 chemical compounds and semi-quantitative determination of another 370 chemicals (See Supplemental Tables 7 and 8, respectively, for list of compounds). A room air blank was collected after the particle size distribution and vapor phase measurements had been collected. Chemicals identified in the room air blanks, regardless of amount, were removed from chemicals reported for the respective vapor samples. All particle size distribution and vapor collections and measurements were conducted in a laboratory with a controlled temperature range of 21–22° C and RH of 56–60%.

The aerosol generation, particle size distribution determination, and volatile chemical sampling were performed by Enthalpy Analytical, Inc. (Durham, NC, USA, division of Montrose Environmental Group) and the volatile chemical analysis was performed by Prism Analytical Technologies (Mt. Pleasant, MI, USA, division of Montrose Environmental Group). The aerosol generation, particle size distribution determination, and volatile chemical analytical methods were validated and were on the respective contract research organizations ISO 17025 scope of accreditation (29) at the time the study was performed.

3.
RESULTS

There were small differences between the ELPI®+ flowrate set by mass flow controller and the ELPI®+ calibration flowrate of 9.930 L/min (Table 1). There were also small differences in the amount of laboratory air that entered the dilution system between puffs for the Mango- and Mint-flavored ENDS products that resulted in a difference of approximately 100-fold in the average puff dilution for the Mango- and Mint-flavored ENDS products (Table 1). ELPI®+ data temporal analysis indicated each ENDS aerosol puff was sampled for approximately 20 s, meaning that the residual ENDS aerosol in the exhaust/room air leg and connection of the two “Y”s of the dilution system at the end of the 8 s puff exhaust, was fully analyzed.

Table 1.

ELPI®+ dilution system measurements and calculations.

ProductReplicateSet by mass flow controllerMeasuredCalculations
ELPI®+ flowrate (L/min)Zero air dilution (L/min)Laboratory air entering dilution system (L/min)Puff volume (cm3)Range of percent of puff sampled by ELPI®+ aPuff dilution ratio b
Mango 5.0% (nicotine by weight)19.8599.8290.03055.45.1 – 12.8329
29.8269.7920.03455.65.8 – 13.4289
39.8179.7940.02355.63.9 – 11.6427
49.9139.8900.02355.63.9 – 11.6431
59.9159.8780.03755.66.3 – 14.0268
69.9399.9110.02855.44.8 – 12.5355
Mean ± SD9.878 ± 0.0519.840 ± 0.0650.029 ± 0.00655.5 ± 0.16.0 ± 1.0 – 12.7 ± 1.0350 ± 68
Mint 5.0% (nicotine by weight)19.9129.9770.03555.46.0 – 13.7283
29.9529.9130.03955.46.7 – 14.3255
39.9629.9280.03455.55.8 – 13.5293
49.9599.9080.05155.28.7 – 16.5195
59.9619.9140.04755.68.0 – 15.6212
610.0249.9840.04055.66.8 – 14.5251
Mean ± SD9.962 ± 0.0369.937 ± 0.0340.041 ± 0.00755.5 ± 0.27.0 ± 1.1 – 14.7 ± 1.2248 ± 39

SD = standard deviation

a

Range of percent of puff sampled = amount of laboratory air sampled over the 8 s exhaust of the pump with and without the 6 cm3 volume from the exhaust arm of dilution system divided by the total puff flow and multiplied by 1 minus aerosol deposition efficiency within the puffing machine piston. Measured aerosol deposition within the puffing machine piston was 29.045% so 1 minus deposition efficiency = 70.955 % of the puff aerosol mass was available for sampling.

b

Ratio of ELPI®+ flowrate set by mass flow controller (column 3) divided by measured amount of laboratory air entering the dilution system (column 5), which is also the amount of the puff that was sampled for the 8 s exhaust from the puffing machine.

There were small differences in the average MMAD and GSD for the Mango- and Mint-flavored ENDS products when measured by the two impactors with greater differences for the Mango-flavored ENDS product (Table 2). This was consistent with puffs from the Mango-flavored ENDS being diluted 100 times more than the puffs from the Mint-flavored ENDS. The ELPI®+ provided consistently larger GSD values, however the biggest difference between the two impactors was the calculated recovery efficiency.

Table 2.

Measurement of Particle Size Distribution and Recovery Efficiency Calculations.

ProductReplicateImpactor
Mini MOUDI™ELPI®+
MMAD (μm)GSDRecovery %MMAD (μm)GSDRecovery range (%) a
Mango 5.0% (nicotine by weight)10.831.4879.11.131.6962.9 – 157.0
20.901.5477.31.131.6959.8 – 138.6
30.781.4677.01.141.6751.3 – 152.2
4NPNPNP1.171.6753.3 – 158.1
5NPNPNP0.981.6659.4 – 132.4
6NPNPNP1.211.7254.6 – 142.6
Mean ± SD0.84 ± 0.061.49 ± 0.0477.8 ± 1.11.13 ± 0.081.68 ± 0.0256.8 ± 4.5 – 146.8 ± 10.5
Mint 5.0% (nicotine by weight)11.041.5079.81.041.6556.2 – 128.7
21.101.4079.30.961.5865.0 – 139.0
31.021.4882.91.091.6662.0 – 144.9
4NPNPNP0.811.6269.4 – 130.6
5NPNPNP0.851.6568.8 – 134.2
6NPNPNP1.071.6763.0 – 134.3
Mean ± SD1.05 ± 0.41.46 ± 0.0580.7 ± 1.90.97 ± 0.121.64 ± 0.0364.1 ± 4.9 – 135.3 ± 5.9

GSD = geometric standard deviation; MMAD = mass median aerodynamic diameter; NP = not performed; SD = standard deviation

a

Based upon percent puff sampled range from Table 1.

There were substantial differences in the number and amount of volatile chemicals identified in the air blank samples from each impactor (Table 3). This trend continued with substantially more volatile chemicals identified in the exhaust from the ELPI®+ impactor compared to the Mini MOUDI™ impactor (Tables 4 and 5).

Table 3.

Chemicals identified in impactor air blank samples.

Mini MOUDI™ELPI®+
QuantitativeSemi-quantitativeQuantitativeSemi-quantitative
ChemicalQuantity (ng)ChemicalQuantity (ng)ChemicalQuantity (μg)ChemicalQuantity (μg)
Acetone77Diethylene glycol100Acetone*5.9Propylene + propane0.3
Methylene chloride22C6-C8 Hydrocarbon230Methylene chloride *14Dichlorodifluoro methane0.1
Hexane11Acetic acid130Hexane1.4Isobutane0.2
Acetonitrile59Acetonitrile0.1Butane0.4
Ethylacetate7.8Ethylacetate0.05Acetaldehyde0.9
Carbon disulfide15Carbon Disulfide0.92-Methylbutane0.7
Ethanol0.4Trichlorofluoro methane0.4
Toluene0.1Pentane0.4
1,1,2-Trichloro-1,2,2-trifluoroethane0.09Isoprene0.4
Isopropanol1.23-Methylpentane0.3
Carbon tetrachloride0.042-Methyl-2-propenal0.2
Isopropylbenzene0.06Methylcyclopentane0.4
Cyclohexane0.02Butanal0.09
Isooctane0.3Methyl ethyl ketone0.1
Chloroform0.05Pentanal0.1
Benzene0.04Propylene glycol0.4
4-Methyl-2-pentanone0.006Hexanal0.2
m,p-Xylene0.03α-Pinene0.1
o-Xylene0.02Octanal0.1
1,2,4-Trimethylbenzene0.009Nonanal0.3
Ethylbenzene0.007Decanal0.4
Naphthalene0.006Undecanal0.1
Diethyl ether*5.9Tetradecane0.2
Dodecanal0.1
Pentadecane0.1
2-Dodecanone0.3
2-Tridecanone0.4
2,6-Ditertbutylphenol0.1
Siloxane0.1
C13-C15 Hydrocarbon0.1
C15-C17 Hydrocarbon0.1
*

Value exceeded calibration range by more than a factor of 10.

Table 4.

Quantitative and semi-quantitative analysis of vapor phase constituents from the Mini MOUDI™ Impactor. Unless otherwise noted chemicals listed were found in all three replicate samples.

ProductChemicalCAS #Amount (Mean ± SD)
Measured (μg/sorbent tube)Calculated for entire sampled volume (μg)Collected aerosol (μg/mg)
Quantitative analysis Mini MOUDI™
Mango 5.0% (nicotine by weight)Total amount aNA65.5 ± 21.11310 ± 42277.6 ± 25.0
Ethanol64-17-538.2 ± 13.3764 ± 26645.2 ± 15.8
Isopropanol67-63-00.076 ± 0.0151.52 ± 0.30.09 ± 0.02
Toluene108-88-30.007 ± 0.0010.14 ± 0.020.008 ± 0.001
Mint 5.0% (nicotine by weight)Total amount aNA34.8 ± 23.6696 ± 47236.2 ± 24.5
Ethanol b64-17-50.051.00.05
Isopropanol c67-63-00.08 ± 0.021.6 ± 0.40.08 ± 0.02
p-Isopropyltoluene c99-87-60.007 ± 0.0010.14 ± 0.020.14 ± 0.02
Benzene b71-43-20.0080.160.008
Toluene b108-88-30.0080.160.008
m,p-Xylene b108-38-3 or 106-42-30.020.40.02
Semi-quantitative analysis Mini MOUDI™
Mango 5.0% (nicotine by weight)Propylene glycol57-55-627.5 ± 11.9550 ± 23832.6 ± 14.1
Nicotine54-11-50.17 ± 0.123.4 ± 2.40.2 ± 0.1
Menthol89-78-1 or 1490-04-60.098 ± 0.0041.96 ± 0.080.116 ± 0.005
Acetaldehyde75-07-00.83 ± 0.1216.6 ± 2.40.98 ± 0.14
Ethylbutyrate105-54-40.35 ± 0.17.0 ± 2.00.41 ± 0.12
Methylcyclopentane c96-37-70.145 ± 0.0782.90 ± 1.560.17 ± 0.09
Limonene or d-limonene138-86-3 or 5989-27-50.182 ± 0.0453.64 ± 0.90.22 ± 0.05
Decanalc112-31-20.15 ± 0.063.0 ± 1.20.18 ± 0.07
3-Methylpentane b96-14-00.24.00.24
Linalool c78-70-60.09 ± 0.011.8 ± 0.20.11 ± 0.01
Nonanal b124-19-60.091.80.11
3-Hexen-1-ol544-12-70.53 ± 0.1810.6 ± 3.60.62 ± 0.21
Allyl hexanoate b123-68-20.088 ± 0.0131.76 ± 0.260.1 ± 0.02
2-Methyl-4-phenyl-1H-imidazole b13739-48-50.043 ± 0.0120.86 ± 0.240.051 ± 0.014
Mint 5.0% (nicotine by weight)*Propylene glycol57-55-620.1 ± 16.8402 ± 33620.9 ± 17.5
Nicotine c54-11-50.15 ± 0.133.0 ± 2.60.16 ± 0.14
Acetol c116-09-60.070.05 – 0.090.003 – 0.005
Acetaldehyde b75-07-00.81.60.083
3-Octanol b589-98-00.051.00.052
Decanal c112-31-20.2 ± 0.14.0 ± 2.00.21 ± 0.1
Eucalyptol c470-82-60.28 ± 0.15.6 ± 2.00.29 ± 0.1
Menthone c89-80-50.67 ± 0.3913.4 ± 7.80.7 ± 0.4
Menthol c89-78-1 or 1490-04-65.4 ± 3.5108 ± 705.6 ± 3.6
Menthyl acetate c89-48-50.11 ± 0.062.2 ± 1.20.11 ± 0.06
Benzyl alcohol c100-51-60.09 ± 0.011.8 ± 0.20.094 ± 0.01
Nonanal b124-19-60.36.00.3
1,1-Difluorethane75-37-64.00.09 – 8.00.005 – 0.42
2-Methyl-4-phenyl-1H-imidazole c13739-48-50.0450.04 – 0.050.002 – 0.003
Nerylacetone b3879-26-30.091.80.09

NA = not available; SD = standard deviation

*

C6-C8 and C9-C11 hydrocarbons found.

a

Calculated based upon internal standard, summation of individual chemical levels will not match due to semi-quantitative analysis.

b

Found in one of three samples.

c

Found in two of three samples (mean ± range).

Table 5.

Quantitative and semi-quantitative analysis of vapor phase constituents from the ELPI®+ Impactor. Unless otherwise noted chemicals listed were found in all six replicate samples.

ProductChemicalCAS #Mean ± SD (μg/collection)Collected aerosol (μg/mg)
Quantitative analysis ELPI®+
Mango 5.0% (nicotine by weight)Total amount aNA9.5 ± 2.90.07 ± 0.03
1,3-Butadiene b106-99-00.020.00014
Tetrahydrofuran109-99-90.0147 ± 0.00060.0001 ± 0.000004
2-Methyl-1-propanol c78-83-10.0009 ± 0.00010.000006 ± 0.0000007
Methyl tertiary butyl ether f1634-04-40.021 ± 0.0270.00015 ± 0.00019
trans 1,2-Dichloroethene c156-60-50.007 ± 0.0010.00005 ± 0.000007
1,2-Dichloroethane b107-06-20.0065 ± 0.00050.000046 ± 0.000035
Bromodichloromethane b75-27-40.0070.00005
Styrene b100-42-50.0050.000035
Bromoform b75-25-20.0090.00006
1,3,5-Trimethylbenzene b108-67-80.0060.000043
p-Isopropyltoluene b99-87-60.0050.000035
Ethylmethacrylate c97-63-20.006 ± 0.0020.0000436 ± 0.000014
Trichloroethene b79-01-20.0070.00005
Mint 5.0% (nicotine by weight)Total amount aNA15.5 ± 7.60.21 ± 0.1
1,3-Butadiene b106-99-00.030.0004
Tetrahydrofuran109-99-90.023 ± 0.0140.0003 ± 0.0002
trans 1,2-Dichloroethene c156-60-50.014 ± 0.0060.0002 ± 0.00008
Styrene b100-42-50.0090.00012
1,3,5-Trimethylbenzene b108-67-80.0050.000067
p-Isopropyltoluene99-87-60.027 ± 0.0050.00036 ± 0.000067
Tetrachloroethene c127-18-40.0135 ± 0.00650.00018 ± 0.000087
Diethyl ether b60-29-70.060.0008
Naphthalene c91-20-30.0075 ± 0.00050.0001 ± 0.000007
1,2-Dichloroethane b107-06-20.0080.00011
Bromodichloromethane b75-27-40.0070.00009
Semi-quantitative analysis ELPI®+
Mango 5.0% (nicotine by weight) f2-Methyl-1-propene b115-11-70.10.0007
1-Propanol b71-23-80.30.002
Menthol89-78-1 or 1490-04-61.08 ± 0.460.008 ± 0.003
Menthone f89-80-50.138 ± 0.0570.001 ± 0.0002
Ethylbutyrate105-54-42.62 ± 1.580.02 ± 0.1
Limonene or d-limonene138-86-3 or 5989-27-50.3 ± 0.170.002 ± 0.001
Decanoic acid e334-48-50.1 ± 0.00.0007 ± 0.0
Noanaoic acid e112-05-00.125 ± 0.050.0009 ± 0.0004
2-Pentanone107-87-90.182 ± 0.0780.001 ± 0.0006
Propanal e123-38-60.15 ± 0.060.001 ± 0.0004
2-Methylpentane e107-83-50.12 ± 0.050.0009 ± 0.0004
tert-Butanol b75-65-00.10.0007
Methyl vinyl ketone d78-94-40.17 ± 0.060.001 ± 0.0004
4-Methyloctane2216-34-40.1 ± 0.00.0007 ± 0.0
Heptanal111-71-70.183 ± 0.0750.0013 ± 0.0005
Ethylhexanoate e123-66-00.125 ± 0.050.0009 ± 0.0004
Acetic acid b64-19-70.10.0007
Diethylene glycol c11-46-60.15 ± 0.050.001 ± 0.0004
3-Methylbutanol b590-86-30.090.0006
Hexamethylchclotrisiloxane c541-05-90.1 ± 0.00.0007 ± 0.0
Mint 5.0% (nicotine by weight) fLimonene or d-limonene c138-86-3 or 5989-27-50.2 ± 0.00.0027 ± 0.0
Heptane b142-82-50.095 ± 0.0050.001 ± 0.0001
Propanal b123-38-60.10.001
Menthol89-78-1 or 1490-04-66.25 ± 3.710.084 ± 0.05
Menthone89-80-51.82 ± 0.780.024 ± 0.01
Eucalyptol470-82-60.72 ± 0.190.01 ± 0.003
Benzaldehyde c100-52-70.15 ± 0.050.002 ± 0.0007
4-methyloctane e2216-34-40.2 ± 0.00.003 ± 0.0
2-Methyl-1-propene d115-11-70.17 ± 0.120.0023 ± 0.0016
Trimethylsilanol d1066-40-60.1 ± 0.00.001 ± 0.0
2-Methylpentane c107-83-50.2 ± 0.140.003 ± 0.002
tert-Butanol e75-65-00.15 ± 0.060.002 ± 0.0008
Methyl vinyl ketone d78-94-40.2 ± 0.10.003 ± 0.001
2-Pentanone107-87-90.2 ± 0.10.003 ± 0.001
Heptanal e111-71-70.3 ± 0.20.004 ± 0.003
Acetic acid e64-19-70.12 ± 0.50.002 ± 0.007
Hexamethylchclotrisiloxane d541-05-90.097 ± 0.0060.013 ± 0.0001
Chlorodifluormethane b75-45-60.10.001
1,1-Difluoroethane b75-37-60.40.005
Nonanoic acid c112-05-00.45 ± 0.150.006 ± 0.002
2-Methyl-2-propenal b78-85-30.40.005
2-Hexanone b591-52-40.10.001
Pentanoic acid b109-52-40.10.001
3-Octanol b589-98-00.20.003
Heptanoic acid b111-14-80.10.001
Triacetin b102-76-10.10.001

NA = not available; SD = standard deviation

a

Calculated based upon internal standard, summation of individual chemical levels will not match due to semi-quantitative analysis.

b

Found in one of six samples.

c

Found in two of six samples (mean ± range).

d

Found in three of six samples.

e

Found in four of six samples.

f

Found in five of six samples.

g

Numerous unidentified hydrocarbons with various number of carbon atoms (C6-C8, C7-C9, C8-C10, C9-C11, C10-C12, C12-C14, C13-C15, C14-C16, and C15-C17) and siloxanes.

The mean total amount of volatile chemicals measured (background corrected) in the exhaust air of the Mini MOUDI™ impactor for JUUL Mango 5.0% ENDS was 65.5 ± 21.1 μg/thermal desorption tube and for Mint 5.0% ENDS was 34.8 ± 23.6 μg/thermal desorption tube. Since the thermal desorption tube flowrate represented only 1/20th of the total Mini MOUDI™ impactor flowrate (100 cm3/min vs. 2000 cm3/min) a mean total of 1.31 mg and 696 μg/Mini MOUDI™ impactor exhaust sample could be used to recalculate the recovery efficiency for JUUL Mango 5.0% and Mint 5.0% ENDS products, respectively (Table 6). For JUUL Mango 5.0%, the range of total device mass loss (mass loss of the pod) for the six samples using the Mini MOUDI™ impactor was from 17.4–23.3 mg. The mean of 1.31 mg of volatile compounds represents 5.6–7.5% of the device mass loss and resulted in an average recovery efficiency of 82.8% for JUUL Mango 5.0% ENDS when volatile chemicals were included. A similar situation occurred with the Mint 5.0% ENDS products where the mean amount of measured volatile compounds was 34.8 μg/thermal desorption tube or 696 μg/Mini MOUDI™ impactor exhaust sample. The range of device mass loss for the three samples of JUUL Mint 5.0% ENDS was 22.20–23.80 mg, so the volatile compounds represented 3.1–3.6% of the collected aerosol mass and resulted in an average recovery efficiency for these samples of 83.1%. The mean total amount of volatile chemicals measured (background corrected) in the ELPI®+ impactor exhaust from JUUL Mango 5.0% ENDS ranged from 40 to 360 μg/replicate and from JUUL Mint 5.0% ENDS from 60 to 430 μg/replicate considering the range in amount of puff sampled in each replicate. With the pod mass loss for the ELPI®+ impactor samples ranging from 34.2–201.2 mg, the measured volatile compounds did not substantively change the recovery efficiency for these samples (< 0.4%; Table 6).

Table 6.

Collection efficiency calculated including the measured vapor phase.

ImpactorParameterMango 5.0% (nicotine by weight)Mint 5.0% (nicotine by weight)
Replicate 1Replicate 2Replicate 3Replicate 1Replicate 2Replicate 3
Mini MOUDIDevice mass loss (mg)20.6023.3017.4023.8022.2023.40
Mass collected (mg)*16.3018.0013.4019.0017.6019.40
Calculated mass of volatiles (mg)1.080.621.240.680.080.96
Recovery efficiency (%) a84.479.984.182.779.687.0
Mango 5.0% (nicotine by weight)
Replicate 1Replicate 2Replicate 3Replicate 4Replicate 5Replicate 6
ELPI®+Device mass loss (mg)116.8122.0185.5201.265.1155.4
Mass collected (mg)*73.4 – 183.472.9 – 169.195.1 – 282.4107.2 – 318.238.1 – 86.284.8 – 221.6
Calculated mass range of volatiles (mg)0.06 – 0.160.04 – 0.10.07 – 0.210.12 – 0.360.09 – 0.190.07 – 0.19
Recovery efficiency range (%)62.9 – 157.259.8 – 138.751.3 – 152.453.3 – 158.358.7 – 132.754.6 – 142.7
Mint 5.0% (nicotine by weight)
Replicate 1Replicate 2Replicate 3Replicate 4Replicate 5Replicate 6
Device mass loss (mg)96.571.9113.734.243.489.1
Mass collected (mg)*54.3 – 124.246.8 – 99.970.5 – 164.723.7 – 45.029.9 – 58.356.1 – 119.7
Calculated mass range of volatiles (mg)0.06 – 0.130.10 – 0.210.14 – 0.330.06 – 0.110.08 – 0.160.20 – 0.43
Recovery efficiency range (%)56.3 – 128.865.1 – 139.267.2 – 145.169.5 – 131.969.1 – 134.763.2 – 134.8
*

Measured on each stage and final filter of Mini MOUDI™ and only on each stage of ELPI®+ based upon percent of puff sampled.

a

Revised mean collection efficiency including calculated volatiles for Mango 5.0% = 84.5 ± 3.2% and for Mint 5.0% = 81.6 ± 7.7%.

4.
DISCUSSION

The differences in total puffing duration, dilution air setups, and the resulting differences in aerosol residence times between the two impactors resulted in:

  • 1)

    Small to no differences in the average MMAD;

  • 2)

    Consistently larger GSD values from the ELPI®+;

  • 3)

    Substantially larger number of vapor phase constituents measured with the ELPI®+;

  • 4)

    Less total vapor phase mass measured by the ELPI®+ from the two JUUL ENDS products used in this study;

  • 5)

    Collection efficiency range consistently greater than 100% for the two JUUL ENDS products from the ELPI®+ measurements.

The difference in residence time of 0.06 s for the Mini MOUDI™ impactor vs. 8 s of exhaust from the puffing machine into the “Y” for the ELPI®+ allowed significantly more time for coagulation due to Brownian motion and probably contributed to consistently larger GSD values from the ELPI®+. The small difference in JUUL ENDS product age is unlikely to explain the difference in number or quantity of volatile substances measured by the ELPI®+ since JUUL ENDS products were tested with the ELPI®+ after the MSP-135-8 Mini Moudi™ impactor and considering the similarity in particle size distribution over 12 months for similar ENDS products (Supplemental Tables 1–6). One difference was the total amount of laboratory air sampled by each impactor. The total Mini MOUDI™ impactor sampling time of 123–125 s resulted in a total of 275 cm3 of laboratory room air (55 cm3/puff × 5 puffs) being pulled through Mini MOUDI™ compared to 3,500–4,523 cm3 of laboratory room air (Laboratory Air Entering Dilution system per second × 22 s/puff × 2 puffs/min × 49.5 min/sampling time + 3.34 cm3/s × 8 s/exhaled puff × 100 puffs) sampled by the ELPI®+ over its 49.5 min sampling time. The substantially greater sampling of laboratory room air by ELPI®+ (> 13 times) is consistent with the difference in the number of volatile compounds detected using the ELPI®+ setup, but not their smaller measured mass. If there were trace levels of volatile compounds in the ultra-zero dilution air (> 122 times), its greater use by the ELPI®+ would also explain the difference in the number of volatile compounds detected, but again not their smaller measured mass. The most likely explanation for the difference in measured volatile compound mass between the two impactors and their setups is that the vacuum oil in the sintered collection plates of the ELPI®+ and/or the pleated HEPA filter collects and retains most of the volatile compounds compared to the uncoated aluminum foils and glass fiber filter used in the Mini MOUDI™ impactor. The collection of vapor phase constituents by the vacuum oil in the sintered collection plates of the ELPI®+ would also explain the larger GSD values provided by the ELPI®+ impactor (See Supplemental Tables 9 and 10). Additionally, some of the unidentified volatile hydrocarbon compounds from sampling with the ELPI®+ (Table 5, Footnote f) may have been volatile components of the vacuum oil used on the ELPI®+ sintered collection plates. The wide recovery efficiency range of the ELPI®+ is primarily due to the variability in estimating the amount of the puff sampled, which is difficult when puffs are diluted between approximately 150–400 to 1. Contrary to what was expected with the ELPI®+ puff dilution ratio and measured MMAD, the ELPI®+ measured MMAD tended to increase with increasing puff dilution ratio suggesting that hygroscopic growth, due to humectants (glycerin and propylene glycol) may have a larger effect than evaporation for the JUUL Products used in this study (Figure 2). This would likely apply to all ENDS product formulations that use propylene glycol and glycerin. To refine the puff sampled amount and reduce the range of calculated recovery efficiency, more work with the existing ELPI®+ dilution system is required. Some of the variability in puff dilution is likely due to the two calibrated mass flow controllers (accurate to ± 1%) that were used to set the ELPI®+ and ultra-zero air dilution resulting in the range of laboratory air entering the dilution system (0.023–0.051 L/min). The wide calculated ELPI®+ collection efficiency could also partially be due to the aerosol density value (1.15 g/cm3) used in the ELPI®+ software, which was based upon the density of the two most abundant primary ingredients, glycerin and propylene glycol. Hygroscopicity of the deposited humectants within the ELPI®+ from the JUUL ENDS products does not explain the range of recovery efficiency since only approximately 60 μg of water vapor was available in 4,523 cm3 of 65% relative humidity air at 21° C (30) and there was reduced humidity in the lower stages (cutoff diameters of 0.17 μm and below) of the ELPI®+ (31). It is also clear that particulate matter in the 3,500–4,523 cm3 of laboratory air (in the μg/m3 range) sampled by the ELPI®+ does not explain the wide calculated recovery efficiency range.

Figure 2.

Correlation of ELPI®+ measured MMAD and puff dilution ratio (values from Tables 1 and 2).

Since some laboratory air was sampled either by puffing or as dilution air in the two impactor setups and these studies were performed in an active environmental chemistry laboratory (Enthalpy Analytical Inc.), it was not surprising to measure some of their commonly used chemicals (e.g., acetone, acetonitrile, carbon disulfide, methylene chloride, hexane, naphthalene, toluene, etc.), especially in the air blank samples, but also in the test samples. This highlights a limitation of this study in that both impactor setups sampled different amounts of normal laboratory room air, instead of exclusively ultra-zero air. This likely led to substantial differences in the number and amounts of measured volatile compounds between the two impactors used in this study. Significant changes to both impactor setups, that may not be feasible within a commercial laboratory situation, would be required to completely eliminate sampling of laboratory room air. Another limitation of this study was use of a standardized analytical assay designed for environmental sampling instead of analytical assays specifically designed for each ENDS product tested. Although some of the volatile chemicals identified in this study have also been identified in the e-liquid for these two JUUL Products (e.g., propylene glycol, nicotine, menthol, menthone, benzyl alcohol, menthyl acetate, eucalyptol, limonene, linalool, etc.) (32), it would have been ideal to develop specific analytical assays for the aerosol from each ENDS product tested. In a commercial laboratory situation, it is not always feasible or cost-effective to develop validated analytical assays for the aerosol from each ENDS product tested or all possible volatile chemicals and have these validated assays on the commercial laboratory's scope of accreditation (29). This is why a standardized analytical assay designed for analysis of environmental air samples was used in this study. Some of the volatile chemicals identified in this study (menthol, menthone, linalool, menthyl acetate, acetic acid, acetaldehyde, etc.) have also been reported in the aerosol from JUUL Menthol (5.0% and 3.0% nicotine) using targeted chemical analysis (33) and more extensive non-targeted analysis (GC/MS and LC-HRMS) (34).

The MMAD values obtained from both impactors used in this study are consistent with the impactor determined MMADs reported for 20 ENDS products (0.9–1.2 μm) while the GSD values from this study are slightly lower than the GSD range (1.7–2.2) reported (3). Bertholon et al. (9) used an ELPI impactor and reported a D50 of 0.65 μm for the Cigarettec® electronic cigarette. Recently a Mini MOUDI™ impactor and an ECAG+ device was used to measure the particle size distribution of JUUL Mango and Mint pods (5% nicotine by weight) and reported MMADs of 1.55 ± 0.11 and 1.41 ± 0.25 with GSDs of 2.65 and 2.64, respectively (35). The larger MMADs and GSDs are likely due to a combination of numerous factors including use of a different puffing regimen (65 cm3 puff volume over 4 s puffs with an interpuff interval of 30 s), use of different relative humidity air for puffing and bias flow, use of 10 priming puffs to preheat the coil, and differential losses within their sampling system. Our work also highlights the importance of appropriate aerosol blank samples, which because of the difference in laboratory air volume sampled between the two impactors were substantially different.

In conclusion, similar MMADs were obtained for the two different JUUL ENDs products tested using two different impactors (Mini MOUDI™ and ELPI®+) that required different setups resulting in ENDS product aerosol dilutions ranging from 1:0.82 to 400:1 dilution and a range of post aerosol generation residence times (0.06 to 20 s). Differences in GSD were likely due to the collection of volatile compounds by the vacuum oil used on the ELPI®+ sintered collection plates and its pleated HEPA filter. Additional work is required to narrow the calculated percent of puff sampled and calculated recovery efficiency range of the ELPI®+ setup that consistently exceeded the aerosol mass generated by the two ENDS products.

Language: English
Page range: 113 - 126
Submitted on: Apr 7, 2023
Accepted on: Aug 17, 2023
Published on: Sep 28, 2023
Published by: Institut für Tabakforschung GmbH
In partnership with: Paradigm Publishing Services
Publication frequency: 4 times per year

© 2023 Michael J. Oldham, Lena Jeong, Adam Ozvald, I. Gene Gillman, published by Institut für Tabakforschung GmbH
This work is licensed under the Creative Commons Attribution-NonCommercial-ShareAlike 4.0 License.