Effect of nanometric-scale roughness on slip at the wall of simple fluids

It is commonly acknowledged that roughness decreases the aptitude of simple liquids to exhibit flow with slip at solid interfaces. Most available studies have, however, been conducted on substrates for which both the surface chemistry and the roughness were varied simultaneously, making it difficult to identify their respective role on wall slip. To overcome this difficulty, we have developed a series of surfaces formed by grafting hyperbranched polymeric nanoparticles on a smooth, dense, self-assembled monolayer of SiH-terminated short poly(dimethylsiloxane) oligomers, allowing us to vary independently the surface density, the height, and the width of the grafted nanoparticles, and thereby the roughness parameters, while keeping similar surface chemistry. On such substrates, the boundary condition for the flow velocity of hexadecane has been characterized through near-field laser velocimetry. We demonstrate that decreasing the wavelength of the roughness at a fixed height strongly decreases slip, while increasing the height of the nanoparticles at a fixed aspect ratio of the roughness also dramatically affects slippage.     

Multi length scale analysis of the microstructure in sticky sphere dispersions during shear flow

The effect of shear flow on the microstructure in a weakly aggregated suspension is investigated. Monodisperse small silica particles with a grafted layer of 1-octadecanol are dispersed in n-tetradecane, yielding a thermoreversible sticky sphere model suspension. A combination of small angle light scattering and ultra small and small-angle X-ray scattering techniques have been used, in situ and time resolved, to study the flow-induced anisotropy of the microstructure. In this manner, the length scales from the single particle size to that of the spatial organization of the aggregates can be covered. Harmonic expansion of the structure factor demonstrates that anisotropy develops in the microstructure on all relevant length scales. Possible real space interpretations of the scattering information are discussed in conjunction with implications for the nonlinear rheological behavior.     

Molecular dynamics simulation study of the dynamics of fluids at solid-liquid interfaces

The dynamics of fluids at solid-liquid interfaces is investigated. In particular, we consider a simple Lennard-Jones fluid as well as a melt of hexadecane chains. For the Lennard-Jones fluid, the numerical results are compared with analytical calculations based on the diffusion equation, which shows that the numerical results can very well by described by the solution of the diffusion equation for reflecting surfaces. The diffusion coefficient is practically independent of the position within the film, although the fluid is inhomogeneous perpendicular to the surface. In contrast, the dynamics of the centers of mass of hexadecane molecules perpendicular to repulsive surfaces is severely slowed down due to their extended and anisotropic nature and cannot be described by a single particle diffusion equation.     

Molecular dynamics study of atomic transport properties in rapidly cooling liquid copper

Based on Mei's embedded atom model molecular dynamics simulations have been performed to investigate the rapidly cooling processes of Cu. The atomic transport property, namely the self-diffusion coefficient, is computed in the liquid state, and the results near the melting point of Cu are in good agreement with experimental data and other computational values. The atom diffusion movements during the long period of relaxation have been also studied around the solidification temperature Tc. To describe the complex microstructural evolutions during the rapidly cooling processes and the long relaxation processes, the pair correlation function and the pair analysis technique are used. It is demonstrated that the crystallization of amorphous Cu is caused by the atomic diffusion.     

Molecular dynamics study of screening at ionic surfaces

Molecular dynamics simulations of NaCl fluid are used to understand the behavior of ionic fluid to screen the field generated by charges on the ionic crystal surfaces in absence of any external electric field. The NaCl fluid in the strongly coupled regime (corresponding to the melt) in contact with the charged octopolar (111) NaCl surface shows that the spatial correlations decay in an oscillatory manner, with a screening length lambdaQ given by the envelope of the damped oscillations. By contrast to the Debye-Huckel theory, in the strongly coupled regime, lambdaQ increases with increasing coupling strength (also seen in bulk ionic simulations). The NaCl fluid confined between neutral (100) NaCl surfaces also shows weak oscillatory charge decay near the surface. Similar oscillatory exponential decay was seen when the NaCl fluid was confined between two analytically smooth neutral walls. The origin of these oscillations was due to the difference in ion sizes. NaCl fluid confined between neutral octopolar (110) and dipolar (110) surface show stronger density oscillations than (100) surface but comparatively very weak charge oscillations. This paper shows that the strength of the charges on the crystal surfaces is enough to induce a characteristic spatial distribution of charges in the contacting fluid and the extent of distribution depends on the type of surface.     

