Comparison of random-walk density functional theory to simulation for bead-spring homopolymer melts

Density profiles for a homopolymer melt near a surface are calculated using a random-walk polymeric density functional theory, and compared to results from molecular dynamics simulations. All interactions are of a Lennard-Jones form, for both monomer-monomer interactions and surface-monomer interactions, rather than the hard core interactions which have been most investigated in the literature. For repulsive systems, the theory somewhat overpredicts the density oscillations near a surface. Nevertheless, near quantitative agreement with simulation can be obtained with an empirical scaling of the direct correlation function. Use of the random phase approximation to treat attractive interactions between polymer chains gives reasonable agreement with simulation of dense liquids near neutral and attractive surfaces.     

Structure of Lennard-Jones fluids confined in square nanoscale channels from density functional theory

The density distribution of Lennard-Jones fluids confined in square nanoscale channels with Lennard-Jones walls has been studied using the nonlocal density functional theory (DFT) based on the Tarazona model. The effect of channel lengths on the density profiles with various chemical potentials was discussed. It was found that there is an apparent layering phenomenon for the confined fluids due to the combining influences of the enhancing solid-fluid interaction and the excluded volume effect. The pronounced density peaks were observed at the corners of square channels due to the strong fluid-solid interactions. The grand canonical ensemble Monte Carlo simulation (GCEMC) was applied to test the nonlocal DFT results. The DFT calculations are in relatively good agreement with the GCEMC simulations. The adsorption isotherms in a series of square channels were evaluated based on the obtained density distributions. The adsorption mechanism within the square pores was investigated. A comparison between the adsorptions of the square pores with those of the corresponding slit-size pores has been given.     

Interfacial properties of Lennard-Jones chains by direct simulation and density gradient theory

We perform a series of molecular dynamics simulations of Lennard-Jones chains systems, up to tetramers, in order to investigate the influence of temperature and chain length on their phase separation and interfacial properties. Simulation results serve as a test to check the accuracy of a statistical associated fluid theory (soft-SAFT) coupled with the density gradient theory. We focus on surface tension and density profiles. The simulations allow us to discuss the success and limitations of the theory and how to estimate the only adjustable parameter, the influence parameter. This parameter is obtained by fitting the surface tension, and then used to obtain the density profiles in a predictive manner. A good agreement is found if the temperature dependence of this parameter is neglected.(c) 2004 American Institute of Physics.     

Study of structural and thermodynamic parameters of adsorption layers using density functional theory. Local density and adsorption isotherms on spherical surfaces

Local density profiles in adsorption layers of Lennard-Jones fluids on two-dimensional adsorbents with spherical geometry and isotherms of excess (Gibbs) adsorption have been calculated using the classical density functional theory (approximations with weighting coefficients). The local density profiles have been found in hydrogen adsorption layers on C₆₀, C₂₄₀, and C₅₄₀ fullerene molecules. The calculations have been performed for both subcritical and supercritical temperature ranges. It has been shown that, at a pressure of 10 MPa and a temperature of 77 K, the gravimetric (mass) hydrogen density on C₆₀ fullerene is 7.6 wt %, which is in good agreement with the results of molecular dynamics simulation and experimental data. It has also been established that the gravimetric hydrogen density on C₆₀ fullerene is higher than that on C₂₄₀ and C₅₄₀ fullerenes, being comparable with its value in a slitlike pore of a carbon adsorbent.     

Application of the extended LJ potential-based equation of state to predict the density of five different classes of refrigerant systems (HCFCs,HFCs, HFEs,PFAs, and PFAAs)

The major goal of this work is to apply the extended Lennard-Jones potential-based equation of state (ELJ-based EoS) to predict the density of five different classes of refrigerant systems including chlorofluorocarbons, hydrochlorofluorocarbons, hydrofluorocarbons, hydrofluoroethers, perfluoroalkanes, and perfluoroalkylalkane. This EoS is based on an effective near-neighbor pair potential of the Lennard-Jones (12,6,3) type. The temperature dependencies of the parameters of the equation of state can be calculated at any temperature for each refrigerant. The calculated parameters along with the ELJ-based EoS have been used to calculate the density and isothermal compressibility coefficient of different refrigerants. A comparison between the predicted results and experimental data shows that the agreement is good. The total absolute average deviation of density for 14,871 data was found to be 0.34 compared with experimental data. Comparisons with the other EoSs show that the ELJ-based EoS is more accurate than other EoSs for most of the studied refrigerants.     

Calculations of van der Waals interaction energies for Al, Cr, Fe, and Ir acetylacetonate crystals

