Computation of the properties of liquid neon, methane, and gas helium at low temperature by the Feynman-Hibbs approach

The properties of liquid methane, liquid neon, and gas helium are calculated at low temperatures over a large range of pressure from the classical molecular-dynamics simulations. The molecular interactions are represented by the Lennard-Jones pair potentials supplemented by quantum corrections following the Feynman-Hibbs approach. The equations of state, diffusion, and shear viscosity coefficients are determined for neon at 45 K, helium at 80 K, and methane at 110 K. A comparison is made with the existing experimental data and for thermodynamical quantities, with results computed from quantum numerical simulations when they are available. The theoretical variation of the viscosity coefficient with pressure is in good agreement with the experimental data when the quantum corrections are taken into account, thus reducing considerably the 60% discrepancy between the simulations and experiments in the absence of these corrections.     

How short-range attractions impact the structural order, self-diffusivity, and viscosity of a fluid

We present molecular simulation data for viscosity, self-diffusivity, and the local structural ordering of (i) a hard-sphere fluid and (ii) a square-well fluid with short-range attractions. The latter fluid exhibits a region of dynamic anomalies in its phase diagram, where its mobility increases upon isochoric cooling, which is found to be a subset of a larger region of structural anomalies, in which its pair correlations strengthen upon isochoric heating. This "cascade of anomalies" qualitatively resembles that found in recent simulations of liquid water. The results for the hard-sphere and square-well systems also show that the breakdown of the Stokes-Einstein relation upon supercooling occurs for conditions where viscosity and self-diffusivity develop different couplings to the degree of pairwise structural ordering of the liquid. We discuss how these couplings reflect dynamic heterogeneities. Finally, we note that the simulation data suggest how repulsive and attractive glasses may generally be characterized by two distinct levels of short-range structural order.     

Modeling spontaneous formation of precursor nanoparticles in clear-solution zeolite synthesis

We present a lattice model describing the formation of silica nanoparticles in the early stages of the clear-solution templated synthesis of silicalite-1 zeolite. Silica condensation/hydrolysis is modeled by a nearest-neighbor attraction, while the electrostatics are represented by an orientation-dependent, short-range interaction. Using this simplified model, we show excellent qualitative agreement with published experimental observations. The nanoparticles are identified as a metastable state, stabilized by electrostatic interactions between the negatively charged silica surface and a layer of organic cations. Nanoparticle size is controlled mainly by the solution pH, through nanoparticle surface charge. The size and concentration of the charge-balancing cation are found to have a negligible effect on nanoparticle size. Increasing the temperature allows for further particle growth by Ostwald ripening. We suggest that this mechanism may play a role in the growth of zeolite crystals.     

The possibility of different time scales in the dynamics of pore imbibition

To model the imbibition of liquids into porous solids, use is often made of the Lucas-Washburn equation, which relates the distance of penetration of a liquid at a given time to the pore radius, the viscosity and surface tension of the liquid, and the effective contact angle between the liquid and the solid. In this paper, we extend previous large-scale molecular dynamics simulations to show how this tool can be used to study the details of liquid imbibition, including the impact of the contact angle on the dynamics of penetration and the evolution of the internal flow field. In particular, we show that the asymptotic behavior of the contact angle versus time for a completely wetting liquid is given by approximately t(-1/4).     

Analysis on wetting and local dynamic properties of single water droplet on a polarized solid surface: a molecular dynamics study

Molecular dynamics simulations of single water droplets on a solid surface were carried out in order to investigate the effects that the Coulomb interaction between liquid and solid molecules has on wetting behavior by appending vertical electric polarization on a solid surface. The water droplet became more wettable both on upward and downward polarized surfaces, although structures of the adsorption layer appearing near the solid surface were clearly different, and the relation between droplet contact angle and surface polarization was also different for upward and downward polarization directions. The probability density distribution of molecular orientation around the adsorption layer indicated that preferable water molecule orientations varied largely by the surface polarization, and the rotational mobility around the preferable orientations was also affected. The dynamic property due to this rotational mobility was clearly captured by means of distribution of rotational diffusion coefficient, which potentially corresponded to local viscosity distribution.     