Nonequilibrium molecular dynamics simulation study on the orientation transition in the amphiphilic lamellar phase under shear flow

By the extensive large-scale nonequilibrium molecular dynamics simulation on an effective generic model-A2B2 tetramer for amphiphiles, we investigate the shear-induced parallel to perpendicular orientation transition in the lamellar phase as a function of segregation degree and shear rate. Under low rate shear flow the evolution of parallel lamellar configurations at different segregation strengths shows a similar kinetic pathway independent of the segregation degree. While under high rate shear flow in which the lifetime of undulation instability exceeds the characteristic time of the applied shear flow, the kinetic pathway of the shear-induced parallel-to-perpendicular orientation transition in lamellar systems is the segregation degree dependent. Comparing the temporal mesoscopic domain morphology, the microscopic chain conformation, and macroscopic observable-viscosity changes with the experimentally proposed mechanisms, we find that the undulation instability, partial breakup of monodomain, grain rotation, and recombination combined with defect migration and annihilation are the kinetic pathway for the parallel-to-perpendicular orientation transition in the lamellar phase in or near the intermediate segregation limit, and that the undulation instability, domain dissolution, and reformation along the preferred direction combined with defect migration and annihilation are the kinetic pathway for the parallel-to-perpendicular orientation transition in the lamellar phase close to the order-to-disorder phase transition point. A detailed underlying microscopic picture of the alignment process illustrates that the orientation transition is driven by the alignment of molecules with shear flow. The orientation diagram that characterizes the steady-state orientations as a function of shear rate and attractive potential depth is built, in which the attractive potential depth takes the role of an inverse temperature, somewhat like the Flory-Huggins interaction parameter. The microscopic mechanism of the critical orientation transition condition is discussed.     

Molecular dynamics study of the methanol effect on the membrane morphology of perfluorosulfonic ionomers

The membranes of a perfluorosulfonic acid polymer swollen in 10-80 wt % methanol solution were investigated to elucidate the methanol effect on their morphologies, such as size of the solvent cluster, solvent location, and polymer structure, by using isothermal-isobaric molecular dynamics simulations. In higher methanol concentrations, we found less-spherical solvent aggregation and a more spread polymer structure because of the ampholytic nature of methanol. The partial radial distribution functions between solvent oxygen and fluorocarbons, which are composed of the main chain, clearly show that methanol is located closer to the polymer matrix than water. On the other hand, water is preferentially located in the vicinity of an acidic headgroup, SO(3)(-), compared with methanol, although both have similar attractive interaction energies to the acidic group. Furthermore, we discussed solvent dynamics and hydrogen bonding between sulfonic oxygen and solvent O-H groups.     

Crosslinking PMMA: Molecular dynamics investigation of the shear response

Crosslinking can fundamentally change the mechanical properties of a linear glassy polymer. It has been experimentally observed that when lightly crosslinked, poly(methyl‐methacrylate) (PMMA) has a characteristically more ductile response to mechanical loading than does linear PMMA despite having a higher glass transition temperature. Here, molecular dynamics (MD) simulations are used to investigate conformational and energetic differences between linear PMMA and lightly crosslinked PMMA under shear deformation. As consistent with experiments, crosslinked PMMA is found to have a reduced yield stress relative to linear PMMA. Using the probing capabilities of our explicit atom MD approach, it is observed that while the crosslinks have a minimal direct energy contribution to the total system, they can alter how the main chains conform to macroscopic loading. In crosslinked PMMA, the backbone aligns more with the direction of external loading, thereby reducing the force applied to (and associated deformation of) the polymer bonds. © 2013 Wiley Periodicals, Inc. J. Polym. Sci., Part B: Polym. Phys. 2014 , 52, 444–449     

Numerical study of electroosmotic slip flow of fractional Oldroyd-B fluids at high zeta potentials

