Molecular simulation of interaction between passivated gold nanoparticles in supercritical CO2

Molecular dynamics simulations have been performed to study the potential of mean force (PMF) between passivated gold nanoparticles (NPs) in supercritical CO(2) (scCO(2)). The nanoparticle model consists of a 140 atom gold nanocore and a surface self-assembled monolayer, in which two kinds of fluorinated alkanethiols were considered. The molecular origin of the thermodynamics interaction and the solvation effect has been comprehensively studied. The simulation results demonstrate that increasing the solvent density and ligand length can enhance the repulsive feature of the free energy between the passivated Au nanoparticles in scCO(2), which is in good agreement with previous experimental results. The interaction forces between the two passivated NPs have been decomposed to reveal various contributions to the free energy. It was revealed that the interaction between capping ligands and the interaction between the capping ligands and scCO(2) solvent molecules cooperatively determine the total PMF. A thermodynamic entropy-energy analysis for each PMF contribution was used to explain the density dependence of PMF in scCO(2) fluid. Our simulation study is expected to provide a novel microscopic understanding of the effect of scCO(2) solvent on the interaction between passivated Au nanoparticles, which is helpful to the dispersion and preparation of functional metal nanoparticles in supercritical fluids.     

Adsorption and self-assembly of surfactant/supercritical CO2 systems in confined pores: a molecular dynamics simulation

A coarse-grained molecular dynamics simulation has been carried out to study the adsorption and self-organization for a model surfactant/supercritical CO2 system confined in the slit-shape nanopores with amorphous silica-like surfaces. The solid surfaces were designed to be CO2-philic and CO2-phobic, respectively. For the CO2-philic surface, obviously surface adsorption is observed for the surfactant molecules. The various energy profiles were used to monitor the lengthy dynamics process of the adsorption and self-assembly for surfactant micelles or monomers in the confined spaces. The equilibrium properties, including the morphologies and micelle-size distributions of absorbed surfactants, were evaluated based on the equilibrium trajectory data. The interaction between the surfactant and the surface produces an obvious effect on the dynamics rate of surfactant adsorption and aggregation, as well as the final self-assembly equilibrium structures of the adsorbed surfactants. However, for the CO2-phobic surfaces, there are scarcely adsorption layers of surfactant molecules, meaning that the CO2-phobic surface repels the surfactant molecules. It seems to conclude that the CO2 solvent depletion near the interfaces determines the surface repellence to the surfactant molecules. The effect of the CO2-phobic surface confinement on the surfactant micelle structure in the supercritical CO2 has also been discussed. In summary, this study on the microscopic behaviors of surfactant/Sc-CO2 in confined pores will help to shed light on the surfactant self-assembly from the Sc-CO2 fluid phase onto solid surfaces and nanoporous media.     

Solvation shell dynamics studied by molecular dynamics simulation in relation to the translational and rotational dynamics of supercritical water and benzene

The solvation shell dynamics of supercritical water is analyzed by molecular dynamics simulation with emphasis on its relationship to the translational and rotational dynamics. The relaxation times of the solvation number (tau S), the velocity autocorrelation function (tau D), the angular momentum correlation function (tau J), and the second-order reorientational correlation function (tau 2R) are studied at a supercritical temperature of 400 degrees C over a wide density region of 0.01-1.5 g cm(-3). The relaxation times are decomposed into those conditioned by the solvation number n, and the effect of the short-ranged structure is examined in terms of its probability Pn of occurrence. In the low to medium-density range of 0.01-0.4 g cm(-3), the time scales of water dynamics are in the following sequence: tau D>tau S approximately or > tau J approximately or > tau 2R. This means that the rotation in supercritical water is of the "in-shell" type while the translational diffusion is not. The comparison to supercritical benzene is also performed and the effect of hydrogen bonding is examined. The water diffusion is not of the in-shell type up to the ambient density of 1.0 g cm(-3), which corresponds to the absence of the transition from the collision to the Brownian picture, whereas such transition is present in the case of benzene. The absence of the transition in water comes from the fast reorganization of the hydrogen bonds and the enhanced mobility of the solvation shell in supercritical conditions.     

