An equation of state for fused hard‐sphere polyatomic molecules

Equation of state and Helmholtz energy functional for fused heterosegmented hard chains

Modern equations of state for real nonspherical molecules are often based on Wertheim's first-order thermodynamic perturbation theory (TPT1). A major drawback of TPT1 is that it assumes tangentially bonded spheres. In this work, we develop a Helmholtz energy functional for systems comprising hard heterosegmented chains with arbitrary bond lengths. This is achieved by using hard-sphere fragments (i.e., hard spheres with spherical caps removed at the intersection to their neighbors) as monomers as opposed to full hard spheres. The model is written as a Helmholtz energy functional for inhomogeneous systems and the equation of state for a homogeneous system is determined as a special case. We thereby obtain an equation of state that can be used as a reference to develop statistical associating fluid theory models that more accurately describe the thermodynamic properties of nonspherical molecules. The model is validated against molecular simulation results of bulk pressures and density profiles in slit pores. For the bulk pressures, we show that the equation of state is in excellent agreement with results from molecular simulation for dimers, trimers, and chains of up to 20 segments. The density profiles of individual segments of the chains are regarded in slit pores. Some deviations of the theory from results of molecular simulations are observed for strongly fused chains. Overall, however, good agreement is found for inhomogeneous systems.

Thermodynamic perturbation theory for fused sphere hard chain fluids using nonadditive interactions

A model is developed for the equation of state of fused chains based on Wertheim thermodynamic perturbation theory and nonadditive size interactions. The model also assumes that the structure (represented by the radial distribution function) of the fused chain fluid is the same as that of the touching hard sphere chain fluid. The model is completely based on spherical additive and nonadditive size interactions. The model has the advantage of offering good agreement with simulation data while at the same time being independent of fitted parameters. The model is most accurate for short chains, small values of Delta (slightly fused spheres) and at intermediate (liquidlike) densities.