Mesoscopic simulation of the crossing dynamics at an entanglement point of surfactant threadlike micelles

Flow-induced scission of wormlike micelles in nonionic surfactant solutions under shear flow

We investigate flow-induced scission of wormlike micelles with dissipative particle dynamics simulations of nonionic surfactant solutions under shear flow. To understand flow-induced scission in terms of micellar timescales, we propose a method to evaluate the longest relaxation time of unentangled surfactant micelles from the rotational relaxation time and the average lifetime at equilibrium. The mean squared displacement of surfactant molecules provides evidence that the longest relaxation time estimated by the proposed method serves as the characteristic timescale at equilibrium. We also demonstrate that the longest relaxation time plays an essential role in flow-induced scission. Using conditional statistics based on the aggregation number of micelles, we examine the statistical properties of the lifetime of wormlike micelles. We then conclude that flow-induced scission occurs when the Weissenberg number defined as the product of the longest relaxation time and the shear rate is larger than a threshold value.

Rhamnolipid Biosurfactants for Oil Recovery: Salt Effects on the Structural Properties Investigated by Mesoscale Simulations

Rhamnolipids (RLs) are biosurfactants produced by Pseudomonas. The biodegradability and the variety of their functionality make them suitable for environmental remediation and oil recovery. We use dissipative particle dynamics simulations to investigate the aggregation behaviors of ionic RL congeners with nonane in various operating conditions. Under zero-salinity conditions, all RL congeners studied here form small ellipsoidal clusters with detectable free surfactants. When salt ions are present, the electrostatic repulsion between the ionized heads is overcome, resulting in the formation of larger aggregates of unique structures. RLs with C10-alkyl tails tend to form elongated wormlike micelles, while RLs with C16-alkyl tails tend to form clusters in spherical symmetry, including vesicles. Di-rhamnolipids (dRLs) require stronger solvation than monorhamnolipids (mRLs) to form clusters, and the resulting size of micelles is decreased. The morphology of the mixed dRL/mRL/oil systems is controlled based on the type of the congeners and the oil contents. In addition, the divalent calcium ions are found to be influential to the structure of the micelles through different mechanisms. For 5 wt % salinity, the ionic RLs can form oil-swollen micelles up to a 1:1 surfactant-to-oil ratio, suggesting that ionic RLs are superb to act as cleaning agents for petroleum hydrocarbons in the marine area. These key findings may guide the design for RL-based washing techniques in enhanced oil recovery.

Dissipative particle dynamics simulations in colloid and Interface science: a review

Dissipative particle dynamics (DPD) is one of the most efficient mesoscale coarse-grained methodologies for modeling soft matter systems. Here, we comprehensively review the progress in theoretical formulations, parametrization strategies, and applications of DPD over the last two decades. DPD bridges the gap between the microscopic atomistic and macroscopic continuum length and time scales. Numerous efforts have been performed to improve the computational efficiency and to develop advanced versions and modifications of the original DPD framework. The progress in the parametrization techniques that can reproduce the engineering properties of experimental systems attracted a lot of interest from the industrial community longing to use DPD to characterize, help design and optimize the practical products. While there are still areas for improvements, DPD has been efficiently applied to numerous colloidal and interfacial phenomena involving phase separations, self-assembly, and transport in polymeric, surfactant, nanoparticle, and biomolecules systems.

Molecular Insight into the Possible Mechanism of Drag Reduction of Surfactant Aqueous Solution in Pipe Flow

The phenomenon of drag reduction (known as the “Toms effect”) has many industrial and engineering applications, but a definitive molecular-level theory has not yet been constructed. This is due both to the multiscale nature of complex fluids and to the difficulty of directly observing self-assembled structures in nonequilibrium states. On the basis of a large-scale coarse-grained molecular simulation that we conducted, we propose a possible mechanism of turbulence suppression in surfactant aqueous solution. We demonstrate that maintaining sufficiently large micellar structures and a homogeneous radial distribution of surfactant molecules is necessary to obtain the drag-reduction effect. This is the first molecular-simulation evidence that a micellar structure is responsible for drag reduction in pipe flow, and should help in understanding the mechanisms underlying drag reduction by surfactant molecules under nonequilibrium conditions.

