Molecular simulations of droplet coalescence in oil/water/surfactant systems

Comprehensive review of the interfacial behavior of water/oil/surfactant systems using dissipative particle dynamics simulation

A comprehensive understanding of interfacial behavior in water/oil/surfactant systems is critical to evaluating the performance of emulsions in various industries, specifically in the oil and gas industry. To gain fundamental knowledge regarding this interfacial behavior, atomistic methods, e.g., molecular dynamics (MD) simulation, can be employed; however, MD simulation cannot handle phenomena that require more than a million atoms. The coarse-grained mesoscale methods were introduced to resolve this issue. One of the most effective mesoscale coarse-grained approaches for simulating colloidal systems is dissipative particle dynamics (DPD), which bridges the gap between macroscopic time and length scales and molecular-scale simulation. This work reviews the fundamentals of DPD simulation and its progress on colloids and interface systems, especially surfactant/water/oil mixtures. The effects of temperature, salt content, a water/oil ratio, a shear rate, and a type of surfactant on the interfacial behavior in water/oil/surfactant systems using DPD simulation are evaluated. In addition, the obtained results are also investigated through the lens of the chemistry of surfactants and emulsions. The outcome of this comprehensive review demonstrates the importance of DPD simulation in various processes with a focus on the colloidal and interfacial behavior of surfactants at water-oil interfaces.

Direct Observation of Emulsion Morphology, Dynamics, and Demulsification

Herein, we present the direct observation and quantification of a water-in-oil (w/o) emulsion, its destabilization, and the effect of additives on such processes at the nanoscale. This is achieved via liquid phase transmission electron microscopy (LPTEM), wherein a small volume of emulsion is encapsulated against vacuum in its liquid state to allow observation of its initial morphology and its evolution over time at excellent spatial and temporal resolution. Emulsions of this class are useful for delivering payloads of materials insoluble in their delivery medium and are currently widely used across food science, pharmaceuticals, and environmental applications. However, their utility is inherently limited by their thermodynamic tendency to demulsify, eventually leading to bulk phase separation. This occurs via several degradation mechanisms, operating at times collectively, and which are difficult to differentiate via traditional ensemble methods (e.g., light scattering), obscuring mechanistic nuances. LPTEM as a characterization technique has the potential to augment our understanding of emulsion behavior and improve performance and formulations. In this work, we also emphasize the importance of the included videographic Supporting Information data in demonstrating the behavior of the studied materials.

Nonane and Hexanol Adsorption in the Lamellar Phase of a Nonionic Surfactant: Molecular Simulations and Comparison to Ideal Adsorbed Solution Theory

Adsorption of n-nonane/1-hexanol (C9/C6OH) mixtures into the lamellar phase formed by a 50/50 w/w triethylene glycol mono-n-decyl ether (C10E3)/water system was studied using configurational-bias Monte Carlo simulations in the osmotic Gibbs ensemble. The interactions were described by the Shinoda-Devane-Klein coarse-grained force field. Prior simulations probing single-component adsorption indicated that C9 molecules preferentially load near the center of the bilayer, increasing the bilayer thickness, whereas C6OH molecules are more likely to be found near the interface of the polar and nonpolar moieties, swelling the bilayer in the lateral dimension. Here, we extend this work to binary C9/C6OH adsorption to probe whether the difference in the spatial preferences may lead to a synergistic effect and enhanced loadings for the mixture. Comparing loading trends and the thermodynamics of binary adsorption to unary adsorption reveals that C9-C9 interactions lead to the largest enhancement, whereas C9-C6OH and C6OH-C6OH interactions are less favorable for this bilayer system. Ideal adsorbed solution theory yields satisfactory predictions of the binary loading.

Responsive Liquid Metal Droplets: From Bulk to Nano

