Analysis of Gas Nanoclusters in Water Using All-Atom Molecular Dynamics.

The Young-Laplace equation suggests that nanosized gas clusters would dissolve under the effects of perturbation. The fact that nanobubbles are observed raises questions as to the mechanism underlying their stability. In the current study, we used all-atom molecular dynamics simulations to investigate the gas-water interfacial properties of gas clusters. We employed the instantaneous coarse-graining method to define the fluctuating boundaries and analyze the deformation of gas clusters. Fourier transform analysis of the cluster morphology revealed that the radius and morphology deformation variations exhibit power law relationships with the vibrational frequency, indicating that the surface energy dissipated through morphology variations. Increasing pressure in the liquid region was found to alter the network of water molecules at the interface, whereas increasing pressure in the gas region did not exhibit this effect. The overall gas concentration was oversaturated and proportional to the gas density inside the clusters. However, the result of comparison with Henry's law reveals that the gas pressure at the interface reduced by the interfacial effects is much lower than that inside the gas region, thus reducing the demanding degree of oversaturation. Originating from the interfacial charge allocation, the magnitude of the electrostatic stress is greater than that of the gas pressure inside the cluster. However, the magnitude of the reversed tension induced by electrostatic stress is far below the value of interfacial tension. The potential of mean force (PMF) profiles revealed that a barrier potential at the interface hindered gas particles from escaping the cluster. Several effects contribute to stabilizing the gas clusters in water, including high-frequency morphological deformation, electrostatic stress, reduced interfacial tension, and gas oversaturation conditions. Our results suggest that gas clusters can exist in water under gas oversaturation conditions in the absence of hydrophobic contaminants or pinning charges at interfaces.

Effects of Gas Adsorption and Surface Conditions on Interfacial Nanobubbles.

Gas aggregation and formation of interfacial nanobubbles (INBs) provide challenges and opportunities in the operation of micro-/nanofluidic devices. In the current study, we used molecular dynamics(MD) simulations to investigate the effects of hydrophobicity and various homogeneous surface conditions on gas aggregation and INB stability with a series of 3D argon-water-solid and water-solid systems. Among various signatures of surface hydrophobicity, the potential of mean force (PMF) minima exhibited the strongest correlation with the water molecular orientation at the liquid-solid interface, compared to the depletion layer width and the droplet contact angle. Our results indicated that argon aggregation on the substrate was a function of hydrophobicity as well as competition between gas-solid and water-solid PMFs. Thus, one precondition for gas aggregation on a surface is that the free energy minima of gas induced by the surface be much lower than that induced by water. We found that although the presence of gas molecules had little effect on the measures of wettability, it enhanced density fluctuations near liquid-solid interfaces. The PMF of gas along the surface tangential plane exhibited a small energy barrier between the epitaxial gas layer (EGL) in the bubble and the gas enrichment layer (GEL) in the liquid, which may benefit nanobubble stability. Much lower PMF in the EGL compared to that in the GEL indicated that gas molecules could migrate from the GEL to the nanobubble basement. However, density fluctuations enhanced by the GEL could reduce the energy barrier, thus reducing the stability of INBs.

Rich phase transitions in strongly confined polymer-nanoparticle mixtures: Nematic ordering, crystallization, and liquid-liquid phase separation.

We investigate the rich phase behavior of strongly confined semi-flexible (SFC) polymer-nanoparticle (NP) systems using the graphics processing unit accelerated Langevin dynamics simulation. Hard nanoparticles (HNP) that repel each other and ideal nanoparticles (INP) that do not interact with the same species are used as model additives to a strongly confined semiflexible polymer fluid. Both types of NPs exclude the monomer beads in the same way, but they have qualitatively different effects on the SFC isotropic-nematic (I-N) transition. For the total volume fraction ϕtot < 0.16, adding a low volume fraction of HNPs (ϕp) disrupts the long range nematic order of the polymers, whereas adding HNPs in a moderately packed system (0.16 < ϕtot < 0.32) facilitates polymer alignment due to the restricted polymer orientational degree of freedom. For dense packing (ϕtot > 0.32), polymers and NPs separate into layers along the slit height and the NPs form crystalline microdomains. In contrast, INP additives always promote inter-polymer alignment for low to moderate monomer volume fractions (ϕm). Furthermore, we found that INPs form a droplet-like fluid domain in dense nematic polymer systems.

