Extending Wertheim’s perturbation theory to the solid phase of Lennard-Jones chains: Determination of the global phase diagram

Phase behavior of active and passive dumbbells

We report phase separation in a mixture of "hot" and "cold" three-dimensional dumbbells which interact by Lennard-Jones potential. We also have studied the effect of asymmetry of dumbbells and the variation of ratio of "hot" and "cold" dumbbells on their phase separation. The ratio of the temperature difference between hot and cold dumbbells to the temperature of cold dumbbells is a measure of the activity χ of the system. From constant density simulations of symmetric dumbbells, we observe that the "hot" and "cold" dumbbells phase separate at higher activity ratio (χ>5.80) compared to that of a mixture of hot and cold Lennard-Jones monomers (χ>3.44). We find that, in the phase-separated system, the hot dumbbells have high effective volume and hence high entropy which is calculated by two-phase thermodynamic method. The high kinetic pressure of hot dumbbells forces the cold dumbbells to form dense clusters such that at the interface the high kinetic pressure of hot dumbbells is balanced by the virial pressure of cold dumbbells. We find that phase separation pushes the cluster of cold dumbbells to have solidlike ordering. Bond orientation order parameters reveal that the cold dumbbells form solidlike ordering consisting of predominantly face-centered cubic and hexagonal-close packing packing, but the individual dumbbells have random orientations. The simulation of the nonequilibrium system of symmetric dumbbells at different ratios of number of hot dumbbells to cold dumbbells reveals that the critical activity of phase separation decreases with increase in fraction of hot dumbbells. The simulation of equal mixture of hot and cold asymmetric dumbbells revealed that the critical activity of phase separation was independent of the asymmetry of dumbbells. We also observed that the clusters of cold asymmetric dumbbells showed both crystalline and noncrystalline order depending on the asymmetry of dumbbells.

A perturbed-chain equation of state based on Wertheim TPT for the fully flexible LJ chains in the fluid and solid phases

In the framework of thermodynamic perturbation theory (TPT), a new perturbed-chain equation of state (EOS) is presented for a fully flexible Lennard-Jones (LJ) chain system. The EOS is the sum of repulsive and perturbation contributions. The reference term of the EOS is derived based on first- and second-order TPT of Wertheim for the chains interacting with each other through the Weeks-Chandler-Anderson potential model. In order to derive the perturbation term, we have used the radial distribution function of the hard-chain system with a chain range of m = 2-10 and packing fraction range of η = 0.10-0.72, which cover the entire density range from vapor to solid phases. The performance of the EOS is tested against simulation data of the compressibility factor, residual internal energy, and phase equilibrium. A close agreement was observed across all cases. The EOS has three pure component parameters and is able to describe the global vapor-liquid-solid phase diagram of the LJ chain.

Equations of state for the fully flexible WCA chains in the fluid and solid phases based on Wertheims-TPT2

Based on Wertheim's second order thermodynamic perturbation theory (TPT2), equations of state (EOSs) are presented for the fluid and solid phases of tangent, freely jointed spheres. It is considered that the spheres interact with each other through the Weeks-Chandler-Anderson (WCA) potential. The developed TPT2 EOS is the sum of a monomeric reference term and a perturbation contribution due to bonding. MC NVT simulations are performed to determine the structural properties of the reference system in the reduced temperature range of 0.6 ≤ T* ≤ 4.0 and the packing fraction range of 0.1 ≤ η ≤ 0.72. Mathematical functions are fitted to the simulation results of the reference system and employed in the framework of Wertheim's theory to develop TPT2 EOSs for the fluid and solid phases. The extended EOSs are compared to the MC NPT simulation results of the compressibility factor and internal energy of the fully flexible chain systems. Simulations are performed for the WCA chain system for chain lengths of up to 15 at T* = 1.0, 1.5, 2.0, 3.0. Across all the reduced temperatures, the agreement between the results of the TPT2 EOS and MC simulations is remarkable. Overall Average Absolute Relative Percent Deviation at T* = 1.0 for the compressibility factor in the entire chain lengths we covered is 0.51 and 0.77 for the solid and fluid phases, respectively. Similar features are observed in the case of residual internal energy.