Diffusion of methanol in zeolite NaY: a molecular dynamics study

Molecular dynamics simulations were performed in order to obtain a detailed understanding of the self-diffusion mechanisms of methanol in the zeolite NaY system. We derived a new force-field term to describe the interactions between the methanol molecules and the extraframework cations. From the simulations, we show that diffusive behavior in the high-temperature range consists of a combination of both short- and long-range motions at low and intermediate loadings. This type of motion is characterized by an activation energy that decreases as the loading increases. At low loadings, we also observe short-range diffusive behavior based on a surface-mediated mechanism. The short-range behavior corresponds to motion only on the length scale of an FAU supercage, whereas the long-range behavior involves intercage diffusion. For the saturation loading corresponding to 96 methanol molecules per unit cell, only short-range motions within the same supercage predominate. Finally, the preferential arrangement of the adsorbate molecules around the extraframework cations are examined and compared with those previously deduced from experimental data.     

From molecular dynamics to hydrodynamics: a novel Galilean invariant thermostat

This paper proposes a novel thermostat applicable to any particle-based dynamic simulation. Each pair of particles is thermostated either (with probability P) with a pairwise Lowe-Andersen thermostat [C. P. Lowe, Europhys. Lett. 47, 145 (1999)] or (with probability 1-P) with a thermostat that is introduced here, which is based on a pairwise interaction similar to the Nosé-Hoover thermostat. When the pairwise Nosé-Hoover thermostat dominates (low P), the liquid has a high diffusion coefficient and low viscosity, but when the Lowe-Andersen thermostat dominates, the diffusion coefficient is low and viscosity is high. This novel Nosé-Hoover-Lowe-Andersen thermostat is Galilean invariant and preserves both total linear and angular momentum of the system, due to the fact that the thermostatic forces between each pair of the particles are pairwise additive and central. We show by simulation that this thermostat also preserves hydrodynamics. For the (noninteracting) ideal gas at P = 0, the diffusion coefficient diverges and viscosity is zero, while for P > 0 it has a finite value. By adjusting probability P, the Schmidt number can be varied by orders of magnitude. The temperature deviation from the required value is at least an order of magnitude smaller than in dissipative particle dynamics (DPD), while the equilibrium properties of the system are very well reproduced. The thermostat is easy to implement and offers a computational efficiency better than (DPD), with better temperature control and greater flexibility in terms of adjusting the diffusion coefficient and viscosity of the simulated system. Applications of this thermostat include all standard molecular dynamic simulations of dense liquids and solids with any type of force field, as well as hydrodynamic simulation of multiphase systems with largely different bulk viscosities, including surface viscosity, and of dilute gases and plasmas.     

Computational electrochemistry: finite element simulation of a disk electrode with ultrasonic acoustic streaming

Computer simulations using the finite element method (FEM) are used to predict the correlation between the transport limiting current (I(lim)) and parameters such as diffusion coefficient, source to electrode separation, source power, and medium viscosity for a sonicated disk electrode in "face on" mode. The fluid dynamics and diffusion layer are modeled directly using FEM and predict that the electrode is uniformly accessible, I(lim) is proportional to the diffusion coefficient to the 2/3 power and I(lim) is proportional to the square root of the source power. Curves are also calculated relating I(lim) to the source to electrode separation and liquid viscosity.     

Microscopic interpretation of the dependence of the contact angle on roughness

