Abstract
Obstructive pulmonary diseases, including asthma and chronic obstructive pulmonary disease are widespread and represent a major global health burden. Despite their impact, effective therapeutic delivery to the small airways using inhaled aerosols remains suboptimal. In this study, we present a novel in vitro airway-on-chip platform that mimics both normal and constricted small bronchial geometries to quantify the deposition charged and neutral polystyrene latex aerosol particles ranging from 0.2 to 2 µm. Analytical and numerical solutions were derived from dimensionless scaling laws to further support the experiments and predict deposition location. Our experiments showcase how electrostatic forces significantly alter deposition patterns across particle sizes in these small airways. For submicron particles, we observe the enhancement of proximal airway deposition due to the coupling of electrostatic-diffusive screening effects. For larger particles, which typically deposit only in the direction of gravity, the inclusion of electrostatic forces significantly extends their deposition footprint, enabling deposition even in orientations where gravitational sedimentation is not feasible. Constricted regions consistently exhibit lower deposition across all cases, the presence of electrostatic forces enhanced overall deposition, offering a potential strategy for targeting bronchioles. Together, these findings suggest that electrostatic attraction may be strategically leveraged to enhance aerosol targeting in the small airways, providing new opportunities for optimizing inhaled drug delivery in obstructive lung diseases.
Introduction
Obstructive pulmonary diseases, most notably asthma and chronic obstructive pulmonary disease (COPD), are highly prevalent across the globe and pose significant public health challenges. Despite the severity of the condition, effective treatment with inhaled aerosols remains limited, in part due to the poor deposition efficiencies of common inhalers, which typically deliver less than 50% of the inhaled dose (; ), with even lower efficiencies observed in pediatric populations (Oakes et al., 2023). The inefficiency in aerosol delivery is particularly pronounced in the presence of airway constriction, where the reduction in airway diameters creates higher resistance to airflow (; Oakes et al., 2023; Virchow et al., 2018). For example, in chronic bronchitis, a condition characterized by excessive mucus production and inflammation in the (small) bronchi, the accumulation of phlegm and swelling of the bronchial walls reduces the effective luminal space (Kim and Criner, 2013). Alternatively, the presence of mucus plugs occluding medium to large airways (i.e., approximately 2–10 mm lumen diameter) has been shown to be significantly associated with higher risk of all-cause mortality in COPD ().
In this landscape, delivering inhaled therapeutics to constricted airways, let alone obstructed ones, remains a vast challenge in pulmonary drug delivery. Recent inhalation therapy techniques have attempted to improve the outcomes of inhalation therapy, including modifying breathing patterns, using breath-actuated inhalers, and adjusting inspiratory flow rates (). Nevertheless, the most widely used strategy remains optimizing the aerodynamic diameter (dp) of inhaled particles, as this parameter largely governs key deposition mechanisms such as impaction, gravitational sedimentation, and Brownian diffusion (Ijsebaert et al., 2001). Aerosol design criteria advocate that aerosols in the size range between approximately 2 and 6 μm hold the best potential to deposit in the central and small airways () as larger particles (>6 μm) tend to mainly deposit in the upper airways due to impaction whereas aerosols <2 μm deposit mainly in the deep lungs (i.e., alveolar regions). Despite such general guidelines, one significant particle property that is often overlooked in addressing aerosol deposition () is the electric charge (q) acquired by aerosols during inhaler-generated formation (Kwok and Chan, 2009).
Briefly, when an inhaled charged particle approaches the (neutral) lung tissue, it induces a localized electric field triggering a dielectric effect (Wilson, 1947; ; ). This effect causes surrounding charges or dipole molecules within the tissue to reorient in response to the particle’s projected field (E). The interaction between a particle’s inherent charge q and the induced dipoles or charges within the tissue results in an electrostatic attraction force Fe; a phenomenon shown to enhance particle deposition in proportion to , where B is the particles mechanical mobility (Melandri et al., 1983; ). While factors such as the carrier fluid (i.e., air), the gravitational field, particle density (which typically does not deviate significantly from that of water), and particle shape (often approximated as spherical, most notably for nebulized suspensions), are often challenging to modify significantly, the charge acquired by an aerosol can be engineered and customized (Kwok and Chan, 2009; Wilson, 1947). Inhalers can generate charged particles from near-neutral values following a Boltzmann distribution to the Rayleigh limit () (i.e., scenarios in which electrostatic repulsion overcomes surface tension of the droplet). For instance, a dp= 1 μm droplet can carry a charge ranging from q ∼ 1 e to approximately 44,000 e. In fact, the ability to manipulate electrostatic charge is well established across the broader aerosol industry (; Osman et al., 2015) and thus warrants further investigation in the context of pulmonary drug delivery (; ).
To date, numerical simulations (Koullapis et al., 2016), in vitro studies () and non-dimensional analysis (; ) have underlined the importance of electrostatic charge towards deposition outcomes. Notably, these studies account for proper scaling based on realistic charge levels and particle sizes generated by commercial inhalers (Kwok and Chan, 2008). Electrostatic forces may play a dominant role in the deeper lung regions where characteristic length scales are small (; ). There, electrostatic attraction, which is inversely proportional to the square of the distance between an aerosol and the luminal wall, can potentially overshadow conventional deposition mechanisms when particle charge is sufficiently elevated relative to size (; ; ). Yet, studies addressing the induced electrostatic forces have largely focused on deposition in large airways (; Xi et al., 2014). This observation results amongst others from historical and technical challenges of exploring the phenomenon at true scale in small bronchioles. Concurrently, while electrostatic forces are expected to play a significant role in the smaller lung regions, electrostatic attraction obeys the principle of superposition in addition to other deposition mechanisms and thus also contributes to deposition in the extra-thoracic and upper airways. Despite some recent interest in the field (), there remains a dearth of available data exploring electrostatic-driven pulmonary deposition at small scale (Sznitman, 2022).
Motivated by the ongoing shortcomings on aerosol electrostatics in the lungs, we investigate in vitro the role of electrostatic forces on aerosol deposition within small airway models. Specifically, we focus on how electrostatic forces influence deposition patterns in normal and constricted bronchioles. To this end, we present a novel in vitro airway-on-chip lined with electrically conductive material to mimic the conductive properties of the luminal airway tissue. We capture the size and branching structure of small to terminal bronchioles (i.e., corresponding to generations 12–15 of the seminal Weibel A model (Weibel, 1963)); the primary sites of airflow obstruction in COPD that contributes to the disease’s characteristic symptoms (Xu et al., 2022). We quantify aerosol deposition patterns spanning aerosol sizes of 0.2 µm–2 μm, representing particles typically influenced by Brownian diffusion and gravitational sedimentation, respectively (). Our experiments examine two specific charge distributions: (i) intrinsically charged atomized particles (>100 e) and (ii) neutralized particles following a Boltzmann distribution. By directly comparing these distributions, in conjunction with particle size, we attempt to shed new quantitative light on how electrostatic charge can overcome traditional deposition mechanisms and ultimately alter local deposition outcomes.
Methods
Aerosol exposure experiment
Aerosol exposure assay
The aerosols are generated using a collision-type atomizer (Model 3076, TSI) with filtered air as the gas source, maintained at bars. To minimize water vapor content while preventing particle loss, the aerosol stream is directed through two consecutive diffusion dryers (Model 3062, TSI). The resulting airflow rate is 1.2 L/min (Figure 1a), measured using a flow meter (Model 4100, TSI). Such airflow rate is about three orders of magnitude higher than the desired bronchiolar flow rate (Q0) under quiet breathing conditions (Pedley, 1977) to feed the first generation of the four microchannel models positioned on a single microscope slide (Figures 1a,b). To achieve the required flow conditions, a stepwise air reduction and splitting mechanism is implemented, as detailed in the table of Figure 1c. The first stage of flow reduction is achieved using an assembled array of tubes (all apparatus dimensions are provided in the Supplementary Material (SM) Supplementary Figure S1a and S1b), which gradually reduces the airflow in a manner analogous to airway branching in the lungs, from the bronchi down to the bronchioles. This system, referred to as the “air reducer” (Figure 1a), expels approximately ∼80% of the excess airflow (Qbleed) through a 90° outlet into a chemical hood. The remaining airflow is symmetrically split to feed two parallel lung-on-chip models (Qchip1 and Qchip2). Each of these flows supplies four microchannels (Qmicro) through their respective inlet ports (Figure 1b).
FIGURE 1
Before entering the first microchannel generation (i.e., Gen. 0, see Figure 1d and top view in Figure 2c) the bulk airflow passes through a small vertical pipe for additional reduction Qpipe. The desired flow rate into Gen 0, .designated as Q0, is then drawn under a constant steady flow into the microchannels via a syringe pump (PHD Ultra, Harvard Apparatus) connected to the suction outlet of each chip (Figure 1b; Supplementary Figure S1a). The constant flow approximation in the model is sufficient to accurately represent physiological breathing, as further justified in the device design section. The final reduction step ensures that Q0 is only 6% of Qmicro and is directed to the model inlet to reach the desired quite breathing rate, while the remaining 94% is expelled. The syringe pump simultaneously withdraws a total airflow of QS = 4Q0. Within each microchannel, the airflow is further subdivided in a manner mimicking bronchiolar bifurcation. The airway-on-chip design features individual airway trees with three bifurcating generations, includes both normal and constricted segments.
FIGURE 2
Fluorescent particle aerosolization
Fluorescent polystyrene latex (PSL) particles are used to investigate size-dependent aerosol deposition (ρp = 1,050 kg/m3, Fluoromax red and green, fluorescent microspheres, 1% solid, Thermo Scientific). The use of PSL particles is well established in aerosol exposure studies involving airway-on-chip models (; ; ) due to their well-characterized physiochemical properties and charge acquisition during aerosolization. To examine the effects of particle size and electrostatic charge on deposition, distinct monodisperse PSL sizes are used: = 0.2, 0.5 and 2 μm. While the smallest particles (i.e. 0.2 μm) are known to be significantly affected by Brownian diffusion, the largest ones (i.e. 2 μm) are subject to considerable sedimentation in small airways. In contrast, 0.5 μm particles sit in a size range where airborne transport is primarily driven by convection such that aerosols are exhaled under normal breathing and deposit in minimal amounts due to weak diffusive or gravitational forces (; Sznitman, 2013). Furthermore, and most relevant to the present study, PSL particles are selected due to their natural tendency to acquire electrostatic charge during the aerosolization process. The “Charged Group” consists of particles exiting the aerosol generator after diffusion drying, with an estimated mean charge of >100 e based on previous experimental data (; Whitby and Liu, 1968). In contrast, the “Neutralized Control Group” is obtained by employing an electrical ionizer (Model 1090 MSP, TSI) and ensuring charge neutralization following a Boltzmann charge equilibrium. Under these conditions, the mean charge per particle is expected to be q ∼ 1 e for 0.2 μm and q ∼ 3 e for 2 μm particles (). Note that these levels are insufficient to induce significant electrostatic deposition, allowing for direct comparison between charged and neutralized aerosols (Osman et al., 2015).
Prior to exposure (see previous section), polystyrene latex (PSL) particles are suspended in deionized (DI) water at a high concentration (∼107–108 particles/mL) to ensure sufficient deposition data for ensemble statistics within reasonable exposure times (∼30–90 min, depending on particle size and charge conditions) (). The particle suspension is prepared at a weight concentration of 5%, carefully selected to maintain an aerosol with <1% aggregate formation. This assumption follows from the expected evaporation of water droplets exiting the atomizer, characterized by a count mean diameter of 0.35 μm and a geometric standard deviation of 2 μm, as specified by the manufacturer. The size distribution of the generated aerosol is assumed to follow a log-normal distribution with monodisperse spherical particles, consistent with previous studies (; Raabe, 1968). To minimize aerosol aggregation due to the high particle concentration, each suspension undergoes 30 min of sonication in an ultrasonic water bath (Elma Elmasonic S10) immediately before experimentation. To prevent particle sedimentation during experiments, a magnetic stirrer continuously operates inside the suspension container.
Deposition quantification
Following each aerosol exposure experiment, the bronchiole microchannels were examined using an inverted fluorescent microscope (Nikon Eclipse Ti) at varying magnifications (×2, ×10 and ×20) depending on particle size (i.e., 0.2 μm, 0.5 μm and 2 μm respectively), see Figures 1e,f. Using ANDOR Zyla sCMOS camera, we captured the entire model by stitching them together with zero overlap, using large free-shape images function. To reduce the image memory space, we used 4 × 4 binning. The quantification of deposition involved identifying the precise planar locations (x, z) of individual deposited fluorescent particles through digital processing of local intensity maxima using ImageJ software. Subsequently, the particles were counted with correspondence to their two-dimensional (2D) location and the sum of deposited particles were obtained for each airway generation. The deposition fraction was then calculated as the percentage of deposited particles in a specific airway generation out of the total number of deposited particles across the entire model. This approach enables us to determine the spatial distribution of particle deposition within the in vitro model.
Airway-on-chip
Device design
The in vitro airway-on-chip platform consists of four planar, symmetric airway path trees spanning four generations 0–3, as schematically shown in Figures 1, 2c. The tree design broadly mimics the deep bronchial airway branches (i.e., bronchioles) in the distal region of the conducting zone, with hydraulic diameters of <1 mm (Figures 1a,b). Each airway segment consists of square cross-sections, with dimensions summarized in Figure 2b. The chosen generation dimensions were based on morphometric measurements representative of typical length and diameter of an adult human lung, corresponding to generations 12–15, following the seminal works of Weibel (1963), Horsfield et al. (1971). These airway dimensions mimic the deep bronchial regions, referring to airways with a diameter smaller than 2 mm, which are associated with constricted diseases (McNulty et al., 2014). While the Weibel A model assumes idealized cylindrical airways, actual distal bronchiolar cross-sections in the human lung are often irregular and deviate from perfect circles (; Wang et al., 2020). In our airway-on-chip design, we employed square cross-sections (Figures 2c,d) as a deliberate engineering choice to enable precise microfabrication (Figure 1b), optical access for high-resolution imaging (Figures 1e,f), and consistent channel alignment across multiple generations. Importantly, flow behavior and particle transport in rectangular microchannels remain comparable to those in circular channels when the hydraulic diameter (defined as Dℎ = 4A/P, where A is the cross-sectional area and P is the wetted perimeter) and flow conditions (Re) are matched. To ensure physiologically relevant flow, we designed each generation as a square channel to have a hydraulic diameter Dℎ equal to that of the corresponding cylindrical segment in the Weibel model (see Figure 2b). While minor deviations in flow may occur near corners due to reduced velocities, these are localized and do not significantly affect the overall transport and deposition behavior of aerosols. Therefore, the underlying physics remain valid despite the geometric simplification. The first airway (Gen. 0, see Figure 1d and see top view Figure 2c) serves as the model’s inlet, ensuring fully-developed flow conditions over a distance exceeding the anticipated entrance length Le by more than two orders of magnitude at such low Reynolds numbers (). These conditions result in fully-developed laminar Poiseuille flow profiles within the model’s bronchioles (see Figures 2c,d).
Approximating a constant inhalation condition may be assumed for modeling airflow in the deep lung regions during quiet breathing, as the Womersley number (Wo = fL2ρf/μf), which is a dimensionless expression of the oscillatory airflow frequency (f ∼ 0.25 Hz for quite breathing) in relation to viscous effects, remains well below unity in our bronchial region of interest (Wo ≪ 1), indicating quasi-steady flow conditions that can be approximated as steady flow conditions (Ménache et al., 2008). However, it is important to note that while the airflow can be considered steady, the absence of exhalation in the model may lead to overestimating particle retention, as particles that would otherwise exit the lungs during exhalation are not accounted for. The Gen. 0 segment terminates in a bifurcation with a 25° split angle. Beyond this point, the airway tree bifurcates symmetrically across all generations, maintaining a consistent 37° branching angle along the symmetry line (Ménache et al., 2008). The model branches into two distinct pathways: one representing normal anatomical bronchioles and the other featuring a 30% constriction, designed to mimic smaller airways or regions affected by obstructive pulmonary diseases (Figure 2). Airflow regulation is achieved using a syringe pump connected to the model outlet port (Figure 1b). The present design enables the investigation of electrostatic deposition by incorporating a conductive coating, while electrical grounding prevents charge buildup ().
Device fabrication
Models were fabricated using 3D-printed molds printed using a PRUSA SL1 3D printer and 3DM-ABS Orange Tough Resin), see Supplementary MaterialSupplementary Figures S1c and S1d. The models were filled with polydimethylsiloxane (PDMS), inspired by recent microfluidic lung airway models (; ). PDMS, a well accepted material for microfluidic devices, was selected for its optical transparency and high flexibility, allowing for easy molding into complex shapes and structures. The printing settings included a layer height of 25 μm and an exposure time of 10 s. Post printing, the molds were cleaned with isopropyl alcohol, and to solidify the resin the molds were cured at 60 °C for 1 h under ultraviolet light (Formlabs Form Cure FH-CU-01 curer). After pouring the PDMS, the PDMS-filled molds were placed in desiccators for 1 hour to release trapped air bubbles. The PDMS models were peeled after 24 h of resting within the mold under ambient conditions. The outlet port for the syringe insertion (Figure 1b) was created using a 2 mm biopsy punch (Miltex, 3331).
Conductive Layering
The electrical conductivity of lung parenchymal tissue closely approximates that of saline water S/m (Osman et al., 2015). This high conductivity allows electrical charges within lung tissue to redistribute rapidly, facilitating the manifestation of image charge phenomena (Melandri et al., 1983). In contrast, glass has a significantly lower conductivity ( 10−11 S/m), which necessitates an alternative approach to achieve a physiologically relevant conductive environment.
To address this limitation, we applied a conductive indium tin oxide (ITO) coating to our microscope slide, following previous methodologies (). The ITO pattern was selectively deposited only in regions corresponding to the airway microchannels (see Supplementary Figure S2) to avoid interfering with the adhesion of the PDMS to the glass. The bonding process was performed via a 1 min plasma treatment (ETP, INC. Model BD-20). The ITO coating, provided by Huizhou Konshen Glass Co., Ltd., had a sheet resistivity of approximately , significantly enhancing the surface conductivity to ∼102 S/m (see Supplementary MaterialSupplementary Table S1 for details on sheet resistivity). This elevated conductivity enables the manifestation of image charge effects, more faithfully replicating the electrostatic interactions expected in lung tissue. Notably, the electrostatic relaxation time scale ( pico s) in real lung tissue is orders of magnitude shorter than the characteristic timescale of a particle ( μs). Thus, achieving a sufficiently conductive environment in our model mimics more faithfully electrostatic behavior. Additional details on the relevant timescales are provided in the Supplementary Material (see Supplementary Material Conductive Layering of the Microchannel). To maintain a consistent electrical setup across experiments and prevent excess charge buildup during exposure assays, the coated conductive glass was electrically grounded using a custom external connector wired to the inner surface of the microfluidic channel (See Supplementary Figure S1a).
Computational fluid dynamics (CFD)
Numerical simulations (see Figure 2; Supplementary Figure S3b) were performed using the commercial solver Fluent 24.2 (ANSYS, Inc.), employing a steady-state laminar flow model with air as the working fluid ( = 1.22 kg/m3, = 1.78 × 10−5 Pa·s). The model’s boundary conditions included a zero-pressure inlet, an outlet flow rate of 2.5 mL/min, and a no-slip condition on the walls. The momentum equations were discretized using a second-order upwind scheme for velocity and a second-order scheme for pressure. Velocity-pressure coupling was handled via the SIMPLE algorithm, with a least-squares-based scheme for gradient calculations. The computational domain was meshed using Ansys Meshing, employing polyhedral cells ranging from 400K to 3M as part of a mesh convergence study.
Analytical and numerical modeling
To gain a deeper understanding of particle deposition within the model, we explore the dynamics of aerosol transport and apply Newton’s second law to charged particles and derive their expected trajectory within the mid-channel xy plane (see Figure 1d; Supplementary Figure S3b). The governing forces considered in this analysis include (i) electrostatic attraction due to induced charge (−keq2/4y2), where ke is the Coulomb constant, (ii) gravitational sedimentation (+mg), (iii) viscous Stokes drag (−vrel/B), where vrel is the relative velocity of the particle with respect to the surrounding streamwise flow velocity uf (whose profile varies along the y-direction) and B is the aerosol mechanical mobility, and (iv) Brownian motion (FB). Assuming no hygroscopic growth (i.e., constant particle mass m), the particle dynamics are described by the coupled equations of motion 1 and 2, presented below:
These equations capture the interplay between electrostatic, drag, gravitational, and diffusive transport mechanisms, providing a framework for predicting particle trajectories within the airway-on-chip model. Note that due to the electrostatic term (∝ y−2), the governing equation in the y-axis (Equation 1) is a nonlinear, second-order, non-homogeneous differential equation, such that a direct analytical solution for the particle height as a function of time is not feasible. Consequently, we extract particle trajectories over time and space using two approaches: (i) an Explicit Euler method, and (ii) non-dimensional analysis.
Explicit Euler method
We propose a canonical model featuring rectangular channels connected in series according to Weibel dimensions mimicking our in vitro model (Supplementary Figure S5). In each channel, the streamwise flow (uf) is approximated using the laminar Hagen–Poiseuille velocity profile (see Supplementary Figures S6 and S7, and derivations for flow in a square cross-section chanel) at low Reynolds number (Re ∼ 3), corresponding approximately to quiet breathing of a healthy person at these deep lung generations. The numerical simulations are initialized with an aerosol starting at a vertical distance of height , chosen to align with the image charge method (. Additionally, this distance from the tissue corresponds to approximately 10% of the annular volume in an ideal cylindrical airway (see Supplementary Figure S8 and derivation in the Supplementary Material), making it a widely accepted threshold for assessing whether electrostatic deposition is significant, as supported by previous studies (). This value will later be chosen as the characteristic length (Lc= y0 = 0.025D). The initial velocity is set based on the corresponding position in the Hagen–Poiseuille flow profile. Following each iteration, based on the time interval Δt, a 1D random Brownian step ( was added to each calculated location, for each axis separately. Here, Ddiff represents the Stokes–Einstein diffusion coefficient. To ensure convergence, the time step was set to 0.1 ms (Δt = 0.1 ms). To validate this choice, we performed simulations with progressively smaller time steps, ranging from 10 ms down to 1 μs, confirming that the solution converges (see Supplementary Figure S8a). Additionally, to account for the inherent error in the Explicit Euler method, where terms of O (Δt2) and higher are neglected, we tracked the local truncation error (LTE) at each iteration. The resulting global truncation errors (GTE) remained around single microns and are summarized in the Supplementary Material (Supplementary Figure S8b; Supplementary Table S3). Additionally, to validate our numerical approach, we compared the results with an analytical solution derived from scaling laws. As the time step (Δt) decreases, both solutions converge well (see Supplementary MaterialSupplementary Figure S8 for comparison). Since our time intervals are much larger compared to the examined particle relaxation time (τp ≪Δt), the particle’s velocity instantaneously follows the surrounding flow velocity ().
Non-dimensional analysis
Non-dimensional analysis is helpful for understanding aerosol transport and its leading physical mechanisms, in particular when electrostatic forces are involved. This approach expresses the dominant mechanisms in terms of characteristic length (Lc∼ D), velocity (Uc∼Q/A) and time (tc∼ Lc/Uc) scales relevant to the bronchioles. We revisit the governing mechanisms of aerosol deposition in our model using dimensionless parameters (Sznitman, 2022; Tsuda et al., 2013). The deposition of small (approximately <0.5 µm) particles in the bronchioles occurs mainly via Brownian motion, characterized by the inverse Péclet number: Pe−1= Ddiff/UcLc. For larger particles (>1 µm), bronchiolar deposition is primarily governed by gravitational sedimentation, quantified by the gravity number H = vter/Uc, where vter is the settling velocity given by vter = τpg, and τp represents the particle relaxation time. Concurrently, particle inertia is captured by the particle Stokes number: Stk = τpUc/Lc.
We then extend the discussion to when electrostatic charge is present. When Coulomb forces arising from particle–wall electrostatic interactions introduce an additional dimensionless group relevant to represent the Induced Charge versus the airflow convection (): Inc = (Bke/Uc) (q/Lc)2. By applying dimensional analysis to Newton’s second law, while omitting the noise introduced by Brownian motion (which can be superimposed subsequently), we obtain the non-dimensional governing equations for particle dynamics, Equations 3, 4, presented below:where is the particle’s normalized horizontal velocity relative to the flow. Non-dimensional analysis in the deep bronchioles reveals that impaction is negligible (Stk ≪ Inc, H), and primarily occurs in the upper lung regions (Stahlhofen et al., 1989). As a result, the governing Equations 3 and 4, can be simplified by neglecting the inertial term, reducing the overall complexity of the system to Equations 5 and 6.
Whereas the horizontal particle trajectory follows a straightforward path, showcasing how the particle velocity aligns with that of the surrounding flow (Equation 6), the non-linear equation in the vertical direction (Equation 5) remains more complex and requires further analysis.
We define the characteristic length scale as the “Equilibrium Distance” from the wall, where the electrostatic force balances gravity, i.e., (see Figure 3A). It is important to note that this length scale is specific to each particle’s size and charge levels and is given by . To better characterize electrostatic deposition, we define the normalized distance from the tissue as: where the initial condition is given by Hence, electrostatic deposition against gravity is expected only within the range . A careful simplification of Equation 6 (detailed in the Supplementary Material, in section “Analytical and Mathematical Derivation”) leads to the following form (Equation 7):
FIGURE 3
This non-dimensional equation is no longer a differential equation, making it significantly easier to solve and hence analyze particle behavior over time. We can write it dimensionally as:
This Equation 8 describes displacement through two distinct terms. The first term on the left represents the classic free-fall sedimentation over time, while the nonlinear logarithmic term captures the opposing electrostatic attraction counteracting gravitational pull. Despite this clear separation, the Equation 8 is implicit in meaning an explicit expression for the particle’s position as a function of time is not straightforward. Yet, for very small particles, highly charged ones, or in a near proxinity to the tissue, gravitational settling is negligible (H≪Inc), the equilibrium distance extends toward infinity in comperison to . Under these conditions , we can expand the logarithmic term in a Taylor series around zero up to the fourth order (), leading to a simplified explicit results presented in Equation 9:
This result offers a simplified yet direct analytical expression for the vertical displacement of the particle trajectory, providing a fundamental framework to understand how electrostatic charge influences deposition dynamics under Stokes drag.
An alternative approach to Equation (8) can be formulated using a non-dimensional framework. By considering the deposition condition the equation reduces to involve only two non-dimensional parameters: , representing the free-fall distance normalized by the equilibrium distance, and (analogous to ). This reformulation leads to Equation 10
Equation 10 provides a generalized expression for the final deposition time , accommodating any combination of charge, particle size and shape across all normalized initial conditions (Figure 3B). Upon reorganizing, the final deposition time is explicitly given in Equation 11:
Statistical analysis
We conducted independent experiments for the combination of each particle size and group (i.e., charged vs neutralized). A three-way analysis of variance (ANOVA) was used to examine the effect of the experiment settings. The independent variables were the aerosol size (i.e., dp ∈ {0.2, 0.5, 1.1} µm), the electrostatic charge group (i.e., q ∈ {Neutralized, Charged}), generation (i.e., Gen ∈ {G1,G2,G3}) and their interactions and state (constricted or Normal). The dependent variable was the Deposition Fraction (DF). In this work, we report both the p-values and the estimated effect sizes (β), where β represents the magnitude of change in DF and is expressed as a percentage to align with the units of DF itself. Additional details on the statistical analysis outcomes, beyond those presented in the Results and Discussion section, can be found in the Supplementary Material.
Results and discussion
In what follows, we present deposition data obtained from quantitative microscopy and numerical simulations, focusing on particle deposition in the microchannels. We extract from microscopy imaging the deposition fractions (DF) for each airway after completion of each exposure assay. A summary of the results is shown in Figure 4, where we present depositions according to particle sizes (Figure 4, rows) and regions: i.e., Constricted and Normal (Figure 4, columns). The histograms are normalized according to the total number of deposited particles. We discriminate between experiments where aerosols were neutralized (i.e., light gray histograms) and when we allow them bypass the neutralizer and maintain highly charged conditions (i.e., dark gray histograms). We note that the effect of lung generation (Gen. 1, 2 and 3) on deposition is not significant (p = 0.095), suggesting that it is not enough to conclude airway branching contributes to deposition in this particular model. This result is reasonable, as the airway generations examined in this study share similar characteristic length scales and flow velocities. However, in a larger model spanning additional airway generations, we would likely expect lung generation to exert a stronger influence on deposition patterns.
FIGURE 4
To ensure that the model reflects physiological flow, we conducted CFD simulations where we observe Hagen–Poiseuille-like flow across all the generations. In Figure 2a, we notice a minor overestimation of flow reduction due to constriction in the theoretical analysis (circles) compared to the CFD results (squares). This discrepancy likely arises from the minor pressure drop added by changes in airway diameter, as well as the bifurcation angles, that are omitted in the theoretical model. Moreover, in our analysis, each airway is modeled as an ideal rectangular resistor (see Supplementary MaterialSupplementary Figure S3a), based on the dimensions shown in Figure 2b. However, further down the airway, both the CFD results and the theoretical analysis converge as anticipated.
We recall that in our experiments, we examine a mild constriction (30%) with a moderate reduction in flow, which approximates the airflow reduction in the constricted regions about two-thirds of that in an open healthy airway and hence less deposition expected there. The influence of lung state (Normal vs Constricted) on deposition is found to be highly significant (p < 0.0001, β = 8%), underscoring that airway narrowing indeed significantly alters deposition patterns. This finding aligns with clinical studies showing that constricted airways are associated with impaired particle transport and reduced deposition efficiency (; Xu et al., 2022). It is important to highlight how the characteristic length scale differs in the constricted region versus normal in the model (Figure 3d); a point we will revisit due to the greater sensitivity of electrostatic forces to length scale compared to other mechanisms. Particle size has a significant impact on deposition in the model (p < 0.0001, β = 14% for 0.2 μm–2 μm), confirming that different particle sizes lead to distinct DF. This finding is consistent with existing literature across the aerosol drug delivery community (Ijsebaert et al., 2001; ; ICRP, 1994).
Role of electrostatics versus diffusion
Upon examining the smallest particle size (0.2 μm), deposition primarily occurs due to Brownian motion (Figures 4a,b) since Pe−1 is eight times larger than H. When neutralized (light gray), deposition in the normal region remains constant, plateauing at ∼20% across the three generations. As expected, in the constricted region, there is less flow with available particles to deposit, leading to a lower deposition rate of around ∼12% per generation. However, when the particles are charged, a different deposition pattern emerges. Charge has a strong and highly significant effect on deposition fraction DF (p < 0.0001, β = 15%), indicating that electrostatic forces play a major role in particle behavior in our model. We observe a screening effect, where deposition is dominated in the first generation, and this effect extends to the rest of the model. Screening refers to particles that are unlikely to reach more distal regions because they first collide with nearby surfaces, effectively “screening” the regions farther downstream. This results in the streamwise depletion of particles. The phenomenon of “diffusional screening” is well known in the respiratory system (; ). Here, we clearly observe that the superposition of Brownian motion and electrostatic forces leads to screening in proximal airways. The interaction between charge and generation is highly significant (β ∼ −18% for Gen 2 and β ∼ −22% for Gen 3 with respect to Gen 1,p < 0.0001 for both), demonstrating how charge-driven screening influences DF. This effect has also been observed in a recent in vitro electrostatic study in the acinar region (). When comparing the charged cases between normal and constricted regions, the DF ratios remain similar, reflecting the mild degree of constriction applied in this particulate study. As a result, no major change in deposition patterns was observed. However, with more severe constriction, we would indeed expect enhanced electrostatic effects, as the reduced characteristic length (Lc) favor the electrostatic attraction (Inc⋅Pe ; ).
To gain further insight, we analyzed particle trajectories using numerical simulations. Figure 5 presents five representative cases that are anticipated to occur in our model, considering the given particle sizes and expected charges (; Whitby and Liu, 1968). Each particle in our numerical simulations experiences the same four main forces (Figure 2a; Supplementary Figure S5b inset): an electrostatic force directed upward toward the top of the model (), gravitational force acting downward (, Brownian random diffusion, and viscous drag, which mitigate it vertical motion and primarily pulls the aerosol downstream in the horizontal direction in a convective manner. It is important to recognize that the characteristic velocity in the bronchioles Uc, which represents the convective nature of airflow, appears in the denominator of each dimensionless term (Inc, H, Pe−1∝ Uc−1). Since Uc is significantly larger than the remaining parameters, these terms become much smaller in magnitude (Inc, H, Pe−1≪1), supporting the phenomenon of poor aerosol deposition in the bronchioles due to the dominance of convection over other mechanisms (ICRP, 1994). Hence, to better understand the interplay between electrostatic forces (Inc) and other deposition mechanisms (H, Pe−1), we compare the relative contributions of electrostatic effects to the dominant transport mechanisms. For the largest particles (2 μm), where Pe−1≪H, the ratio Inc/H compares electrostatic attraction versus gravitational settling (Figures 5b,d). Conversely, for the smallest particles (0.2 μm), where H ≪ Pe−1, the term Inc⋅Pe contrasts electrostatic with random diffusive motion (Figures 5a,c). These ratios effectively eliminate the dependence on Uc, allowing us to isolate mechanism interplay from overwhelming convection effect. This approach provides a clearer framework for analyzing the physical mechanisms governing aerosol deposition and enhances our discussion of electrostatic effects in the deep lung.
FIGURE 5
Figure 5a highlights the dominance of electrostatic forces over diffusion (Inc⋅Pe ∼ 10 and inset Inc⋅Pe ∼ 100), whereas Figure 5c provides a comparison by illustrating deposition mechanisms when electrostatic effects are negligible in comparison to diffusion (Inc⋅Pe = 0.1). Note that all particles were released from a height y0, corresponding to 5% of the radius of generation 0, coinciding with the characteristic length scale (Lc ∼ y0). This choice satisfies the image charge assumption (y0 << D0) and follows the conventional approach for examining electrostatic effects through dimensionless analysis (). This distance y0 presents ∼10% of the total airway volume (or cross-section) populated by aerosols prone to deposition. We immediately observe a clear difference between the two cases for 0.2 µm particles. When electrostatic forces are superimposed with Brownian motion (Figure 5a), the particle trajectory is not entirely random; rather, it oscillates slightly around a dominant, sharply directed path toward the tissue, dictated by electrostatic attraction (dashed green line). This sharp path can be predicted using our analytic solution (Equation 9) which closely matches numerical simulations, with even greater agreement observed in the inset (Inc⋅Pe ∼ 100). This equation closely approximated the trajectories of small-mass particles (e.g., 0.2 ) where gravity is almost neglected (H < Pe−1) or for highly charged particles in close proximity to the tissue, where other forces become insignificant.
The charged particles exhibit very short airborne times, with deposition occurring within tenths of a second (e.g., 0.134 s). In this regime (Inc·Pe ∼ 10), diffusion can either hinder or enhance deposition time when combined with electrostatic effects (Figure 5a). Our non-dimensional plot (Figure 2b), which also accounts for the minor influence of gravity, predicts a deposition time of tf = 0.138 s, closely matching with the numerical simulations. As implied by Equation 9, the higher charge presented in Figure 5a inset, results in even shorter deposition times, reaching approximately 75 ms. In contrast, neutralized particles (Figure 5c) remained airborne for over a second and ultimately escaped the model, despite reaching 10 µm away from the tissue due to Brownian motion. A slightly higher charge would likely have ensured deposition and close this gap. These results highlight how electrostatics effects accelerate screening behavior and enhance deposition efficiency in the proximal generations of the bronchioles.
Equation (9) further highlights the critical role of mechanical mobility B, derived from the particle’s size and shape. Mechanical mobility playing a role as significant as the charge squared itself. This finding aligns well with clinical in vivo observations in humans, where DF for similar sizes (i.e. 0.3, 0.6 and 1 µm) has been found to follow a similar scaling relationship, (Melandri et al., 1983). This analytical result concurs with our finding of highly significant interaction between particle size and charge (p < 0.0001, for 0.2 β∼15% for Gen 1, see Supplementary MaterialSupplementary Table S6 for more) since B is linear with the particle diameter dp. In the future, more complex CFD models that include real lung anatomies (Koullapis et al., 2016) can leverage our approach to predict deposition location in the lungs based on the airborne time t of these two particle dynamic approaches.
We acknowledge that our airway-on-chip platform is based on a planar geometry, which does not fully capture the complex three-dimensional architecture of the human lung. However, the two primary mechanisms examined in this study, electrostatic attraction and diffusion, are largely independent of spatial orientation. Electrostatic forces follow electric field lines toward the nearest surface, while diffusion is isotropic and random. Therefore, in clinically relevant scenarios involving small, charged particles under low to moderate flow conditions, the deposition patterns observed in our planar model are anticipated to approximate sufficiently closely those in a realistic 3D airway structure. Nevertheless, we recognize that in regions with higher airflow velocities (e.g., the proximal bronchi), inertial effects may become more pronounced and sensitive to airway geometry. Additionally, in the distal bronchioles, gravitational settling becomes increasingly important for heavier particles, and the 3D orientation and curvature of the airways are likely to influence deposition. This is addressed below.
Role of electrostatic versus sedimentation
We now turn our attention to the largest particles, i.e. 2 μm. Recalling that H·Pe ∼ 235, these particles, when neutralized, are governed almost entirely by gravitational pull. In this regime, the planar orientation of our model becomes particularly relevant, as gravity acts in a fixed direction, and the absence of 3D airway curvature may influence the observed deposition patterns for clinical interpretation of the results. In order to distinguish the electrostatic forces from gravitational effects in our experiments, we conducted experiments with the model turned upside-down such that the ITO-coated surface is at the top, hence electrostatic may be interpreted as a “lifting force” against the downward pull of gravity. As expected, in the neutralized case, no deposition occurred on the top even when we extended the exposure time (see missing bars in Figures 4c,d). To validate this, as discussed in the Supplementary Material, we repeated the experiment with the model flipped to align with gravity and observed deposition biased by gravity direction, both with and without charge yields significant deposition which is hard to distinguish (see Supplementary MaterialSupplementary Figure S4c and S4d). This confirms the significant role of gravity in the bronchiolar region, as suggested previously (Ma and Darquenne, 2012). As expected, gravitational deposition occurs only in the direction of the gravitational pull. This can be visually shown in our numerical analysis of particle trajectories (see Figure 5d, Inc/H ∼ 0.01) where particles do not reach deposition at the top of the model when electrostatic forces and diffusion are absent. The particle follows its free-fall trajectory (see Figure 5d, dashed blue line) away from the tissue, with only slight deviations caused by weak diffusion. Its airborne time is approximately half a second (0.48 s), which is shorter compared to the 0.2 µm diffusive particle (1.200 s). This difference arises because the free-falling particle enters closer to the channel centerline, leading it to gain horizontal velocity as it moves into the high-speed region of the parabolic Poiseuille flow profile, resulting in a faster outward convection.
When these same 2 μm particles are charged, they have a chance to deposit. A distinct electrostatic-driven screening effect emerges leading to deposition within the bronchial model (Figure 4b). This suggests that electrostatic forces are enough to counteract gravity, encouraging greater bronchiolar deposition, even in regions where gravitational settling is unlikely to occur. We recall that the characteristic length scale represents the “Equilibrium Distance” from the wall at which the electrostatic force balances gravity. This length scale is crucial for understanding particle dynamics under the influence of both electrostatic forces and gravitational effects. Mathematically, it serves as a so-called separatrix, distinguishing between two distinct scenarios: (i) deposition due to electrostatic attraction versus (ii) convective transport through the channel to the next-generation. Most likely, in our case, this results in (; Whitby and Liu, 1968). Based on our in-vitro deposition findings (Figures 4c,d), we can infer that regions where are populated with aerosols that are most likely deposited due to electrostatic contributions. In comparison to the screening observed for the 0.2 µm particles, the gradient of the screening is less intense here. This can be attributed to the fact that while Inc⋅Pe is not dependent on particle size, the ratio Inc/H ∝ dp−3 is highly size-dependent. For the larger 2 µm particles, this strong dependence reduces the effectiveness of electrostatic deposition, as their increased dimensions hinders the influence of electrostatic forces. This is clearly shown and aligns with an in vitro study () and an additional upper airway CFD study, where Inc/Stk ∝ dp−3 follows the same dependency (Koullapis et al., 2016). As a result, more substantial deposition in proximal airways occurs for smaller particles, reinforcing the weight of electrostatic forces for smaller aerosols (Koullapis et al., 2016). Accordingly, a highly significant interaction between particle size and charge (p < 0.0001, for 2 β∼21% for Gen 1, see Supplementary MaterialSupplementary Table S6 for more), indicating that the effect of particle size on deposition depends on charge, and vice versa. Hence, charge does not influence all particle sizes uniformly, certain sizes are more affected by electrostatic forces than others.
Although electrostatics generally increased overall targeted bronchiolar deposition, it is important to note that across all tested combinations of charge states (neutralized and charged) and particle sizes (0.2, 0.5, and 2 µm), no significant increase in deposition was observed within the constricted regions compared to the normal airways. This finding indicates that, within the limited range of particle sizes and electrostatic charge levels explored in this study, we did not identify conditions that preferentially enhance deposition in constricted regions. Specifically, for 2 µm particles, deposition in the constricted airway was consistently about half of that observed in the normal airway, regardless of charge. These results underscore the ongoing challenge of achieving efficient, targeted aerosol delivery to the bronchioles, especially in distinguishing deposition behavior between constricted and non-constricted geometries. Future studies should explore a broader range of particle sizes, electrostatic charge levels, and airway constriction severities, as this work scope focused mainly on the effect of electrostatics in the bronchiole region and consider only one case of a mild constriction. Given that electrostatic forces become increasingly dominant at smaller length scales, we hypothesize that a promising window of opportunity remains to be uncovered—one that could enable selective targeting of diseased airways through optimized charge and size combinations.
In Figure 5b, we examine the particle trajectory under conditions where electrostatic and gravitational forces are of comparable magnitude (Inc ∼ H), as suggested by Equation (8). The simulation reveals that charged particles can ascend against gravity and deposit on airway surfaces that would otherwise remain unaffected. This behavior aligns with our in vitro findings (Figure 4a), where charged particles exhibited enhanced deposition compared to their neutralized counterparts. Notably, since the release height is close to the equilibrium distance (), the particle spends most of its airborne time hovering near this equilibrium position. However, beyond a critical threshold, the electrostatic force becomes strong enough to overcome gravitational pull entirely, and the particle trajectory converges with the analytical prediction of Equation (9) (green dashed line). As a result, the final velocity of the particle moving upward against gravity surpasses its free-fall terminal velocity by three orders of magnitude (). In summary, with clinical relevance, electrostatic forces expand the effective bronchial tissue area available for aerosol delivery by enabling particles to reach the top surfaces throughout the airway cross-section, and not only those who favor gravity direction. Hence, electrostatics demonstrate substantial potential for enhancing drug delivery efficiency, particularly in regions where conventional mechanisms fall short.
As a final step, we demonstrated the non-dimensional dependency (Figure 2b) that applies to all deposition scenarios involving electrostatic forces versus gravity. The result of our theoretical non-dimensional plot perfectly aligns with the deterministic trajectory and the one predicted by the numerical simulation, yielding a deposition time of = 0.239 s in both approaches. This framework serves as a valuable and comprehensive summary, capturing the interplay between particle shape, size, charge, and initial conditions through two dimensionless parameters: y0/yE and vtertf/yE. This relationship provides an effective tool for estimating the deposition time tf, which can be further related to the convective timescale. This representation also accurately predicts deposition times for our specific cases, such as 0.2 and 2 µm particles, showing strong agreement with our numerical solutions. Notably, deposition time approaches infinity when particles are uncharged or if they are initially positioned near the equilibrium stagnation point . The initial height plays a crucial role in determining deposition time. In constricted airways, particles are more likely to start from a lower initial height, smaller distance leading to shorter deposition periods. Additionally, as observed, the particle’s terminal velocity , influenced by its shape, remains a key factor in deposition dynamics. While our simulations assume spherical particles, it is important to note that alternative particle shapes could exhibit significantly different behaviors (Kleinstreuer et al., 2008; Kleinstreuer and Zhang, 2010). For example, long, straight fibers interact differently with the drug medium (; ; Shachar-Berman et al., 2018). Asbestos fibers, for instance, with lengths of 50–200 μm, are known to deposit extensively in the lungs and cause severe health conditions (Timbrell, 1965; ). In the context of electrostatic, in vivo studies in rodents have found much higher acinar deposition in charged asbestos fibers compared to neutralized ones (). This suggests that combining non-spherical shaped particles with electrostatic calls for future investigation (). While most accepted approaches neglect the role of electrostatic effects, a growing body of evidence suggests that they play a significant role in altering deposition (Kwok and Chan, 2009).
Conclusion
Using a novel in vitro airway-on-chip model, we systematically explored how electrostatic charge alters deposition patterns compared to classical mechanisms such as Brownian motion and gravitational sedimentation. Our results support that charge plays a dominant role in aerosol deposition at the bronchiole length scale, particularly for smaller aerosols, where electrostatic attraction accelerates deposition and causes diffusion screening to occur earlier. For neutralized larger particles, gravity primarily dictates a downward trajectory, but when charged, electrostatic forces influence their motion by defying gravity and promoting deposition in wider areas such as the “ceilings” of the airway wall. Hence, an increase in charge enhances the effective deposition region in bronchioles for 2 µm particles, while for smaller particles, it accelerates their deposition by electrostatic screening in this region. Additionally, our numerical simulations and analytical derivations offer a valuable step for overcoming the complexity of introducing electrostatics into particle dynamics. They provide a visual representation of the transport mechanisms and highlight key findings, such how the aerosols’ mechanical mobility can counteract electrostatic attraction, potentially allowing aerosols to penetrate deeper into the distal lung regions. These findings have relevant implications for inhaled drug delivery, suggesting that electrostatic manipulation along with size could serve as a tunable parameters for optimizing deposition in targeted airway regions, particularly where traditional aerodynamic strategies fail.
This suggests the need for further investigation into a broader variety of constriction severity and charge-diameter combinations that may selectively favor deposition in diseased airways, offering a potential avenue that may improve therapeutic targeting in obstructive lung diseases beyond the limited cases examined in this study (). An open question remains regarding the spatial distribution and the accumulative effect of electrostatics throughout the lungs. Some studies suggest that electrostatic forces may enhance deposition in the upper airways, raising concerns that charged particles may be lost before reaching critical therapeutic regions like the bronchioles and alveoli (Melandri et al., 1983; Koullapis et al., 2016; ). Alternatively, the finding here may suggest that the bronchioles themselves may act as a screening mechanism, limiting deposition in the deeper acinar regions. To fully address this question, multiscale models are needed to identify the window of opportunity in which charged aerosols can provide therapeutic benefit. Although the effect of electrical charge has been recognized for nearly a century (Wilson, 1947), the role of electrostatic forces in pulmonary deposition remains fundamentally underexplored. Overshadowed by more widely studied mechanisms such as impaction, sedimentation, and diffusion, electrostatics has often received limited attention, potentially hindering its broader application in aerosol drug delivery. To fully harness its potential, a deeper understanding is needed of the conditions under which electrostatic forces can meaningfully influence particle behavior throughout the respiratory tract. Addressing these overarching questions may require the development of multiscale lung models capable of capturing the cumulative effects of electrostatic interactions across multiple airway generations (Kassinos and Sznitman, 2024). Future studies should focus on expanding the range of charge distributions, incorporating asymmetric airway geometries, and exploring alternative particle morphologies, including fibers and aggregates, to better capture real-world deposition dynamics. Ultimately, this work reinforces the importance of electrostatics as a fundamental mechanism in pulmonary aerosol transport.
Statements
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding authors.
Author contributions
RB: Writing – review and editing, Writing – original draft, Conceptualization, Project administration, Software, Supervision, Data curation, Methodology, Visualization, Investigation. JS: Writing – review and editing, Writing – original draft, Supervision. TM: Writing – review and editing, Software. RF: Writing – review and editing. OF: Writing – review and editing, Visualization, Software. SD: Writing – review and editing, Investigation, Data curation. TD: Writing – review and editing, Investigation, Data curation. BK: Writing – review and editing, Software. KK: Writing – review and editing. NK: Writing – review and editing.
Funding
The author(s) declare that financial support was received for the research and/or publication of this article. This research was supported by the Israel Science Foundation (Grant no. 1840/21) and the Nichia Corporation via a Tokushima University–Technion grant. The funder was not involved in the study design, collection, analysis, interpretation of data, the writing of this article, or the decision to submit it for publication.
Acknowledgments
The authors thank Morane Elbaz, Leo Ibghi, Daniel Malka and Amjad Okasha, for helpful support with experiments.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The author(s) declare that no Generative AI was used in the creation of this manuscript.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphys.2025.1621177/full#supplementary-material
References
1
AriA.FinkJ. B. (2020). Recent advances in aerosol devices for the delivery of inhaled medications. Expert Opin. Drug Deliv.17 (2), 133–144. 10.1080/17425247.2020.1712356
2
AsgharianB.YuC. P. (1988). Deposition of inhaled fibrous particles in the human lung. J. Aerosol Med. Depos. Clear. Eff. Lung1 (1), 37–50. 10.1089/jam.1988.1.37
3
AzhdarzadehM.OlfertJ. S.VehringR.FinlayW. H. (2014). Effect of electrostatic charge on oral-extrathoracic deposition for uniformly charged monodisperse aerosols. J. Aerosol Sci.68, 38–45. 10.1016/j.jaerosci.2013.11.002
4
BaileyA. G.HashishA. H.WilliamsT. J. (1998). Drug delivery by inhalation of charged particles. J. Electrostat.44 (1–2), 3–10. 10.1016/S0304-3886(98)00017-5
5
BalachandranW.MachowskiW.GauraE.HudsonC.SystemsE. (1997). Control of drug aerosol in human airways using electrostatic forces. J. Electrostat.3886 (97), 579–584. 10.1016/s0304-3886(97)00106-x
6
BesslerR.SznitmanJ. (2024). The potential of leveraging electrostatics for improved inhaled drug delivery to the lungs. Front. Med. Eng.1 (January), 1298251–1298258. 10.3389/fmede.2023.1298251
7
BesslerR.BhardwajS.MalkaD.FishlerR.SznitmanJ. (2023). Exploring the role of electrostatic deposition on inhaled aerosols in alveolated microchannels. Sci. Rep.13 (1), 23069–14. 10.1038/s41598-023-49946-w
8
ChoiJ.LeBlancL. J.ChoiS.HaghighiB.HoffmanE. A.O'ShaughnessyP.et al (2019). Differences in particle deposition between members of imaging-based asthma clusters. J. Aerosol Med. Pulm. Drug Deliv.32 (4), 213–223. 10.1089/jamp.2018.1487
9
CohenB. S.XiongJ. Q.AsgharianB.AyresL. (1995). Deposition of inhaled charged ultrafine particles in a simple tracheal model. J. Aerosol Sci.26 (7), 1149–1160. 10.1016/0021-8502(95)00039-f
10
CohenB. S.XiongJ. Q.FangC.-P.LiW. (1998). Deposition of charged particles on lung airways. Health Phys.74 (5), 554–560. 10.1097/00004032-199805000-00002
11
DarquenneC. (2012). Aerosol deposition in health and disease. J. Aerosol Med. Pulm. Drug Deliv.25 (3), 140–147. 10.1089/jamp.2011.0916
12
DasP.NofE.AmiravI.KassinosS. C.SznitmanJ. (2018). Targeting inhaled aerosol delivery to upper airways in children: insight from computational fluid dynamics (CFD). PLoS One13 (11), e0207711–e0207720. 10.1371/journal.pone.0207711
13
DavisJ. M. G.BoltonR. E.DouglasA. N.JonesA. D.SmithT. (1988). Effects of electrostatic charge on the pathogenicity of chrysotile asbestos. Br. J. Ind. Med.45 (5), 292–299. 10.1136/oem.45.5.292
14
de BoerA. H.HagedoornP.HoppentochtM.ButtiniF.GrasmeijerF.FrijlinkH. W. (2017). Dry powder inhalation: past, present and future. Expert Opin. Drug Deliv.14 (4), 499–512. 10.1080/17425247.2016.1224846
15
DiazA. A.OrejasJ. L.GrumleyS.NathH. P.WangW.DolliverW. R.et al (2023). Airway-occluding mucus plugs and mortality in patients with chronic obstructive pulmonary disease. JAMA329 (21), 1832–1839. 10.1001/jama.2023.2065
16
DonaldsonK.BrownG. M.BrownD. M.BoltonR. E.DavisJ. M. G. (1989). Inflammation generating potential of long and short fibre amosite asbestos samples. Br. J. Ind. Med.46 (4), 271–276. 10.1136/oem.46.4.271
17
Elias-KirmaS.Artzy-SchnirmanA.DasP.Heller-AlgaziM.KorinN.SznitmanJ. (2020). In situ-like aerosol inhalation exposure for cytotoxicity assessment using airway-on-chips platforms. Front. Bioeng. Biotechnol.8 (February), 91–13. 10.3389/fbioe.2020.00091
18
FeliciM.FilocheM.SapovalB. (2003). Diffusional screening in the human pulmonary acinus. J. Appl. Physiol.94 (5), 2010–2016. 10.1152/japplphysiol.00913.2002
19
FinlayW. H. (2001). The mechanics of inhaled pharmaceutical aerosols: an introduction. San Diego: Academic Press.
20
FinlayW. H. (2021). Deposition of aerosols in the lungs: particle characteristics. J. Aerosol Med. Pulm. Drug Deliv.34 (4), 213–216. 10.1089/jamp.2021.29040.whf
21
FishlerR.HofemeierP.EtzionY.DubowskiY.SznitmanJ. (2015). Particle dynamics and deposition in true-scale pulmonary acinar models. Sci. Rep.5, 14071–11. 10.1038/srep14071
22
ForsythB.LiuB. Y. H.RomayF. J. (1998). Particle charge distribution measurement for commonly generated laboratory aerosols. Aerosol Sci. Technol.28 (6), 489–501. 10.1080/02786829808965540
23
HarrisR.TimbrellV. (1975). The influence of fibre shape in lung deposition-mathematical estimates. Inhaled Part.4 Pt 1, 75–89.
24
HindsW. C.ZhuY. (2022). “Electrical properties,” in Aerosol technology: properties, behavior, and measurement of airborne particles. (Hoboken: John Wiley and Sons), 277–304.
25
HofemeierP.SznitmanJ. (2015). Revisiting pulmonary acinar particle transport: convection, sedimentation, diffusion, and their interplay. J. Appl. Physiol.118 (11), 1375–1385. 10.1152/japplphysiol.01117.2014
26
HofemeierP.Shachar-BermanL.Tenenbaum-KatanJ.FilocheM.SznitmanJ. (2016). Unsteady diffusional screening in 3D pulmonary acinar structures: from infancy to adulthood. J. Biomechanics49 (11), 2193–2200. 10.1016/j.jbiomech.2015.11.039
27
HoggJ. C. (2004). Pathophysiology of airflow limitation in chronic obstructive pulmonary disease. Lancet364 (9435), 709–721. 10.1016/S0140-6736(04)16900-6
28
HorsfieldK.DartG.OlsonD. E.FilleyG. F.CummingG. (1971). Models of the human bronchial tree. J. Appl. Physiol.31 (2), 207–217. 10.1152/jappl.1971.31.2.207
29
ICRP (1994). Human respiratory tract model for radiological protection. Ann. ICRP66. 1–3.
30
IjsebaertJ. C.GeerseK. B.MarijnissenJ. C. M.LammersJ.-W. J.ZanenP. (2001). Electro-hydrodynamic atomization of drug solutions for inhalation purposes. J. Appl. Physiol.91 (6), 2735–2741. 10.1152/jappl.2001.91.6.2735
31
KassinosS. C.SznitmanJ. (2024). Multiscale modeling of respiratory transport phenomena and intersubject variability. Annu. Rev. Fluid Mech.57, 141–165. 10.1146/annurev-fluid-031424-103721
32
KimV.CrinerG. J. (2013). Chronic bronchitis and chronic obstructive pulmonary disease. Am. J. Respir. Crit. Care Med.187 (3), 228–237. 10.1164/rccm.201210-1843CI
33
KleinstreuerC.ZhangZ. (2010). Airflow and particle transport in the human respiratory system. Annu. Rev. Fluid Mech.42, 301–334. 10.1146/annurev-fluid-121108-145453
34
KleinstreuerC.ZhangZ.LiZ. (2008). Modeling airflow and particle transport/deposition in pulmonary airways. Respir. Physiol. Neurobiol.163 (1–3), 128–138. 10.1016/j.resp.2008.07.002
35
KoullapisP. G.KassinosS. C.BivolarovaM. P.MelikovA. K. (2016). Particle deposition in a realistic geometry of the human conducting airways: effects of inlet velocity profile, inhalation flowrate and electrostatic charge. J. Biomech.49 (11), 2201–2212. 10.1016/j.jbiomech.2015.11.029
36
KwokP. C. L.ChanH. K. (2008). Effect of relative humidity on the electrostatic charge properties of dry powder inhaler aerosols. Pharm. Res.25 (2), 277–288. 10.1007/s11095-007-9377-2
37
KwokP. C. L.ChanH.-K. (2009). Electrostatics of pharmaceutical inhalation aerosols. J. Pharm. Pharmacol.61 (12), 1587–1599. 10.1211/jpp/61.12.0002
38
MaB.DarquenneC. (2012). Aerosol bolus dispersion in acinar airways - influence of gravity and airway asymmetry. J. Appl. Physiol.113 (3), 442–450. 10.1152/JAPPLPHYSIOL.01549.2011
39
McNultyW.UsmaniO. S. (2014). Techniques of assessing small airways dysfunction. Eur. Clin. Respir. J.1 (1), 25898. 10.3402/ecrj.v1.25898
40
MelandriC.TarroniG.ProdiV.De ZaiacomoT.FormignaniM.LombardiC. C. (1983). Deposition of charged particles in the human airways. J. Aerosol Sci.14 (5), 657–669. 10.1016/0021-8502(83)90070-8
41
MénacheM. G.HofmannW.AshgarianB.MillerF. J. (2008). Airway geometry models of children’s lungs for use in dosimetry modeling. Inhal. Toxicol.20 (2), 101–126. 10.1080/08958370701821433
42
OakesJ. M.AmiravI.SznitmanJ. (2023). Pediatric inhalation therapy and the aerodynamic rationale for age-based aerosol sizes. Expert Opin. Drug Deliv.00 (00), 1037–1040. 10.1080/17425247.2023.2209314
43
OsmanH.CastleG. S. P.AdamiakK. (2015). Numerical study of particle deposition in electrostatic painting near a protrusion or indentation on a planar surface. J. Electrostat.77, 58–68. 10.1016/j.elstat.2015.07.005
44
PedleyT. J. (1977). Pulmonary fluid dynamics. Annu. Rev. Fluid Mech.9, 229–274. 10.1146/annurev.fl.09.010177.001305
45
RaabeO. G. (1968). The dilution of monodisperse suspensions for aerosolization. Am. Ind. Hyg. Assoc. J.29 (5), 439–443. 10.1080/00028896809343031
46
Shachar-BermanL.OstrovskiY.De RosisA.KassinosS.SznitmanJ. (2018). Transport of ellipsoid fibers in oscillatory shear flows: implications for aerosol deposition in deep airways. Eur. J. Pharm. Sci.113, 145–151. 10.1016/j.ejps.2017.09.023
47
StahlhofenW.RudolfG.JamesA. C. (1989). Intercomparison of experimental regional aerosol deposition data. J. Aerosol Med. Depos. Clear. Eff. Lung2 (3), 285–308. 10.1089/jam.1989.2.285
48
SznitmanJ. (2013). Respiratory microflows in the pulmonary acinus. J. Biomech.46 (2), 284–298. 10.1016/j.jbiomech.2012.10.028
49
SznitmanJ. (2022). Revisiting airflow and aerosol transport phenomena in the deep lungs with microfluidics. Chem. Rev.122 (7), 7182–7204. 10.1021/acs.chemrev.1c00621
50
TimbrellV. (1965). Human exposure to asbestos: dust controls and standards. The inhalation of fibrous dusts. Ann. N. Y. Acad. Sci.132 (1), 255–273. 10.1111/j.1749-6632.1965.tb41107.x
51
TsudaA.HenryF. S. S.ButlerJ. P. P. (2013). Particle transport and deposition: basic physics of particle kinetics. Compr. Physiol.3 (4), 1437–1471. 10.1002/cphy.c100085
52
VirchowJ. C.PoliG.HerpichC.KietzigC.EhlichH.BraeutigamD.et al (2018). Lung deposition of the dry powder fixed combination beclometasone dipropionate plus formoterol fumarate using NEXThaler ® device in healthy subjects, asthmatic patients, and COPD patients. J. Aerosol Med. Pulm. Drug Deliv.31 (5), 269–280. 10.1089/jamp.2016.1359
53
WangC.ZhouJ.WangJ.LiS.FukunagaA.YodoiJ.et al (2020). Progress in the mechanism and targeted drug therapy for COPD. Signal Transduct. Target. Ther.5 (1), 248. 10.1038/s41392-020-00345-x
54
WeibelE. R. (1963). Morphometry of the human lung with a foreword. Berlin, Heidelberg: Springer.
55
WhitbyK. T.LiuB. Y. H. (1968). Polystyrene aerosols-electrical charge and residue size distribution. Atmos. Environ.2 (2), 103–116. 10.1016/0004-6981(68)90026-7
56
WilsonI. B. (1947). The deposition of charged particles in tubes, with reference to the retention of therapeutic aerosols in the human lung. J. Colloid Sci.2 (2), 271–276. 10.1016/0095-8522(47)90028-7
57
XiJ.SiX.LongestW. (2014). Electrostatic charge effects on pharmaceutical aerosol deposition in human nasal-laryngeal airways. Pharmaceutics6 (1), 26–35. 10.3390/pharmaceutics6010026
58
XuF.VasilescuD. M.KinoseD.TanabeN.NgK. W.CoxsonH. O.et al (2022). The molecular and cellular mechanisms associated with the destruction of terminal bronchioles in COPD. Eur. Respir. J.59 (5), 2101411. 10.1183/13993003.01411-2021
Summary
Keywords
airway-on-chip, pulmonary drug delivery, obstructive pulmonary disease, non-dimensional analysis, aerosol inhalation, electrostatics
Citation
Bessler R, Mekler T, Fishler R, Farhana O, Dhatavkar S, Daniel T, Kalifa B, Koshiyama K, Korin N and Sznitman J (2025) Small scale electrostatically-driven aerosol deposition in airway-on-chip models of bronchial constriction. Front. Physiol. 16:1621177. doi: 10.3389/fphys.2025.1621177
Received
30 April 2025
Accepted
14 August 2025
Published
16 September 2025
Volume
16 - 2025
Edited by
Aranyak Chakravarty, Jadavpur University, India
Reviewed by
Jingliang Dong, Victoria University, Australia
Eknath Kole, Kavayitri Bahinabai Chaudhari North Maharashtra University, India
Updates
Copyright
© 2025 Bessler, Mekler, Fishler, Farhana, Dhatavkar, Daniel, Kalifa, Koshiyama, Korin and Sznitman.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Ron Bessler, sronbess@campus.technion.ac.il; Josué Sznitman, sznitman@bm.technion.ac.il
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.