Droplets exist widely in nature and play an extremely important role in a broad variety of industrial processes. Typical droplets, including water and oil droplets, have received extensive attention and research, however their single properties still cannot meet diverse needs. Fortunately, liquid metal droplets emerging in recent years possess outstanding properties, including large surface tension, excellent electrical and thermal conductivity, convenient chemical processing, easy transition between liquid and solid phase state, and large-scale deformability, etc. More interestingly, liquid metal droplets with unique features can respond to external factors, including the electronic field, magnetic field, acoustic field, chemical field, temperature, and light, exhibiting extraordinary intelligent response characteristics. Their development over the past decade has brought substantial breakthroughs and progress. To better promote the advancement of this field, the present article is devoted to systematically summarizing and analyzing the recent fundamental progress of responsive liquid metal droplets, not only involving droplet characteristics and preparation methods, but also focusing on their diverse response behaviors and mechanisms. On this basis, the challenges and prospects related to the following development of liquid metal droplets are also proposed. In the future, responsive liquid metal droplets with a rapid development trend are expected to play a key role in soft robots, biomedicine, smart matter, and a variety of other fields.

Viscoelastic necking dynamics between attractive microgels

HYPOTHESIS

Microgels can deform and interpenetrate and display colloid/polymer duality. The effective interaction of microgels in the collapsed state is governed by the interplay of polymer-solvent interfacial tension and bulk elasticity. A connecting neck is shown to mediate microgel interaction, but its temporal evolution has not been addressed. We hypothesize that the necking dynamics of attractive microgels exhibits liquid-like or solid-like behavior over different time and length scales.

EXPERIMENTS

We simulate the merging and pinching of attractive microgels with different crosslinking densities in explicit solvent using dissipative particle dynamics. The temporal coalescence dynamics of microgels is investigated and compared with simple liquid and polymeric droplets. We model the neck growth on long time scales using Maxwell model of polymer relaxation and compare the theoretical prediction with simulation data. The mechanical strength of the neck is characterized systematically via simulated pinch-off of microgels by steered molecular dynamics.

FINDINGS

We evidence a crossover in the coalescence dynamics reflecting the viscoelastic signature of microgels. In contrast to the common knowledge that viscoelastic materials respond elastically on short time scales, the early expansion of the microgel neck exhibits a linear behavior, similar to the viscous coalescence of liquid droplets. However, the late regime with arrested dynamics resembles sintering of solid particles. Through an analytical model relating microgel dynamics to neck growth, we show that the long-term behavior is governed by stress relaxation of the polymers in the neck region and predict an exponential decay in the rate of growth, which agrees favorably with the simulation. Different from coalescence, the thread thinning in microgel breakup primarily highlights its polymeric characteristics.

Effect of a Surfactant Mixture on Coalescence Occurring in Concentrated Emulsions: The Hole Nucleation Theory Revisited

By conducting both a bottle test and isolate drop-drop experiments, we determine the coalescence rates of water droplets within water-in-oil emulsions stabilized by a large amount of Span 80 in the presence of Tween 20, a surfactant that acts as a demulsifier. Using a microscopic model based on a theory of hole nucleation, we establish an analytical formula that quantitatively predicts the coalescence frequency per unit area of droplets whose interfaces are fully covered by surfactant molecules. Despite its simplicity and the strong assumptions made for its derivation, this formula captures our experimental findings on Span 80-stabilized emulsions as well as other results, found in the literature, remarkably well on a wide range of water-in-crude oil systems.

Exploring the effects of approach velocity on depletion force and coalescence in oil-in-water emulsions

An emulsion is a thermodynamically unstable system consisting of at least two immiscible liquid phases, one of which is dispersed in the other in the form of droplets of varying size. Most studies on emulsions have focused on the behaviour of emulsion droplets with diameter from ∼50 μm and upwards. However, the properties of smaller droplets may be highly relevant in order to understand the behaviour of emulsions, including their performance in numerous applications within the fields of food, industry, and medical science. The relatively long life-time and small size of these droplets compared to other emulsion droplets, make them suited for optical trapping and micromanipulation technologies. Optical tweezers have previously shown potential in the study of stabilized emulsions. Here we employ optical tweezers to examine unstable oil-in-water emulsions to determine the effects of system parameters on depletion force and coalescence times.

