Population genetics in microchannels.

SignificanceMany microbial populations proliferate in small channels. In such environments, reproducing cells organize in parallel lanes. Reproducing cells shift these lanes, potentially expelling other cells from the channel. In this paper, we combine theory and experiments to understand how these dynamics affects the diversity of a microbial population. We theoretically predict that genetic diversity is quickly lost along lanes of cells. Our experiments confirm that a population of proliferating Escherichia coli in a microchannel organizes into lanes of genetically identical cells within a few generations. Our findings elucidate the effect of lane formation on populations evolution, with potential applications ranging from microbial ecology in soil to dynamics of epithelial tissues in higher organisms.

Chemotactic Interactions Drive Migration of Membraneless Active Droplets.

In nature, chemotactic interactions are ubiquitous and play a critical role in driving the collective behavior of living organisms. Reproducing these interactions in vitro is still a paramount challenge due to the complexity of mimicking and controlling cellular features, such as tangled metabolic networks, cytosolic macromolecular crowding, and cellular migration, on a microorganism size scale. Here, we generate enzymatically active cell-sized droplets able to move freely, and by following a chemical gradient, able to interact with the surrounding droplets in a collective manner. The enzyme within the droplets generates a pH gradient that extends outside the edge of the droplets. We discovered that the external pH gradient triggers droplet migration and controls its directionality, which is selectively toward the neighboring droplets. Hence, by changing the enzyme activity inside the droplet, we tuned the droplet migration speed. Furthermore, we showed that these cellular-like features can facilitate the reconstitution of a simple and linear protometabolic pathway and increase the final reaction product generation. Our work suggests that simple and stable membraneless droplets can reproduce complex biological phenomena, opening new perspectives as bioinspired materials and synthetic biology tools.

Localized Surface Plasmon Resonance Sensing and its Interplay with Fluidics.

In this Feature Article, we discuss the interplay between fluidics and the localized surface plasmon resonance (LSPR) sensing technique, primarily focusing on its applications in the realm of bio/chemical sensing within fluidic environments. Commencing with a foundational overview of LSPR principles from a sensing perspective, we subsequently showcase the development of a streamlined LSPR chip integrated with microfluidic structures. This integration opens the doors to advanced experiments involving fluid dynamics, greatly expanding the scope of LSPR-based research. Our discussions then turn to the practical implementation of LSPR and microfluidics in real-time biosensing, with a specific emphasis on monitoring DNA polymerase activity. Additionally, we illustrate the direct sensing of biological fluids, exemplified by the analysis of urine, while also shedding light on a unique particle assembly process that occurs on LSPR chips. We not only discuss the significance of LSPR sensing but also explore its potential to investigate a plethora of phenomena at liquid-liquid and solid-liquid interfaces. This is particularly noteworthy, as existing transduction methods and sensors fall short in fully comprehending these interfacial phenomena. Concluding our discussion, we present a futuristic perspective that provides insights into potential opportunities emerging at the intersection of fluidics and LSPR sensing.

Understanding the Role of Loss Modulus of Viscoelastic Substrates in the Evaporation Dynamics of Sessile Drops.