In this paper, an investigation of the electroosmotic flow of fractional Oldroyd-B fluids in a narrow circular tube with high zeta potential is presented. The Navier linear slip law at the walls is considered. The potential field is applied along the walls described by the nonlinear Poisson–Boltzmann equation. It's worth noting here that the linear Debye–Hückel approximation can't be used at the condition of high zeta potential and the exact solution of potential in cylindrical coordinates can't be obtained. Therefore, the Matlab bvp4c solver method and the finite difference method are employed to numerically solve the nonlinear Poisson–Boltzmann equation and the governing equations of the velocity distribution, respectively. To verify the validity of our numerical approach, a comparison has been made with the previous work in the case of low zeta potential and the excellent agreement between the solutions is clear. Then, in view of the obtained numerical solution for the velocity distribution, the numerical solutions of the flow rate and the shear stress are derived. Furthermore, based on numerical analysis, the influence of pertinent parameters on the potential distribution and the generation of flow is presented graphically.     

The effect of shear flow on the conformation of polyelectrolyte brushes

The behavior of flat polyelectrolyte brushes under the action of a lateral force or flow was studied. Special attention was focused on the case when a lateral force acts on a brush that occurs near the point of phase transition from the swollen state to the collapsed state. The difference between phase transitions in a brush induced by isotropic and anisotropic interactions is analyzed. As examples of such transitions, the collapse of a polyelectrolyte brush upon cooling and the nematic collapse of an anisotropic brush are considered. It was shown that lateral force (flow), exerting a marked effect on the nematic collapse of an anisotropic brush, has no practical effect on the collapse of a brush with isotropic interactions.     

Molecular dynamics study of the hydration of the hydroxyl radical at body temperature

Classical molecular dynamics (MD) simulation of ˙OH in liquid water at 37 °C has been performed using flexible models of the solute and solvent molecules. We derived the Morse function describing the bond stretching of the radical and the potential for ˙OH-H(2)O interactions, including short-range interactions of hydrogen atoms. Scans of the potential energy surface of the ˙OH-H(2)O complex have been performed using the DFT method with the B3LYP functional and the 6-311G(d,p) basis set. The DFT-derived partial charges, ±0.375e, and the equilibrium bond-length, 0.975 ?, of ˙OH resulted in the dipole moment of 1.76 D. The radical-water radial distribution functions revealed that ˙OH is not built into the solvent structure but it rather occupies distortions or cavities in the hydrogen-bonded network. The solvent structure at 37 °C has been found to be the same as that of pure water. The hydration cage of the radical comprises 13-14 water molecules. The estimated hydration enthalpy -42 ± 5 kJ mol(-1) is comparable with the experimental value -39 ± 6 kJ mol(-1) for 25 °C. Inspection of hydrogen bonds showed the importance of short-range interaction of hydrogen atoms and indicated that neglect of the angular condition greatly overestimates the number of the H-acceptor radical-water bonds. The mean number ?n = 0.85 of radical-water H-bonds has been calculated using geometric definition of H-bond and ?n = 0.62 has been obtained when the energetic condition, E(da)≤-8 kJ mol(-1), was additionally considered. The continuous lifetimes of 0.033 ps and 0.024 ps have been estimated for the radical H-donor and the H-acceptor bonds, respectively. Within statistical uncertainty the radical self-diffusion coefficient, (2.9 ± 0.6) × 10(-9) m(2) s(-1), is the same as (3.1 ± 0.5) × 10(-9) m(2) s(-1) calculated for water in solution and in pure solvent. To the best of our knowledge, this is the first study of the ˙OH(aq) properties at a biologically relevant body temperature.     

Effect of wall slip on the shear flow of a polymer in a channel with a wavy bottom

The effect of stick and wall slip boundary conditions on the specific features of the shear flow of viscous polymers in a confined two-dimensional channel with a wavy bottom is studied. The distribution of flow-rate disturbances across the transverse cross section of the channel is calculated by the numerical simulation of the Navier-Stokes equation for an incompressible fluid at arbitrary amplitudes and an arbitrary wave number of the wall. The wall slip is modeled by the introduction of a thin layer of a low-viscosity fluid at the bottom face. Slippage leads to a marked enhancement of flow rate disturbances including inertial advection. The results agree with the known analytical solutions for the low-amplitude wall wave.     

Molecular dynamics simulation of planar elongational flow at constant pressure and constant temperature

Molecular dynamics simulations of liquid systems under planar elongational flow have mainly been performed in the NVT ensemble. However, in most material processing techniques and common experimental settings, at least one surface of the fluid is kept in contact with the atmosphere, thus maintaining the sample in the NpT ensemble. For this reason, an implementation of the Nose-Hoover integral-feedback mechanism for constant pressure is presented, implemented via the SLLOD algorithm for elongational flow. The authors test their procedure for an atomic liquid and compare the viscosity obtained with that in the NVT ensemble. The scheme is easy to implement, self-starting and reliable, and can be a useful tool for the simulation of more complex liquid systems, such as polymer melts and solutions.     