Coalescence in concentrated emulsions: theoretical predictions and comparison with experimental bottle test behaviour

Fusion between emulsion drops, also called coalescence, may be undesirable for storage or sought after depending on the desired application. In this latter case, a complete separation of the two liquids composing the emulsion is required. The same objective may be applicable to foams. We have performed bottle test experiments on a model system of water in oil (w/o) emulsion stabilized by high amounts of hydrophobic surfactant Span 80. We observe two regimes for emulsion separation: the first regime, which is fast and includes sedimentation of the water droplets, and the second regime, which exhibits a very dense and stable emulsion zone. We predict the initial thickness of the dense zone as a simple function of surfactant concentration and mean droplet size. From the assumption that the coalescence rate depends only on the area of the thin film between two contacted droplets, we quantitatively model the separation kinetics of the dense emulsion zone. Our results give rise to a simple method that allows measuring the coalescence frequency per unit area, only by monitoring bottle test experiments.

Temporally Coarse-Grained All-Atom Molecular Dynamics Achieved via Stochastic Padé Approximants

A Padé approximant scheme for realizing the discrete-time evolution of the state of a many-atom system is introduced. This temporal coarse-graining scheme accounts for the underlying Newtonian physics and avoids the need for construction of spatially coarse-grained variables. Newtonian physics is incorporated through short molecular dynamics simulations at the beginning of each of the large coarse-grained timesteps. The balance between stochastic and coherent dynamics expressed by many-atom systems is captured via incorporation of the Ito formula into a Padé approximant for the time dependence of individual atom positions over large timesteps. Since the time for a many-atom system to express a characteristic ensemble of atomic velocity fluctuations is typically short relative to the characteristic time of large-scale atomic displacements, a computationally efficient and accurate temporal coarse-graining of the atom-resolved Newtonian dynamics is formulated, denoted all-atom Padé-Ito molecular dynamics (APIMD). Evolution of the system over a time step much longer than that required for standard molecular dynamics (MD) is achieved via incorporation of information from the short MD simulations into a Padé approximant extrapolation in time. The extrapolated atomic configuration is subjected to energy minimization and, when needed, thermal equilibration so as to avoid occasional unphysical close encounters deriving from the Padé approximant extrapolation and to represent configurations appropriate for the temperature of interest. APIMD is implemented and tested via comparison with traditional MD simulations of five phenomena: (1) pertussis toxin subunit deformation, (2) structural transition in a T = 1 capsid-like structure of HPV16 L1 protein, (3) coalescence of argon nanodroplets, and structural transitions in dialanine in (4) vacuum, and (5) water. Accuracy of APIMD is demonstrated using semimicroscopic descriptors (rmsd, radius of gyration, residue-residue contact maps, and densities) and the free energy. Significant computational acceleration relative to traditional molecular dynamics is illustrated.

Quantification of Ostwald Ripening in Emulsions via Coarse-Grained Simulations