Solvation of calcium ions in methanol-water mixtures: molecular dynamics simulation

Molecular dynamics simulations of CaCl2 solutions in water and methanol-water mixtures, with methanol concentrations of 5, 10, 50, and 90 mol %, at room temperature, have been performed. The methanol and water molecules have been modeled as flexible three-site bodies. Solvation of the calcium ions has been discussed on the basis of the radial and angular distribution functions, the orientation of the solvent molecules, and their geometrical arrangement in the coordination shells. Analysis of the H-bonds of the solvent molecules coordinated by Ca2+ has been done. Residence time of the solvent molecules in the coordination shell has been calculated. The preferential hydration of the calcium ions has been found over the whole range of the mixture composition. The water concentration in the first and second coordination shells of Ca2+ significantly exceeds the water content in the solution, despite the very similar interaction energy of the calcium ion with water and methanol. In aqueous solution and methanol-water mixtures, the first coordination shell of Ca2+ is irregular and long-living. The solvent molecules prefer the anti-dipole arrangement, but, in aqueous solutions and water-rich mixtures, the water molecules in the primary shell have only one H-bonded neighbor.     

Molecular dynamics simulation of a reverse micelle self assembly in supercritical CO2

In this communication we report on molecular dynamics computer simulations of self-assembly of reverse micelles in supercritical carbon dioxide. The reverse micelles contain perfluoropolyether ammonium carboxylate surfactants and an aqueous core. We observed a quick self-assembly of these micelles over time periods of approximately 5 ns, irrespective of initial conditions. In most cases, the self-assembled perfluorinated reverse micelles have a nice spherical shape and properties consistent with experiments. When the fluorinated surfactant is replaced by its hydrogenated analogue, the assembled aggregate contains a region of direct contact between water and carbon dioxide, indicating that hydrogenated surfactant is not a good agent for creation of microemulsions in water/carbon dioxide mixtures.     

A molecular dynamics simulation of the adsorption of water molecules surrounding an Au nanoparticle

This study uses molecular dynamics simulations performed in a parallel computing environment to investigate the adsorption of water molecules surrounding Au nanoparticles of various sizes. An observation of the oxygen and hydrogen atom distributions reveals that the adsorption of the water molecules creates two shell-like formations of water in close vicinity to the Au nanoparticle surface. These shell-like formations are found to be more pronounced around smaller Au nanoparticles. The rearrangement of water molecules in this region reduces the local hydrogen bond strength to below that which is observed in the bulk region. Finally, the simulation results indicate that the absolute value of the interaction energy between the water molecules and the Au nanoparticle is reduced when the water molecules surround a nanoparticle of larger diameter. This observation implies that a stronger adsorption effect exists between smaller Au nanoparticles and water molecules. Hence, the value of the adsorption constant increases for smaller Au nanoparticles.     

Structural features of binary mixtures of supercritical CO2 with polar entrainers by molecular dynamics simulation

Computer simulations of supercritical carbon dioxide and its mixtures with polar cosolvents: water, methanol, and ethanol (concentration, 0.125 mole fractions) at T = 318 K and ρ = 0.7 g/cm³ are performed. Atom-atom radial distribution functions are calculated by classical molecular dynamics, while the probability distributions of relative orientation of CO₂ molecules in the first and second coordination spheres describing the geometry of the nearest environment of CO₂ molecules and the trajectories of cosolvent molecules are found using Car-Parrinello molecular dynamics. Based on the latter, the conclusions regarding structure and interactions of polar entrainers in their mixtures with supercritical CO₂ are made. It is shown that the microstructure of carbon dioxide varies only slightly upon the introduction of cosolvents.     

Solvation of phenylglycine- and leucine-derived chiral stationary phases: molecular dynamics simulation study