Molecular dynamics study of nanoparticle stability at liquid interfaces: effect of nanoparticle-solvent interaction and capillary waves

While the interaction of colloidal particles (sizes in excess of 100 nm) with liquid interfaces may be understood in terms of continuum models, which are grounded in macroscopic properties such as surface and line tensions, the behaviour of nanoparticles at liquid interfaces may be more complex. Recent simulations [D. L. Cheung and S. A. F. Bon, Phys. Rev. Lett. 102, 066103 (2009)] of nanoparticles at an idealised liquid-liquid interface showed that the nanoparticle-interface interaction range was larger than expected due, in part, to the action of thermal capillary waves. In this paper, molecular dynamics simulations of a Lennard-Jones nanoparticle in a binary Lennard-Jones mixture are used to confirm that these previous results hold for more realistic models. Furthermore by including attractive interactions between the nanoparticle and the solvent, it is found that the detachment energy decreases as the nanoparticle-solvent attraction increases. Comparison between the simulation results and recent theoretical predictions [H. Lehle and M. Oettel, J. Phys. Condens. Matter 20, 404224 (2008)] shows that for small particles the incorporation of capillary waves into the predicted effective nanoparticle-interface interaction improves agreement between simulation and theory.     

Molecular dynamics study of polymeric chemical gels: Effect of persistence length and crosslinker density

Using molecular dynamics, we study the formation of chemical gels from an initial solution of reactive polymers that undergo a crosslinking reaction. We study the effect of the polymer persistence length and different densities of crosslinkers along the chains. As the reaction progresses, different structural features are identified in the system leading to the development of a percolated cluster. These features are (a) single strands, (b) double strands, and (c) bridges. We found that the total numbers of these three kinds of features are roughly independent of the persistence length; however, the average lengths of single and double strands grow with this variable. The average length of double strands strongly increases with increasing crosslinker density and the amount of single strands sharply falls as crosslinker density grows. We also found that general structural features of polymer networks are highly dependent on chain persistence length and crosslinker density. Fully flexible chains with high density of crosslinkers result in inhomogeneous network structures with large voids. In contrast, precursor chains with high rigidity and scarce number of crosslinkers result in homogeneous networks having small cavities. © 2019 Wiley Periodicals, Inc. J. Polym. Sci., Part B: Polym. Phys. 2019 , 57, 1343–1350     

Molecular dynamics study of the effect of calcium ions on the monolayer of SDC and SDSn surfactants at the vapor/liquid interface

The effect of Ca(2+) ions on the hydration shell of sodium dodecyl carboxylate (SDC) and sodium dodecyl sulfonate (SDSn) monolayer at vapor/liquid interfaces was studied using molecular dynamics simulations. For each surfactant, two different surface concentrations were used to perform the simulations, and the aggregation morphologies and structural details have been reported. The results showed that the aggregation structures relate to both the surface coverage and the calcium ions. The divalent ions can screen the interaction between the polar head and Na(+) ions. Thus, Ca(2+) ions locate near the vapor/liquid interface to bind to the headgroup, making the aggregations much more compact via the salt bridge. The potential of mean force (PMF) between Ca(2+) and the headgroups shows that the interaction is decided by a stabilizing solvent-separated minimum in the PMF. To bind to the headgroup, Ca(2+) should overcome the energy barrier. Among contributions to the PMF, the major repulsive interaction was due to the rearrangement of the hydration shell after the calcium ions entered into the hydration shell of the headgroup. The PMFs between the headgroup and Ca(2+) in the SDSn systems showed higher energy barriers than those in the SDC systems. This result indicated that SDSn binds the divalent ions with more difficulty compared with SDC, so the ions have a strong effect on the hydration shell of SDC. That is why sulfonate surfactants have better efficiency in salt solutions with Ca(2+) ions for enhanced oil recovery.     

Molecular dynamics study on evaporation and condensation of n-dodecane at liquid-vapor phase equilibria