Clusters of circulating tumor cells traverse capillary-sized vessels.

Multicellular aggregates of circulating tumor cells (CTC clusters) are potent initiators of distant organ metastasis. However, it is currently assumed that CTC clusters are too large to pass through narrow vessels to reach these organs. Here, we present evidence that challenges this assumption through the use of microfluidic devices designed to mimic human capillary constrictions and CTC clusters obtained from patient and cancer cell origins. Over 90% of clusters containing up to 20 cells successfully traversed 5- to 10-µm constrictions even in whole blood. Clusters rapidly and reversibly reorganized into single-file chain-like geometries that substantially reduced their hydrodynamic resistances. Xenotransplantation of human CTC clusters into zebrafish showed similar reorganization and transit through capillary-sized vessels in vivo. Preliminary experiments demonstrated that clusters could be disrupted during transit using drugs that affected cellular interaction energies. These findings suggest that CTC clusters may contribute a greater role to tumor dissemination than previously believed and may point to strategies for combating CTC cluster-initiated metastasis.

Electrofluidic Circuit-Based Microfluidic Viscometer for Analysis of Newtonian and Non-Newtonian Liquids under Different Temperatures.

This paper reports a microfluidic viscometer with an integrated pressure sensor based on electrofluidic circuits, which are electrical circuits constructed by ionic liquid-filled microfluidic channels. The electrofluidic circuit provides a pressure-sensing scheme with great long-term and thermal stability. The viscosity of the tested fluidic sample is estimated by its flow resistance, which is a function of pressure drop, flow rate, and the geometry of the microfluidic channel. The viscometer can be exploited to measure viscosity of either Newtonian or non-Newtonian power-law fluid under various shear rates (3-500 1/s) and temperatures (4-70 °C) with small sample volume (less than 400 µL). The developed sensor-integrated microfluidic viscometer is made of poly(dimethylsiloxane) (PDMS) with transparent electrofluidic circuit, which makes it feasible to simultaneously image samples under tests. In addition, the entire device is disposable to prevent cross-contamination between samples, which is desired for various chemical and biomedical applications. In the experiments, viscosities of Newtonian fluids, glycerol water solutions with different concentrations and a mixture of pyrogallol and sodium hydroxide (NaOH), and non-Newtonian fluids, xanthan gum solutions and human blood samples, have been characterized. The results demonstrate that the developed microfluidic viscometer provides a convenient and useful platform for practical viscosity characterization of fluidic samples for a wide variety of applications.

Emerging Roles of Air Gases in Lipid Bilayers.

Recent studies indicate that changing the physical properties of lipid bilayers may profoundly change the function of membrane proteins. Here, the effects of dissolved nitrogen and oxygen molecules on the mechanical properties and stability of lipid bilayers are investigated using differential confocal microscopy, atomic force microscopy, and molecular dynamics simulations. All experiments evidence the presence of dissolved air gas in lipid bilayers prepared without gas control. The lipid bilayers in degassed solutions are softer and less stable than those in ambient solutions. High concentrations of nitrogen increase the bending moduli and stability of the lipid bilayers and impede phase separation in ternary lipid bilayers. The effect of oxygen is less prominent. Molecular dynamics simulations indicate that higher nitrogen affinity accounts for increased rigidity. These findings have fundamental and wide implications for phenomena related to lipid bilayers and cell membranes, including the origin of life.

Investigating Interfacial Effects on Surface Nanobubbles without Pinning Using Molecular Dynamics Simulation.

We investigated how the stability of aqueous argon surface nanobubbles on hydrophobic surfaces depends on gas adsorption, solid-gas interaction energy, and the bulk gas concentration using molecular dynamics simulation with the SPC/E water solvent. We observed stable surface nanobubbles without surface pinning sites for longer than 160 ns, contrary to previous findings using monoatomic Lennard-Jones solvent. In addition, the hydrophobicity of a substrate has an effect to reduce the requirement degree of oversaturation on water bulk. We found that the gas enrichment layer, gas adsorption monolayer on the hydrophobic substrate, and water hydrogen bonding near the interface are likely necessary conditions for nanobubble stability. We concluded that gas nanobubble stability does not necessarily require three-phase pinning sites.