A theoretical study of the solvation of ( R)- N-(3,5-dinitrobenzoyl)phenylglycine- and ( R)- N-(3,5-dinitrobenzoyl)leucine-derived chiral stationary phases (CSPs) is presented. Semiflexible models of the chiral selectors are prepared from B3LYP/6-311G** calculations, and these are used in the molecular dynamics simulations of the corresponding interface. The chiral interface is examined for four solvents: 100% hexane, 90:10 hexane:2-propanol, 80:20 hexane:2-propanol, and 100% 2-propanol. Despite the similarities between phenylglycine and leucine, the interfaces are distinct both in terms of the selector orientations at the surface and in the number of hydrogen bonds formed with 2-propanol. We also find that an increase in alcohol concentration alters the preferred orientations of the selectors.     

Local density inhomogeneities and dynamics in supercritical water: A molecular dynamics simulation approach

Molecular dynamics atomistic simulations in the canonical ensemble (NVT-MD) have been used to investigate the "Local Density Inhomogeneities and their Dynamics" in pure supercritical water. The simulations were carried out along a near-critical isotherm (Tr = T/Tc = 1.03) and for a wide range of densities below and above the critical one (0.2 rho(c) - 2.0 rho(c)). The results obtained reveal the existence of significant local density augmentation effects, which are found to be sufficiently larger in comparison to those reported for nonassociated fluids. The time evolution of the local density distribution around each molecule was studied in terms of the appropriate time correlation functions C(Delta)rhol(t). It is found that the shape of these functions changes significantly by increasing the density of the fluid. Finally, the local density reorganization times for the first and second coordination shell derived from these correlations exhibit a decreasing behavior by increasing the density of the system, signifying the density effect upon the dynamics of the local environment around each molecule.     

Systematic coarse-graining of nanoparticle interactions in molecular dynamics simulation

A recently developed multiscale coarse-graining procedure [Izvekov, S.; Voth, G. A. J. Phys. Chem. B 2005, 109, 2469] is extended to derive coarse-grained models for nanoparticles. The methodology is applied to C(60) and to carbonaceous nanoparticles produced in combustion environments. The coarse-graining of the interparticle force field is accomplished applying a force-matching procedure to data obtained from trajectories and forces from all-atom MD simulations. The CG models are shown to reproduce accurately the structural properties of the nanoparticle systems studied, while allowing for MD simulations of much larger self-assembled nanoparticle systems.     

Molecular simulation of swelling and interlayer structure for organoclay in supercritical CO(2)

In this work, Monte Carlo simulations have been carried out to investigate the swelling stability and interlayer structures of alkylammonium-modified montmorillonite both in vacuum and in supercritical CO(2) (scCO(2)) fluid. In the vacuum (dry) condition, the stable spacing for this kind of organoclay was determined based on the energy minimum. In the stable spacing, the corresponding interlayer structure of dry organoclay is the monolayer arrangement with the intercalated surfactant chains lying parallel to the silicate surface. In scCO(2) fluid medium, the normal pressures within the organoclay gallery and the swelling free energy have been obtained from Gibbs ensemble Monte Carlo simulation. The mechanically and thermodynamically stable spacings of the organoclay have been determined. As compared with the case in vacuum, the simulation shows that the swelling of the organoclay is thermodynamically favorable in the environment of scCO(2) fluid. The interlayer structure and conformation have been used to analyze the mechanism of swelling. The headgroups of surfactant cations are distributed close to the clay surfaces. The presence of CO(2) molecules within the clay gallery can cause a specific steric arrangement of the long-chain alkylammonium cations.     

Solvation of xyloglucan in water/alcohol systems by molecular dynamics simulation

Xyloglucan in water solution turns into a gel with addition of alcohol such as methanol and ethanol. In regard to this phenomenon, we investigated the adhesive property of alcohol to xyloglucan and proposed the mechanism of the gelation by molecular dynamics (MD) simulation of a xyloglucan in water, water/methanol, and water/ethanol solution for 10 ns. The alcohol molecules showed its adhesive property to the xyloglucan and made the swelling-shrinking motion of the xyloglucan slow. Alcohol molecules solvated to the xyloglucan mainly in hydrophobic way so as to fill the void of water hydration shell, resulting in reformation of the hydrogen-bond network of water molecules around the solute. We also found that alcohol molecules have strong tendency to hydrogen-bond on xylose O3 in xyloglucan. According to these results, we proposed the gelation mechanism of xyloglucan in water/alcohol solution.     

