Interfacial colloidal sedimentation equilibrium. II. Closure-based density functional theory

In Part I [R. E. Beckham and M. A. Bevan, J. Chem. Phys. 127, 164708 (2007)], results were presented for the sedimentation equilibrium of concentrated colloidal dispersions using confocal scanning laser microscopy experiments, Monte Carlo (MC) simulations, and a local density approximation perturbation theory. In this paper, we extended the modeling effort on those systems to include nonlocal density functional theory (DFT), which is capable of predicting the microstructure of the sediment at length scales comparable to the colloidal particle dimension. Specifically, we use a closure-based DFT formulation to predict interfacial colloidal sedimentation equilibrium density profiles. The colloid-colloid and colloid-surface interactions were modeled with DLVO screened electrostatic potentials using parameters taken directly from the experimental work. The DFT profiles were compared to the experimental and MC results from Part I. Good agreement was found for relatively dilute interfacial colloidal fluids, but agreement was less satisfactory as interfacial layering became more pronounced for conditions approaching the onset of interfacial crystallization. We also applied DFT in an inverse sense using the measured colloid density profile to extract the underlying colloid-surface potential; this can be thought of as a microscopic analog to the well-known procedure of using the macroscopic (coarse-grained) density profile to extract the osmotic equation of state. For the dilute interfacial fluid, the inverse DFT calculations reproduced the true colloid-surface potential to within 0.5kT at all elevations.     

Inverse density-functional theory as an interpretive tool for measuring colloid-surface interactions in dense systems

Recent advances in optical microscopy, such as total internal reflection and confocal scanning laser techniques, now permit the direct three-dimensional tracking of large numbers of colloidal particles both near and far from interfaces. A novel application of this technology, currently being developed by one of the authors under the name of diffusing colloidal probe microscopy (DCPM), is to use colloidal particles as probes of the energetic characteristics of a surface. A major theoretical challenge in implementing DCPM is to obtain the potential energy of a single particle in the external field created by the surface, from the measured particle trajectories in a dense colloidal system. In this paper we develop an approach based on an inversion of density-functional theory (DFT), where we calculate the single-particle-surface potential from the experimentally measured equilibrium density profile in a nondilute colloidal fluid. The underlying DFT formulation is based on the recent work of Zhou and Ruckenstein [Zhou and Ruckenstein, J. Chem. Phys. 112, 8079 (2000)]. For model hard-sphere and Lennard-Jones systems, using Monte Carlo simulation to provide the "experimental" density profiles, we found that the inversion procedure reproduces the true particle-surface-potential energy to an accuracy within typical DCPM experimental limitations (approximately 0.1 kT) at low to moderate colloidal densities. The choice of DFT closures also significantly affects the accuracy.     

Structure and adsorption of a hard-core multi-Yukawa fluid confined in a slitlike pore: grand canonical Monte Carlo simulation and density functional study

Because of the increasing interest in studying the phenomenon exhibited by charge-stabilized colloidal suspensions in confining geometry, we present a density functional theory (DFT) for a hard-core multi-Yukawa fluid. The excess Helmholtz free-energy functional is constructed by using the modified fundamental measure theory and Rosenfeld's perturbative method, in which the bulk direct correlation function is obtained from the first-order mean spherical approximation. To validate the established theory, grand canonical ensemble Monte Carlo (GCMC) simulations are carried out to determine the density profiles and surface excesses of multi-Yukawa fluid in a slitlike pore. Comparisons of the theoretical results with the GCMC data suggest that the present DFT gives very accurate density profiles and surface excesses of multi-Yukawa fluid in the slitlike pore as well as the radial distribution functions of the bulk fluid. Both the DFT and the GCMC simulations predict the depletion of the multi-Yukawa fluid near a nonattractive wall, while the mean-field theory fails to describe this depletion in some cases. Because the simple form of the direct correlation function is used, the present DFT is computationally as efficient as the mean-field theory, but reproduces the simulation data much better than the mean-field theory.     

Further test of third order + second-order perturbation DFT approach: hard core repulsive Yukawa fluid subjected to diverse external fields

