RESUMO
We demonstrate that a cavitation bubble initiated by a Nd:YAG laser pulse below breakdown threshold induces crystallization from supersaturated aqueous solutions with supersaturation and laser-energy-dependent nucleation kinetics. Combining high-speed video microscopy and simulations, we argue that a competition between the dissipation of absorbed laser energy as latent and sensible heat dictates the solvent evaporation rate and creates a momentary supersaturation peak at the vapor-liquid interface. The number and morphology of crystals correlate to the characteristics of the simulated supersaturation peak.
RESUMO
Surface conductivity in the electrical double layer (EDL) is known to be affected by proton hopping and diffusion at solid-liquid interfaces. Yet, the role of surface protolysis and its kinetics on the thermodynamic and transport properties of the EDL are usually ignored as physical models consider static surfaces. Here, using a novel molecular dynamics method mimicking surface protolysis, we unveil the impact of such chemical events on the system's response. Protolysis is found to strongly affect the EDL and electrokinetic aspects with major changes in ζ potential and electro-osmotic flow.
RESUMO
In heterogeneous catalysis, reactivity and selectivity are not only influenced by chemical processes occurring on catalytic surfaces but also by physical transport phenomena in the bulk fluid and fluid near the reactive surfaces. Because these processes take place at a large range of time and length scales, it is a challenge to model catalytic reactors, especially when dealing with complex surface reactions that cannot be reduced to simple mean-field boundary conditions. As a particle-based mesoscale method, Stochastic Rotation Dynamics (SRD) is well suited for studying problems that include both microscale effects on surfaces and transport phenomena in fluids. In this work, we demonstrate how to simulate heterogeneous catalytic reactors by coupling an SRD fluid with a catalytic surface on which complex surface reactions are explicitly modeled. We provide a theoretical background for modeling different stages of heterogeneous surface reactions. After validating the simulation method for surface reactions with mean-field assumptions, we apply the method to non-mean-field reactions in which surface species interact with each other through a Monte Carlo scheme, leading to island formation on the catalytic surface. We show the potential of the method by simulating a more complex three-step reaction mechanism with reactant dissociation.
RESUMO
Surfactant molecules, known as organic friction modifiers (OFMs), are routinely added to lubricants to reduce friction and wear between sliding surfaces. In macroscale experiments, friction generally decreases as the coverage of OFM molecules on the sliding surfaces increases; however, recent nanoscale experiments with sharp atomic force microscopy (AFM) tips have shown increasing friction. To elucidate the origin of these opposite trends, we use nonequilibrium molecular dynamics (NEMD) simulations and study kinetic friction between OFM monolayers and an indenting nanoscale asperity. For this purpose, we investigate various coverages of stearamide OFMs on iron oxide surfaces and silica AFM tips with different radii of curvature. We show that the differences between the friction-coverage relations from macroscale and nanoscale experiments are due to molecular plowing in the latter. For our small tip radii, the friction coefficient and indentation depth both have a nonmonotonic dependence on OFM surface coverage, with maxima occurring at intermediate coverage. We rationalize the nonmonotonic relations through a competition of two effects (confinement and packing density) that varying the surface coverage has on the effective stiffness of the OFM monolayers. We also show that kinetic friction is not very sensitive to the sliding velocity in the range studied, indicating that it originates from instabilities. Indeed, we find that friction predominately originates from plowing of the monolayers by the leading edge of the tip, where gauche defects are created, while thermal dissipation is mostly localized in molecules toward the trailing edge of the tip, where the chains return to a more extended conformation.
RESUMO
Molecular Dynamics (MD) simulations are uniquely suitable for providing molecular-level insights into the Electric Double Layer (EDL) that forms when a charged surface is in contact with an aqueous solution. However, simulations are only as accurate in predicting EDL properties as permitted by the atomic interaction models. Experimental ζ-potential values and surface charges could provide a potentially suitable reference to validate and tune the interaction models, if not for the fact that they themselves are a product of imperfect models used to interpret the raw measurement data. Here, we present an approach to tune an interaction model by comparing Electro-Osmotic Flow (EOF) MD simulations against experimental Streaming Current (SC) measurements while minimizing potential modeling errors arising from both approaches. The point that is least susceptible to interpretation and modeling errors is argued to be at the concentration for which zero flow velocity is observed in EOF simulations and a net zero electric current is measured in SC experiments. At this concentration, the ζ-potential is also zero. We were able to match the experimental concentration at which ζ = 0 in MD simulations for a CaCl2 solution at pH 7.5 in contact with fused silica by tuning the ion-surface Lennard-Jones cross interactions. These interactions were found to greatly affect the ion distribution within the EDL and particularly the formation of inner-sphere surface-complexes, which, in turn, affects the electrokinetic flow. With the ion distribution determined explicitly, a series of properties can be calculated unambiguously, such as the capacitance needed for surface complexation models.
RESUMO
Countless molecular dynamics studies have relied on available ion and water force field parameters to model aqueous electrolyte solutions. The TIP4P/2005 model has proven itself to be among the best rigid water force fields, whereas many of the most successful ion parameters were optimized in combination with SPC/E, TIP3P, or TIP4P/Ew water. Many researchers have combined these ions with TIP4P/2005, hoping to leverage the strengths of both parameter sets. To assess if this widely used approach is justified and to provide a guide in selecting ion parameters, we investigated the transferability of various commonly used monovalent and multivalent ion parameters to the TIP4P/2005 water model. The transferability is evaluated in terms of ion hydration free energy, hydration radius, coordination number, and self-diffusion coefficient at infinite dilution. For selected ion parameters, we also investigated density, ion pairing, chemical potential, and mean ionic activity coefficients at finite concentrations. We found that not all ions are equally transferable to TIP4P/2005 without compromising their performance. In particular, ions optimized for TIP3P water were found to be poorly transferable to TIP4P/2005, whereas ions optimized for TIP4P/Ew water provided nearly perfect transferability. The latter ions also showed good overall agreement with experimental values. The one exception is that no combination of ion parameters and water model considered here was found to accurately reproduce experimental self-diffusion coefficients. Additionally, we found that cations optimized for SPC/E and TIP3P water displayed consistent underpredictions in the hydration free energy, whereas anions consistently overpredicted the hydration free energy.
RESUMO
We present a new plugin for LAMMPS for on-the-fly computation of transport properties (OCTP) in equilibrium molecular dynamics. OCTP computes the self- and Maxwell-Stefan diffusivities, bulk and shear viscosities, and thermal conductivities of pure fluids and mixtures in a single simulation. OCTP is the first implementation in LAMMPS that uses the Einstein relations combined with the order- n algorithm for the efficient sampling of dynamic variables. OCTP has low computational requirements and is easy to use because it follows the native input file format of LAMMPS. A tool for calculating the radial distribution function (RDF) of the fluid beyond the cutoff radius, while taking into account the system size effects, is also part of the new plugin. The RDFs computed from OCTP are needed to obtain the thermodynamic factor, which relates Maxwell-Stefan and Fick diffusivities. To demonstrate the efficiency of the new plugin, the transport properties of an equimolar mixture of water-methanol were computed at 298 K and 1 bar.
Assuntos
Algoritmos , Hidrodinâmica , Simulação de Dinâmica Molecular , Transporte Biológico , Difusão , ViscosidadeRESUMO
Strongly confined fluids exhibit inhomogeneous properties due to atomistic structuring in close proximity to a solid surface. State variables and transport coefficients at a solid-fluid interface vary locally and become dependent on the properties of the confining walls. However, the precise mechanisms for these effects are not known as of yet. Here, we make use of nonequilibrium molecular dynamics simulations to scrutinize the local fluid properties at the solid-fluid interface for a range of surface conditions and temperatures. We also derive microscopic relations connecting fluid viscosity and density profiles for dense fluids. Moreover, we propose empirical ready-to-use relations to express the average density and viscosity in the channel as a function of temperature, wall interaction strength, and bulk density or viscosity. Such relations are key to technological applications such as micro-/nanofluidics and tribology but also natural phenomena.
RESUMO
Slip lengths reported from molecular dynamics (MD) simulations of water flow in graphene nanochannels show significant scatter in the literature. These discrepancies are in part due to the used water models. We demonstrate self-consistent comparisons of slip characteristics between the SPC, SPC/E, SPC/Fw, TIP3P, TIP4P, and TIP4P/2005 water models. The slip lengths are inferred using an analytical model that employs the shear viscosity of water and channel average velocities obtained from nonequilibrium MD simulations. First, viscosities for each water model are quantified using MD simulations of counterflowing, force-driven flows in periodic domains in the absence of physical walls. While the TIP4P/2005 model predicts water viscosity at the specified thermodynamic state with 1.7% error, the predictions of SPC/Fw and SPC/E models exhibit 13.9% and 23.1% deviations, respectively. Water viscosities obtained from SPC, TIP4P, and TIP3P models show larger deviations. Next, force-driven water flows in rigid (cold) and thermally vibrating (thermal) graphene nanochannels are simulated, resulting in pluglike velocity profiles. Large differences in the flow velocities are observed depending on the used water model and to a lesser extent on the choice of rigid vs thermal walls. Depending on the water model, the slip length of water on cold graphene walls varied between 34.2 nm and 62.9 nm, while the slip lengths of water on thermal graphene walls varied in the range of 38.1 nm-84.3 nm.
RESUMO
Lipid bilayers composed of non-hydroxy sphingosine ceramide (CER NS), cholesterol (CHOL), and free fatty acids (FFAs), which are components of the human skin barrier, are studied via molecular dynamics simulations. Since mixtures of these lipids exist in dense gel phases with little molecular mobility at physiological conditions, care must be taken to ensure that the simulations become decorrelated from the initial conditions. Thus, we propose and validate an equilibration protocol based on simulated tempering, in which the simulation takes a random walk through temperature space, allowing the system to break out of metastable configurations and hence become decorrelated from its initial configuration. After validating the equilibration protocol, which we refer to as random-walk molecular dynamics, the effects of the lipid composition and ceramide tail length on bilayer properties are studied. Systems containing pure CER NS, CER NS + CHOL, and CER NS + CHOL + FFA, with the CER NS fatty acid tail length varied within each CER NS-CHOL-FFA composition, are simulated. The bilayer thickness is found to depend on the structure of the center of the bilayer, which arises as a result of the tail-length asymmetry between the lipids studied. The hydrogen bonding between the lipid headgroups and with water is found to change with the overall lipid composition, but is mostly independent of the CER fatty acid tail length. Subtle differences in the lateral packing of the lipid tails are also found as a function of CER tail length. Overall, these results provide insight into the experimentally observed trend of altered barrier properties in skin systems where there are more CERs with shorter tails present.
Assuntos
Ceramidas/química , Células Epidérmicas/citologia , Bicamadas Lipídicas/química , Ligação de Hidrogênio , Conformação Molecular , Simulação de Dinâmica MolecularRESUMO
We introduce an analytical method to predict the slip length (L s) in cylindrical nanopores using equilibrium molecular dynamics (EMD) simulations, following the approach proposed by Sokhan and Quirke for planar channels [39]. Using this approach, we determined the slip length of water in carbon nanotubes (CNTs) of various diameters. The slip length predicted from our method shows excellent agreement with the results obtained from nonequilibrium molecular dynamics (NEMD) simulations. The data show a monotonically decreasing slip length with an increasing nanotube diameter. The proposed EMD method can be used to precisely estimate slip length in high slip cylindrical systems, whereas, L s calculated from NEMD is highly sensitive to the velocity profile and may cause large statistical errors due to large velocity slip at the channel surface. We also demonstrated the validity of the EMD method in a BNNT-water system, where the slip length is very small compared to that in a CNT pore of similar diameter. The developed method enables us to calculate the interfacial friction coefficient directly from EMD simulations, while friction can be estimated using NEMD by performing simulations at various external driving forces, thereby increasing the overall computational time. The EMD analysis revealed a curvature dependence in the friction coefficient, which induces the slip length dependency on the tube diameter. Conversely, in flat graphene nanopores, both L s and friction coefficient show no strong dependency on the channel width.
RESUMO
We investigate thermally driven water droplet transport on graphene and hexagonal boron nitride (h-BN) surfaces using molecular dynamics simulations. The two surfaces considered here have different wettabilities with a significant difference in the mode of droplet transport. The water droplet travels along a straighter path on the h-BN sheet than on graphene. The h-BN surface produced a higher driving force on the droplet than the graphene surface. The water droplet is found to move faster on h-BN surface compared to graphene surface. The instantaneous contact angle was monitored as a measure of droplet deformation during thermal transport. The characteristics of the droplet motion on both surfaces is determined through the moment scaling spectrum. The water droplet on h-BN surface showed the attributes of the super-diffusive process, whereas it was sub-diffusive on the graphene surface.
RESUMO
Molecular dynamics simulations are widely employed to analyze water and ion permeation through nanoporous membranes for reverse osmosis applications. In such simulations, water models play an important role in accurately reproducing the properties of water. We investigated the water and ion transport across a hydroxyl (OH) functionalized graphene nanopore using six water models: SPC, SPC/E, SPC/Fw, TIP3P, TIP4P, and TIP4P/2005. The water flux thus obtained varied up to 84% between the models. The water and ion flux showed a correlation with the bulk transport properties of the models such as the diffusion coefficient and shear viscosity. We found that the hydrogen-bond lifetime, resulting from the partial charges of the model, influenced the flux. Our results are useful in the selection of a water model for computer simulations of desalination using nanomembranes. Our findings also suggest that lowering the hydrogen-bond lifetime and enhancing the rate of diffusion of water would lead to enhanced water/ion flux.
RESUMO
Thermal-gradient induced transport of ionic liquid (IL) and water droplets through a carbon nanotube (CNT) is investigated in this study using molecular dynamics simulations. Energetic analysis indicates that IL transport through a CNT is driven primarily by the fluid-solid interaction, while fluid-fluid interactions dominate in water-CNT systems. Droplet diffusion analysis via the moment scaling spectrum reveals sub-diffusive motion of the IL droplet, in contrast to the self-diffusive motion of the water droplet. The Soret coefficient and energetic analysis of the systems suggest that the CNT shows more affinity for interaction with IL than with the water droplet. Thermophoretic transport of IL is shown to be feasible, which can create new opportunities in nanofluidic applications.
RESUMO
Although the importance of temperature control in nonequilibrium molecular dynamics simulations is widely accepted, the consequences of the thermostatting approach in the case of strongly confined fluids are underappreciated. We show the strong influence of the thermostatting method on the water transport in carbon nanotubes (CNTs) by considering simulations in which the system temperature is controlled via the walls or via the fluid. Streaming velocities and mass flow rates are found to depend on the tube flexibility and on the thermostatting algorithm, with flow rates up to 20% larger when the walls are flexible. The larger flow rates in flexible CNTs are explained by a lower friction coefficient between water and the wall. Despite the lower friction, a larger solid-fluid interaction energy is found for flexible CNTs than for rigid ones. Furthermore, a comparison of thermostat schemes has shown that the Berendsen and Nosé-Hoover thermostats result in very similar transport rates, while lower flow rates are found under the influence of the Langevin thermostat. These findings illustrate the significant influence of the thermostatting methods on the simulated confined fluid transport.
RESUMO
Single- and multicomponent lipid bilayers of 1,2-dipalmitoyl-sn-glycero-3-phosphatidylcholine (DPPC), 1,2-distearoyl-sn-glycero-3-phosphatidylcholine (DSPC), isostearyl isostearate, and heptadecanoyl heptadecanoate in the gel phase are studied via molecular dynamics simulations. It is shown that the structural properties of multicomponent bilayers can deviate strongly from the structures of their single-component counterparts. Specifically, the lipid mixtures are shown to adopt a compact packing by offsetting the positioning depths at which different lipid species are located in the bilayer. This packing mechanism affects the area per lipid, the bilayer height, and the chain tilt angles and has important consequences for other bilayer properties, such as interfacial hydrogen bonding and bilayer permeability. In particular, the simulations suggest that bilayers containing isostearyl isostearate or heptadecanoyl heptadecanoate are less permeable than pure 1,2-dipalmitoyl-sn-glycero-3-phosphatidylcholine or DSPC bilayers. Furthermore, hydrogen-bond analysis shows that the residence times of lipid-water hydrogen bonds depend strongly on the bilayer composition, with longer residence times for bilayers that have a higher DSPC content. The findings illustrate and explain the fundamental differences between the properties of single- and multicomponent bilayers.
Assuntos
Bicamadas Lipídicas/química , Géis , Ligação de Hidrogênio , Modelos Moleculares , Conformação Molecular , Água/químicaRESUMO
Monovalent and divalent aqueous electrolytes confined in negatively charged porous silica are studied by means of molecular simulations including free energy calculations. Owing to the strong cation adsorption at the surface, surface charge overcompensation (overscreening) occurs which leads to an effective positive surface next to the Stern layer, followed by a negatively charged diffuse layer. A simple Poisson-Boltzmann model in which the single-ion potential of mean force is introduced is shown to capture the most prominent features of ion density profiles near an amorphous silica surface. Nevertheless, due to its mean-field nature, which fails to account for correlations, this simple model does not predict overscreening corresponding to charge inversion at the surface. Such an overscreening drastically affects the transport of confined electrolytes as it leads to flow reversal when subjected to an electric field. A simple continuum theory is shown to capture how the electro-osmotic flow is affected by overscreening and by the apparent enhanced viscosity of the confined electrolytes. Comparison with available experimental data is discussed, as well as the implications of these phenomena for ζ-potential measurements.
RESUMO
Molecular simulations are used to compare the structure and dynamics of conventional and radioactive aqueous electrolytes: chloride solutions with sodium, potassium, cesium, calcium, and strontium. The study of Cs(+) and Sr(2+) is important because these radioactive ions can be extremely harmful and are often confused by living organisms for K(+) and Ca(2+), respectively. Na(+), Ca(2+), and Sr(2+) are strongly bonded to their hydration shell because of their large charge density. We find that the water molecules in the first hydration shell around Na(+) form hydrogen bonds between each other, whereas molecules in the first hydration shell around Ca(2+) and Sr(2+) predominantly form hydrogen bonds with water molecules in the second shell. In contrast to these three ions, K(+) and Cs(+) have low charge densities so that they are weakly bonded to their hydration shell. Overall, the structural differences between Ca(2+) and Sr(2+) are small, but the difference between their coordination numbers relative to their surface areas could potentially be used to separate these ions. Moreover, the different decays of the velocity-autocorrelation functions corresponding to these ions indicates that the difference in mass could be used to separate these cations. In this work, we also propose a new definition of the pairing time that is easy to calculate and of physical significance regardless of the problem at hand.
Assuntos
Íons/química , Água/química , Cálcio/química , Césio/química , Cloretos/química , Simulação por Computador , Ligação de Hidrogênio , Modelos Moleculares , Potássio/química , Radioisótopos/química , Sais/química , Sódio/química , Estrôncio/química , ViscosidadeRESUMO
Oscillatory flow in confined spaces is central to understanding physiological flows and rational design of synthetic periodic-actuation based micromachines. Using theory and experiments on oscillating flows generated through a laser-induced cavitation bubble, we associate the dynamic bubble size (fluid velocity) and bubble lifetime to the laser energy supplied-a control parameter in experiments. Employing different channel cross-section shapes, sizes and lengths, we demonstrate the characteristic scales for velocity, time and energy to depend solely on the channel geometry. Contrary to the generally assumed absence of instability in low Reynolds number flows ( < 1000 ), we report a momentary flow distortion that originates due to the boundary layer separation near channel walls during flow deceleration. The emergence of distorted laminar states is characterized using two stages. First the conditions for the onset of instabilities is analyzed using the Reynolds number and Womersley number for oscillating flows. Second the growth and the ability of an instability to prevail is analyzed using the convective time scale of the flow. Our findings inform rational design of microsystems leveraging pulsatile flows via cavitation-powered microactuation.
RESUMO
Non-equilibrium molecular dynamics simulations of an atomic fluid under shear flow, planar elongational flow, and a combination of shear and elongational flow are unified consistently with a tensorial model over a wide range of strain rates. A model is presented that predicts the pressure tensor for a non-Newtonian bulk fluid under a homogeneous planar flow field. The model provides a quantitative description of the strain-thinning viscosity, pressure dilatancy, deviatoric viscoelastic lagging, and out-of-flow-plane pressure anisotropy. The non-equilibrium pressure tensor is completely described through these four quantities and can be calculated as a function of the equilibrium material constants and the velocity gradient. This constitutive framework in terms of invariants of the pressure tensor departs from the conventional description that deals with an orientation-dependent description of shear stresses and normal stresses. The present model makes it possible to predict the full pressure tensor for a simple fluid under various types of flows without having to produce these flow types explicitly in a simulation or experiment.