Effect of the Janus Amphiphilic Wall on the Viscosity Behavior of Aqueous Surfactant Solutions

The effects of the chemical nature of an interface are one of the key parameters which can affect self-assembly and rheological behavior. To date, several studies have reported self-assembled structures and rheological behaviors in the development of various functional materials. In this study, we investigated the self-assembly and viscosity behavior of aqueous surfactant solutions confined in three types of Janus amphiphilic nanotubes (JANTs), which have two, four, and eight sequential domains, respectively, using molecular simulation. We found that the viscosity behavior depends on the surfactant concentration and the chemical nature of the wall surface. For instance, although the concentration levels of the surfactants are the same (c = 10%), completely different viscosity behaviors were observed in the two sequential domains (Newtonian-like) and the four and eight sequential domains (strong shear-thinning) of the JANTs. Our simulations demonstrated how the rheological properties of aqueous surfactant solutions, including viscosity and velocity profiles, can be controlled by the chemical nature of the JANT wall surface, effect of confinement, and their self-assembly structures. Considering the foregoing, we hope that our study offers new knowledge on nanofluid systems.

Replica exchange dissipative particle dynamics method on threadlike micellar aqueous solutions

The self-assembly of surfactant molecules can spontaneously result in a variety of micelle morphologies, such as spherical micelles, threadlike micelles, and vesicles, and it is therefore crucial to predict and control the self-assembly to achieve a helpful process in the fields of materials chemistry and engineering. A dissipative particle dynamics (DPD) method used in a coarse-grained molecular simulation is applied to simulate various self-assembling soft matter systems because it can handle greater length and time scales than a typical molecular dynamics simulation (MD). It should be noted that the thorough sampling of a system is not assured at low temperatures because of large complex systems with coarse-grained representations. In this article, we demonstrate that the replica exchange method (REM) is very effective for even a DPD in which the energy barrier is comparatively lower than that of a MD. A replica exchange on DPD (REDPD) simulation for threadlike micellar aqueous solutions was conducted, and the values of the potential energy and the mean aggregation number were compared. As a result, the correct values and a self-assembled structure within a low-temperature range can only be obtained through the REDPD.

Multiscale prediction of functional self-assembled materials using machine learning: high-performance surfactant molecules

Various physical properties of functional materials can be induced by controlling their chemical molecular structures. Therefore, molecular design is crucial in the fields of engineering and materials science. With its remarkable development in various fields, machine learning combined with molecular simulation has recently been found to be effective at predicting the electronic structure of materials (Nat. Commun., 2017, 8, 872 and Nat. Commun., 2017, 8, 13890). However, previous studies have used similar microscale information as input and output data for machine learning, i.e., molecular structures and electronic structures. In this study, we determined whether multiscale data can be predicted using machine learning via a self-assembly functional material system. In particular, we investigated whether machine learning can be used to predict dispersion and viscosity, as the representative physical properties of a self-assembled surfactant solution, from the chemical molecular structures of a surfactant. The results showed that relatively accurate information on these physical properties can be predicted from the molecular structure, suggesting that machine learning can be used to predict multiscale systems, such as surfactant molecules, self-assembled micelle structures, and physical properties of solutions. The results of this study will aid in further development of the application of machine learning to materials science and molecular design.

Self-assembly of peptide amphiphiles by vapor pressure osmometry and dissipative particle dynamics

Peptide amphiphiles are one of the most promising materials in the biomedical field, so much effort has been devoted to characterizing the mechanism of their self-assembly and thermosensitive gelation. In this work, vapor pressure osmometry measurements were carried out to parameterize the thermosensitivity of interactions between peptide amphiphiles in an aqueous solution. The osmometry measurement verified that the peptides became more hydrophobic as temperature increased, which was quantitatively described with the Flory-Huggins χ parameter. Thereafter, a coarse-grained molecular model was used to simulate peptide amphiphiles dissolved in an aqueous solution. The temperature sensitive coarse-grained parameter a HW, which is the repulsive force between the hydrophilic head of the peptide amphiphile and water was estimated from the aforementioned experimentally obtained χ. Furthermore, the effects of concentration and temperature on the self-assembly behavior of peptide amphiphiles were quantitatively studied by dissipative particle dynamics. The simulation results revealed that a HW plays an important role in self-assembly characteristics and in the resulting microstructure of the peptide amphiphiles, which coincides with previous experimental and computational findings. The methodology in quantitatively linking the coarse-grained parameter from experiment and theory provides a sensible foundation for bridging future simulation studies with experimental work on macromolecules.

