Modeling Gas-Liquid Interfaces by Dissipative Particle Dynamics: Adsorption and Surface Tension of Cetyl Trimethyl Ammonium Bromide at the Air-Water Interface

Dissipative Particle Dynamics: Simulation of Chitosan-Citral Microcapsules

In this paper, the dissipative particle dynamics (DPD) method is used to simulate the self-assembly process, appearance, mesoscopic structure, and wrapping properties of microcapsules formed with citral as the core material and chitosan and sodium alginate as the single-wall materials, and with citral as the core material and chitosan-sodium alginate, chitosan-methylcellulose, sodium alginate-chitosan, and sodium alginate-methylcellulose as the double-wall materials. The effects of chitosan content and wall material composition on the structure, morphology, encapsulation performance, and stability of microcapsules are compared and analyzed. In addition, the microcapsules are deeply analyzed by using the mesoscopic structure, radial distribution function, and diffusion coefficient. This study provides a new idea and method for the preparation of citral microcapsules, and is of great significance for the design and development of new composite wall microcapsules.

Dissipative particle dynamics parametrisation using infinite dilution activity coefficients: the impact of bonding

Dissipative particle dynamics (DPD) simulations have proven to be a valuable coarse-grained simulation technique for studying complex systems such as surfactant and polymer solutions. However, the best method to use in parametrising DPD systems is not universally agreed. One common approach is to map infinite dilution activity coefficients to the DPD simulation 'beads' that represent molecular fragments. However, we show that here that this approach can lead to serious errors when bonding beads together to create molecules. We show errors arise from the verlaps between bonded beads, which alters their solubility. In this article, we demonstrate how these bonding errors can be accounted for when defining DPD force fields using simple theoretical methods to account for the overlapping volumes, and we demonstrate the validity of our approach by calculating the partition coefficients for a series of solutes into two immiscible solvents.

Effect mechanism of metal cations on the interface interaction of cell-collector-bubble for microalgal foam flotation

Foam flotation is generally recognized as a low-cost and efficient technology for the harvesting of microalgae for food, feed and fuel production, as well as environmental remediation. However, the harvesting efficiency of microalgae using foam flotation is restricted by the residual metal cations in the medium, and the corresponding inhibition mechanism has not yet been revealed. This study investigated the effects of metal cations in the medium on the harvesting efficiency and concentration factor during the foam flotation of Scenedesmus acuminatus. The interface interaction of cell-collector-bubble effected by metal cations was revealed by quantifying the amount of collector (cetyl trimethylammonium bromide, CTAB) between cells and bubbles, as well as the response of bubble interface characteristics. Results showed that the harvesting efficiency dropped linearly as the increase of cationic concentrations. Under the CTAB dose of 20 mg L-1, the harvesting efficiency decreased from 98.65% to 56.77% with a decrease of concentration factor from 25.41 to 9.05 in the presence of metal cations. The Na+ and Mg2+ in the medium were the major inhibitors. The inhibitory mechanisms revealed that metal cations obviously impeded the adsorption of CTAB onto the cells by competing adsorption site, resulting in a low harvesting efficiency. The presence of metal cations also inhibited the bubble coalescence and slowed down drainage velocity in the plateau channel of foam layer, forming foam with higher water content, thus reducing the concentration factor. A schematic illustration is proposed to better understand the effect mechanism of metal cations on microalgal foam flotation. This study might facilitate the process development in an effort to overcome the inhibition of cations during microalgal foam flotation.

A many-body dissipative particle dynamics parametrisation scheme to study behaviour at air-water interfaces

In this article, we present a general parametrisation scheme for many-body dissipative particle dynamics (MDPD). The scheme is based on matching model components to experimental surface tensions and chemical potentials. This allows us to obtain the correct surface and mixing behaviours of complex, multicomponent systems. The methodology is tested by modelling the behaviour of nonionic polyoxyethylene alkyl ether surfactants at an air/water interface. In particular, the influence of the number of ethylene oxide units in the surfactant head group is investigated. We find good agreement with many experimentally obtained parameters, such as minimum surface area per molecule; and a decrease in the surface tension with increasing surfactant surface density. Moreover, we observe an orientational transition, from surfactants lying directly on the water surface at low surface coverage, to surfactants lying parallel or tilted with respect to the surface normal at high surface coverage. The parametrisation scheme is also extended to cover the zwitterionic surfactant lauryldimethylamine oxide (LDAO), where we provide good predictions for the surface tension at maximum surface coverage. Here, if we exceed this coverage, we are able to demonstrate the spontaneous production of micelles from the surface surfactant layer.