The macroscopic contact angle theta(m) of a liquid drop on a rough solid surface in the presence of a gas is calculated microscopically on the basis of a variational minimization of the total potential energy of the drop. Two limiting cases are considered: the liquid penetrates into the space between asperities (Wenzel regime) and the liquid resides on the top of asperities (Cassie-Baxter regime). Long-range as well as short-range interactions between the molecules of liquid, solid, and gas are taken into account. It was also assumed that small portions of insoluble gas are accumulated near the edges of the asperities during the formation of the droplet. The contact angle depends on several parameters involved in the microscopic interactions as well as on the fractions of solid surface between asperities, and of the surface of the asperities themselves, that are in contact with the liquid. It is shown that the theory can explain the nonlinear dependence of cos theta(m) on roughness observed by Krupenkin et al. [Krupenkin, T. N.; Taylor, J. A.; Schneider, T. M.; Yang, S. Langmuir 2004, 20, 3824].(1).     

Simulation analysis of contact angles and retention forces of liquid drops on inclined surfaces

A simulation study of liquid drops on inclined surfaces is performed in order to understand the evolution of drop shapes, contact angles, and retention forces with the tilt angle. The simulations are made by means of a method recently developed for dealing with contact angle hysteresis in the public-domain Surface Evolver software. The results of our simulations are highly dependent on the initial contact angle of the drop. For a drop with an initial contact angle equal to the advancing angle, we obtain results similar to those of experiments in which a drop is placed on a horizontal surface that is slowly tilted. For drops with an initial contact angle equal to the mean between the advancing and the receding contact angles, we recover previous results of finite element studies of drops on inclined surfaces. Comparison with experimental results for molten Sn-Ag-Cu on a tilted Cu substrate shows excellent agreement.     

Thermodiffusion in model nanofluids by molecular dynamics simulations

In this work, a new algorithm is proposed to compute single particle (infinite dilution) thermodiffusion using nonequilibrium molecular dynamics simulations through the estimation of the thermophoretic force that applies on a solute particle. This scheme is shown to provide consistent results for model nanofluids in the liquid state (spherical nonmetallic nanoparticles+Lennard-Jones fluid) where it appears that thermodiffusion amplitude, as well as thermal conductivity, decreases with nanoparticle concentration. Then, by changing the nature of the nanoparticle (size, mass, and internal stiffness) and that of the solvent (quality and viscosity), various trends are exhibited. In all cases, the single particle thermodiffusion is positive, i.e., the nanoparticle tends to migrate toward the cold area. The single particle thermal diffusion coefficient is shown to be independent of the size of the nanoparticle (diameter of 0.8-4 nm), whereas it increases with the quality of the solvent and is inversely proportional to the viscosity of the fluid. In addition, this coefficient is shown to be independent of the mass of the nanoparticle and to increase with the stiffness of the nanoparticle internal bonds. Besides, for these configurations, the mass diffusion coefficient behavior appears to be consistent with a Stokes-Einstein-like law.     

Molecular alignment in a liquid induced by a nonresonant laser field: Molecular dynamics simulation

We carried out molecular dynamics (MD) simulations for a dilute aqueous solution of pyrimidine in order to investigate the mechanisms of field-induced molecular alignment in a liquid phase. An anisotopically polarizable molecule can be aligned in a liquid phase by the interaction with a nonresonant intense laser field. We derived the effective forces induced by a nonresonant field on the basis of the concept of the average of the total potential over one optical cycle. The results of MD simulations show that a pyrimidine molecule is aligned in an aqueous solution by a linearly polarized field of light intensity I approximately 10(13) W/cm2 and wavelength lambda = 800 nm. The temporal behavior of field-induced alignment is adequately reproduced by the solution of the Fokker-Planck equation for a model system in which environmental fluctuations are represented by Gaussian white noise. From this analysis, we have revealed that the time required for alignment in a liquid phase is in the order of the reciprocals of rotational diffusion coefficients of a solute molecule. The degree of alignment is determined by the anisotropy of the polarizability of a molecule, light intensity, and temperature. We also discuss differences between the mechanisms of optical alignment in a gas phase and a liquid phase.     