Grand canonical Monte Carlo simulation is used to investigate density profiles of hard-core repulsive Yukawa (HCRY) model fluid under the influence of various external fields and radial distribution function (RDF) of the bulk HCRY system. The aim of these extensive simulations is to provide exact data for purely repulsive interaction potential against which the validity of a third order + second-order perturbation DFT approach can be tested. It is found that a semiempirical parametrized bridge function due to Malijevsky and Labik performs very well for the RDF of the bulk HCRY fluid. Incorporation of a bulk second-order direct correlation function (DCF) of the HCRY fluid based on the Malijevsky-Labik bridge function into the third order + second-order perturbation DFT approach yields the resulting theoretical predictions for the density profiles of inhomogeneous HCRY fluid that are in a very good agreement with the simulation data, an exception being somewhat larger deviations appearing for the structure of the fluid around the center of a hard spherical cavity. Both theory and simulation predict layering transition and gas-liquid coexistence phenomena occurring with the HCRY model fluid under confined conditions. For the case of an inverse sixth-power repulsive potential under the influence of a flat stationary wall defined by an inverse twelfth-power repulsive potential, the present third order + second-order perturbation DFT approach is found to be superior to several existing weighted density approximations (WDA) and partitioned WDA.     

Is perturbation DFT approach applicable to purely repulsive fluids?

A recently proposed third order + second order perturbation density functional theory (DFT) approach is tested for the validity and applicability to purely repulsive model fluids subjected to various external fields. Hard core repulsive Yukawa potential, point particle Yukawa potential, and inverse power potential are employed as sample models. Theoretical DFT results are compared with the corresponding simulation data obtained by grand canonical ensemble Monte Carlo simulation. This comparison indicates that the third order + second order perturbation DFT approach is suitable for these purely repulsive fluids only on condition of high accuracy of the imported bulk second order direct correlation function (DCF). However, in this case the origin of the successful performance somewhat differs from that observed for the mean field approximation applied to van der Waals fluids. In the present case it originates from the observation that the bulk second order DCF is strongly dependent on the density argument for the hard-core part, while for the distances exceeding the core dimension this dependence is considerably weaker.     

Vertical motion of a charged colloidal particle near an AC polarized electrode with a nonuniform potential distribution: theory and experimental evidence

Electroosmotic flow in the vicinity of a colloidal particle suspended over an electrode accounts for observed changes in the average height of the particle when the electrode passes alternating current at 100 Hz. The main findings are (1) electroosmotic flow provides sufficient force to move the particle and (2) a phase shift between the purely electrical force on the particle and the particle's motion provides evidence of an E2 force acting on the particle. The electroosmotic force in this case arises from the boundary condition applied when faradaic reactions occur on the electrode. The presence of a potential-dependent electrode reaction moves the likely distribution of electrical current at the electrode surface toward uniform current density around the particle. In the presence of a particle the uniform current density is associated with a nonuniform potential; thus, the electric field around the particle has a nonzero radial component along the electrode surface, which interacts with unbalanced charge in the diffuse double layer on the electrode to create a flow pattern and impose an electroosmotic-flow-based force on the particle. Numerical solutions are presented for these additional height-dependent forces on the particle as a function of the current distribution on the electrode and for the time-dependent probability density of a charged colloidal particle near a planar electrode with a nonuniform electrical potential boundary condition. The electrical potential distribution on the electrode, combined with a phase difference between the electric field in solution and the electrode potential, can account for the experimentally observed motion of particles in ac electric fields in the frequency range from approximately 10 to 200 Hz.     

Atomic scale insights on chlorinated gamma-alumina surfaces

