Open Access Article
This Open Access Article is licensed under a Creative Commons Attribution-Non Commercial 3.0 Unported Licence

Characteristics of nanoparticle formation and hazardous air pollutants emitted by 3D printer operations: from emission to inhalation

Jong-Sang Youna, Jeong-Won Seob, Sehyun Hana and Ki-Joon Jeon*a
aDepartment of Environmental Engineering, Inha University, Michuhol-gu, Incheon 22212, Korea. E-mail: kjjeon@inha.ac.kr
bDepartment of Ophthalmology, Hallym University, Dongtan Sacred Heart Hospital 7, Hwaseong-si, Gyeonggi-do, Republic of Korea

Received 1st May 2019 , Accepted 10th June 2019

First published on 24th June 2019


Abstract

This study examined the emissions of nanoparticles and hazardous air pollutants (HAPs) by 3D printer operations and evaluated nanoparticle deposition behavior using a prediction model. Nanoparticles and HAPs were sampled at the Inha University 3D printing center with five fused filament fabrication (FFF)-type 3D printers. The number size distribution of the nanoparticles exhibited a bimodal distribution with dominant peaks over a large size range between 70 and 100 nm and a smaller size range between 10 and 20 nm. With increasing 3D printer operation, the number concentration of 10 nm particles increased, and the final number concentration was 3.6 times higher than that of the background concentration. Nanoparticle formation and agglomeration were characterized by transmission electron microscopy (TEM) and scanning electron microscopy (SEM). Model calculations revealed that a large number of nanoparticles between 10 and 30 nm in size are deposited in the lower human respiratory tract (generation number: 16–22). A total of 14 HAPs species were detected, among which hexane, acrylonitrile, and benzene concentrations were the highest.


1. Introduction

The popularity of desktop three dimensional (3D) printers is rapidly increasing.1,2 The original purpose of 3D printers was to create prototypes before launching new products.3 However, recently the scope of 3D printers has been expanded to the medical and education fields, and industry.4–6 Various 3D printing methods have been developed, such as fused filament fabrication (FFF), selective laser sintering (SLS), and stereo lithographic apparatus (SLA).1,7 Among these methods, the FFF method is most frequently used because it is inexpensive and easy to use.8

FFF-type 3D printers use the extrusion of polymer base thermoplastic materials as a filament under high temperature conditions (∼200 °C).8 It is well-known that the thermal processing of polymer materials release toxic air pollutants, including carcinogens and chemicals that cause nervous system disorders.9–11 The composition of the toxic air pollutants produced during the thermal heating process is very complex because of the many types of thermoplastics.12 Furthermore, FFF-type 3D printers are known to emit particulate matter with sizes of <100 nm and volatile organic compounds (VOCs) because of filament extrusion.13 Moreover, this study evaluated nanoparticles and VOCs emissions from 3D printers by type of 3D printer and filament and used it to evaluate nanoparticle and total VOC emission factors. Particulate matter emission rate of 3D printers have been studied by work place conditions,14 type of 3D printers,13,15 and printing parameters.16 Kim et al. (2015) and Vance et al. (2017) measured the aerosol and VOCs emissions from acrylonitrile butadiene styrene (ABS) and polylactic acid (PLA) filaments and showed large difference by types of filaments.17,18 While significant attention has been focused on particulate matter and total VOC emissions from 3D printers in previous studies, no study has been investigated nanoparticle formation, emissions of hazardous air pollutants (HAPs), and emissions upon inhalation in the human respiratory tract.

Although most of the particles size larger than 0.1 μm effectively collect in the upper respiratory tract because of inertial collision, nanoparticles (Dp < 100 nm) penetrate deeper into the lower respiratory tract, infiltrating the extra-thoracic, conducting, and pulmonary systems upon inhalation.19–21 These nanoparticles can then pass through the bloodstream and contaminants contained in the nanoparticles can diffuse in bloodstream and can be consumed by macrophages.22 Therefore, understanding deposition behavior of nanoparticle emitted from 3D printer is crucial factor in terms of human health perspective. Moreover, HAPs are toxic air pollutants that can cause cancer and other severe reproductive and birth defects.23,24 Hence, the US environmental protection agency (EPA) manages 187 species of toxic air pollutants classified as HAPs which include some VOCs, polycyclic aromatic hydrocarbons (PAHs), pesticides, and heavy metals.25