Investigation of nematic to smectic phase transition and dynamical properties of strongly confined semiflexible polymers using Langevin dynamics.

We investigated the nematic to smectic phase transition for strongly confined semiflexible polymer solutions in slit-like confinements using GPU-accelerated Langevin dynamics. We characterized the phase transitions from the nematic to smectic phases for semi-flexible polymer solutions as the polymer density increased. The dependence for the lyotropic nematic to smectic transition can be collapsed by scaling exponents between 0.2 and 0.3. The smectic C phase is found for all the cases with the polymer orientation director tilted with respect to smectic layer lateral alignment. As the chain rigidity increases, the transition density decreases for systems in which the polymer persistence length (P) to slit height (H) ratios are 1.25, 2.5, 3.75, 5 and 25. We also characterized the polymer dynamics for the isotropic-nematic-smectic transitions. The overall polymer diffusivity decreased steadily as the polymer density increased. We observed anomalous polymer diffusion along the nematic director near the isotropic-nematic transition, similar to previously reported behavior for nematic-forming ellipsoids. Polymer diffusivity decreased sharply by two orders of magnitude upon the nematic-smectic transition.

Confinement, curvature, and attractive interaction effects on polymer surface adsorption.

We investigate the conformation and dynamics of a semi-flexible polymer near an attractive plane or a cylindrical post using Langevin dynamics. We characterize the transition from the desorbed to absorbed state and quantify how absorption depends on the attraction interaction, polymer molecular weight, polymer flexibility, intra-polymer interaction, and micro-confinement. We find that the critical point of adsorption for ideal flexible polymers only weakly depends on confinement. However, the critical point of adsorption increases significantly for self-avoiding flexible polymers and under confinement, deviating from scaling theory predictions. These findings provide insights into DNA surface adsorption in nanoslits and nanochannels.

Abnormal polymer transport in crowded attractive micropost arrays.

We investigate polymer diffusion in a quasi-two-dimensional environment decorated with attractive cylindrical posts using Langevin dynamics simulation. We find that the polymer diffusivity has non-monotonic dependence on the post array density. This diffusive behavior strongly depends on the adsorption-desorption transition and the critical adsorption strength εc. For ε < εc, the polymer undergoes normal diffusion and the diffusivity decreases as the post density increases due to the reduction of the void volume. For ε > εc, polymer dynamics is strongly mediated by post adsorption, and we observe a regime where the polymer diffusivity increases as the post density increases. The polymer diffusivity reaches a maximum, which can be attributed to cross-post translation enabled by large polymer conformation fluctuations. We find both cross-post transport and polymer conformation fluctuations strongly depend on the post absorption strength and the chain length.

Entropic attraction: Polymer compaction and expansion induced by nano-particles in confinement.

We investigated nanoparticle (NP)-induced coil-to-globule transition of a semi-flexible polymer in a confined suspension of ideal NP using Langevin dynamics. DNA molecules are often found to be highly compact, bound with oppositely charged proteins in a crowded environment within cells and viruses. Recent studies found that high concentration of electrostatically neutral NP also condenses DNA due to entropically induced depletion attraction between DNA segments. Langevin dynamics simulations with a semi-flexible chain under strong confinement were performed to investigate the competition between NP-induced monomer-monomer and monomer-wall attraction under different confinement heights and NP volume fractions. We found that whether NP induce polymer segments to adsorb to the walls and swell or to attract one another and compact strongly depends on the relative strength of the monomer-wall and the NP-wall interactions.

Entropy-driven single molecule tug-of-war of DNA at micro-nanofluidic interfaces.

Entropy-driven polymer dynamics at the nanoscale is fundamentally important in biological systems but the dependence of the entropic force on the nanoconfinement remains elusive. Here, we established an entropy-driven single molecule tug-of-war (TOW) at two micro-nanofluidic interfaces bridged by a nanoslit, performed the force analysis from a modified wormlike chain in the TOW scenario and the entropic recoiling process, and determined the associated scalings on the nanoconfinement. Our results provide a direct experimental evidence that the entropic forces in these two regimes, though unequal, are essentially constant at defined slit heights, irrespective of the slit lengths and the DNA segments within. Our findings have the implications to polymer transport at the nanoscale, device design for single molecule analysis, and biotechnological applications.