Ostwald ripening is a diffusional mass transfer process that occurs in polydisperse emulsions, often with the result of threatening the emulsion stability. In this work, we design a simulation protocol that is capable of quantifying the process of Ostwald ripening at the molecular level. To achieve experimentally relevant time scales, the dissipative particle dynamics (DPD) simulation protocol is implemented. The simulation parameters are tuned to represent two benzene droplets dispersed in water. The coalescence between the two droplets is prevented via the introduction of membranes, which allow diffusion of benzene from one droplet to the other. The simulation results are quantified in terms of the changes in the droplet volume as a function of time. The results are in qualitative agreement with experiments. The agreement with the Lifshitz-Slyozov-Wagner theory becomes quantitative when the simulated solubility and diffusion coefficient of benzene-in-water are considered. The effect of two different surfactants was also investigated. In agreement with both experimental observations and theory, the addition of surfactants at moderate concentrations decreased the Ostwald ripening rate because of the reduction in the interfacial tension between benzene and water; as the surfactant film becomes dense, other phenomena are likely to further delay the Ostwald ripening. In fact, the results suggest that the surfactant that yields higher density at the benzene-water interface delayed more effectively Ostwald ripening. The formation of micelles can also affect the ripening rate, in qualitative agreement with experiments, although our simulations are not conclusive on such effects. Our simulations show that the coarse-grained DPD formalism is able to capture the molecular phenomena related to Ostwald ripening and reveal molecular level features that could help to understand experimental observations. The results could be useful for predicting and eventually controlling the long-term stability of emulsions.

Molecular Dynamics Simulations of the Electrocoalescence Behaviors of Two Unequally Sized Conducting Droplets

The electrocoalescence of droplets plays a crucial role in various fields. However, studies on the effects of droplet radius on the electrocoalescence behaviors of droplets have not been conducted until now. In this work, the electrocoalescence behaviors of two unequally sized conducting droplets are investigated via molecular dynamics (MD) simulations. The influences of electric field strength and droplet radius on the electrocoalescence behaviors of two unequally sized droplets are investigated. When the electric field strength increases, the contact cone angle between the droplets increases, and the two droplets are more likely to partially coalesce and bounce. When the radius of the smaller droplet between the two droplets increases at the same electric field strength, the contact cone angle, daughter droplet size, and ions in the daughter droplet increase, whereas the critical electric field strength ( E_d) for generating the daughter droplet decreases. Furthermore, the daughter droplet is ejected from the smaller droplet when the two droplets have different radii.

Probing Additive Loading in the Lamellar Phase of a Nonionic Surfactant: Gibbs Ensemble Monte Carlo Simulations Using the SDK Force Field

Understanding solute uptake into soft microstructured materials, such as bilayers and worm-like and spherical micelles, is of interest in the pharmaceutical, agricultural, and personal care industries. To obtain molecular-level insight on the effects of solutes loading into a lamellar phase, we utilize the Shinoda-Devane-Klein (SDK) coarse-grained force field in conjunction with configurational-bias Monte Carlo simulations in the osmotic Gibbs ensemble. The lamellar phase is comprised of a bilayer formed by triethylene glycol mono- n-decyl ether (C10E3) surfactants surrounded by water with a 50:50 surfactant/water weight ratio. We study both the unary adsorption isotherm and the effects on bilayer structure and stability caused by n-nonane, 1-hexanol, and ethyl butyrate at several different reduced reservoir pressures. The nonpolar n-nonane molecules load near the center of the bilayer. In contrast, the polar 1-hexanol and ethyl butyrate molecules both load with their polar bead close to the surfactant head groups. Near the center of the bilayer, none of the solute molecules exhibits a significant orientational preference. Solute molecules adsorbed near the polar groups of the surfactant chains show a preference for orientations perpendicular to the interface, and this alignment with the long axis of the surfactant molecules is most pronounced for 1-hexanol. Loading of n-nonane leads to an increase of the bilayer thickness, but does not affect the surface area per surfactant. Loading of polar additives leads to both lateral and transverse swelling. The reduced Henry's law constants of adsorption (expressed as a molar ratio of additive to surfactant per reduced pressure) are 0.23, 1.4, and 14 for n-nonane, 1-hexanol, and ethyl butyrate, respectively, and it appears that the SDK force field significantly overestimates the ethyl butyrate-surfactant interactions.

Drop coalescence in technical liquid/liquid applications: a review on experimental techniques and modeling approaches