Ruthenium nanoparticles in ionic liquids: structural and stability effects of polar solutes

Ionic liquids are a stabilizing medium for the in situ synthesis of ruthenium nanoparticles. Herein we show that the addition of molecular polar solutes to the ionic liquid, even in low concentrations, eliminates the role of the ionic liquid 3D structure in controlling the size of ruthenium nanoparticles, and can induce their aggregation. We have performed the synthesis of ruthenium nanoparticles by decomposition of [Ru(COD)(COT)] in 1-butyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide, [C(1)C(4)Im][NTf(2)], under H(2) in the presence of varying amounts of water or 1-octylamine. For water added during the synthesis of metallic nanoparticles, a decrease of the solubility in the ionic liquid was observed, showed by nanoparticles located at the interface between aqueous and ionic phases. When 1-octylamine is present during the synthesis, stable nanoparticles of a constant size are obtained. When 1-octylamine is added after the synthesis, aggregation of the ruthenium nanoparticles is observed. In order to explain these phenomena, we have explored the molecular interactions between the different species using (13)C-NMR and DOSY (Diffusional Order Spectroscopy) experiments, mixing calorimetry, surface tension measurements and molecular simulations. We conclude that the behaviour of the ruthenium nanoparticles in [C(1)C(4)Im][NTf(2)] in the presence of 1-octylamine depends on the interaction between the ligand and the nanoparticles in terms of the energetics but also of the structural arrangement of the amine at the nanoparticle's surface.     

Molecular dynamics simulation of the ionic liquid N-ethyl-N,N-dimethyl-N-(2-methoxyethyl)ammonium bis(trifluoromethanesulfonyl)imide

Thermodynamics, structure, and dynamics of an ionic liquid based on a quaternary ammonium salt with ether side chain, namely, N-ethyl-N,N-dimethyl-N-(2-methoxyethyl)ammonium bis(trifluoromethanesulfonyl)imide, MOENM2E TFSI, are investigated by molecular dynamics (MD) simulations. Average density and configurational energy of simulated MOENM2E TFSI are interpreted with models that take into account empirical ionic volumes. A throughout comparison of the equilibrium structure of MOENM2E TFSI with previous results for the more common ionic liquids based on imidazolium cations is provided. Several time correlation functions are used to reveal the microscopic dynamics of MOENM2E TFSI. Structural relaxation is discussed by the calculation of simultaneous space-time correlation functions. Temperature effects on transport coefficients (diffusion, conductivity, and viscosity) are investigated. The ratio between the actual conductivity and the estimate from ionic diffusion by the Nernst-Einstein equation indicates that correlated motion of neighboring ions in MOENM2E TFSI is similar to imidazolium ionic liquids. In line with experiment, Walden plot of conductivity and viscosity indicates that simulated MOENM2E TFSI should be classified as a poor ionic liquid.     

Molecular simulation of fluid-solid interfaces at nanoscale

The equilibrium states of vapor and liquid coexistent phases in contact with a solid surface are studied at the nanoscale by molecular dynamics simulations for a temperature close to the fluid triple point. The characteristics of the solid-fluid interfaces are determined when the interaction strength between the fluid and the solid varies in order to go from a situation of complete drying to that of complete wetting. From the vapor-liquid density profiles of liquid drops lying on the substrate surface or menisci of liquid films confined in slit pores, the contact angles made by the vapor-liquid interface with the solid are computed. The angle values are similar for the drops and the films. They are also in good qualitative agreement with the estimates obtained through the Young's relation from the surface tensions associated with the vapor-solid, liquid-solid, and vapor-liquid interfaces. However, at this scale, the uncertainties inherent to the angle computation and, to a lesser extent, to size effects seem to preclude that the quantitative agreement between the angle estimates obtained from the interface geometry and calculated from the Young's relation can be better than few degrees.     

Effect of two different size chiral ligand-capped gold nanoparticle dopants on the electro-optic and dielectric dynamics of a ferroelectric liquid crystal mixture