Flow-induced "waltzing" red blood cells: Microstructural reorganization and the corresponding rheological response.

We investigate flow-induced structural organization in a dilute suspension of tumbling red blood cells (RBCs) under confined shear flow. For small Reynolds (Re = 0.1) and capillary numbers (Ca), with fully coupled hydrodynamic interaction (HI) and without interparticle adhesion, we find that HI between the biconcave discoid particles prompts the formation of layered RBC chains and synchronized rotating RBC pairs, referred here as "waltzing doublets." As the volume fraction ϕ increases, more waltzing doublets appear in RBC files. Stronger shear stress disrupts structural arrangements at higher Ca. We find that the flow-induced organization of waltzing doublets changes how the suspension viscosity varies with ϕ qualitatively. The intrinsic viscosity is particularly sensitive to microstructural rearrangement, increasing (decreasing) with ϕ at low (high) Ca that correlates with the change in the fraction of doublets. We verified flow-induced collective motion with comparison to two-cell simulations in which the cell volume fraction is controlled by varying the domain volume.

Probing the Stoichiometry Dependence of Enzyme-Catalyzed Junction Zone Network Formation in Aiyu Pectin Gel via a Reaction Kinetics Model.

We investigate the enzymatic self-catalyzed gelation process in aiyu gel, a natural ion crosslinked polysaccharide gel. The gelation process depends on the concentration ratio (Rmax) of the crosslinking calcium ions and all galacturonic acid binding sites. The physical gel network formation relies on the assembly of calcium-polysaccharide crosslink bonds. The crosslinks are initially transient and through break-up/rebinding gradually re-organizing into long, stable junction zones. Our previous study formulated a reaction kinetics model to describe enzymatic activation, crosslinker binding, and crosslink microstructural reorganization, in order to model the complex growth of elasticity. In this study, we extend the theory for the time-dependent profile of complex moduli and examine the interplay of enzyme conversion, crosslink formation, and crosslink re-organization. The adjusted model captures how the gelation and structural rearrangement characteristic times vary with the polymer and calcium concentrations. Furthermore, we find that calcium ions act as both crosslinkers and dopants in the excess calcium ion scenario and the binding dynamics is determined by Rmax. This study provides perspectives on the dynamic binding behaviors of aiyu pectin gel system and the theoretical approach can be generalized to enzyme-catalyzed ionic gel systems.

Vibrational and electrochemical studies of pectin-a candidate towards environmental friendly lithium-ion battery development.

Pectin polymers are considered for lithium-ion battery electrodes. To understand the performance of pectin as an applied buffer layer, the electrical, magnetic, and optical properties of pectin films are investigated. This work describes a methodology for creating pectin films, including both pristine pectin and Fe-doped pectin, which are optically translucent, and explores their potential for lithium-ion battery application. The transmission response is found extended in optimally Fe-doped pectin, and prominent modes for cation bonding are identified. Fe doping enhances the conductivity observed in electrochemical impedance spectroscopy, and from the magnetic response of pectin evidence for Fe3+ is identified. The Li-ion half-cell prepared with pectin as binder for anode materials such as graphite shows stable charge capacity over long cycle life, and with slightly higher specific capacity compare with the cell prepared using polyvinylidene fluoride (PVDF) as binder. A novel enhanced charging specific capacity at a high C-rate is observed in cells with pectin binder, suggesting that within a certain rate (∼5 C), pectin has higher capacity at faster charge rates. The pectin system is found as a viable base material for organic-inorganic synthesis studies.

Conformation and diffusion behavior of ring polymers in solution: a comparison between molecular dynamics, multiparticle collision dynamics, and lattice Boltzmann simulations.