Conformations and dynamic behaviors of confined wormlike chains in a pressure-driven flow

The conformations and dynamic behaviors of wormlike chains confined by a slit in a pressure-driven flow were investigated using dissipative particle dynamics method. The wormlike chains exhibit varying conformations due to the varying shear stresses across the slit. The wormlike chain solution can be well described by the power-law fluid, and the power-law index decreases with the increase in chain rigidity. We also presented that the wormlike chain undergoes tumbling motion in the vicinity of the wall in the presence of pressure-driven flow. We also found that the wormlike chains can migrate both away from the wall and slightly away from the slit center, and the migration away from the slit center increases as the chain rigidity is increased because of hydrodynamic interactions induced in a more rigid wormlike chain.

Coupling Effects of Ionic Surfactants and Electrolytes on the Stability of Bulk Nanobubbles

As interest in the extensive application of bulk nanobubbles increases, it is becoming progressively important to understand the key factors affecting their anomalous stability. The scientific intrigue over nanobubbles originates from the discrepancy between the Epstein-Plesset prediction and experimental observations. Herein, the coupling effects of ionic surfactants and electrolytes on the stability of bulk nanobubbles is studied. Experimental results show that ionic surfactants not only reduce the surface tension but also promote the accumulation of net charges, which facilitate the nucleation and stabilization of bulk nanobubbles. The addition of an electrolyte in a surfactant solution further results in a decrease in the zeta potential and the number concentration of nanobubbles due to the ion shielding effect, essentially colloidal stability. An adsorption model for the coexistence of ionic surfactants and electrolytes in solution, that specifically considers the effect of the adsorption layer thickness within the framework of the modified Poisson-Boltzmann equation, is developed. A quantitative agreement between the predicted and experimental surface tension is found in a wide range of bulk concentrations. The spatial distribution of the surface potential, surfactant ions and counterions in the vicinity of the interface of bulk nanobubbles are described. Our study intrinsically paves a route to investigate the stability of bulk nanobubbles.

Self-Assembly and Micellar Transition in CTAB Solutions Triggered by 1-Octanol

This study exploits higher-order micellar transition ranging from ellipsoidal to rodlike to wormlike induced by 1-octanol (C₈OH) in an aqueous solution of cetyltrimethylammonium bromide (CTAB), characterizing phase behavior, rheology, and small-angle neutron scattering (SANS). The phase diagram for the ternary system CTAB-C₈OH-water was constructed, which depicted the varied solution behavior. Such performance was further inferred from the rheology study (oscillatory-shear frequency sweep (ω) and viscosity (η)) that displayed an interesting solution behavior of CTAB solutions as a function of C₈OH. It was observed that at low C₈OH concentrations, the solutions appeared viscous/viscoelastic fluids that changed to an elastic gel with an infinite relaxation time at higher concentrations of C₈OH, thereby confirming the existence of distinct micelle morphologies. Small-angle neutron scattering (SANS) provided various micellar parameters such as aggregation numbers (N_agg) and micellar size/shape. The experimental results were further validated with a computational simulation approach. The molecular dynamic (MD) study offered an insight into the molecular interactions and aggregation behavior through different analyses, including radial distribution function (RDF), radius of gyration (R_g), and solvent-accessible surface area (SASA).

Mesoscale Modeling of Agglomeration of Molecular Bottlebrushes: Focus on Conformations and Clustering Criteria

Using dissipative particle dynamics, we characterize dynamics of aggregation of molecular bottlebrushes in solvents of various qualities by tracking the number of clusters, the size of the largest cluster, and an average aggregation number. We focus on a low volume fraction of bottlebrushes in a range of solvents and probe three different cutoff criteria to identify bottlebrushes belonging to the same cluster. We demonstrate that the cutoff criteria which depend on both the coordination number and the length of the side chain allows one to correlate the agglomeration status with the structural characteristics of bottlebrushes in solvents of various qualities. We characterize conformational changes of the bottlebrush within the agglomerates with respect to those of an isolated bottlebrush in the same solvents. The characterization of bottlebrush conformations within the agglomerates is an important step in understanding the relationship between the bottlebrush architecture and material properties. An analysis of three distinct cutoff criteria to identify bottlebrushes belonging to the same cluster introduces a framework to identify both short-lived transient and long-lived agglomerates; the same approach could be further extended to characterize agglomerates of various macromolecules with complex architectures beyond the specific bottlebrush architecture considered herein.