Experimental data have been analyzed and interpreted for four volatile acetylacetonates of trivalent metals Al, Cr, Fe, and Ir. The crystal lattice energies were calculated by the atom-atom potential method. The lattice energies obtained by using the Buckingham potential are in better agreement with the sublimation heats of these metal complexes than those calculated from the Lennard-Jones potential. The experimental dependences of vapor pressure for the complexes are in satisfactory agreement with the values obtained from the calculated lattice energies and entropies of crystal-gas transitions.     

Interatomic potential-based semiclassical theory for Lennard-Jones fluids

An interatomic potential based semiclassical theory is proposed to predict the concentration and potential profiles of a Lennard-Jones (LJ) fluid confined in a channel. The inputs to the semiclassical formulation are the LJ parameters of the fluid and the wall, the density of channel wall atoms, and the average concentration of the fluid inside the channel. Using the semiclassical formulation, fluid confinement in channel with widths ranging from 2sigma ff to 100sigma ff, where sigma ff is the fluid-fluid LJ distance parameter, is investigated. The concentration and potential predicted by the semiclassical formulation are found to be in good agreement with those from equilibrium molecular dynamics simulations. While atomistic simulations in large channels are computationally expensive, the proposed semiclassical formulation can rapidly and accurately predict the concentration and potential profiles. The proposed semiclassical theory is thus a robust and fast method to predict the interfacial and "bulk" fluid phenomena in channels with widths ranging from the macroscale down to the scale of a few atomic diameters.     

Comparing the density of states of binary Lennard-Jones glasses in bulk and film

We used Wang-Landau density of states Monte Carlo to study a binary Lennard-Jones glass-forming mixture in bulk and films between noninteracting walls. Thermodynamic properties are calculated using two different ensembles and film data are compared with the bulk. Bulk properties are in good agreement with previous simulations. We confirm the formation of a glass using various properties, e.g., energy, heat capacity, and pressure with temperature. We find a change in slope in the energy per particle and pressure as a function of temperature. We do not find any defined crystal structure. A higher glass transition temperature is found for the film.     

Perturbed-chain equation of state for the solid phase

A perturbed chain equation of state for the solid phase has been derived. Although the equation is general with respect to intermolecular potential, we incorporate the Lennard-Jones potential in this work in order to compare results from the model with available Monte Carlo simulation data. Two forms of the radial distribution function for the hard-sphere solid chain reference state are used in the model. First, a theoretically rigorous approach is taken by using a correlation of actual solid-phase Monte Carlo hard-sphere chain data for the radial distribution function. This results in good agreement with the Monte Carlo data only at high density. Second, a simple extended-density approximation was used for the radial distribution function. This second approach was found to work well across the entire density range including the vicinity of the solid-fluid equilibrium.     

Criticality of a liquid-vapor interface from an inhomogeneous integral equation theory

A microscopic theory is developed to study the liquid-vapor interfacial properties of simple fluids with ab initio treatment of the inhomogeneous two-body correlation functions, without any interpolation. It consists of the inhomogeneous Ornstein-Zernike equation coupled with the Duh-Henderson-Verlet closure and the Lovett-Mou-Buff-Wertheim equation. For the liquid-vapor interface of the Lennard-Jones fluid, we obtained the density profile and the surface tension, as well as their critical behaviour. In particular, we identified non-classical critical exponents. The theory accurately predicts the phase diagram and the interfacial properties in a very good agreement with simulations. We also showed that the method leads to true capillary-wave asymptotics in the macroscopic limit.     

Mean-field kinetic nucleation theory

A new semiphenomenological model of homogeneous vapor-liquid nucleation is proposed in which the cluster kinetics follows the "kinetic approach to nucleation" and the thermodynamic part is based on the revised Fisher droplet model with the mean-field argument for the cluster configuration integral. The theory is nonperturbative in a cluster size and as such is valid for all clusters down to monomers. It contains two surface tensions: macroscopic (planar) and microscopic. The latter is a temperature dependent quantity related to the vapor compressibility factor at saturation. For Lennard-Jones fluids the microscopic surface tension possesses a universal behavior with the parameters found from the mean-field density functional calculations. The theory is verified against nucleation experiments for argon, nitrogen, water, and mercury, demonstrating very good agreement with experimental data. Classical nucleation theory fails to predict experimental results when a critical cluster becomes small.     

Fundamental measure density functional theory studies on the freezing of binary hard-sphere and Lennard-Jones mixtures