The coalescence phenomenon of drops in liquid/liquid systems is reviewed with particular focus on its technical relevance and application. Due to the complexity of coalescence, a comprehensive survey of the coalescence process and the numerous influencing factors is given. Subsequently, available experimental techniques with different levels of detail are summarized and compared. These techniques can be divided in simple settling tests for qualitative coalescence behavior investigations and gravity settler design, single-drop coalescence studies at flat interfaces as well as between droplets, and detailed film drainage analysis. To model the coalescence rate in liquid/liquid systems on a technical scale, the generic population balance framework is introduced. Additionally, different coalescence modeling approaches are reviewed with ascending level of detail from empirical correlations to comprehensive film drainage models and detailed computational fluid and particle dynamics.

Computational investigation of structure, dynamics and nucleation kinetics of a family of modified Stillinger-Weber model fluids in bulk and free-standing thin films

In recent years, computer simulations have found increasingly widespread use as powerful tools for studying phase transitions in wide variety of systems. In the particular and very important case of aqueous systems, the commonly used force-fields tend to offer quite different predictions with respect to a wide range of thermodynamic and kinetic properties, including the ease of ice nucleation, the propensity to freeze at a vapor-liquid interface, and the existence of a liquid-liquid phase transition. It is thus of fundamental and practical interest to understand how different features of a given water model affect its thermodynamic and kinetic properties. In this work, we use the forward-flux sampling technique to study the crystallization kinetics of a family of modified Stillinger-Weber (SW) potentials with energy (ε) and length (σ) scales taken from the monoatomic water (mW) model, but with different tetrahedrality parameters (λ). By increasing λ from 21 to 24, we observe the nucleation rate increases by 48 orders of magnitude at a supercooling of ζ = T/Tm = 0.845. Using classical nucleation theory, we are able to demonstrate that this change can largely be accounted for by the increase in |Δμ|, the thermodynamic driving force. We also perform rate calculations in freestanding thin films of the supercooled liquid, and observe a crossover from surface-enhanced crystallization at λ = 21 to bulk-dominated crystallization for λ ≥ 22.

Molecular Dynamics Simulations on Coalescence and Non-coalescence of Conducting Droplets

When an electric field with various strengths is applied to two adjacent conducting droplets, the droplets may completely coalesce, partially coalesce, or bounce off one another. To reveal an atom-scale mechanism of coalescence or non-coalescence, dynamic behaviors of two conducting nanodroplets at a homogeneous electric field are studied via molecular dynamics simulations in this work. The results show that there is a critical field strength and a critical cone angle above which the two droplets partially coalesce or bounce off. Charge transfer between the two droplets is observed when the droplets are brought into contact. The partial coalescence and the bounce-off of the two droplets at strong field strengths are found to be due to the high charge transfer rate, which leads to the breakup of the coalescing droplet at different locations.

Investigation of demulsification efficiency in water-in-crude oil emulsions using dissipative particle dynamics

Demulsification efficiency with alternating hydrophobic blocks of the polyether is investigated by dissipative particle dynamics.

Mesoscale models of dispersions stabilized by surfactants and colloids

In this paper we discuss and give an outlook on numerical models describing dispersions, stabilized by surfactants and colloidal particles. Examples of these dispersions are foams and emulsions. In particular, we focus on the potential of the diffuse interface models based on a free energy approach, which describe dispersions with the surface-active agent soluble in one of the bulk phases. The free energy approach renders thermodynamic consistent models with realistic sorption isotherms and adsorption kinetics. The free energy approach is attractive because of its ability to describe highly complex dispersions, such as emulsions stabilized by ionic surfactants, or surfactant mixtures and dispersions with surfactant micelles. We have classified existing numerical methods into classes, using either a Eulerian or a Lagrangian representation for fluid and for the surfactant/colloid. A Eulerian representation gives a more coarse-grained, mean field description of the surface-active agent, while a Lagrangian representation can deal with steric effects and larger complexity concerning geometry and (amphiphilic) wetting properties of colloids and surfactants. However, the similarity between the description of wetting properties of both Eulerian and Lagrangian models allows for the development of hybrid Eulerian/Lagrangian models having advantages of both representations.