Osmotic Method for Calculating Surface Pressure of Monolayers in Molecular Dynamics Simulations

Surface pressure is a fundamental thermodynamic property related to the activity of molecules at interfaces. In molecular simulations, it is typically calculated from its definition: the difference between the surface tension of the air-water and air-surfactant interfaces. In this Letter, we show how to connect the surface pressure with a two-dimensional osmotic pressure and how to take advantage of this analogy to obtain a practical method of calculating surface pressure-area isotherms in molecular simulation. As a proof-of-concept, compression curves of zwitterionic and ionic surfactant monolayers were obtained using the osmotic approach and the curves were compared with the ones from the traditional pressure tensor-based scheme. The results shown an excellent agreement between both alternatives. Advantageously, the osmotic approach is simple to use and allows to obtain the surface pressure-area isotherm on the fly with a single simulation using equilibration stages.

DPD Simulation on the Transformation and Stability of O/W and W/O Microemulsions

The dissipative particle dynamics simulation method is adopted to investigate the microemulsion systems prepared with surfactant (H1T1), oil (O) and water (W), which are expressed by coarse-grained models. Two topologies of O/W and W/O microemulsions are simulated with various oil and water ratios. Inverse W/O microemulsion transform to O/W microemulsion by decreasing the ratio of oil-water from 3:1 to 1:3. The stability of O/W and W/O microemulsion is controlled by shear rate, inorganic salt and the temperature, and the corresponding results are analyzed by the translucent three-dimensional structure, the mean interfacial tension and end-to-end distance of H1T1. The results show that W/O microemulsion is more stable than O/W microemulsion to resist higher inorganic salt concentration, shear rate and temperature. This investigation provides a powerful tool to predict the structure and the stability of various microemulsion systems, which is of great importance to developing new multifunctional microemulsions for multiple applications.

Controlling Degradation and Erosion of Polymer Networks: Insights from Mesoscale Modeling

Understanding and controlling degradation of polymer networks on the mesoscale is critical for a range of applications. We utilize dissipative particle dynamics to capture photocontrolled degradation and erosion processes in hydrogels formed by end-linking of four-arm polyethylene glycol precursors. We demonstrate that the polydispersity and the fraction of broken-off fragments scale with the relative extent of reaction. The reverse gel point measured is close to the value predicted by the bond percolation theory on a diamond lattice. We characterize the erosion process via tracking the mass loss that accounts for the fragments remaining in contact with the percolated network. We quantify the dependence of the mass loss on the extent of reaction and on the properties of the film prior to degradation. These results elucidate the main features of degradation and erosion on the mesoscale and could provide guidelines for future design of degrading materials with dynamically controlled properties.

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.

Tuning Cationic Micelle Properties with an Antioxidant Additive: A Molecular Perspective

In this work, we characterize the micellization and morphology transition induced in aqueous cetyltrimethylammonium bromide (CTAB) solution by the addition of the antioxidant propyl gallate (PG) using tensiometry, rheology, and small-angle neutron scattering (SANS) techniques combined with the molecular dynamics (MD) simulation approach. The adsorption of CTAB at the air-water interface in the presence of varying [PG] revealed a progressive decrease in the critical micelle concentration (CMC), while the changes in different interfacial parameters indicated enhancement of the hydrophobicity induced by PG in the CTAB micellar system. The dynamic rheology behavior indicated an increase in the flow viscosity (η) as a function of [PG]. Moreover, the rheological components (storage modulus, G', and loss modulus, G″) depicted the viscoelastic features. SANS measurements depicted the existence of ellipsoidal micelles with varying sizes and aggregation number (Nagg) as a function of [PG] and temperature. Computational simulation performed using density functional theory (DFT) calculations and molecular dynamics (MD) provided an insight into the atomic composition of the examined system. The molecular electrostatic potential (MEP) analysis depicted a close proximity of CTAB, i.e., emphasized favorable interactions between the quaternary nitrogen of CTAB and the hydroxyl group of the PG monomer, further validated by the two-dimensional nuclear Overhauser enhancement spectroscopy (2D-NOESY), which showed the penetration of PG inside the CTAB micelles. In addition, various dynamic properties, viz., the radial distribution function (RDF), the radius of gyration (Rg), and solvent-accessible surface area (SASA), showed a significant microstructural evolution of the ellipsoidal micelles in the examined CTAB-PG system, where the changes in the micellar morphology with a more elongated hydrophobic chain and the increased Rg and SASA values indicated the notable intercalation of PG in the CTAB micelles.