The thermochemistry of chlorinated gamma-alumina surfaces is explored by means of density functional calculations as a function of relevant reaction conditions used in experiments and in high-octane fuel production in the refining industry such as hydrocarbon isomerization and reforming. The role of chlorine as a dope of the Br?nsted acidity of gamma-alumina surfaces is investigated at an atomic scale. Combining infrared spectroscopy and density functional theory calculations, the most favorable location of chlorine atoms on the (110), (100) and (111) surfaces of gamma-alumina is found to result either from direct adsorption or from the exchange of basic hydroxyl groups. Moreover, the modification of the hydrogen bond network upon chlorine adsorption is put forward as a key parameter for changing the Br?nsted acidity. In a second step, we use a thermodynamic approach based on DFT total energy calculations corrected by the chemical potentials of HCl and H2O to determine the adsorption isotherms of chlorine and the relative surface concentration of hydroxyl groups and chlorine species on the gamma-alumina surfaces. The determination of chlorine content as a function of temperature and partial pressures of H2O and HCl offers new quantitative data required for optimizing the state of the support surface in industrial conditions. The mechanisms of chlorination are also discussed as a function of reaction conditions.     

Photodesorption of diatomic molecules from surfaces: A theoretical approach based on first principles

Photodesorption of small molecules from surfaces is one of the most fundamental processes in surface photochemistry. Despite its apparent simplicity, a microscopic understanding beyond a qualitative picture still poses a true challenge for theory. While the dynamics of nuclear motion can be treated on various levels of sophistication, all approaches suffer from the lack of sufficiently accurate potential energy surfaces, in particular for electronically excited states involved in the desorption scenario.In the last decade, a systematic and accurate methodology has been developed which allows a reliable calculation of accurate ground and excited state potential energy surfaces (PES) for different adsorbate–substrate systems. These potential energy surfaces serve as a prerequisite for subsequent quantum dynamical wave packet calculations, which allow for a direct simulation of experimentally observable quantities such as quantum state resolved velocity distributions.In the first part of this review, we will focus on scalar properties of desorbing diatomic molecules from insulating surfaces, where we also present a recently developed strategy of obtaining accurate potential energy surfaces using quantum chemical approaches. In general, diatomic molecules on large band gap materials such as oxide surfaces are studied which allows the use of sufficiently large cluster models and accurate ab initio methods beyond density functional theory (DFT). In the second part, we will focus on the vectorial aspects of the dynamics of nuclear motion and present simulations of experimentally accessible observables such as velocity distributions, Doppler profiles and alignment parameters. For each system, the microscopic mechanism of photodesorption is elucidated. We will demonstrate that the driving force of surface photochemistry is strongly dependent on details of the electronic structure of the adsorbate–substrate systems. This implies that great caution is advisable if experimental results are interpreted using empirical or semi-empirical models.     

A theoretical study for nanoparticle partitioning in the lamellae of diblock copolymers

Morphology control is important for practical applications of composite materials that consist of functional polymers and nanoparticles. Toward that end, block copolymers provide useful templates to arrange nanoparticles in the scaffold of self-organized polymer microdomains. This paper reports theoretical predictions for the distribution of nanoparticles in the lamellar structures of symmetric diblock copolymers on the basis of a polymer density functional theory (DFT) and the potential distribution theorem (PDT). The DFT predicts periodic spacing of lamellar structures in good agreement with molecular dynamics simulations. With the polymer structure from DFT as the input, the PDT is used to examine the effects of particle size, surface energy, polymer chain length, and compressibility on the distribution of nanoparticles in the limit of low particle density. It is found that the nanoparticle distribution depends not only on the particle size and surface energy but also on the local structure of the microdomain interface, polymer chain length, and compressibility. The theoretical predictions are compared well with experiments and simulations.     

Anomalous potentials from inverse analyses of interfacial polydisperse attractive colloidal fluids

This paper investigates effects of using monodisperse inverse analyses to extract particle-particle and particle-surface potentials from simulated interfacial colloidal fluids of polydisperse attractive particles. Effects of polydispersity are investigated as functions of particle concentration and attractive well depth and range for van der Waals and depletion potentials. Forward Monte Carlo simulations are used to generate particle distribution functions for polydisperse interfacial colloidal fluids from which inverted potentials are obtained using an inverse Ornstein-Zernike analysis and an inverse Monte Carlo simulation method. Attractive potentials are successfully recovered for monodisperse colloidal fluids, but polydispersity that is unaccounted for in inverse analyses produces (1) apparent softening of strong forces, (2) anomalous repulsive and attractive interactions, and (3) aphysical particle overlaps. This investigation provides insights into the role of polydispersity in altering the equilibrium structure and corresponding inverted potentials of attractive colloidal fluids near surfaces. These findings should assist the design and interpretation of optical microscopy experiments involving interfacial colloidal fluids similar to the simulated experiments reported here.     