Cluster formation of anchored proteins induced by membrane-mediated interaction

Computer simulations were used to study the cluster formation of anchored proteins in a membrane. The rate and extent of clustering was found to be dependent upon the hydrophobic length of the anchored proteins embedded in the membrane. The cluster formation mechanism of anchored proteins in our work was ascribed to the different local perturbations on the upper and lower monolayers of the membrane and the intermonolayer coupling. Simulation results demonstrated that only when the penetration depth of anchored proteins was larger than half the membrane thickness, could the structure of the lower monolayer be significantly deformed. Additionally, studies on the local structures of membranes indicated weak perturbation of bilayer thickness for a shallowly inserted protein, while there was significant perturbation for a more deeply inserted protein. The origin of membrane-mediated protein-protein interaction is therefore due to the local perturbation of the membrane thickness, and the entropy loss-both of which are caused by the conformation restriction on the lipid chains and the enhanced intermonolayer coupling for a deeply inserted protein. Finally, in this study we addressed the difference of cluster formation mechanisms between anchored proteins and transmembrane proteins.

Coarse-grained model for mechanosensitive ion channels

In this work, a coarse-grained model for a kind of membrane protein, the mechanosensitive channel of small conductance (MscS), is proposed. The basic structure of the MscS is preserved when the protein is coarse-grained. For the coarse-grained model, the channels show two different states, namely the open and closed states, depending on the model parameters. Under the same membrane tension, the state of the ion channel is found to be critically determined by the protein structure, especially the length of the transmembrane alpha-helix. It is also found that for the protein with certain size, the gating transition occurs when the membrane tension is applied, resembling in a real mechanosensitive channel.

Forward flux sampling for rare event simulations

Rare events are ubiquitous in many different fields, yet they are notoriously difficult to simulate because few, if any, events are observed in a conventional simulation run. Over the past several decades, specialized simulation methods have been developed to overcome this problem. We review one recently developed class of such methods, known as forward flux sampling. Forward flux sampling uses a series of interfaces between the initial and final states to calculate rate constants and generate transition paths for rare events in equilibrium or nonequilibrium systems with stochastic dynamics. This review draws together a number of recent advances, summarizes several applications of the method and highlights challenges that remain to be overcome.

Stability and rupture of archaebacterial cell membrane: a model study

It is known that the thermoacidophilic archaebacterium Sulfolobus acidocaldarius can grow in hot springs at 65-80 degrees C and live in acidic environments (pH 2-3); however, the origin of its unusual thermal stability remains unclear. In this work, using a vesicle as a model, we study the thermal stability and rupture of archaebacterial cell membrane. We perform a simulation investigation of the structure-property relationship of monolayer membrane formed by bolaform lipids and compare it with that of bilayer membrane formed by monopolar lipids. The origin of the unusually thermal stability of archaebacterial cell and the mechanism for its rupture are presented in molecular details.

Dynamics of coalescence of plugs with a hydrophilic wetting layer induced by flow in a microfluidic chemistrode