We have studied the effect of chain topology on the structural properties and diffusion of polymers in a dilute solution in a good solvent. Specifically, we have used three different simulation techniques to compare the chain size and diffusion coefficient of linear and ring polymers in solution. The polymer chain is modeled using a bead-spring representation. The solvent is modeled using three different techniques: molecular dynamics (MD) simulations with a particulate solvent in which hydrodynamic interactions are accounted through the intermolecular interactions, multiparticle collision dynamics (MPCD) with a point particle solvent which has stochastic interactions with the polymer, and the lattice Boltzmann method in which the polymer chains are coupled to the lattice fluid through friction. Our results show that the three methods give quantitatively similar results for the effect of chain topology on the conformation and diffusion behavior of the polymer chain in a good solvent. The ratio of diffusivities of ring and linear polymers is observed to be close to that predicted by perturbation calculations based on the Kirkwood hydrodynamic theory.

Migration and fractionation of deformable particles in microchannel.

The complexity of the coupling between soft particle deformation and fluid perturbation has limited studies of soft particle hydrodynamics to dilute suspensions. A hybrid Brownian dynamics-lattice Boltzmann method is presented that models nondilute soft spherical deformable particle (DP) suspensions in flow. Dependences on particle size and density are investigated for suspensions with over 100 DP. Multi-DP interactions lead to complex dependence of particle distributions on concentration and flow rate. Flow-induced DP migration toward channel center for DP in narrow channels is found. In wide channels, off-center peaks in the center of mass distribution for DP are found. The migration of DP leads to faster average speed of DP than the flow, which can be exploited for fractionating DPs of different sizes.

Elongation and migration of single DNA molecules in microchannels using oscillatory shear flows.

Much of modern biology relies on the strategic manipulation of molecules for creating ordered arrays prior to high throughput molecular analysis. Normally, DNA arrays involve deposition on surfaces, or confinement in nanochannels; however, we show that microfluidic devices can present stretched molecules within a controlled flow in ways complementing surface modalities, or extreme confinement conditions. Here we utilize pressure-driven oscillatory shear flows generated in microchannels as a new way of stretching DNA molecules for imaging "arrays" of individual DNA molecules. Fluid shear effects both stretch DNA molecules and cause them to migrate away from the walls becoming focused in the centerline of a channel. We show experimental findings confirming simulations using Brownian dynamics accounting for hydrodynamic interactions between molecules and channel-flow boundary conditions. Our findings characterize DNA elongation and migration phenomena as a function of molecular size, shear rate, oscillatory frequency with comparisons to computer simulation studies.

Shear-induced non-monotonic viscosity dependence for model red blood cell suspensions in microvessels.

The cell-free layer thickness of an aggregating red blood cell (RBC) suspension in a rectangular microchannel is investigated by hybrid fluid-particle numerical modeling. Several factors affect the suspension viscosity, cell-free layer thickness, and the cell aggregate distribution. These include the hematocrit, vessel size, red cell stiffness, aggregation interaction, and shear rate. In particular, the effect of the shear rate on the cell-free layer thickness is controversial. We found that the suspension viscosity increases along with a decrease in the cell-free layer thickness as the shear rate increases for aggregating model RBCs at low shear rates. At moderate to high shear rates, the cell-free layer thickness increases with the increasing shear rate from medium to strong shear flow for both 10% and 20% red blood cell suspensions.

Depletion-induced surface alignment of asymmetric diblock copolymer in selective solvents.

Phase separation of asymmetric diblock copolymer near surfaces in selective solvents is theoretically investigated by using the real-space version of self-consistent field theory (SCFT). Several morphologies are predicted and the phase diagram is constructed by varying the distance between two parallel hard surfaces (or the film thickness) W and the block copolymer concentration f(P). Morphologies of the diblock copolymer in dilute solution are found to change significantly with different film thicknesses. In confined systems, stable morphologies found in the bulk solution become unstable due to the loss of polymer conformation entropy. The vesicle phase region contracts when the repulsive interaction between the blocks is strong (strong segregation regime). The mixture of vesicles, rodlike and spherelike micelles and the mixture of vesicles and sphere-like micelles disappear in contrast to the weakly segregating regime. The walls strongly affect the phase separation of block copolymer in selective solvent, and the depletion layer near the surface contributes much to the micelle formation of the block copolymer. Interestingly, the self-assembled morphologies stay near the walls with the distance on the order of the radius of gyration of the block copolymer. The oscillation of the polymer distribution near the walls allows the surface phase separation to be observed due to the strong repulsion between the blocks A and B.