Hydrogen-bonded clusters and solvate structures in the supercritical CO2-water-o-hydroxybenzoic acid system: the car-parrinello molecular dynamics

Hydrogen-bonded clusters and solvate structures formed by o-hydroxybenzoic acid (o-HBA) and water in supercritical CO₂ were studied (T = 318 K, 348 K, ρ = 0.7 g/cm³). The atom-atom radial distribution functions, coordination numbers, average numbers of hydrogen bonds for individual atomic groups, and power spectrum were calculated by the Car-Parrinello molecular dynamics. Despite the high polarity of the cosolvent, the hydroxyl group of o-HBA predominantly forms intramolecular hydrogen bond, while hydrogen bonds with water involve only the atoms of carboxyl groups. The temperature effect on the stability of these bonds showed itself in different ways. The intermolecular interactions of o-HBA with carbon dioxide were found to be weaker than those with water. It was established that the Lewis acid-Lewis base interactions between CO₂ and the hydroxyl group of the solute increase with increasing temperature. Instantaneous configurations illustrating the temperature effects on the molecular structures were obtained.     

Density dependence of the entropy and the solvation shell structure in supercritical water via molecular dynamics simulation

We perform molecular dynamics simulations of supercritical water (SCW) with a wide range of densities along a near critical isotherm using the simple point charge extended (SPC/E) pair potential in order to study the entropy and the solvation shell structure around a central water molecule. It is shown that both the translational and orientational two-particle correlation entropy terms can serve as the metrics of the translational-orientational structural orders in water and it is revealed that the translational structural order is very sensitive to the density variation in the gas-like and liquid-like region, while the orientational structural order is much more dependent upon compression in the medium-density SCW region. The comparison of the magnitudes of the full thermodynamic excess entropy and two-particle correlation entropy confirms the recent findings that the many-body terms other than two-body ones also present significant and non-neglectable contributions to the full excess entropy for the highly anomalous fluids like water. The analysis of entropy terms as a function of intermolecular distance and the orientational distribution functions as well as the three-dimensional spatial distribution functions indicate that the structural order occurs only in a much more diffused first solvation shell due to the elongated hydrogen bonds under supercritical conditions. It is revealed that no obvious second or higher neighbor shells occur in SCW, in contrast with the feature of normal liquid water that the anomalous decrease of translational order upon compression occurs mainly in the second shell.     

IBIsCO: a molecular dynamics simulation package for coarse-grained simulation

IBIsCO is a parallel molecular dynamics simulation package developed specially for coarse-grained simulations with numerical potentials derived by the iterative Boltzmann inversion (IBI) method (Reith et al., J Comput Chem 2003, 24, 1624). In addition to common features of molecular dynamics programs, the techniques of dissipative particle dynamics (Groot and Warren, J Chem Phys 1997, 107, 4423) and Lowe-Andersen dynamics (Lowe, Europhys Lett 1999, 47, 145) are implemented, which can be used both as thermostats and as sources of friction to compensate the loss of degrees of freedom by coarse-graining. The reverse nonequilibrium molecular dynamics simulation method (Müller-Plathe, Phys Rev E 1999, 59, 4894) for the calculation of viscosities is also implemented. Details of the algorithms, functionalities, implementation, user interfaces, and file formats are described. The code is parallelized using PE_MPI on PowerPC architecture. The execution time scales satisfactorily with the number of processors.     

Electron transport dynamics in a room-temperature Au nanoparticle molten salt