Sedimentation equilibrium of colloidal suspensions in a planar pore based on density functional theory and the hard-core attractive Yukawa model

The sedimentation equilibrium of colloidal suspensions modeled by hard-core attractive Yukawa (HCAY) fluids in a planar pore is studied. The density profile of the HCAY fluid in a gravitational field and its distribution between the pore and uniform phases are investigated by a density functional theory (DFT) approach, which results from employing a recently proposed parameter-free version of the Lagrangian theorem-based density functional approximation (Zhou, S. Phys. Lett. A 2003, 319, 279) for hard-sphere fluids to the hard-core part of the HCAY fluid, and the second-order functional perturbation expansion approximation to the tail part as was done in a recent partitioned density functional approximation (Zhou, S. Phys. Rev. E: Stat. Phys., Plasmas, Fluids, Relat. Interdiscip. Top. 2003, 68, 061201). The resultant DFT approach is, thus, the first adjustable parameter-free DFT for HCAY fluids. The validity of the present DFT for HCAY fluids of reduced range parameter z(red) = 1.8 under various external potentials is established in the first of the papers cited previously. The present DFT for HCAY fluids can predict the radial distribution function for the bulk HCAY fluid accurately in the colloidal limit (large value of z(red)), and in the hard-sphere limit, its prediction for the density profile of the hard-sphere fluid in a gravitational field is in very good agreement with the existing simulation data. The dependence of the density profile and distribution coefficient on the magnitude of the interparticle attraction, gravitational field, and degree of confinement is investigated in detail by the present DFT approach. Intuitive and qualitative analyses are also compared with the quantitative DFT calculational results.     

Effective interactions in colloid-semipermeable membrane systems

We investigate effective interactions between a colloidal particle, immersed in a binary mixture of smaller spheres, and a semipermeable membrane. The colloid is modeled as a big hard sphere, and the membrane is represented as an infinitely thin surface, which is fully permeable to one of the smaller spheres and impermeable to the other one. Within the framework of the density functional theory, we evaluate depletion potentials and we consider two different approximate theories: the simple Asakura-Oosawa approximation and the accurate White-Bear version of the fundamental measure theory. The effective potentials are compared with the corresponding potentials for the hard, nonpermeable wall. Using statistical-mechanical sum rules, we argue that the contact value of the depletion potential between a colloid and a semipermeable membrane is smaller in magnitude than the potential between a colloid and a hard wall. A heuristic argument is provided that the colloid-semipermeable membrane effective interactions are generally weaker than these near a hard nonpermeable wall. These predictions are confirmed by explicit calculations, and the effect is more pronounced for smaller osmotic pressures. The depletion potential for a colloidal particle inside a semipermeable vesicle is stronger than the potential for the colloidal particle located outside of a vesicle. We find that the asymptotic decay of the depletion potential for the semipermeable membrane is similar to that for the nonpermeable wall and reflects the asymptotics of the total correlation function of the corresponding binary mixture of smaller spheres. Our results demonstrate that the ability of the membrane to change its shape as well as specific interactions constitute an important factor in determining the effective interactions between the semipermeable membrane and the colloidal macroparticle.     

A combined density functional and x-ray diffraction study of Pt nanoparticle structure