The goal of this study was to measure the nanoparticles and HAPs produced during 3D printing to obtain the following results: (1) nanoparticle size distribution as a function of time; (2) nanoparticle morphology and elemental composition; (3) nanoparticle deposition behavior in the human airway; and (4) characteristics of HAP emissions.

2. Materials and methods

2.1 FFF-type 3D printers

Nanoparticle and HAPs emissions from 5 FFF-type 3D printers were measured. Although Yi et al. (2016) reported that the air pollutant emission characteristics vary depending on filament color, so we used common white filaments to minimize the influence of color pigments.26 Filaments were extruded and injected from the printer nozzle at 230 °C onto the sample bed which heated at 80 °C. The detailed operation conditions are summarized in Table 1. In this system, the nozzle moved on the X and Y axis, and the sample bed was downed from the top or fixed at the bottom as the Z axis. We used commercially available FFF-type 3D printers (FINTBOT Z420, TPC Mechatronics Corp., South Korea). To confirm that nanoparticle formation was enhanced by the increased number of 3D printer operating, we designated five measurement periods over a total of 3 h sampling time: (1) background (22 min), (2) a single 3D printer operation (phase 1; 23 min), (3) two 3D printer operation (phase 2; 21 min), (4) five 3D printer operation (phase 3; 40 min), and (5) background (22 min). A 10 min break was taken between second and third phases to confirm that the nanoparticle formation was inhibited.
Table 1 The 3D printer operating conditions
Printer type FFF
Filament type PLA
Hot-end temperature (°C) 230
Bed temperature (°C) 60
Layer height (mm) 0.2
Bottom/top thickness (mm) 1.0
Fill density (%) 20
Filament diameter (mm) 1.75
Nozzle size (mm) 1.0


2.2 Sampling and analytical methods

Measurements were performed in a 126 m3 size 3D printing center in Inha University. Room layout and experimental setup is provided in Fig. S1. With the exception of the office furniture, instruments, and 3D printers, the room was empty during sampling and the ventilation system was turned off to reduce the interruption of background emissions. Real-time particle number size distribution was measured using a scanning mobility particle sizer (SMPS; TSI, USA) at 80 s scan time. The sampling ports were placed 1 m from each of the 3D printers and 1 m from the office floor. To analyze nanoparticle morphology and elemental composition, the samples were characterized by transmission electron microscopy (TEM; JEM-2100F, JEOL, Japan) at 200 kV with energy dispersive spectroscopy (EDS). A TEM grid (Ted Pella, Inc., USA) was placed in the mini particle sampler (MPS; Ecomesure Inc, France) connected to a mini pump (MP-∑30, SIBATA, Japan). TEM grid sampling was conducted in phase 3. Moreover, we analyzed the 3D printer filaments to confirm the emission source of the nanoparticles. The morphology and elemental composition of the filaments were characterized using a field emission scanning electron microscope (FE-SEM; S-4300SE, Hitachi, Japan) equipped with EDS.

2.3 Aerosol deposition model

Aerosol particle deposition modeling for the human respiratory system was conducted using particle dosimetry27–29 and pulmonary morphometric models30,31 ranging in complexity. In the simplest human respiratory models, various semi-empirical models can provide accurate results using a small number of parameters. In this study, we used the multiple-path particle dosimetry (MPPD, v3.04) model,32,33 which evaluated deposition of 10 nm to 20 μm particles over the entire respiratory system. The MPPD model has been used in diverse research fields including pure inorganic aerosol, bioaerosol, indoor sources, cigarettes, and biomass burning.34–38 To calculate the number dose of particle deposition in each airway generation, the deposition rate (DR) during 3D printing operation was calculated by the eqn (1):
 