Self-Assembly of Janus Oligomers into Onion-like Vesicles with Layer-by-Layer Water Discharging Capability: A Minimalist Model

A vesicle in a cell is an enclosed structure in which the interior fluid is encompassed by a lipid bilayer. Synthetic vesicles are known as the liposomes. Liposomes with a single phospholipid bilayer are called unilamellar liposomes; otherwise, they are called multilamellar liposomes or onion-like liposomes (vesicles). One prototype synthetic onion-like vesicle, namely, onion-like dendrimersomes, have been recently produced via the self-assembly of amphiphilic Janus dendrimers (Proc. Natl. Acad. Sci. U.S.A. 2016, 113, 1162). Herein, we show computer simulation evidence of another type of onion-like vesicle, namely, onion-like oligomersomes, via the self-assembly of amphiphilic Janus oligomers in water. Specifically, we investigate the minimum-sized oligomers (or minimalist model) that can give rise to the onion-like oligomersomes as well as the composition-dependent phase diagrams. Insights into the formation condition and formation process of the onion-like oligomersomes are obtained. We demonstrate that the discharge of the in-vesicle water is through the remarkable "peeling-one-onion-layer-at-a-time" fashion, a feature that can be utilized for a clinical dosing regimen. The ability to control the formation of onion-like oligomersomes by design can be exploited for applications in drug and gene delivery.

Self-consistent description of electrokinetic phenomena in particle-based simulations

A new computational method is presented for study suspensions of charged particles undergoing fluctuating hydrodynamic and electrostatic interactions. The proposed model is appropriate for polymers, proteins, and porous particles embedded in a continuum electrolyte. A self-consistent Langevin description of the particles is adopted in which hydrodynamic and electrostatic interactions are included through a Green's function formalism. An Ewald-like split is adopted in order to satisfy arbitrary boundary conditions for the Stokeslet and Poisson Green functions, thereby providing a formalism that is applicable to any geometry and that can be extended to deformable objects. The convection-diffusion equation for the continuum ions is solved simultaneously considering Nernst-Planck diffusion. The method can be applied to systems at equilibrium and far from equilibrium. Its applicability is demonstrated in the context of electrokinetic motion, where it is shown that the ionic clouds associated with individual particles can be severely altered by the flow and concentration, leading to intriguing cooperative effects.

Turning up the heat on wormlike micelles with a hydrotopic salt in microfluidics

In equilibrium, wormlike micelles can transition from entangled to branched structures with increasing surfactant concentrations and ionic strength. Under flow conditions, structural transition of micellar solutions can follow very different trajectories. In this study we consider the flow of a semi-dilute wormlike micellar solution through an array of microposts, with focus on its rheological and microstructural evolutions. Specifically, the micellar solution (precursor) contains cationic surfactant cetyltrimethylammonium bromide (CTAB) and hydrotropic organic salt 3-hydroxynaphthalene-2-carboxylate (SHNC). We report the formation of a flow induced structured phase (FISP), with entangled, branched, and multi-connected micellar bundles, evidenced by electron microscopy and small-angle neutron scattering (SANS). By integrating gold-etched microheaters with the micropost design in a microfluidic device, we investigate the localized temperature effect on both the precursor and FISP, with complementary investigations from SANS. We observe that the FISP does not completely disintegrate at high temperatures, whereas, the precursor exhibits shortening of wormlike micelles as temperature increases. We also correlate the microstructure of both FISP and precursor with two point passive microrheology and bulk rheology characterizations.

Microstructure and rheology of a flow-induced structured phase in wormlike micellar solutions

