A diffusion model for the fluids confined in micropores

The self-diffusion coefficients were calculated by molecular dynamics simulations and the effects of pore width, temperature, and fluid density on diffusion behavior of simple fluid argon and polar fluid water confined in micropores were analyzed and studied. A mathematical model describing diffusion behavior of fluids confined in micropores was proposed from the theories of molecular dynamics and molecular kinematics, and validated on the basis of the simulation results at various conditions. The model indicates that the diffusion coefficient is proportional to the square root of the pore width and to the temperature divided by the density squared. It is applicable to either liquid or gas states and only two parameters are required.     

Molecular dynamics simulation of the diffusion of nanoconfined fluids

A new molecular dynamics simulation technique for simulating fluids in confinement [H. Eslami, F. Mozaffari, J. Moghadasi, F. Müller-Plathe, J. Chem. Phys. 129 (2008) 194702] is employed to simulate the diffusion coefficient of nanoconfined Lennard-Jones fluid. The diffusing fluid is liquid Ar and the confining surfaces are solid Ar fcc (100) surfaces, which are kept frozen during the simulation. In this simulation just the fluid in confinement is simulated at a constant temperature and a constant parallel component of pressure, which is assumed to be equal to the bulk pressure. It is shown that the calculated parallel (to the surfaces) component of the diffusion coefficients depends on the distance between the surfaces (pore size) and shows oscillatory behavior with respect to the intersurface separations. Our results show that on formation of well-organized layers between the surfaces, the parallel diffusion coefficients decrease considerably with respect to the bulk fluid. The effect of pressure on the parallel diffusion coefficients has also been studied. Better organized layers, and hence, lower diffusion coefficients are observed with increasing the pressure.     

Density of He adsorbed in micropores at 4.2K

The density of He adsorbed in the cylindrical micropores of zeolites NaY and KL has been determined by He adsorption at 4.2K. He adsorption isotherms were then compared with N₂ adsorption isotherms at 77K. Crystallographic considerations of the micropore volumes gave the density of the He adsorbed layer, which is necessary for assessment of ultramicroporosity of less-crystalline microporous solids, such as activated carbons. The determined density of He adsorbed in the cylindrical micropores of the zeolite was in the range 0.22 to 0.26 gml^–1, greater than that of He adsorbed on a flat surface (0.202 gml^–1). A value for the density of He between 0.20 to 0.22 gml^–1 is recommended for evaluation of ultramicroporosity of a slit-shaped microporous system such as activated carbon.     

Thermodynamics of d-dimensional hard sphere fluids confined to micropores

We derive an analytical expression of the second virial coefficient of d-dimensional hard sphere fluids confined to slit pores by applying Speedy and Reiss' interpretation of cavity space. We confirm that this coefficient is identical to the one obtained from the Mayer cluster expansion up to second order with respect to fugacity. The key step of both approaches is to evaluate either the surface area or the volume of the d-dimensional exclusion sphere confined to a slit pore. We, further, present an analytical form of thermodynamic functions such as entropy and pressure tensor as a function of the size of the slit pore. Molecular dynamics simulations are performed for d = 2 and d = 3, and the results are compared with analytically obtained equations of state. They agree satisfactorily in the low density regime, and, for given density, the agreement of the results becomes excellent as the width of the slit pore gets smaller, because the higher order virial coefficients become unimportant.     

Effect of micropores on the effective diffusion coefficient

The effective diffusion coefficient for catalysts differing in their porous structure has been derived from experimental data on H₂S conversion in the Claus reaction. The effective diffusion coefficient increases under conditions of catalyst deactivation due to sulfur condensation in micropores. A mathematical model is suggested to describe the micropore effect on the effective diffusion coefficient.     

Measure of accuracy for multicanonical molecular-dynamics simulation

A measure of the flatness of the energy probability distribution for multicanonical molecular-dynamics (MMD) simulation is presented. It (the flatness measure) can be introduced by a slight change in the renewing scheme of the MMD potential energy. Our proposed measure is applied to liquid Ar with a Lennard-Jones potential system in order to investigate the influence of flatness on the simulation results such as internal energy and specific heat at constant volume. We find that the accuracy of MMD simulation is influenced not only by the flatness of the energy probability distribution but also by the width of the energy region that is accessible during the MMD simulation.     

Classical molecular-dynamics simulation of the hydroxyl radical in water

We have studied the hydration and diffusion of the hydroxyl radical OH0 in water using classical molecular dynamics. We report the atomic radial distribution functions, hydrogen-bond distributions, angular distribution functions, and lifetimes of the hydration structures. The most frequent hydration structure in the OH0 has one water molecule bound to the OH0 oxygen (57% of the time), and one water molecule bound to the OH0 hydrogen (88% of the time). In the hydrogen bonds between the OH0 and the water that surrounds it the OH0 acts mainly as proton donor. These hydrogen bonds take place in a low percentage, indicating little adaptability of the molecule to the structure of the solvent. All hydration structures of the OH0 have shorter lifetimes than those corresponding to the hydration structures of the water molecule. The value of the diffusion coefficient of the OH0 obtained from the simulation was 7.1x10(-9) m2 s(-1), which is higher than those of the water and the OH-.     