This manuscript analyzes the dynamics of coalescence of an incoming aqueous plug with a wetting layer above a hydrophilic surface in the chemistrode. The chemistrode is a recently described (Chen, D.; Du, W.; Liu, Y.; Liu, W.; Kuznetsov, A.; Mendez, F. E.; Philipson, L. H.; Ismagilov, R. F. Proc. Natl. Acad. Sci. U.S.A. 2008, 105, 16843-16848) microfluidic analogue of an electrode, but operating at the chemical rather than electrical level, developed with the aim of capturing local stimulus-response processes in chemistry and biology. The chemistrode consists of open-ended V-shaped microfluidic channels that can be brought into contact with a chemical or biological hydrophilic substrate. The chemistrode relies on multiphase aqueous/fluorous flow and uses plugs to achieve high temporal resolution of stimulation and sampling. Coalescence of the incoming plugs, containing the stimuli, with the liquid in the wetting layer is required for chemical exchange to take place in the chemistrode. Here, we investigate the system with triethyleneglycol mono[1H,1H-perfluorooctyl]ether RfOEG as the surfactant. This surfactant was necessary to prevent nonspecific absorption of proteins to the aqueous fluorous interface and to ensure biocompatibility of the system, but too much surfactant increased the barrier for coalescence. In this system, coalescence was controlled by the capillary number. At a higher value of the capillary number, coalescence took more time, and deformation of the interface of the incoming plug and the wetting layer was more significant. Above a critical capillary number, coalescence did not occur between the incoming plug and the wetting layer. The critical capillary number was an increasing function of surface tension but was independent of viscosity ratio. Coalescence was surprisingly reproducible, presumably because film rupture during coalescence was reliably initiated at the hydrophilic substrate. These results are useful in rational operation of the chemistrode and also provide an experimental description of deformation, film drainage, and coalescence of surfactant-coated droplets in an external flow field.

The chemistrode: a droplet-based microfluidic device for stimulation and recording with high temporal, spatial, and chemical resolution

Microelectrodes enable localized electrical stimulation and recording, and they have revolutionized our understanding of the spatiotemporal dynamics of systems that generate or respond to electrical signals. However, such comprehensive understanding of systems that rely on molecular signals-e.g., chemical communication in multicellular neural, developmental, or immune systems-remains elusive because of the inability to deliver, capture, and interpret complex chemical information. To overcome this challenge, we developed the "chemistrode," a plug-based microfluidic device that enables stimulation, recording, and analysis of molecular signals with high spatial and temporal resolution. Stimulation with and recording of pulses as short as 50 ms was demonstrated. A pair of chemistrodes fabricated by multilayer soft lithography recorded independent signals from 2 locations separated by 15 mum. Like an electrode, the chemistrode does not need to be built into an experimental system-it is simply brought into contact with a chemical or biological substrate, and, instead of electrical signals, molecular signals are exchanged. Recorded molecular signals can be injected with additional reagents and analyzed off-line by multiple, independent techniques in parallel (e.g., fluorescence correlation spectroscopy, MALDI-MS, and fluorescence microscopy). When recombined, these analyses provide a time-resolved chemical record of a system's response to stimulation. Insulin secretion from a single murine islet of Langerhans was measured at a frequency of 0.67 Hz by using the chemistrode. This article characterizes and tests the physical principles that govern the operation of the chemistrode to enable its application to probing local dynamics of chemically responsive matter in chemistry and biology.

Computer simulations of solute exchange using micelles by a collision-driven fusion process

In this work, the kinetic process of collision-driven solute exchange in an aqueous phase in which micelles are used as solute carriers is investigated by dissipative particle dynamics simulations. Here, we try to answer two questions about the exchange process of hydrophobic solute molecules: How the solute molecules are exchanged and what factors affect the process. For the first question, the simulation results indicate that, after a stage of intermittent collision between two neighboring aggregates, there are roughly three sequential events in a coalescence stage: (1) molecular contact, (2) neck formation, and (3) neck growth. The coalescence stage is followed by a stage of solute transfer and diffusion. It is found that there are two rate-limiting steps in the whole process of solute exchange, i.e., the break of the water film between two neighboring aggregates and the nucleation of a pore between two surfactant films. For the second question, the effects of the collision velocity, the surface tension, the repulsive interaction between the surfactant films of the colliding aggregates, as well as the steric repulsion are examined. For example, the simulation results show that the depletion force plays an important role during the coalescence stage, while the initial collision velocity basically does not change the fusion ratio. The results also demonstrate that the surface tension and interaction show different effects on the different stages of a solute exchange process.