Surfactant molecules can self-assemble into various morphologies under proper combinations of ionic strength, temperature, and flow conditions. At equilibrium, wormlike micelles can transition from entangled to branched and multiconnected structures with increasing salt concentration. Under certain flow conditions, micellar structural transitions follow different trajectories. In this work, we consider the flow of two semidilute wormlike micellar solutions through microposts, focusing on their microstructural and rheological evolutions. Both solutions contain cetyltrimethylammonium bromide and sodium salicylate. One is weakly viscoelastic and shear thickening, whereas the other is strongly viscoelastic and shear thinning. When subjected to strain rates of ∼10(3) s(-1) and strains of ∼10(3), we observe the formation of a stable flow-induced structured phase (FISP), with entangled, branched, and multiconnected micellar bundles, as evidenced by electron microscopy. The high stretching and flow alignment in the microposts enhance the flexibility and lower the bending modulus of the wormlike micelles. As flexible micelles flow through the microposts, it becomes energetically favorable to minimize the number of end caps while concurrently promoting the formation of cross-links. The presence of spatial confinement and extensional flow also enhances entropic fluctuations, lowering the energy barrier between states, thus increasing transition frequencies between states and enabling FISP formation. Whereas the rheological properties (zero-shear viscosity, plateau modulus, and stress relaxation time) of the shear-thickening precursor are smaller than those of the FISP, those of the shear-thinning precursor are several times larger than those of the FISP. This rheological property variation stems from differences in the structural evolution from the precursor to the FISP.

Tuning the viscoelastic properties of bis(urea)-based supramolecular polymer solutions by adding cosolutes

Polymers formed by the self-assembly of a bis(urea)-based polymer, 2,4-bis(2-ethylhexylureido)toluene (EHUT), in organic solvents such as octane are promising systems with remarkable rheological properties. This is the first self-assembled polymer recently reported as a hydrodynamic drag reducer for hydrocarbons. The rheology of diluted and semidiluted EHUT solutions can be tuned by specific interactions between the chains, modulated by the nature of the solvent and the presence of additives. In this article, rheological, thermal and SANS measurements were performed in order to investigate the competition between EHUT self-assembly and its interaction with specific molecules (benzene, benzyl alcohol, and ethanol) that can interact with EHUT unimers via hydrogen bonds and π-π interactions. No substantial rheological, thermal, or structural effect is observed when benzene is added to the systems. However, ethanol and benzyl alcohol interact with EHUT unimers through hydrogen bonds, drastically decreasing the viscoelasticity of the solutions. In addition, benzyl alcohol can interact with EHUT polymers by π-stacking interactions, playing an important role in tuning the rheological properties of the systems.

Nanochannel with uniform and Janus surfaces: shear thinning and thickening in surfactant solution

On basis of molecular simulation of confined surfactant solutions, we show that by adding chemical patterns on the inner surface of nanochannels dynamical properties of the confined surfactant solutions could be modified from shear thinning to shear thickening. To this end, we select uniformly hydrophobic and hydrophilic surfaces as well as a stripe-patterned Janus surface as three prototype confining surfaces of nanochannels. In all three nanochannels, when the surfactant solution is under relatively low shear rates, it shears thin. Under moderate shear rates, a sharp decrease in the shear viscosity could occur due to surfactant morphology transition. Under relatively high shear rates, a shear-thinning-to-thickening transition can emerge due to the tendency of stratification normal to the confining surface. Our simulation study offers a guide to steering dynamic properties of surfactant fluids in nanofluidic devices through engineering surfaces of nanochannels by design.

Asymmetric Brownian motor driven by bubble formation in a hydrophobic channel