Test of classical nucleation theory via molecular-dynamics simulation

A direct test of classical nucleation theory (CNT) is made using molecular-dynamics simulations. The relation between critical nucleus size and undercooling temperature is extracted and the result yields the solid-liquid interfacial energy. It is shown that the CNT, within the assumptions made for spherical nucleus in supercooled liquid, is valid in the critical regime of nucleation for a large range of undercooling and nucleus size.     

Density functional theory based molecular-dynamics study of aqueous iodide solvation

We study the solvation of iodide in water using density functional theory based molecular-dynamics simulations. Detailed analysis of the structural and dynamical properties of the first solvation shell is presented, showing a disruptive influence of the ion on the local water structure. Iodide-water hydrogen bonding is weak, compared to water-water hydrogen bonds. This effective repulsive ion-water interaction leads to the formation of a quite unstructured solvation shell. The dynamics of water molecules surrounding the iodide is relatively fast. The intramolecular structural and electronical properties of water molecules around the ion are not affected.     

Density functional theory based molecular-dynamics study of aqueous fluoride solvation

We use density functional theory based molecular-dynamics simulations to study the aqueous solvation of the fluoride anion. Our studies are focused on the first solvation shell and have resulted in detailed information on its structural and dynamical properties. The fluoride ion leads to the formation of a rigid solvation shell, qualitatively consistent with simulation and experimental studies, classifying fluoride as a "structure making" particle. However, quantitatively we find the solvation shell to be less structured and more mobile than predicted from empirical force-field simulation. The influence on the intramolecular electronical and structural properties of water is minimal, as observed for other halogens. We propose two distinct mechanisms for the exchange of bulk and first solvation shell water molecules.     

Large-scale molecular-dynamics simulation of nanoscale hydrophobic interaction and nanobubble formation

We performed large-scale molecular-dynamics simulation of nanoscale hydrophobic interaction manifested by the formation of nanobubble between nanometer-sized hydrophobic clusters at constrained equilibrium. Particular attention is placed on the tendency of formation and stability of nanobubbles in between model nanoassemblies which are composed of hydrophobic clusters (or patches) embedded in a hydrophilic substrate. On the basis of physical behavior of nanobubble formation, we observed a change from short-range molecular hydrophobic interaction to midrange nanoscopic interaction when the length scale of hydrophobe approaches to about 1 nm. We investigated the behavior of nanobubble formation with several different patterns of nonpolar-site distribution on the nanoassemblies but always keeping a constant ratio of nonpolar to polar monomer sites. Dynamical properties of confined water molecules in between nanoassemblies are also calculated.     

Adsorption of gas-water binary systems in carbon micropores: Computer simulation

The adsorption of gas-water mixture in micropores of carbon materials at 298 K has been studied using computer simulation. Methane, nitrogen, ammonia, carbon dioxide, and hydrogen sulfide were considered as gas components. In the grand canonical ensemble Monte-Carlo simulation of adsorption, the displacement of a gas component from a pore as a result of the formation of water microclusters was observed for all systems studied. Cluster growth conditions on graphite-like and activated surfaces differ significantly. The comparative stability of adsorbed gas-water mixtures has been determined for all gases.     

Ab initio molecular-dynamics simulation of aqueous proton solvation and transport revisited

The solvation and transport of the hydrated excess proton is studied using the Car-Parrinello molecular-dynamics (CPMD) simulation method. The simulations were performed using BLYP and HCTH gradient-corrected exchange-correlation energy functionals. The fictitious electronic mass was chosen to be small enough so that the underlying water structural and dynamical properties were converged with respect to this important CPMD simulation parameter. An unphysical overstructuring of liquid water in the CPMD simulations using the BLYP functional resulted in the formation of long-lived hydrogen-bonding structures involving the excess proton and a particular (special) water oxygen. The excess proton was observed to be attracted to the special oxygen through the entire length of the BLYP CPMD simulations. Consequently, the excess proton diffusion was limited by the mobility of the special oxygen in the slowly diffusing water network and, in turn, the excess proton self-diffusion coefficient was found to be significantly below the experimental value. On the other hand, the structural properties of liquid water in the HCTH CPMD simulation were seen to be in better agreement with experiment, although the water and excess proton diffusions were still well below the experimental value.     

Phase behavior of attractive and repulsive ramp fluids: integral equation and computer simulation studies

Using computer simulations and a thermodynamically self-consistent integral equation we investigate the phase behavior and thermodynamic anomalies of a fluid composed of spherical particles interacting via a two-scale ramp potential (a hard core plus a repulsive and an attractive ramp) and the corresponding purely repulsive model. Both simulation and integral equation results predict a liquid-liquid demixing when attractive forces are present, in addition to a gas-liquid transition. Furthermore, a fluid-solid transition emerges in the neighborhood of the liquid-liquid transition region, leading to a phase diagram with a somewhat complicated topology. This solidification at moderate densities is also present in the repulsive ramp fluid, but in this case inhibits the fluid-fluid separation.     

