Perturbation density functional theory and Monte Carlo simulations for the structure of hard triatomic fluids in slitlike pores

Microscopic density functional theory for monolayers of diblock copolymers

We propose density functional theory for diblock copolymers in two dimensions. Our theoretical framework is based on Wertheim's first order thermodynamic perturbation theory. Using the proposed approach, we investigate the structure and phase behavior of monolayers of symmetric diblock copolymers. We find that the phase behavior of symmetric diblock copolymer monolayers is similar to that in 3D. This includes the scaling of the equilibrium lamellar width with chain length. We find that the topology of the resulting phase diagrams depends on the chain length and the unlike segment interaction incompatibility and involves either one, two, or three triple points (one of them being the peritectic point). We expect that a similar phase behavior could be obtained for monolayers of colloidal suspensions with carefully tuned interparticle interactions.

Structure of room temperature ionic liquids

The structure of room temperature ionic liquids is studied using molecular dynamics simulations and integral equation theory. Three ionic liquids 1-alkyl-3-methylimidazolium hexfluorophosphate, [C n MIM] [PF6], for n  =  1, 4, and 8, are studied using a united atom model of the ions. The primary interest is a study of the pair correlation functions and a test of the reference interaction site model theory. There is liquid-like ordering in the liquid that arises from electrostatic attractions and steric packing considerations. The theory is not in quantitative agreement with the simulation results and underestimates the degree of liquid-like order. A pre-peak in the static structure factor is seen in both simulations and theory, suggesting that this is a geometric effect arising from a packing of the alkyl chains.

Density functional theory for polymeric systems in 2D

We propose density functional theory for polymeric fluids in two dimensions. The approach is based on Wertheim's first order thermodynamic perturbation theory (TPT) and closely follows density functional theory for polymers proposed by Yu and Wu (2002 J. Chem. Phys. 117 2368). As a simple application we evaluate the density profiles of tangent hard-disk polymers at hard walls. The theoretical predictions are compared against the results of the Monte Carlo simulations. We find that for short chain lengths the theoretical density profiles are in an excellent agreement with the Monte Carlo data. The agreement is less satisfactory for longer chains. The performance of the theory can be improved by recasting the approach using the self-consistent field theory formalism. When the self-avoiding chain statistics is used, the theory yields a marked improvement in the low density limit. Further improvements for long chains could be reached by going beyond the first order of TPT.

A hybrid method for predicting the microstructure of polymers with complex architecture: Combination of single-chain simulation with density functional theory

A hybrid method is proposed to investigate the microstructure of various polymeric fluids confined between two parallel surfaces. The hybrid method combines a single-chain Monte Carlo (MC) simulation for the ideal-gas part of the Helmholtz energy and a density functional theory (DFT) for the excess part that arises from nonbonded intersegment interactions. The latter consists of a modified fundamental measure theory for excluded-volume effect, the first-order thermodynamics perturbation theory for chain connectivity, and a mean-field approximation for the van der Waals attraction. In comparison with a conventional DFT, the hybrid method avoids calculation of the time-consuming recursive functions and is directly applicable to polymers with arbitrary molecular architecture. Its numerical performance has been validated by extensive comparisons with MC data for the density distributions of totally flexible, semiflexible, or rigid polymers and those with starlike architecture. Special attention is also given to the formation of a nematic monolayer by rigid molecules laying perpendicular to a planar surface. The hybrid method predicts the surface pressure versus surface coverage in good agreement with experiment.

Effects of potential models on the adsorption of ethane and ethylene on graphitized thermal carbon black. Study of two-dimensional critical temperature and isosteric heat versus loading

Adsorption of ethylene and ethane on graphitized thermal carbon black and in slit pores whose walls are composed of graphene layers is studied in detail to investigate the packing efficiency, the two-dimensional critical temperature, and the variation of the isosteric heat of adsorption with loading and temperature. Here we used a Monte Carlo simulation method with a grand canonical Monte Carlo ensemble. A number of two-center Lennard-Jones (LJ) potential models are investigated to study the impact of the choice of potential models in the description of adsorption behavior. We chose two 2C-LJ potential models in our investigation of the (i) UA-TraPPE-LJ model of Martin and Siepmann for ethane and Wick et al. for ethylene and (ii) AUA4-LJ model of Ungerer et al. for ethane and Bourasseau et al. for ethylene. These models are used to study the adsorption of ethane and ethylene on graphitized thermal carbon black. It is found that the solid-fluid binary interaction parameter is a function of adsorbate and temperature, and the adsorption isotherms and heat of adsorption are well described by both the UA-TraPPE and AUA models, although the UA-TraPPE model performs slightly better. However, the local distributions predicted by these two models are slightly different. These two models are used to explore the two-dimensional condensation for the graphitized thermal carbon black, and these values are 110 K for ethylene and 120 K for ethane.

Density functional theory for semiflexible and cyclic polyatomic fluids