The structure of 1.7 nm Pt nanoparticles is investigated using x-ray diffraction (XRD) measurements and density functional theory (DFT) calculations. Two types of particles are compared, those made by solution chemistry which are capped either by thiol or amine ligands, and dendrimer encapsulated particles (DENs) which do not have capping ligands. All particles were dried before analyzing their structure. Pair distribution function (PDF) data from XRD measurements show that the ligand-capped particles are more disordered than the DENs. To determine the structure of the particles and the nature of the ligand-induced disorder, we use a hybrid reverse Monte Carlo approach. A weighted average of the calculated binding energy of the particles and a goodness-of-fit parameter to the PDF data is taken as the object function, which is minimized to determine the optimal structure. A scan over different weights gives the set of pareto optimal structures, which show how well simultaneous agreement can be reached to both experiment and theory. Using an embedded atom potential to sample configuration space and DFT to refine the optimal structures, we show that the DEN structure is most consistent with a face centered cubic lattice of truncated octahedral shape. The disorder induced by the capping ligands is consistent with surface relaxation of the particle rather than disorder of the crystal structure.     

Transient electrophoresis in a suspension of charged particles with arbitrary electric double layers

The startup of electrophoretic motion in a suspension of spherical colloidal particles, which may be charged with constant zeta potential or constant surface charge density, due to the sudden application of an electric field is analytically examined. The unsteady modified Stokes equation governing the fluid velocity field is solved with unit cell models. Explicit formulas for the transient electrophoretic velocity of the particle in a cell in the Laplace transforms are obtained as functions of relevant parameters. The transient electrophoretic mobility is a monotonic decreasing function of the particle-to-fluid density ratio and in general a decreasing function of the particle volume fraction, but it increases and decreases with a raise in the ratio of the particle radius to the Debye length for the particles with constant zeta potential and constant surface charge density, respectively. On the other hand, the relaxation time in the growth of the electrophoretic mobility increases substantially with an increase in the particle-to-fluid density ratio and with a decrease in the particle volume fraction but is not a sensitive function of the ratio of the particle radius to the Debye length. For specified values of the particle volume fraction and particle-to-fluid density ratio in a suspension, the relaxation times in the growth of the particle mobility in transient electrophoresis and transient sedimentation are equivalent.     

A Cu/Pt near-surface alloy for water-gas shift catalysis

The primary route to hydrogen production from fossil fuels involves the water-gas shift (WGS) reaction, and an improvement in the efficiency of WGS catalysts could therefore lead to a major leap forward in the realization of hydrogen economy. On the basis of a combination of high-resolution scanning tunneling microscopy, X-ray photoelectron spectroscopy, and density functional theory (DFT) calculations, we suggest the existence of a new thermodynamically stable Cu/Pt near-surface alloy (NSA). Temperature-programmed desorption and DFT reveal that this Cu/Pt NSA binds CO significantly more weakly than does Pt alone, thereby implying a considerable reduction in the potential for CO poisoning of the Cu/Pt NSA surface as compared to that of pure Pt. In addition, DFT calculations show that this Cu/Pt NSA is able to activate H2O easily, which is the rate-determining step for the WGS on several metal surfaces, and, at the same time, to bind the products of that reaction and formate intermediates rather weakly, thus avoiding possible poisoning of the catalyst surface. The Cu/Pt NSA is thus a promising candidate for an improved WGS catalyst.     

The colloidal force of bead-spring chains in a good solvent

A recently developed density functional theory (DFT) for tethered bead-spring chains is used to investigate colloidal forces for the good solvent case. A planar surface of tethered chains is opposed to a bare, hard wall and the force exerted on the bare wall is calculated by way of the contact density. Previously, the case of large wall separation was investigated. The density profiles of the unperturbed chains, in that case, were found to be neither stepfunctions nor parabolas and were shown to accurately predict computer simulation results. In the present paper, the surface forces that result from the distortion of these density profiles at finite wall separation is studied. The resulting force function is analyzed for varying surface coverages, wall separations, and chain lengths. The results are found to be in near quantitative agreement with the scaling predictions of Alexander [S. Alexander, J. Phys. (Paris) 38, 983 (1977)] when the layer thickness is "correctly" defined. Finally, a hybrid Alexander-DFT theory is suggested for the analysis of experimental results.     

Selective hydrogenation of the C=O bond in acrolein through the architecture of bimetallic surface structures