Molecular simulation of the phase behavior of fluids and fluid mixtures using the synthetic method

A new molecular simulation procedure is reported for determining the phase behavior of fluids and fluid mixtures, which closely follows the experimental synthetic method. The simulation procedure can be implemented using Monte Calro or molecular dynamics in either the microcanonical or canonical statistical ensembles. Microcanonical molecular dynamics simulations are reported for the phase behavior of both the pure Lennard-Jones fluid and a Lennard-Jones mixture. The vapor pressures for the pure fluid are in good agreement with Monte Carlo Gibbs ensemble and Gibbs-Duhem calculations. The Lennard-Jones mixture is composed of equal size particles, with dissimilar energy parameters (?(2)∕?(1) = 1∕2, ?(12)∕?(1) = 1∕2). The binary Lennard-Jones mixture exhibits liquid-liquid equilibria at high pressures and the simulation procedure allows us to estimate the coordinates of the high-pressure branch of the critical curve.     

The behavior of fluids near solutes: a density functional theory and computer simulation study

The density distribution of solvent near a solute particle is studied using density functional theory and Monte Carlo simulation. The fluid atoms interact with each other via a hard sphere plus Yukawa potential, and interact with the solute via a hard sphere potential. For small solute sizes, the solvent displays liquidlike ordering near the particle. When the solute become larger, a drying transition is observed at state points near the coexistence conditions of the solvent. These predictions are similar to those of a recent theory for the hydrophobic effect by Lum, Chandler, and Weeks [J. Phys. Chem. 103, 4570 (1999)], although a comparison with simulations shows that the theory of this work is quantitatively more accurate. The connection between density functional methods and the LCW approach is also established.     

Non-Newtonian behavior in simple fluids

Using nonequilibrium molecular dynamics simulations, we study the non-Newtonian rheology of a microscopic sample of simple fluid. The calculations were performed using a configurational thermostat which unlike previous nonequilibrium molecular dynamics or nonequilibrium Brownian dynamics methods does not exert any additional constraint on the flow profile. Our findings are in agreement with experimental results on concentrated "hard sphere"-like colloidal suspensions. We observe: (i) a shear thickening regime under steady shear; (ii) a strain thickening regime under oscillatory shear at low frequencies; and (iii) shear-induced ordering under oscillatory shear at higher frequencies. These results significantly differ from previous simulation results which showed systematically a strong ordering for all frequencies. They also indicate that shear thickening can occur even in the absence of a solvent.     

A molecular-dynamics simulation study of solvent-induced repulsion between C60 fullerenes in water

Molecular-dynamics simulations of a single C(60) fullerene and pairs of C(60) fullerenes in aqueous solution have been performed for the purpose of obtaining improved understanding of the nature of solvent-induced interactions between C(60) fullerenes in water. Our simulations reveal repulsive solvent-induced interactions between two C(60) fullerenes in aqueous solution in contrast to the associative effects observed for conventional nonpolar solutes. A decomposition of the solvent-induced potential of mean force between fullerenes into entropy and energy (enthalpy) contributions reveals that the water-induced repulsion between fullerenes is energetic in origin, contrasting strongly to entropy-driven association observed for conventional nonpolar solutes. The dominance of energy in the solvent-induced interactions between C(60) fullerenes arises primarily from the high atomic density of the C(60) molecule, resulting in strong C(60)-water van der Waals attraction that is reduced upon association of the fullerenes. The water-induced repulsion is found to decrease with increasing temperature due largely to an increasing contribution from a relatively weak entropy-driven association.     

Geometric structure of a molecule of Epitalon tetrapeptide: A molecular-dynamics simulation

Variation of the geometric parameters of a molecule of Epitalon tetrapeptide (Ala-Glu-Asp-Gly) over a period of 1500 ps was simulated by the method of molecular dynamics using AMBER force field. The structure of the molecule is stabilized by two salt bridges formed by the N-terminal nitrogen atom and oxygen atoms of Asp and Glu side chains. The biological effect of Epitalon was attributed to formation of salt or hydrogen bonds involving one or several ionizable functional groups of the molecules.     

Density functional theory of fluids in the isothermal-isobaric ensemble

We present a density functional theory for inhomogeneous fluids at constant external pressure. The theory is formulated for a volume-dependent density, n(r,V), defined as the conjugate variable of a generalized external potential, nu(r,V), that conveys the information on the pressure. An exact expression for the isothermal-isobaric free-energy density functional is obtained in terms of the corresponding canonical ensemble functional. As an application we consider a hard-sphere system in a spherical pore with fluctuating radius. In general we obtain very good agreement with simulation. However, in some situations a peak develops in the center of the cavity and the agreement between theory and simulation becomes worse. This happens for systems where the number of particles is close to the magic numbers N=13, 55, and 147.