The effects of bond angle and chain stiffness on the structures of semiflexible polyatomic fluids are investigated by incorporating the bending potential into a density functional theory [Y. X. Yu and J. Z. Wu, J. Chem. Phys. 117, 2368 (2002)] that combines a modified fundamental measure theory for the excluded-volume effects and the first-order thermodynamics perturbation theory for the chain connectivity. The refined density functional theory faithfully reproduces the density profiles and conformational properties of a variety of triatomic fluids near a hard wall in which extensive Monte Carlo simulation data are available. In particular, the theory is able to capture the structures of rigid cyclic trimers where all segments are identical. The variation of local density profiles with respect to the chain length of confined polyatomic fluids is also explored. For quadratomic fluids confined in slit pores, the density profile of the middle segments exhibits novel double peaks that are absent in a fully flexible chain model. In addition, the density functional theory is applied to predicting the conformational properties and adsorption behavior of heterogeneous triatomic fluids of type "ABB" mimicking surfactant molecules. The competition between surface adsorption and self-association of trimers consisting of surface active and self-binding "A" segments and neutral "B" segment is explored.

An accurate density functional theory for the vapor-liquid interface of associating chain molecules based on the statistical associating fluid theory for potentials of variable range

A Helmholtz free energy density functional is developed to describe the vapor-liquid interface of associating chain molecules. The functional is based on the statistical associating fluid theory with attractive potentials of variable range (SAFT-VR) for the homogenous fluid [A. Gil-Villegas, A. Galindo, P. J. Whitehead, S. J. Mills, G. Jackson, and A. N. Burgess, J. Chem. Phys. 106, 4168 (1997)]. A standard perturbative density functional theory (DFT) is constructed by partitioning the free energy density into a reference term (which incorporates all of the short-range interactions, and is treated locally) and an attractive perturbation (which incorporates the long-range dispersion interactions). In our previous work [F. J. Blas, E. Martin del Rio, E. de Miguel, and G. Jackson, Mol. Phys. 99, 1851 (2001); G. J. Gloor, F. J. Blas, E. Martin del Rio, E. de Miguel, and G. Jackson, Fluid Phase Equil. 194, 521 (2002)] we used a mean-field version of the theory (SAFT-HS) in which the pair correlations were neglected in the attractive term. This provides only a qualitative description of the vapor-liquid interface, due to the inadequate mean-field treatment of the vapor-liquid equilibria. Two different approaches are used to include the correlations in the attractive term: in the first, the free energy of the homogeneous fluid is partitioned such that the effect of correlations are incorporated in the local reference term; in the second, a density averaged correlation function is incorporated into the perturbative term in a similar way to that proposed by Toxvaerd [S. Toxvaerd, J. Chem. Phys. 64, 2863 (1976)]. The latter is found to provide the most accurate description of the vapor-liquid surface tension on comparison with new simulation data for a square-well fluid of variable range. The SAFT-VR DFT is used to examine the effect of molecular chain length and association on the surface tension. Different association schemes (dimerization, straight and branched chain formation, and network structures) are examined separately. The surface tension of the associating fluid is found to be bounded between the nonassociating and fully associated limits (both of which correspond to equivalent nonassociating systems). The temperature dependence of the surface tension is found to depend strongly on the balance between the strength and range of the association, and on the particular association scheme. In the case of a system with a strong but very localized association interaction, the surface tension exhibits the characteristic "s shaped" behavior with temperature observed in fluids such as water and alkanols. The various types of curves observed in real substances can be reproduced by the theory. It is very gratifying that a DFT based on SAFT-VR free energy can provide an accurate quantitative description of the surface tension of both the model and experimental systems.

Microscopic approach for the site distribution and thermodynamic properties of a single-component polymer subjected to an external field

A microscopic statistical mechanics approach is proposed for a nonuniform single-component freely jointed tangential hard-sphere polymer in the framework of density-functional theory. The present approach avoids the use of single-chain simulation in the theory by treating bonding interaction on the basis of the properties of the Dirac delta function. The present excess free energy includes all terms of functional perturbative expansion around the uniform bulk fluid in the form of the Verlet-modified bridge function. The second-order direct correlation function of a uniform polymer melt as the input parameter is obtained by solving numerically the polymer-reference-interaction-site-model integral equation with the Percus-Yevick closure. Predictions of the present approach for such microscopic structural and thermodynamics properties as the site density distribution, the partition coefficient, and the adsorption isotherm near a hard wall or between two hard walls are compared with computer-simulation results and with those of previous theories. The comparison indicates that the present approach is more accurate than the previous integral equation theory and the most accurate Monte Carlo density-functional theories. The predicted oscillations of the medium-induced force between two hard walls immersed in polymer melts are consistent with the experimental results available in the literature. The relation of the present approach with self-consistent-field theory, as well as the differences between the two, are discussed.

Density functional theories and molecular simulations of adsorption and phase transitions in nanopores

The nonlocal density functional theory (NLDFT) of confined fluids is tested against Monte Carlo simulations by using the example of Lennard-Jones (LJ) fluid sorption in slit-shaped and cylindrical nanopores ranging from 0.3 to 10 nm in width. The fluid-fluid and solid-fluid parameters of the LJ potentials were chosen to represent several experimentally important adsorption systems: nitrogen and carbon dioxide in activated carbons, zeolites, and mesoporous molecular sieves of the MCM-41 type. Freezing in nanopores is discussed using the example of methane sorption in carbon at 111 K. Comparison with reference experiments is given when available. Two versions of NLDFT, the smoothed density approximation and the fundamental measure theory, are considered. It is shown that NLDFT approaches with properly chosen parameters provide quantitative agreement with the results of Monte Carlo simulations and reference experiments. Appreciable deviations are found in extremely narrow pores of less than two molecular diameters in width. In wider pores, NLDFT models can be used for quantitative predictions of reversible and hysteretic adsorption isotherms and analyses of the specifics of phase transformations, including the equilibrium and spinodal phase transitions.