In the current study we have performed experimental studies and density functional theory (DFT) modeling to investigate the selective hydrogenation of the C=O bond in acrolein on two bimetallic surface structures, the subsurface Pt-Ni-Pt(111) and surface Ni-Pt-Pt(111). We have observed for the first time the production of the desirable unsaturated alcohol (2-propenol) on Pt-Ni-Pt(111) under ultra-high vacuum conditions. Furthermore, our DFT modeling revealed a general trend in the binding energy and bonding configuration of acrolein with the surface d-band center of Pt-Ni-Pt(111), Ni-Pt-Pt(111), and Pt(111), suggesting the possibility of using the value of the surface d-band center as a parameter to predict other bimetallic surfaces for the selective hydrogenation of acrolein.     

Motion of a colloidal sphere with interfacial self-electrochemical reactions induced by a magnetic field

The motion of a spherical colloidal particle with spontaneous electrochemical reactions occurring on its surface in an ionic solution subjected to an applied magnetic field is analyzed for an arbitrary zeta potential distribution. The thickness of the electric double layer adjacent to the particle surface is assumed to be much less than the particle radius. The solutions of the Laplace equations governing the magnetic scalar potential and electric potential, respectively, lead to the magnetic flux and electric current density distributions in the particle and fluid phases of arbitrary magnetic permeabilities and electric conductivities. The Stokes equations modified with the Lorentz force contribution for the fluid motion are dealt by using a generalized reciprocal theorem, and closed-form formulas for the translational and angular velocities of the colloidal sphere induced by the magnetohydrodynamic effect are obtained. The dipole and quadrupole moments of the zeta potential distribution over the particle surface cause the particle translation and rotation, respectively. The induced velocities of the particle are unexpectedly significant, and their dependence on the characteristics of the particle-fluid system is physically different from that for electromagnetophoretic particles or phoretic swimmers.     

Interaction of micrometer-scale particles with nanotextured surfaces in shear flow

Dynamic particle adhesion from flow over collecting surfaces with nanoscale heterogeneity occurs in important natural systems and current technologies. Accurate modeling and prediction of the dynamics of particles interacting with such surfaces will facilitate their use in applications for sensing, separating, and sorting colloidal-scale objects. In this paper, the interaction of micrometer-scale particles with electrostatically heterogeneous surfaces is analyzed. The deposited polymeric patches that provide the charge heterogeneity in experiments are modeled as 11-nm disks randomly distributed on a planar surface. A novel technique based on surface discretization is introduced to facilitate computation of the colloidal interactions between a particle and the heterogeneous surface based on expressions for parallel plates. Combining these interactions with hydrodynamic forces and torques on a particle in a low Reynolds number shear flow allows particle dynamics to be computed for varying net surface coverage. Spatial fluctuations in the local surface density of the deposited patches are shown responsible for the dynamic adhesion phenomena observed experimentally, including particle capture on a net-repulsive surface.     

Relevance of electrokinetic theory for "soft" particles to bacterial cells: implications for bacterial adhesion

Bacterial cells and other biological particles carry charged macromolecules on their surface that form a "soft" ion-permeable layer. In this paper, we test the applicability of an electrokinetic theory for soft particles to characterize the electrophoretic mobility (EPM) and adhesion kinetics of bacterial cells. The theory allows the calculation of two parameters--the electrophoretic softness and the fixed charged density--that define the characteristics of the polyelectrolyte layer at the soft particle surface. The theory also allows the calculation of an outer-surface potential that may better predict the electrostatic interaction of soft particles with solid surfaces. To verify its relevance for bacterial cells, the theory was applied to EPM measurements of two well-characterized Escherichia coli K12 mutants having lipopolysaccharide (LPS) layers of different lengths and molecular compositions. Results showed that the obtained softness and fixed charge density were not directly related to the known characteristics of the LPS of the selected strains. Interaction energy profiles calculated from Derjaguin-Landau-Verwey-Overbeek (DLVO) theory were used to interpret bacterial deposition (adhesion) rates on a pure quartz surface. The outer surface potential failed to predict the low attachment efficiencies of the two bacterial strains. The lack of success in the application of the theory for soft particles to bacterial cells is attributed to chemical and physical heterogeneities of the polyelectrolyte layer at the cell surface.