The "asymmetric brownian ratchet model" is a variation of Feynman's ratchet and pawl system proposed. In this model, a system consisting of a motor and a rail has two binding states. One is the random brownian state, and the other is the asymmetric potential state. When the system is alternatively switched between these states, the motor can be driven in one direction. This model is believed to explain nanomotor behavior in biological systems. The feasibility of the model has been demonstrated using electrical and magnetic forces; however, switching of these forces is unlikely to be found in biological systems. In this paper, we propose an original mechanism of transition between states by bubble formation in a nanosized channel surrounded by hydrophobic atoms. This amounts to a nanoscale motor system using bubble propulsion. The motor system consists of a hydrophobic motor and a rail on which hydrophobic patterns are printed. Potential asymmetry can be produced by using a left-right asymmetric pattern shape. Hydrophobic interactions are believed to play an important role in the binding of biomolecules and molecular recognition. The bubble formation is controlled by changing the width of the channel by an atomic distance (∼0.1 nm). Therefore, the motor is potentially more efficient than systems controlled by other forces, in which a much larger change in the motor position is necessary. We have simulated the bubble-powered motor using dissipative particle dynamics and found behavior in good agreement with that of motor proteins. Energy efficiency is as high as 60%.

Wrapping nanotubes with micelles, hemimicelles, and cylindrical micelles

This work uses a simple model based on hydrophobic and hydrophilic forces to investigate the molecular dynamics that lead to the supramolecular self-assembly of surfactants around carbon nanotubes (CNTs). The effects of the concentration and the structure of surfactants are explored. The bead-based mesoscopic description spontaneously develops the several micellar morphologies that are known to wrap and solvate CNTs.

Multicompartment micelles and vesicles from pi-shaped ABC block copolymers: a dissipative particle dynamics study

Dissipative particle dynamics simulations were performed on the morphology and structure of multicompartment micelles and vesicles formed from pi-shaped ABC block copolymers in water. By varying the chain architecture and the composition of copolymers, a rich variety of morphologies were observed, such as oblate vesicles, trumpet vesicles, layered ribbon-like micelles and Y-shaped micelles. The simulations show that the hydrophilic block length and the distance between the two grafts play important roles in the control of the morphology. Since pi-shaped ABC block copolymers can reduce to linear ABC and star ABC block copolymers, they are good model copolymers for studying the self-assembly of complex block copolymers into micelles or vesicles. Thus, the knowledge obtained in this work as well as the new morphologies identified provides useful information for future rational design and synthesis of novel multicompartment micelles and vesicles.

Self-assembly of surfactants and polymorphic transition in nanotubes

We study self-assembly and polymorphic transitions of surfactant molecules in water within a nanotube and the effect of water-nanotube interactions on the self-assembly morphologies. We present a simulation evidence of a cornucopia of polymorphic structures of surfactant assemblies--many of which have not been observed in bulk solutions--through adjusting the water-nanotube chemical interactions which range from hydrophilic to hydroneutral and to hydrophobic. The ability to control the morphologies of surfactant assemblies within nanoscale confinement can be used for patterning the interior surface of nanochannels for application in nanofluidics and nanomedical devices.

Spontaneous self-assembly process for threadlike micelles

More than 100 micros dissipative particle dynamics simulations were carried out to investigate the spontaneous formation process of threadlike micelles from the random configuration for surfactant molecules. Stable spherical micelles were formed during the earlier stage. These spherical micelles fused to each other and grew into rodlike and threadlike micelles during the later stage. The length and radius of a micelle were estimated by tracing the backbone positions and the distance between the head group particles and the backbone of the micelles, respectively. The ratio of the largest to the smallest principal moments of inertia for each micelle was calculated as the micelle shape.

Budding Dynamics of Individual Domains in Multicomponent Membranes Simulated by N-Varied Dissipative Particle Dynamics

We study the budding dynamics of individual domains in flat, multicomponent membranes using dissipative particle dynamics (DPD) simulations with varied bead number N, in which addition and deletion of beads based on their density at the membrane boundary is introduced. The budding process of a tubular bud, accompanied by a dynamical transition reflected in the energy and morphology evolutions, is investigated. The simulations show that budding duration is shortened with increasing line tension and depends on the domain size quadratically. At low line tension, increasing bending modulus accelerates budding at first, but suppresses the process as it increases further. In addition, the controlling role of the surface tension in the budding process is also explored. Finally, we use the N-varied DPD to simulate the experimentally observed multicomponent tubular vesicles, and the three bud growth modes are confirmed.