Molecular dynamics simulations are performed to study the evaporation and condensation of n-dodecane (C(12)H(26)) at temperatures in the range 400-600 K. A modified optimized potential for liquid simulation model is applied to take into account the Lennard-Jones, bond bending and torsion potentials with the bond length constrained. The equilibrium liquid-vapor n-dodecane interface thickness is predicted to be ~1.2-2.0 nm. It is shown that the molecular chains lie preferentially parallel to the interface in the liquid-vapor transition region. The predicted evaporation/condensation coefficient decreased from 0.9 to 0.3 when temperature increased from 400 to 600 K. These values can be used for the formulation of boundary conditions in the kinetic modeling of droplet heating and evaporation processes; they are noticeably different from those predicted by the transition state theory. We also present the typical molecular behaviors in the evaporation and condensation processes. The molecular exchange in condensation, typical for simple molecules, has never been observed for n-dodecane molecular chains.     

Molecular dynamics study of the ribosomal A-site

Many aminoglycosidic antibiotics target the A-site of 16S RNA in the small ribosomal subunit and affect the fidelity of protein translation in bacteria. Upon binding, aminoglycosides displace two adenines (A1492 and A1493 for E. coli numbering) that are involved in tRNA anticodon loop recognition. The major difference in the aminoglycosidic binding site between the prokaryota and eukaryota is an adenine into guanine substitution in the position 1408. This mutation likely affects the dynamics of near A1492 and A1493 and hinders the binding of aminoglycosides to eukaryotic ribosomes. With multiple 20 ns long all-atom molecular dynamics simulations, we study the flexibility of a 22 nucleotide RNA fragment which mimics the aminoglycosidic binding site. Simulations are carried out for both native and A1408G mutated RNA as well as for their complexes with aminoglycosidic representative paromomycin. We observe intra- and extrahelical configurations of A1492 and A1493, which differ between the prokaryotic and the mutated structure. We obtain configurations of the A-site that are also observed in the NMR and crystal structures. Our studies show the differences in the internal mobility of the A-site, as well as that in ion and water density distributions inside of the binding cleft, between the prokaryotic and mutated RNA. We also compare the performance of two force field parameters for RNA, Amber and Charmm.     

Molecular dynamics simulation study on the transient response of solvation structure during the translational diffusion of solute

The transient response function of the density profile of the solvent around a solute during the translational diffusion of the solute is formulated based on the generalized Langevin formalism. The resultant theory is applied to both neat Lennard-Jones fluids and cations in liquid water, and the response functions are obtained from the analysis of the molecular dynamics simulations. In the case of the self-diffusion of Lennard-Jones fluids, the responses of the solvation structures are in harmony with conventional pictures based on the mode-coupling theory, that is, the binary collision in the low-density fluids, the backflow effect from medium to high density fluids, and the backscatter effect in the liquids near the triple point. In the case of cations in water, the qualitative behavior is strongly dependent on the size of cations. The pictures similar to simple dense liquids are obtained for the large ion and the neutral molecule, while the solvent waters within the first solvation shell of small ions show an oscillatory response in the short-time region. In particular, the oscillation is remarkably underdumped for lithium ion. The origin of the oscillation is discussed in relation to the theoretical treatment of the translational diffusion of ions in water.     

Molecular dynamics simulations of ionic liquid-vapour interfaces: effect of cation symmetry on structure at the interface

The influence of alkyl chain symmetry of the imidazolium cation on the structure and properties of the ionic liquid-vapour interface has been addressed through molecular dynamics simulations. The anion chosen is bis(trifluoromethylsulfonyl)imide (NTf(2)). Profiles of number densities, orientation of cations, charge density, electrostatic potential, and surface tension have been obtained. At the interface, both cations and anions were present, and the alkyl chains of the former preferred to orient out into the vapour phase. A large fraction of cations preferred to be oriented with their ring-normal parallel to the surface and alkyl chains perpendicular to it. These orientational preferences are reduced in ionic liquids with symmetric cations. Although the charge densities at the interface were largely negative, an additional small positive charge density has been observed for systems with longer alkyl chains. The electrostatic potential difference developed between the liquid and the vapour phases were positive and decreased with increasing length of the alkyl group. The calculated surface tension of the liquids also decreased with increasing alkyl chain length, in agreement with experiment. The surface tension of an ionic liquid with symmetric cation was marginally higher than that of one with an asymmetric, isomeric cation.