Experimental Study and Molecular Simulation of the Effect of Temperature on the Stability of Surfactant Foam

Temperature changes in CO2 foam-fracturing construction can easily affect surfactant foam stability. To investigate the effect of temperature on the foam stability of different types of surfactants, this study measured the foam half-life and viscosity of four typical surfactants, CTAB, LAS-30, HSB1214, and TX-10, using a novel self-designed and built foam performance measurement device. The effects of temperature on foam half-life and viscosity were studied. The results show that as the temperature increased, the half-life shortened, and the viscosity of the liquid phase decreased, which led to a decrease in foam stability. Moreover, using Materials Studio, a type of molecular simulation software, an interfacial model of the foam film was constructed to calculate the IFE and the self-diffusion coefficient of water molecules at 300 ps after the equilibrium of the foam system to investigate the mechanism of temperature influence on the stability of the foam. The results show that, for CTAB, LAS-30, HSB1214, and TX-10, the temperature increases from 15 °C to 45 °C, the IFE is enhanced by −50.05%, −59.10%, −64.21%, and −44.26%, respectively, the interfacial system changes from a low-energy state to a high-energy state, and the interfacial stability decreases. Meanwhile, Dwater increased 1.10-fold, 0.78-fold, 1.43-fold, and 0.64-fold, respectively, which accelerated the diffusion and migration of water molecules, weakened the intermolecular forces, and accelerated the instability of the foam system.

Effect of Wetting and Dewetting Dynamics on Atomic Force Microscopy Measurements

Water bridge dynamics between an atomic force microscopy (AFM) tip and a flat substrate is studied by using a multibody dissipative particle dynamics (MDPD) model. First, the numerical model is validated by comparing the present results of droplet contact angles and liquid bridges with those reported in the literature. Then, the ability of MDPD to capture the meniscus shape and behavior for different operating conditions and geometric parameters is examined for both static and dynamic cases. Hence, several parametric studies and analyses of the AFM tip configuration and its operating conditions are reported. It is found that a critical capillary number of about 0.001 is calculated based on 5% change on the force measurements between the static and dynamic results. It is also demonstrated that the hysteresis behavior in the capillary force exerted on the AFM tip can be successfully predicted by using the MDPD model when the tip approaches or retracts from the substrate. Moreover, there is an excellent agreement in the results of breakup distance for different water bridge volumes between the predictions of the MDPD model and the theory. Also, the hysteresis of capillary force exerted on an AFM tip composed of multibody design is studied. The prediction on the transition of the capillary force vs distance between the AFM tip and the substrate is in good agreement with the experimental results. Therefore, we demonstrate a validated MDPD model which can successfully capture liquid bridge dynamics. This model can be used as a powerful design tool for meniscus manipulation technology, such as dip-pen nanolithography, as well as for studying dynamic, e.g., tapping mode AFM tip, interactions with a liquid bridge.

Atomically detailed simulations of concentrated protein solutions: the effects of salt, pH, point mutations, and protein concentration in simulations of 1000-molecule systems

An ability to accurately simulate the dynamic behavior of concentrated macromolecular solutions would be of considerable utility in studies of a wide range of biological systems. With this goal in mind, a Brownian dynamics (BD) simulation method is reported here that allows systems to be modeled that comprise in excess of 1000 protein molecules, all of which are treated in atomic detail. Intermolecular forces are described in the method using an energy function that incorporates electrostatic and hydrophobic interactions and that is calibrated to reproduce experimental thermodynamic information with a single adjustable parameter. Using the method, BD simulations have been performed over a wide range of pH and ionic strengths for three proteins: hen egg white lysozyme (HEWL), chymotrypsinogen, and T4 lysozyme. The simulations reproduce experimental trends in second virial coefficients (B(22)) and translational diffusion coefficients, correctly capture changes in B(22) values due to single amino acid substitutions, and reveal a new explanation for the difficulties reported previously in the literature in reproducing B(22) values for protein solutions of very low ionic strength. In addition, a strong correlation is found between a residue's probability of being involved in a protein-protein contact in the simulations and its probability of being involved in an experimental crystal contact. Finally, exploratory simulations of HEWL indicate that the simulation model also gives a promising description of behavior at very high protein concentrations (approximately 250 g/L), suggesting that it may provide a suitable computational framework for modeling the complex behavior exhibited by macromolecules in cellular conditions.

Modulation of electroosmotic flow by neutral polymers

Polymer coating is widely used to modulate the fluid flow in micro- and nanometer pores and flows that are sensitive to surface properties such as electroosmotic flow. Here we report on the dissipative particle dynamics simulations of the modulation of electroosmotic flow by neutral polymers. In these coarse-grained simulations, fluid and polymers are resolved at a scale comparable to polymer size and the two-way coupling between polymer conformation and fluid flow are explicitly accounted for. The simulations indicate that, in the parameter space explored, the screening of electroosmotic flow by polymers decreases nonlinearly as the external electric field increases. Such an observation is understood by analyzing the surface coverage by polymers, height and orientation of the grafted polymers, and the two different modes of flow screening by polymer segments as a function of the external electric field. Understanding the effects and interplay of these physical processes is crucial for the rational design of polymer coating for flow control in microfluidic and nanofluidic systems.

Modeling 2D and 3D diffusion

Modeling obstructed diffusion is essential to the understanding of diffusion-mediated processes in the crowded cellular environment. Simple Monte Carlo techniques for modeling obstructed random walks are explained and related to Brownian dynamics and more complicated Monte Carlo methods. Random number generation is reviewed in the context of random walk simulations. Programming techniques and event-driven algorithms are discussed as ways to speed simulations.