Vapor-liquid interfacial properties of fully flexible Lennard-Jones chains

We consider the computation of the interfacial properties of molecular chains from direct simulation of the vapor-liquid interface. The molecules are modeled as fully flexible chains formed from tangentially bonded monomers with truncated Lennard-Jones interactions. Four different model systems comprising of 4, 8, 12, and 16 monomers per molecule are considered. The simulations are performed in the canonical ensemble, and the vapor-liquid interfacial tension is evaluated using the test area and the wandering interface methods. In addition to the surface tension, we also obtain density profiles, coexistence densities, critical temperature and density, and interfacial thickness as functions of temperature, paying particular attention to the effect of the chain length on these properties. According to our results, the main effect of increasing the chain length (at fixed temperature) is to sharpen the vapor-liquid interface and to increase the width of the biphasic coexistence region. As a result, the interfacial thickness decreases and the surface tension increases as the molecular chains get longer. The interfacial thickness and surface tension appear to exhibit an asymptotic limiting behavior for long chains. A similar behavior is also observed for the coexistence densities and critical properties. Our simulation results indicate that the asymptotic regime is reached for Lennard-Jones chains formed from eight monomer segments. We also include a preliminary study on the effect of the cutoff distance on the interfacial properties. Our results indicate that all of the properties exhibit a dependence with the distance at which the interactions are truncated, though the relative effect varies from one property to the other. The interfacial thickness and, more particularly, the interfacial tension are found to be strongly dependent on the particular choice of cutoff, whereas the density profiles and coexistence densities are, in general, less sensitive to the truncation.

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.

Radial distribution function of freely jointed hard-sphere chains in the solid phase

Monte Carlo simulation is used to generate the radial distribution function of freely jointed tangent-bonded hard-sphere chains in the disordered solid phase for chain lengths of three, four, six, and eight segments. The data are used to create an accurate analytical expression of the total radial distribution function of the hard-sphere chains that covers a density range from the solidification point up to a packing fraction of 0.71. It is envisioned that the correlation will help further progress toward molecular thermodynamic treatment of the solid phase in general and toward perturbed chain theories for the solid phase, in particular.

Phase Transition of Short Linear Molecules Adsorbed on Solid Surfaces from a Density Functional Approach

A microscopic density functional theory is used to investigate the adsorption of short chains on strongly attractive solid surfaces. We analyze the structure of the adsorbed fluid and investigate how the layering transitions change with the change of the chain length and with relative strength of the fluid-solid interaction. The critical temperature of the first layering transition, rescaled by the bulk critical temperature, increases slightly with an increase of the chain length. We have found that for longer chains the layering transitions within consecutive layers are shifted toward very low temperatures and that their sequence is finally replaced by a single transition.

Computer simulation study of the global phase behavior of linear rigid Lennard-Jones chain molecules: Comparison with flexible models

The global phase behavior (i.e., vapor-liquid and fluid-solid equilibria) of rigid linear Lennard-Jones (LJ) chain molecules is studied. The phase diagrams for three-center and five-center rigid model molecules are obtained by computer simulation. The segment-segment bond lengths are L = sigma, so that models of tangent monomers are considered in this study. The vapor-liquid equilibrium conditions are obtained using the Gibbs ensemble Monte Carlo method and by performing isobaric-isothermal NPT calculations at zero pressure. The phase envelopes and critical conditions are compared with those of flexible LJ molecules of tangent segments. An increase in the critical temperature of linear rigid chains with respect to their flexible counterparts is observed. In the limit of infinitely long chains the critical temperature of linear rigid LJ chains of tangent segments seems to be higher than that of flexible LJ chains. The solid-fluid equilibrium is obtained by Gibbs-Duhem integration, and by performing NPT simulations at zero pressure. A stabilization of the solid phase, an increase in the triple-point temperature, and a widening of the transition region are observed for linear rigid chains when compared to flexible chains with the same number of segments. The triple-point temperature of linear rigid LJ chains increases dramatically with chain length. The results of this work suggest that the fluid-vapor transition could be metastable with respect to the fluid-solid transition for chains with more than six LJ monomer units.