A room-temperature Au38 nanoparticle polyether melt has been prepared by exchanging poly(ethylene glycol) (PEG) thiolate ligands, HS-C6-PEG163, into the organic protecting monolayer of Au38(PhC2)24 nanoparticles. Spectral and electrochemical properties verify that the Au38 core size is preserved during the exchange. Adding LiClO4 electrolyte, free PEG plasticizer, and/or partitioned CO2 leads to an ionically conductive nanoparticle melt, on which voltammetric, chronoamperometric, and impedance measurements have been made, respectively, of the rates of electron and ion transport in the melt. Electron transport occurs by electron self-exchange reactions, or electron hopping, between diffusively relatively immobile Au38(0) and Au38(1+) nanoparticles. The rates of physical diffusion of electrolyte ions (diffusion coefficients DCION) are obtained from ionic conductivities. The measured rates of electron and of electrolyte ion transport are very similar, as are their thermal activation energy barriers, observations that are consistent with a recently introduced ion atmosphere relaxation model describing control of electron transfer in semisolid ion and electron-conductive media. The model has been previously demonstrated using a variety of metal complex polyether melts; the present results extend it to electron transfers between Au nanoparticles. In ion atmosphere relaxation control, measured rates and energy barriers for electron transfer are not intrinsic values but are instead characteristic of competition between back-electron transfer caused by a Coulombic disequilibrium resulting from an electron transfer and relaxation of counterions around donor-acceptor reaction partners so as to reachieve local electroneutrality.     

Local density augmentation and dynamic properties of hydrogen-and non-hydrogen-bonded supercritical fluids: a molecular dynamics study

The local density inhomogeneities in neat supercritical fluids were investigated via canonical molecular dynamics simulations. The selected systems under investigation were the polar and hydrogen-bonded fluid methanol as well as the quadrupolar non-hydrogen-bonded carbon dioxide one. Effective local densities, local density augmentation, and enhancement factors were calculated at state points along an isotherm close to the critical temperature of each system (T(r)=1.03). The results obtained reveal strong influence of the polarity and hydrogen bonding upon the intensity of the local density augmentation. It is found that this effect is sufficiently larger in the case of the polar and associated methanol in comparison to those predicted for carbon dioxide. For both fluids the local density augmentation values are maximized in the bulk density region near 0.7rho(c), a result that is in agreement with experiment. In addition, the local density dynamics of each fluid were investigated in terms of the appropriate time correlation functions. The behavior of these functions reveals that the bulk density dependence of the local density reorganization times is very sensitive to the specific intermolecular interactions and to the size of the local region. Also, the estimated local density reorganization time as a function of bulk density of each fluid was further analyzed and successfully related to two different time-scale relaxation mechanisms. Finally, the results obtained indicate a possible relationship between the single-molecule reorientational dynamics and the local density reorganization ones.     

Solvent density effects on the solvation behavior and configurational structure of bare and passivated 38-atom gold nanoparticle in supercritical ethane

In exploring the effects of solvent density on the mode and the degree of solvation of the bare and passivated 38-atom gold particle in supercritical ethane, we have extended the molecular dynamics simulations of the system, reported previously,(34) to cover a range of isotherms in the T > T(c) regime, where T(c) is the critical temperature of the solvent. Consonant with our previous observations, the modes of solvation of the bare and the passivated particle, deduced from the radial distribution of the solvent about the metal core center of mass, are found to be vastly different from each other at all solvent densities: while the molecules solvating the bare particle form a well-defined, two-region layer around it, those solvating the passivated particle are loosely dispersed in the passivating layer. For the bare particle, the degree of solvation (vartheta) as a function of solvent density passes through a maximum occurring in the close vicinity of the critical point, consistent with our previous results and in agreement with Debenedetti's theoretical analysis,(22,23) which predicts a solvation enhancement effect in the critical region for systems where the unlike solvent/solute interaction is much stronger than the solvent/solvent interaction. Taking the degree of solvation (vartheta) as a measure of solvent quality, we have investigated how the solvent quality would vary along the solvent-density isotherms. In the solvent-density regime rho > rho(c), the solvent quality is found to be a decreasing function of the density as a result of progressive dominance of the excluded volume effect over the attractive particle/solvent interactions. The particle/solvent affinity is greatly reduced in the presence of the passivating layer, resulting in considerable shrinkage of the good-solvent-quality domain in the supercritical regime. The solvent environment and the presence of the passivating chains produce significant disorder in the equilibrium structure assumed by the nanoparticle core.