Viscoelastic properties of soft substrates play a crucial role in the evaporation dynamics of sessile drops. Recent studies have revealed that the modification of the viscoelastic properties of substrates changes the dynamics of the three-phase contact line, consequently affecting the evaporation behavior of sessile drops. Notably, these modifications occur without any noticeable changes to the substrate's wetting characteristics or surface topography. However, the individual role of storage (G') and loss (G″) moduli of substrates on drop evaporation dynamics remains unexplored. In this study, we investigate the evaporation dynamics of water drops on two groups of poly(dimethylsiloxane)-based viscoelastic substrates possessing either identical G' with varying G″ or identical G″ with varying G'. Our study reveals that on a substrate with constant shear modulus (G'), a reduction of an order of magnitude in loss modulus shifts the evaporation process from the constant contact radius mode to the constant contact angle mode. We hypothesize that this observed shift in behavior stems from the varying viscoelastic dissipation influenced by the plateau modulus and characteristic relaxation time of polymer gels. Our hypothesis is further supported from the observation that the evaporation process persists on the substrate with constant loss modulus (G″). Our study advances the current understanding of drop evaporation on soft substrates that may find potential applications involving soft composites, biological entities, tissue engineering, and wearable electronics.

Dilute polyelectrolyte solutions: recent progress and open questions.

Polyelectrolytes are a class of polymers possessing ionic groups on their repeating units. Since counterions can dissociate from the polymer backbone, polyelectrolyte chains are strongly influenced by electrostatic interactions. As a result, the physical properties of polyelectrolyte solutions are significantly different from those of electrically neutral polymers. The aim of this article is to highlight key results and some outstanding questions in the polyelectrolyte research from recent literature. We focus on the influence of electrostatics on conformational and hydrodynamic properties of polyelectrolyte chains. A compilation of experimental results from the literature reveals significant disparities with theoretical predictions. We also discuss a new class of polyelectrolytes called poly(ionic liquid)s that exhibit unique physical properties in comparison to ordinary polyelectrolytes. We conclude this review by listing some key research challenges in order to fully understand the conformation and dynamics of polyelectrolytes in solutions.

Stagnation points control chaotic fluctuations in viscoelastic porous media flow.

Viscoelastic flows through porous media become unstable and chaotic beyond critical flow conditions, impacting widespread industrial and biological processes such as enhanced oil recovery and drug delivery. Understanding the influence of the pore structure or geometry on the onset of flow instability can lead to fundamental insights into these processes and, potentially, to their optimization. Recently, for viscoelastic flows through porous media modeled by arrays of microscopic posts, Walkama et al. [D. M. Walkama, N. Waisbord, J. S. Guasto, Phys. Rev. Lett 124, 164501 (2020)] demonstrated that geometric disorder greatly suppressed the strength of the chaotic fluctuations that arose as the flow rate was increased. However, in that work, disorder was only applied to one originally ordered configuration of posts. Here, we demonstrate experimentally that, given a slightly modified ordered array of posts, introducing disorder can also promote chaotic fluctuations. We provide a unifying explanation for these contrasting results by considering the effect of disorder on the occurrence of stagnation points exposed to the flow field, which depends on the nature of the originally ordered post array. This work provides a general understanding of how pore geometry affects the stability of viscoelastic porous media flows.

Torsional fracture of viscoelastic liquid bridges.

Short liquid bridges are stable under the action of surface tension. In applications like electronic packaging, food engineering, and additive manufacturing, this poses challenges to the clean and fast dispensing of viscoelastic fluids. Here, we investigate how viscoelastic liquid bridges can be destabilized by torsion. By combining high-speed imaging and numerical simulation, we show that concave surfaces of liquid bridges can localize shear, in turn localizing normal stresses and making the surface more concave. Such positive feedback creates an indent, which propagates toward the center and leads to breakup of the liquid bridge. The indent formation mechanism closely resembles edge fracture, an often undesired viscoelastic flow instability characterized by the sudden indentation of the fluid's free surface when the fluid is subjected to shear. By applying torsion, even short, capillary stable liquid bridges can be broken in the order of 1 s. This may lead to the development of dispensing protocols that reduce substrate contamination by the satellite droplets and long capillary tails formed by capillary retraction, which is the current mainstream industrial method for destabilizing viscoelastic liquid bridges.

Duplex Electrochemical Microfluidic Sensor for COVID-19 Antibody Detection: Natural versus Vaccine-Induced Humoral Response.

The rapid transmission and resilience of coronavirus disease 2019 (COVID-19) have led to urgent demands in monitoring humoral response for effective vaccine development, thus a multiplex co-detection platform to discriminate infection-induced from vaccine-induced antibodies is needed. Here a duplex electrochemical immunosensor for co-detection of anti-nucleocapsid IgG (N-IgG) and anti-spike IgG (S-IgG) is developed by using a two-working electrode system, via an indirect immunoassay, with antibody quantification obtained by differential pulse voltammetry. The screen-printed electrodes (SPEs) are modified by carbon black and electrodeposited gold nanoflowers for maximized surface areas, enabling the construction of an immunological chain for S-IgG and N-IgG electrochemical detection with enhanced performance. Using an optimized immunoassay protocol, a wide linear range between 30-750 and 20-1000 ng mL-1 , and a limit of detection of 28 and 15 ng mL-1 are achieved to detect N-IgG and S-IgG simultaneously in serum samples. This duplex immunosensor is then integrated in a microfluidic device to obtain significantly reduced detection time (≤ 7 min) while maintaining its analytical performance. The duplex microfluidic immunosensor can be easily expanded into multiplex format to achieve high throughput screening for the sero-surveillance of COVID-19  and other infectious diseases.

Alignment-Rheology Relationship of Biosourced Rod-Like Colloids and Polymers under Flow.

Fluids composed of biosourced rod-like colloids (RC) and rod-like polymers (RP) have been extensively studied due to various promising applications relying on their flow-induced orientation (e.g., fiber spinning). However, the relationship between RC and RP alignment and the resulting rheological properties is unclear due to experimental challenges. We investigate the alignment-rheology relationship for a variety of biosourced RC and RP, including cellulose-based particles, filamentous viruses, and xanthan gum, by simultaneous measurements of the shear viscosity and fluid anisotropy under rheometric shear flows. For each system, the RC and RP contribution to the fluid viscosity, captured by the specific viscosity ηsp, follows a universal trend with the extent of the RC and RP alignment independent of concentration. We further exploit this unique rheological-structural link to retrieve a dimensionless parameter (ß) directly proportional to ηsp at zero shear rate (η0,sp), a parameter often difficult to access from experimental rheometry for RC and RP with relatively long contour lengths. Our results highlight the unique link between the flow-induced structural and rheological changes occurring in RC and RP fluids. We envision that our findings will be relevant in building and testing microstructural constitutive models to predict the flow-induced structural and rheological evolution of fluids containing RC and RP.

Rheological effects on purely-elastic flow asymmetries in the cross-slot geometry.

Viscoelastic flows in the cross-slot geometry can undergo a transition from a steady symmetric to a steady asymmetric flow state, ostensibly due to purely-elastic effects arising beyond a critical flow rate, or Weissenberg number Wi. However, some reports suggest that shear thinning of the fluid's viscosity may also play an important role in this transition. We employ a series of polymer solutions of varying rheological properties to investigate in detail how the interplay between fluid elasticity and shear thinning affects the onset and development of asymmetric flows in the cross-slot. Flow velocimetry is performed on each of the polymer solutions, and is used to assess the degree of flow asymmetry I in the cross-slot as a function of both Wi and a dimensionless parameter S quantifying the flow-rate-dependent extent of shear thinning. Typically, the flow field breaks symmetry as Wi is increased beyond a critical value, but the magnitude of I is found to also be dependent on S. For a few specific polymer solutions, the flow field recovers symmetry above a second, higher critical Wi as S becomes small. The experimental results are summarized in a flow state diagram in Wi-S space, showing the relationship between flow asymmetry and fluid rheology. Finally, to gain a deeper understanding of the effects of shear thinning, numerical simulations are performed using the linear simplified Phan-Thien-Tanner model. We demonstrate that the degree of both shear thinning and elasticity of the fluid, and their interplay, are important factors controlling elastic instabilities in the cross-slot geometry.

Transition between solid and liquid state of yield-stress fluids under purely extensional deformations.

We report experimental microfluidic measurements and theoretical modeling of elastoviscoplastic materials under steady, planar elongation. Employing a theory that allows the solid state to deform, we predict the yielding and flow dynamics of such complex materials in pure extensional flows. We find a significant deviation of the ratio of the elongational to the shear yield stress from the standard value predicted by ideal viscoplastic theory, which is attributed to the normal stresses that develop in the solid state prior to yielding. Our results show that the yield strain of the material governs the transition dynamics from the solid state to the liquid state. Finally, given the difficulties of quantifying the stress field in such materials under elongational flow conditions, we identify a simple scaling law that enables the determination of the elongational yield stress from experimentally measured velocity fields.

Rheological scaling of ionic-liquid-based polyelectrolytes in ionic liquid solutions: the effect of the ion diameter of ionic liquids.

We investigate the effect of the ion diameter a of ionic liquids (ILs) on the shear viscosity of polymerized ionic liquids (PILs) in IL solutions. When both the PIL and IL contain large PFSI anions (a ≈ 0.57 nm), the specific viscosity ηsp first decreases with increasing IL concentration cIL in the low cIL regime, reaches a minimum and then increases with increasing cIL in the high cIL regime. By comparing the measured ηsp with the modified charge screening model proposed in our previous study [Matsumoto et al., Macromolecules, 2021, 54, 5648-5661], we attribute the observed non-monotonic trend of ηsp against cIL to the charge underscreening phenomenon, i.e., an increase of the screening length at high cIL leads to the upturn of ηsp. On the other hand, when the PIL and IL contain small BF4 anions (a ≈ 0.34 nm), the ηsp decreases asymptotically with increasing cIL, because the charge on the PIL chain is likely screened fully in the entire cIL regime. Our results demonstrate that the ion diameter of ILs plays an important role in governing the charge screening mechanism of PILs in IL solutions, and thus influencing the viscoelastic properties of PIL solutions.

Upstream wall vortices in viscoelastic flow past a cylinder.

We report a novel inertia-less, elastic flow instability for a viscoelastic, shear-thinning wormlike micellar solution flowing past a microcylinder in a channel with blockage ratio BR = 2R/W = 0.5 and aspect ratio α = H/W ≈ 5, where R ≈ 100 µm is the cylinder radius, W is the channel width, and H is the channel height. The instability manifests upstream of the cylinder and changes form with increasing Weissenberg number over the range 0.5 ≲ Wi = Uλ/R ≲ 900, where U is the average flow velocity and λ is the terminal relaxation time of the fluid. Beyond a first critical Wi, the instability begins as a bending of the streamlines near the upstream pole of the cylinder that breaks the symmetry of the flow. Beyond a second critical Wi, small, time-steady, and approximately symmetric wall-attached vortices form upstream of the cylinder. Beyond a third critical Wi, the flow becomes time dependent and pulses with a characteristic frequency commensurate with the breakage timescale of the wormlike micelles. This is accompanied by a breaking of the symmetry of the wall-attached vortices, where one vortex becomes considerably larger than the other. Finally, beyond a fourth critical Wi, a vortex forms attached to the upstream pole of the cylinder whose length fluctuates in time. The flow is highly time dependent, and the cylinder-attached vortex and wall-attached vortices compete dynamically for space and time in the channel. Our results add to the rapidly growing understanding of viscoelastic flow instabilities in microfluidic geometries.

Torsional instability of constant viscosity elastic liquid bridges.

By experiment and simulation, we report that viscoelastic liquid bridges made of constant viscosity elastic liquids, a.k.a. Boger fluids, can be effectively destabilized by torsion. Under torsion, the deformation of the liquid bridge depends on the competition between elastocapillarity and torsion-induced normal stress effects. When the elastocapillary effect dominates, the liquid bridge undergoes elastocapillary instability and thins into a cylindrical thread, whose length increases and whose radius decays exponentially over time. When the torsion-induced normal stress effect dominates, the liquid bridge deforms in a way similar to edge fracture, a flow instability characterized by the sudden indentation of the fluid's free surface when a viscoelastic fluid is sheared at above a critical deformation rate. The vertical component of the normal stress causes the upper and lower portions of the liquid bridge to approach each other, and the radial component of the normal stress results in the liquid bridge thinning more quickly than under elastocapillarity. Whether such quick thinning continues until the bridge breaks depends on both the liquid bridge configuration and the level of torsion applied.

Effect of blockage ratio on flow of a viscoelastic wormlike micellar solution past a cylinder in a microchannel.

We present experiments on the flow of a viscoelastic wormlike micellar solution around cylinders (radius R) confined in straight microchannels (width W). Thirteen flow geometries are tested where the blockage ratio is varied over a wide range 0.055 ≤ BR = 2R/W ≤ 0.63. Experiments are performed at negligible Reynolds number, and for Weissenberg numbers Wi = λU/R up to 1000, where U is the average flow speed and λ is the relaxation time of the fluid. Micro-particle image velocimetry is used to characterise the flow state at each BR and Wi. In all of the geometries, a first critical Weissenberg number marks a transition from symmetric flow to an asymmetric but time-steady flow state, while a second higher critical Weissenberg number marks the onset of time-dependent flows. However, we report a clear shift in behaviour over a narrow intermediate range of 0.33 ≲ BR ≲ 0.41. Channels with BR ≤ 0.33 fall in a 'low' BR regime, with instabilities that originate from the downstream stagnation point, while those with BR ≥ 0.44 fall in a 'high' BR regime, with instabilities developing at the upstream stagnation point. Behaviour within the newly-identified intermediate BR regime is complex due to the competing influence of the two stagnation points. We summarise all our results in a flow state diagram covering Wi-BR parameter space, clearly defining the different regimes of blockage ratio for the first time. Our results contribute to the understanding of the complexities of viscoelastic flow in this benchmark geometry.

Tristability in Viscoelastic Flow Past Side-by-Side Microcylinders.

Viscoelastic flows through microscale porous arrays exhibit complex path selection and switching phenomena. However, understanding this process is limited by a lack of studies linking between a single object and large arrays. Here, we report experiments on viscoelastic flow past side-by-side microcylinders with variable intercylinder gap. With increasing flow rate, a sequence of two imperfect symmetry-breaking bifurcations forces selection of either one or two of the three possible flow paths around the cylinders. Tuning the gap length through the value where the first bifurcation becomes perfect reveals regions of bistability and tristability in a dimensionless flow rate-gap length phase diagram.

Shear thickening behavior in dense repulsive and attractive suspensions of hard spheres.

Shear thickening in stable dense colloidal suspensions is a reversible phenomenon and no hysteresis is observed in the flow curve measurements. However, a reduction in the stability of colloids promotes particle aggregation and introduces a time dependent rheological response. In this work, by using a model colloidal system of hard spherical silica particles (average diameter of 415 nm) with varying particle volume fractions 0.2 ≤ ϕ ≤ 0.56, we study the effect of particle stability on the hysteresis of the shear thickening behavior of these suspensions. The particle stability is manipulated by adding a simple monovalent salt (sodium chloride) in the silica suspension with varying concentrations α ∈ [0,0.5] M. For repulsive and weakly attractive suspensions, the flow behavior is history independent and the shear thickening behavior does not exhibit hysteresis. However, significant hysteresis is observed in rheological measurements for strongly attractive suspensions, with shear history playing a critical role due to the dynamic nature of particle clusters, resulting in time dependent hysteresis behavior. By performing numerical simulations, we find that this hysteresis behavior arises due to the competition among shear, electrostatic repulsive, van der Waals attractive, and frictional contact forces. The critical shear stress (i.e., the onset of shear thickening) decreases with increasing salt concentrations, which can be captured by a scaling relationship based on the force balance between particle-particle contact force and electrostatic repulsive force. Our combined experimental and simulation results imply the formation of particle contacts in our sheared suspensions.

Rheology of the Electric Double Layer in Electrolyte Solutions.

Electric double layers (EDLs) are ionic structures formed on charged surfaces and play an important role in various biological and industrial processes. An extensive study in the past decade has revealed the structure of the EDL in concentrated electrolyte solutions of both ordinary salts and ionic liquids. However, how the EDL structure affects their material properties remains a challenging topic due to technical difficulties of these measurements at nanoscale. In this work, we report the first detailed characterization of the viscoelasticity of the EDL formed over a wide range of ion concentrations, including concentrated electrolyte solutions. Specifically, we investigate the complex shear modulus of the EDL by measuring the resonant frequency and the energy dissipation of a quartz crystal microbalance (QCM), a surface-sensitive device, immersed in aqueous solutions containing three types of solutes: an ionic liquid, 1-butyl-3-methylimidazolium chloride (BmimCl); an ordinary salt, sodium chloride (NaCl); and a nonelectrolyte, ethylene glycol (EG). For the two electrolyte solutions, we observe a monotonic decrease in the resonant frequency and a monotonic increase in the energy dissipation with increasing ion concentrations due to the presence of the EDL. The complex shear modulus of the EDL is estimated through a wave propagation model in which the density and shear modulus of the EDL decay exponentially toward those of the bulk solution. Our results show that both the storage and the loss modulus of the EDL increase rapidly with increasing ion concentrations in the low ion concentration regime (<1 M) but reach saturation values with similar magnitude at a sufficiently high ion concentration. The shear viscosity of the EDL near the charged QCM surface is approximately 50 times for NaCl solutions and 500 times for BmimCl solutions of the bulk solution value at the saturation concentration. We also demonstrate that QCM can be utilized for analyzing the rheological properties of the EDL, thus providing a complementary, low-cost, and portable alternative to conventional laboratory instruments such as the surface force apparatus. Our results elucidate new perspectives on the viscoelastic properties of the EDL and can potentially guide device optimization for applications such as biosensing and fast charging of batteries.

Purely Elastic Fluid-Structure Interactions in Microfluidics: Implications for Mucociliary Flows.

Fluid-structure interactions lie at the heart of the complex, and often highly coordinated, motions of actively driven microscale biological systems (e.g., translating cilia, flagella, and motile cells). Due to the highly viscoelastic nature of most relevant biological fluids and the small length scales involved, the viscous and inertial forces in such flows are dominated by elasticity. However, elastic effects are often overlooked in studies seeking to address phenomena like the synchronization of beating cilia. In this study, unique microfluidic experiments are presented to demonstrate that inertia-free viscoelastic flows can lead to highly regular beating of an immersed (passive) flexible structure, herein named "purely-elastic" fluid-structure interaction. It is also shown how two such flexible structures can achieve an extraordinary degree of synchronization, with a correlation coefficient approaching unity. The synchronization is a result of the generation of localized elastic stresses in the fluid that effectively link the two objects. These purely elastic interactions may be important to consider toward developing a complete understanding of the motions of microscale biological systems.

Optimized Immobilization of Biomolecules on Nonspherical Gold Nanostructures for Efficient Localized Surface Plasmon Resonance Biosensing.

Plasmonic biosensing techniques employ metal nanostructures, commonly gold (Au), often with biomolecules attached to their surfaces either directly or via other linkers. Various surface chemistry methods based on dispersion and covalent interactions are used to attach biomolecules to Au. As a result, when immobilizing a molecule on a metal surface, quantitative estimates of binding efficiency and stability of these surface chemistry methods are needed. Most prior work to compare such methods deals with bulk/thin film configurations or spherical nanoparticles, and very little is known about immobilization of biomolecules on plasmonic nanostructures of different shapes. Besides, due to rapid advancement of modern nanofabrication techniques, there is a growing need to determine an efficient surface chemistry method for immobilization of biomolecules on nonspherical plasmonic nanostructures. Previous comparison of immobilization methods on spherical Au nanoparticles has shown that physical adsorption resulted in the highest concentration of immobilized antibodies. In our work, we conducted a similar study and compared four representative Au surface functionalization methods as well as estimated how efficient these methods are at attaching biomolecules to nonspherical plasmonic Au nanostructures. We estimated the concentration of immobilized antibody that is specific to human C-reactive protein (anti-hCRP) by measuring the localized surface plasmon resonance (LSPR) shifts after exposing the surface of Au nanostructures to the antibody. Our results differ from the previously reported ones since the highest concentration of anti-hCRP was immobilized using 11-mercaptoundecanoic acid (MUA) chemistry. We demonstrated that immobilized antibodies retained their stability and specificity toward hCRP throughout the immunoassay when diluted hCRP or hCRP-spiked human serum samples were used. These findings have important implications for the fields of biosensing and diagnostics that employ nonspherical plasmonic nanostructures since an overall performance of these devices depends on efficient biomolecule immobilization.