Nanocolloids consisting of a ferroelectric liquid crystal (FLC) doped with different concentrations (0.10 and 0.50 wt.%) of surface treated gold nanoparticles (GNPs) differing in size (1.77, 5.5 nm) are prepared and characterised. The effects of doping on the clearing temperatures as well as electro-optic and dielectric parameters of a FLC mixture are presented. The clearing temperatures remain invariant with doping. A remarkable increase in the spontaneous polarisation is noticed due to the addition of the GNPs with chiral monolayer capping. Tilt angle and switching time, at least in their tendency, become slightly reduced and increased, respectively. Depending on the size of the nanoparticles, surface plasmon resonance is observed to be slightly increased by increasing the surface. In addition, a small change in localised electric field is found upon doping. The increase in the dielectric permittivity and the dielectric strength is observed and attributed to the parallel coupling between the dipoles of functionalised GNPs, induced by external electric field, and the vector of the spontaneous polarisation of the FLC matrix. A decrease in relaxation frequency is observed. A substantial increment of one order in the dc conductivity is also observed for the nanocolloids.     

A comparison of implicit- and explicit-solvent simulations of self-assembly in block copolymer and solute systems

We have developed explicit- and implicit-solvent models for the flash nanoprecipitation process, which involves rapid coprecipitation of block copolymers and solutes by changing solvent quality. The explicit-solvent model uses the dissipative particle dynamics (DPD) method and the implicit-solvent model uses the Brownian dynamics (BD) method. Each of the two models was parameterized to match key properties of the diblock copolymer (specifically, critical micelle concentration, diffusion coefficient, polystyrene melt density, and polyethylene glycol radius of gyration) and the hydrophobic solute (aqueous solubility, diffusion coefficient, and solid density). The models were simulated in the limit of instantaneous mixing of solvent with antisolvent. Despite the significant differences in the potentials employed in the implicit- and explicit-solvent models, the polymer-stabilized nanoparticles formed in both sets of simulations are similar in size and structure; however, the dynamic evolution of the two simulations is quite different. Nanoparticles in the BD simulations have diffusion coefficients that follow Rouse behavior (D ∝ M(-1)), whereas those in the DPD simulations have diffusion coefficients that are close to the values predicted by the Stokes-Einstein relation (D ∝ R(-1)). As the nanoparticles become larger, the discrepancy between diffusion coefficients grows. As a consequence, BD simulations produce increasingly slower aggregation dynamics with respect to real time and result in an unphysical evolution of the nanoparticle size distribution. Surface area per polymer of the stable explicit-solvent nanoparticles agrees well with experimental values, whereas the implicit-solvent nanoparticles are stable when the surface area per particle is roughly two to four times larger. We conclude that implicit-solvent models may produce questionable results when simulating nonequilibrium processes in which hydrodynamics play a critical role.     

Dynamics of dewetting at the nanoscale using molecular dynamics

Large-scale molecular dynamics simulations are used to model the dewetting of solid surfaces by partially wetting thin liquid films. Two levels of solid-liquid interaction are considered that give rise to large equilibrium contact angles. The initial length and thickness of the films are varied over a wide range at the nanoscale. Spontaneous dewetting is initiated by removing a band of molecules either from each end of the film or from its center. As observed experimentally and in previous simulations, the films recede at an initially constant speed, creating a growing rim of liquid with a constant receding dynamic contact angle. Consistent with the current understanding of wetting dynamics, film recession is faster on the more poorly wetted surface to an extent that cannot be explained solely by the increase in the surface tension driving force. In addition, the rates of recession of the thinnest films are found to increase with decreasing film thickness. These new results imply not only that the mobility of the liquid molecules adjacent to the solid increases with decreasing solid-liquid interactions, but also that the mobility adjacent to the free surface of the film is higher than in the bulk, so that the effective viscosity of the film decreases with thickness.