DR (# mim−1) = DF × C × TV × BF (1)
here, DF is deposition fraction, C is particle number concentration (# cm−3) of the ith SMPS channel obtained from experimental work, TV is tidal volume (cm3), and BF is breathing frequency (min−1). Then number dose (#) of each airway generation is calculated using DF value by eqn (2):
 
image file: c9ra03248g-t1.tif(2)

The deposition density (# cm−2) was calculated by applying the surface area of each generation provided by the MPPD model. Further details regarding the model input parameters are provided in Tables S1 and S2. The human lung geometry used in the MPPD model of this study is anatomically divided into airway generation numbers.39 Generation number 1–20 and 21–28 represent the tracheobronchial (TB) region and the acinar (AC) region, respectively.

2.4 HAPs sampling and analysis

Gas phase sampling was conducted using a personal mini pump connected to 280 mg of Teneax TA (40/60 mesh, Markers Inc., UK). Sampling was conducted three times for 5 min each in phase 3 and blank samples were collected in first background stage. The blank sample concentrations were subtracted. The samples were then analyzed by gas chromatography-mass spectrometry (GC-MS; HP 6890/5973, Hewlett Packard, USA; DB-1 column, 0.32 mm × 60 m × 3 μm) coupled with a thermal desorption apparatus (TDA). We analyzed total 34 gas species and confirmed that 14 of them were HAPs species.

3. Results and discussion

3.1 Nanoparticle formation

The nanoparticle number concentration and size distribution were measured using SMPS in the size range from 11 to 350 nm (Fig. 1). At the beginning of sampling, the size resolved particle number concentrations of the room was measured as a background reading for 20 min, showing a unimodal distribution with a peak at 90 nm. Subsequently, the 3D printers were turned on in order, and a bimodal distribution emerged with peaks at 15 and 90 nm. In phase 1, the formation of <20 nm nanoparticles began to be detected and 70–100 nm nanoparticles increased only slightly. During phase 2, the nanoparticle concentration with sizes ranging between 10 and 30 nm increased. The number concentration of 70–100 nm particles at the beginning of phase 2 was similar to that observed in phase 1. However, it gradually increased as a function of time because the 10–30 nm nanoparticles presumably agglomerated and 70–100 nm large single nanoparticles were also emitted. The nanoparticle characteristics at the beginning of phase 2 were similar to those observed in phase 1. During the break time in between phase 1 and 2, the 10–30 nm particle number concentration decreased but that of the 70–100 nm particles remained similar to that measured in phase 2. Formation of the 70–100 nm particles occurred via agglomeration40,41 of the smaller nanoparticles emitted during phases 1 and 2. In phase 3, five 3D printers were operated simultaneously. During this phase, it was confirmed that the 10–30 nm particle concentration dramatically increased initially with the highest concentration of all sampling periods. After phase 3, nanoparticle sampling was conducted without 3D printer operation to obtain an after-background reading. Inhibition of 10–30 nm particle emissions from 3D printer was confirmed, whereas that of the 70–100 nm particles remained at similar number concentrations during the after-background measurement due to nanoparticle agglomeration. This indicates that the emitted particles lingered in the room at size of 70–100 nm, which is most inhalable particle in the human respiratory system.42
image file: c9ra03248g-f1.tif
Fig. 1 Real-time measurement of the particle number concentrations as a function of particle diameter and time.

Fig. 2 shows the average nanoparticle size distributions with phase change. The unimodal distribution with a dominant peak at 88 nm was observed in the background reading. During phase 1, the 88 nm peak did not clearly increase, but the 10 nm peak appeared because of the nanoparticle formation originating from the 3D printers. During phases 2 and 3, the prominent peak at 10 nm noticeably increased in intensity and the secondary peak at 88 nm also increased. This indicates nanoparticle formation from the operation of the 3D printers, confirming that nanoparticle formation was accelerated by the use of an increasing number of 3D printers. With increasing duration, the concentration of the 88 nm nanoparticles increased, which indicates that the 10 nm nanoparticles became agglomerated or larger nanoparticles were directly emitted.


image file: c9ra03248g-f2.tif
Fig. 2 Average number size distribution of the nanoparticles during the background measurement and phases 1, 2, and 3.

3.2 Nanoparticle morphology and elemental composition

The morphology and elemental characteristics of the emitted nanoparticles were analyzed using TEM-EDS (Fig. 3). Most of the nanoparticles emitted from 3D printers are agglomerated near-spherical nanoparticles (Fig. 3(a)). In the 50–150 nm range, single spherical nanoparticles were analyzed (Fig. 3(b)). These particles were likely generated when the filament was melted and injected into the 3D printer sampling bed. The majority of the 10–20 nm nanoparticle emissions from the 3D printers agglomerate to form larger particles approximately 100 nm in size. Moreover, the direct emission of the approximately 100 nm sized particles contributed to the increase in the number concentration of the nanoparticle size ranging between 70 and 100 nm. This indicates the emission of gaseous species, which will be further discussed in the last section. The primary particle size is critical, since the deposition fraction and deposited areas in the human respiratory system differ depending on particle size.43 The median primary particle size was 19.85 nm, as shown in Fig. 3(c), which corresponds to the primary particle size shown in Fig. 1 and 2.
image file: c9ra03248g-f3.tif
Fig. 3 Nanoparticle morphologies (a and b), primary particle size distribution (c), and elemental composition (d).

In addition to nanoparticle morphology, the elemental composition of the nanoparticles was analyzed using TEM-EDS (Fig. 3(d)). C and Cu were detected in both the particle and background measurements because holey C coated Cu grids were used. Si, S, and K were mainly detected in the primary particles (20–30 nm). Si and K were expected and originated from the additive material used in the filament manufacturing process. Fig. 4 shows the SEM-EDS results of the filament. Most of the filament surface was smooth, except for some convex protrusions (Fig. 4(a)). The EDS analysis of the filament showed C, Si, and Mg on the surface (Fig. 4(b)). The Si in the filament agreed with the TEM-EDS results. Therefore, it was confirmed that Si and K originate from the additive materials in the polymer used for filament manufacturing. Moreover, the S in the nanoparticles was presumably originated from the heating the nozzle material at the high temperature.


image file: c9ra03248g-f4.tif
Fig. 4 PLA type filament (a) elemental composition and (b) morphology.

3.3 Nanoparticle deposition model results

To study the potential health effects of the nanoparticles emitted from 3D printers, it is important to determine the number of nanoparticles that could be deposited in human respiratory tract as a function of airway generation number. The particle dose in the human respiratory system during 3D printing is shown in Fig. 5. In Fig. 5(a), the maximum particle dose occurred at the 20th airway generation with 4 × 106 particles. Moreover, approximately 70% of the nanoparticles are deposited in the human respiratory system ranging in generation number from 16 to 28. These results indicate that the nanoparticles deposited in the pulmonary region which consists primarily of alveolar airways.44 Fig. 5(b) shows the contour plot of the particle number dose as a function of both particle size and airway generation number. The highest particle number dose (6.6 × 106) was observed at the nanoparticle size of 15.1 nm and the 19th airway generation. The 10–40 nm nanoparticles showed the highest particle number dose range, located between the airway generation number 16 and 22. This indicates that a large number of nanoparticles are deposited in the lower respiratory tract compared to that of the upper respiratory tract.
image file: c9ra03248g-f5.tif
Fig. 5 (a) Total particle number dose (#) per human air generation number and (b) particle number doses per human airway generation number and nanoparticle diameter.

Aerosol particle inhalation research typically considers deposition density and particle number dose per airway generation surface area as important parameters because each generation is characterized by different anatomical surface areas.44 Fig. 6(a) shows the nanoparticle deposition density by airway generation number, showing that maximum deposition density occurred in the third airways (3.7 × 104 cm−2). Approximately 60% of the deposition density appeared in the airway number between 1 and 5, at the starting point of the tracheobronchial region in the upper respiratory tract. This is likely because the surface areas of the airway generation number from 1 to 5 are approximately 2 to 530 times smaller than that of the surface areas of the airway generation ranging from 6 to 28 (Table S2).33 Fig. 6(b) shows the deposition densities according to both nanoparticle diameter and airway generation number. The highest deposition densities were observed for particles with diameters between 10–40 nm.


image file: c9ra03248g-f6.tif
Fig. 6 (a) Total particle deposition density (# cm−2) per human air generation number and (b) particle deposition density per human airway generation number and nanoparticle diameter.

3.4 HAPs in gaseous emissions

Previously, it was shown that VOCs were emitted along with nanoparticles from 3D printer operation, and they focused on total VOCs (TVOCs) emissions and associated emission factors.13,17 However, understanding the species present in the gaseous emissions in terms of exposure assessment for each gas species is important.45 Of the 34 gaseous species analyzed, 14 were considered classified as HAPs, which are known to be carcinogens or cause serious impact on human health and the environment.25 The analyzed HAPs are shown in Fig. 7, and the descriptive concentrations of the 14 HAPs and their associated health effects are summarized in Table S3. The International Agency for Research on Cancer (IARC) classifies chemicals in terms of groups of suspected carcinogens.46 From the HAPs results, benzene and trichlorethylene are designated as group 1 carcinogens which are obviously carcinogenic to humans. The benzene and trichlorethylene concentrations were determined to be 0.52 and 0.09 ppb, respectively. Even though these concentrations are low, it has been reported that long-term exposure to these chemicals at low concentrations can cause cancer.47 Group 2A and 2B chemicals are defined as probably and possibly carcinogenic to humans. Of the HAPs classified as group 2A and 2B, acrylonitrile, methylene chloride, chloroform, tetrachloroethylene, and styrene were detected at concentrations of 3.36, 1.64, 0.32, 0.02, and 0.33 ppb, respectively. The concentrations of hexane, toluene, xylene, and 1,2,4-trichlorobenzene were 8.46, 1.65, 1.28, and 0.26 ppb, respectively. These substances are not listed as carcinogens but are known to effects the central nervous system (CNS) after inhalation.25 Chlorobenzene (0.16 ppb), 1,4-dichlorobenzene (0.30 ppb), and hexachloro-1,3-butadiene (0.27 ppb) are not classified as carcinogen or toxic to the CNS but are hazardous pollutants that can cause narcosis, irritation, and cardiac disease.25
image file: c9ra03248g-f7.tif
Fig. 7 Concentrations of the 14 species classified as hazardous air pollutants (HAPs).

4. Conclusions

This study assessed nanoparticle formation and the emission of HAPs from the operation of 3D printers. Nanoparticle deposition behavior in the human respiratory system was also studied. The main findings are as follows, in order of the research objectives presented at the end of Introduction: (1) nanoparticle emissions exhibited a bimodal size distribution, with dominant peaks at 10 nm and 88 nm. Once the 3D printer was started, the generation of 10 nm size nanoparticles increased, and after a few minutes the concentration of 88 nm size nanoparticles increased; (2) most single nanoparticles were near spherical in shape, with a median primary particle size of 19.85 nm. The single spherical nanoparticles agglomerated and existed in the form of larger particles; (3) 10–40 nm size nanoparticles were mostly deposited in the lower respiratory system with a generation number between 16 and 22, whereas particle deposition density was higher in the upper respiratory system due to a smaller surface area; and (4) 14 HAPs species generated by 3D printer operation were analyzed, indicating that they might cause cancer or other adverse health effects once they are inhaled into the human respiratory system. As an extension of this research, future studies should focus on the evaluation of nanoparticle emission using closed camber study and quantifying the doses of each HAPs species in the human respiratory tract to inform exposure assessments.

Conflicts of interest

There are no conflicts to declare.

Acknowledgements

This work was funded by the Korea Ministry of Environment (MOE) as “the Environmental Health Action Program (2016001360005) and funded by the Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Science and ICT (NRF-2017R1C1B1008811).

Notes and references

  1. H. Lipson and K. Melba, Fabricated the new world of 3D printing, 2013, vol. 53 Search PubMed.
  2. C. Schubert, M. C. Van Langeveld and L. A. Donoso, Br. J. Ophthalmol., 2014, 98, 159–161 CrossRef PubMed.
  3. B. Berman, Bus. Horiz., 2012, 55, 155–162 CrossRef.
  4. F. Rengier, A. Mehndiratta, H. Von Tengg-Kobligk, C. M. Zechmann, R. Unterhinninghofen, H. U. Kauczor and F. L. Giesel, Int. J. Comput. Assist. Radiol. Surg., 2010, 5, 335–341 CrossRef CAS PubMed.
  5. R. Bogue, Assem. Autom., 2013, 33, 307–311 CrossRef.
  6. E. Buehler, S. K. Kane and A. Hurst, in Proceedings of the 16th International ACM SIGACCESS Conference on Computers & Accessibility, 2014, pp. 107–114 Search PubMed.
  7. J. W. Stansbury and M. J. Idacavage, Dent. Mater., 2016, 32, 54–64 CrossRef CAS PubMed.
  8. J. Skowyra, K. Pietrzak and M. A. Alhnan, Eur. J. Pharm. Sci., 2015, 68, 11–17 CrossRef CAS PubMed.
  9. D. Lithner, A. Larsson and G. Dave, Sci. Total Environ., 2011, 409, 3309–3324 CrossRef CAS PubMed.
  10. L. Mendes, A. Kangas, K. Kukko, B. Mølgaard, A. Säämänen, T. Kanerva, I. Flores Ituarte, M. Huhtiniemi, H. Stockmann-Juvala, J. Partanen, K. Hämeri, K. Eleftheriadis and A. K. Viitanen, J. Ind. Ecol., 2017, 21, S94–S106 CrossRef CAS.
  11. C. H. Shin, Y. Yoon, J. H. Park and B. C. Ma, Korean J. Chem. Eng., 2018, 35, 1225–1230 CrossRef CAS.
  12. J. Unwin, M. R. Coldwell, C. Keen and J. J. McAlinden, Ann. Occup. Hyg., 2013, 57, 399–406 CAS.
  13. P. Azimi, D. Zhao, C. Pouzet, N. E. Crain and B. Stephens, Environ. Sci. Technol., 2016, 50, 1260–1268 CrossRef CAS PubMed.
  14. B. Stephens, P. Azimi, Z. El Orch and T. Ramos, Atmos. Environ., 2013, 79, 334–339 CrossRef CAS.
  15. N. Afshar-Mohajer, C. Y. Wu, T. Ladun, D. A. Rajon and Y. Huang, Build. Environ., 2015, 93, 293–301 CrossRef.
  16. O. Kwon, C. Yoon, S. Ham, J. Park, J. Lee, D. Yoo and Y. Kim, Environ. Sci. Technol., 2017, 51, 10357–10368 CrossRef CAS PubMed.
  17. Y. Kim, C. Yoon, S. Ham, J. Park, S. Kim, O. Kwon and P. J. Tsai, Environ. Sci. Technol., 2015, 49, 12044–12053 CrossRef CAS PubMed.
  18. M. E. Vance, V. Pegues, S. Van Montfrans, W. Leng and L. C. Marr, Environ. Sci. Technol., 2017, 51, 9516–9523 CrossRef CAS PubMed.
  19. J. Heyder, J. Gebhart, G. Rudolf, C. F. Schiller and W. Stahlhofen, J. Aerosol Sci., 1986, 17, 811–825 CrossRef.
  20. R. B. Schlesinger, Concepts in Inhalation Toxicology, 1989, 163–192 Search PubMed.
  21. S. S. Park and A. S. Wexler, J. Aerosol Sci., 2008, 39, 266–276 CrossRef CAS.
  22. C. Saraiva, C. Praça, R. Ferreira, T. Santos, L. Ferreira and L. Bernardino, J. Controlled Release, 2016, 235, 34–47 CrossRef CAS PubMed.
  23. K. Sexton, J. L. Adgate, G. Ramachandran, G. C. Pratt, S. J. Mongin, T. H. Stock and M. T. Morandi, Environ. Sci. Technol., 2004, 38, 423–430 CrossRef CAS PubMed.
  24. W. W. Nazaroff and B. C. Singer, J. Exposure Anal. Environ. Epidemiol., 2004, 14, S71–S77 CrossRef CAS PubMed.
  25. US EPA, Health Effects Notebook for Hazardous Air Pollutants, 2017, https://www.epa.gov/haps/health-effects-notebook-hazardous-air-pollutants, accessed 2019 Search PubMed.
  26. J. Yi, R. F. LeBouf, M. G. Duling, T. Nurkiewicz, B. T. Chen, D. Schwegler-Berry, M. A. Virji and A. B. Stefaniak, J. Toxicol. Environ. Health, Part A, 2016, 79, 453–465 CrossRef CAS PubMed.
  27. International Agency for Research on Cancer, Trichloroethylene, Tetrachloroethylene, and Some Other Chlorinated Agents, WHO, Geneva, 2014 Search PubMed.
  28. A. Voutilainen, J. P. Kaipio, J. Pekkanen, K. L. Timonen and J. Ruuskanen, Scand. J. Work, Environ. Health, 2004, 30, 73–79 Search PubMed.
  29. M. Manigrasso, L. Stabile, P. Avino and G. Buonanno, Atmos. Environ., 2013, 67, 130–142 CrossRef CAS.
  30. H. C. Yeh, G. M. Schum and M. T. Duggan, Anat. Rec., 1979, 195, 483–492 CrossRef CAS PubMed.
  31. B. Asgharian, W. Hofmann and R. Bergmann, Aerosol Sci. Technol., 2001, 34, 332–339 CrossRef CAS.
  32. S. Anjilvel and B. Asgharian, Fundam. Appl. Toxicol., 1995, 28, 41–50 CrossRef CAS PubMed.
  33. National Institute for Public Health and the Environment (RIVM), Multiple Path Particle Dosimetry Model (MPPD v 1.0): A Model for Human and Rat Airway Particle Dosimetry, 2002 Search PubMed.
  34. J. Löndahl, A. Massling, E. Swietlicki, E. V. Bräuner, M. Ketzel, J. Pagels and S. Loft, Environ. Sci. Technol., 2009, 43, 4659–4664 CrossRef.
  35. K. Oravisjärvi, M. Pietikäinen, J. Ruuskanen, A. Rautio, A. Voutilainen and R. L. Keiski, Sci. Total Environ., 2011, 409, 4511–4518 CrossRef PubMed.
  36. Y. Zhang, W. Sumner and D. R. Chen, Nicotine Tob. Res., 2013, 15, 501–508 CrossRef CAS PubMed.
  37. S. Guha, P. Hariharan and M. R. Myers, Aerosol Sci. Technol., 2014, 48, 1226–1235 CrossRef CAS.
  38. A. E. Haddrell, J. F. Davies and J. P. Reid, Environ. Sci. Technol., 2015, 49, 14512–14521 CrossRef CAS PubMed.
  39. H. C. Yeh and G. M. Schum, Bull. Math. Biol., 1980, 42, 461–480 CrossRef CAS PubMed.
  40. S. Froeschke, S. Kohler, A. P. Weber and G. Kasper, J. Aerosol Sci., 2003, 34, 275–287 CrossRef CAS.
  41. J. S. Youn, S. J. Park, H. Cho, Y. W. Jung and K. J. Jeon, Korean J. Chem. Eng., 2018, 9, 1948–1953 CrossRef.
  42. W. C. Hinds, J. Aerosol Sci., 1999, 14, 175 Search PubMed.
  43. J. S. Youn, J. Csavina, K. P. Rine, T. Shingler, M. P. Taylor, A. E. Sáez, E. A. Betterton and A. Sorooshian, Environ. Sci. Technol., 2016, 50, 11706–11713 CrossRef CAS PubMed.
  44. H. Salem and S. A. Katz, Inhalation Toxicology, CRC Press, Boca Raton, 2014 Search PubMed.
  45. D. J. Paustenbach, J. Toxicol. Environ. Health, Part B, 2000, 3, 179–291 CAS.
  46. International Agency for Research on Cancer, List of classifications by cancer site, 2015 Search PubMed.
  47. M. C. Rousseau, K. Straif and J. Siemiatycki, Environ. Health Perspect., 2005, 113, A580–A583 CrossRef.

Footnote

Electronic supplementary information (ESI) available. See DOI: 10.1039/c9ra03248g

This journal is © The Royal Society of Chemistry 2019