Free energies and correlation functions of liquid and solid hard-sphere (HS) mixtures are calculated using the fundamental measure density functional theory. Using the thermodynamic perturbation theory the free energies of solid and liquid Lennard-Jones (LJ) mixtures are obtained from correlation functions of HS systems within a single theoretical approach. The resulting azeotrope- and spindle-type solid-liquid phase diagrams of HS and LJ binary mixtures are in good agreement with the corresponding ones from computer simulations.     

Study of structural and thermodynamic characteristics of adsorption layers by the density functional method: Local density in adsorption layer on planar solid surface and adsorption isotherms

The local density profiles in Lennard-Jones adsorption layers, as well as the excess (Gibbs) and absolute adsorption values, are calculated by the density functional method (weighting factor approximations). The substrate is described using the single-particle potential corresponding to the Lennard-Jones potential integrated over the half space occupied by the substrate. The Steele potential is used as a single-particle potential to consider methane adsorption on the surface of nonporous graphite as a specific system. The calculations are performed for both sub- and supercritical temperature regions. It is established that the density profiles are characterized by the existence of one to three maxima, which reflect the positional order of molecules in adsorption layers, i.e., the layered structure of an adsorbate.     

Mercedes-Benz water molecules near hydrophobic wall: integral equation theories vs Monte Carlo simulations

Associative version of Henderson-Abraham-Barker theory is applied for the study of Mercedes-Benz model of water near hydrophobic surface. We calculated density profiles and adsorption coefficients using Percus-Yevick and soft mean spherical associative approximations. The results are compared with Monte Carlo simulation data. It is shown that at higher temperatures both approximations satisfactory reproduce the simulation data. For lower temperatures, soft mean spherical approximation gives good agreement at low and at high densities while in at mid range densities, the prediction is only qualitative. The formation of a depletion layer between water and hydrophobic surface was also demonstrated and studied.     

A mode-coupling theory treatment of the transport coefficients of the Lennard-Jones fluid

We apply mode-coupling theory to study shear viscosity and self-diffusion coefficient of the Lennard-Jones fluid throughout the entire fluid region of the phase diagram. Theoretical results are compared with the extensive simulation data and good agreement is found. In addition, theory is compared to the experimental data on the transport coefficients of inert gas fluids.     

Density functional theory of size-dependent surface tension of Lennard-Jones fluid droplets using a double well type Helmholtz free energy functional

A double well type Helmholtz free energy density functional and a model density profile for a two phase vapor-liquid system are used to obtain the size-dependent interfacial properties of the vapor-liquid interface at coexistence condition along the lines of van der Waals and Cahn and Hilliard density functional formalism of the interface. The surface tension, temperature-density curve, density profile, and thickness of the interface of Lennard-Jones fluid droplet-vapor equilibrium, as predicted in this work are reported. The planar interfacial properties, obtained from consideration of large radius of the liquid drop, are in good agreement with the results of other earlier theories and experiments. The same free energy model has been tested by solving the equations numerically, and the results compare well with those from the use of model density profile.     

Analysis of reduced density matrices in the co-ordinate representation. III. Electron density and correlation in the ground states of H2 and the H6 ring system within some approximations of the simple LCAO type

The electron density and spatial correlation as given by the MO , VB and AMO methods for H₂ and H₆ are studied by means of diagrams. For comparison, diagrams representing accurate wave functions for H₂ are also given. The study of model functions representing localized bonds leads to results concerning the nature of localized electron pairs in agreement with those of Lennard-Jones.     

Thermal conductivity of solid argon from molecular dynamics simulations

The thermal conductivity of solid argon in the classical limit has been calculated by equilibrium molecular dynamic simulations using the Green-Kubo formalism and a Lennard-Jones interatomic potential. Contrary to previous theoretical reports, we find that the computed thermal conductivities are in good agreement with experimental data. The computed values are also in agreement with the high-temperature limit of the three-phonon scattering contribution to the thermal conductivity. We find that finite-size effects are negligible and that phonon lifetimes have two characteristic time scales, so that agreement with kinetic theory is obtained only after appropriate averaging of the calculated phonon lifetimes.     

Calculations of crystal-melt interfacial free energies by nonequilibrium work measurements

We developed a multistep thermodynamic perturbation method to compute the interfacial free energies by nonequilibrium work measurements with cleaving potential procedure. Using this method, we calculated the interfacial free energies of different crystal orientations for the Lennard-Jones system. Our results are in good agreement with the results by thermodynamic integration method. Compared with thermodynamic integration method, the multistep thermodynamic perturbation method is more efficient. For each stage of the cleaving process, only a few thermodynamic perturbation steps are needed, and there is no requirement on the reversibility of the path.